当配对设计资料的差值不满足正态分布时,用Wilcoxon配对符号秩和检验进行计算,发现SPSS结果和SAS结果不一致,如何调整SAS程序,使得与SPSS的结果一致呢?
如举个例子:某研究者研究保健食品对小鼠抗疲劳作用,将同种属的小鼠按性别和年龄相同、体重相近配成对子,共10对,将每对中的两只小鼠随机分到保健食品两个不同的剂量组,过一定时期将小鼠杀死,测得其肝糖原含量(mg/100g),
见下表
问:保健食品是否对小鼠的肝糖原含量有影响?
SAS程序如下: 1./*建立数据集*/
data demo; input mid high; d=high-mid; cards; 620.16 958.47 866.50 838.42 641.22 788.90 812.91 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 ; run; 2.分析:差值d的正态性检验和符号秩和检验
proc univariate data=demo normal; var d; run;
SAS的结果P=0.0313,而SPSS中“分析” -- “非参数检验” -- “旧对话框” -- “2个相关样本” -- “将中剂量组和高剂量组均放入变量列表中“ -- 勾选wilcoxon(W)检验 -- 点击”确定“ -- 得到检验结果。P=0.032;
如何让SAS的结果也等于0.032???谢谢各位大神 |