题目
某研究者研究保健食品对小鼠抗疲劳作用,将同种属的小鼠按性别和年龄相,体重相近配成对子,共10对,将每对中的两只小鼠随机分到保健食品两个不同的剂量组,过一定时期将小鼠杀死,测得其肝糖原含量,见表,问不同剂量组的小鼠肝糖原含量有无差别?
附表
小鼠对号 | 中剂量组 | 高剂量组 |
1 | 620.16 | 958.47 |
2 | 866.50 | 838.42 |
3 | 641.22 | 788.90 |
4 | 812.91 | 815.20 |
5 | 738.96 | 783.17 |
6 | 899.38 | 910.92 |
7 | 760.78 | 758.49 |
8 | 694.95 | 870.80 |
9 | 749.92 | 862.26 |
10 | 794.94 | 805.48 |
答案 n=10, T=6.5.查表的双侧0.02<P<0.05,按α=0.05检验标准,拒绝H0。可以认为该保健食品的不同剂量对小鼠肝糖原含量的作用不同。
用sas编程,运行结果和上述答案大不一样。请教是什么原因?
data prg1;
input x1 x2@@;
d=x1-x2;
cards;
620.16 958.47 866.50 838.42 641.22 788.90 812.90 815.20 738.96 783.17 899.38 910.92 760.78 758.49 694.95 870.80 749.92.862.26 793.94 805.48
;
proc univariate;
var d;
run;
运行结果
Tests for Location: Mu0=0
Test | Statistic | p Value |
Student's t | t 1.93905 | Pr > |t| 0.0937 |
Sign Signed Rank | M -2 S -13 | Pr >= |M| 0.2891 Pr >= |S| 0.0781 |
结果明显不对,程序错在哪里?请高手指教,