请高手指教sas编程非参数检验问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

题目

某研究者研究保健食品对小鼠抗疲劳作用，将同种属的小鼠按性别和年龄相，体重相近配成对子，共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

 

结果明显不对，程序错在哪里？请高手指教,

 

 

 

[此贴子已经被作者于2008-6-8 11:51:43编辑过]

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：非参数检验 SAS编程 参数检验 非参数 保健食品 检验 高手 SAS 指教 参数

data prg1;
input x c @@;
cards;
620.16 1
958.47 2
866.50 1
838.42 2
641.22 1
788.90 2
812.90 1
815.20 2
738.96 1
783.17 2
899.38 1
910.92 2
760.78 1
758.49 2
694.95 1
870.80 2
749.92 1
862.26 2
793.94 1
805.48 2
;

proc npar1way wilcoxon;
var x;
class c;
run;

楼主数据集输入有误，你可以看看prg1数据集中的数据，最后一个观测缺失。

对edragon1983给出的方法做了一下，还是不行？

另外，这是配对资料的秩和检验，用npar1way可以吗？

[此贴子已经被作者于2008-6-8 22:37:46编辑过]

univariate过程名是单一样本的描述性统计过程，能用来检验Ｕ是否等于Ｕ０、检验样本数据是否来自正态总体等。　　对于楼主提出的问题应该用npar1way过程加median或者wilcoxon选项进行检验．还可以加上

exact median 或（exact wilcoxon)利用蒙可卡洛方法作出精确检验及相应的Ｐ值！　　

                                     Tests for Location: Mu0=0

                          Test           -Statistic-    -----p Value------

                          Student's t    t  -2.25921    Pr > |t|    0.0502
                          Sign           M        -3    Pr >= |M|   0.1094
                          Signed Rank    S     -21.5    Pr >= |S|   0.0254

看这个结果

程序：

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 mu0=0;
var d;
run;

多谢各位知道。问题解决。