多个实验组与一个对照组均值间两两比较程序的编写！！！

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

有6组数据（没写具体测量结果），例如分别为0，1，2，3，4，5；
后面的5组为实验组，0组为对照组。
我希望计算后面5组分别和对照组的统计分析结果，能够有精确P值。
由于需要具体的P值，查到可能需要SAS中proc glm语句才能实现，而且和对照组两两比较要用到option:adiust=dunnett
但是不知道具体如何编写代码，请大家帮帮忙！在线等！

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：对照组 实验组 Option 统计分析 Dunn 程序

data prg5_1;
do group=1 to 4;
  input x@@;
  output;
end;
cards;
3.6  3.0  0.4  3.3
4.5  2.3  1.8  1.2
4.2  2.4  2.1  1.3
4.4  1.1  4.5  2.5
3.7  4.0  3.6  3.1
5.6  3.7  1.3  3.2
6.0  2.8  3.2  0.6
4.1  1.9  2.1  1.4
5.0  2.6  2.6  1.3
4.5  1.3  2.3  2.1
;
run;
proc sort data=prg5_1;
by group;
run;
proc univariate data=prg5_1 normal;
var x;
by group;
ods html;
run;
proc glm data=prg5_1;
class group;
model x=group;
means group/hovtest dunnett('1') alpha=0.01;
ods html;
run;

自己再根据实际数据情况简单修改就可以了。另，这个用SAS JMP可以很简单实现。

以上自胡良平教授三型理论定理与定性分册，对他表示感谢！

shenweis 发表于 2015-4-19 13:15
data prg5_1;
do group=1 to 4;
  input x@@;

首先表示感谢。另外，你的alpha取值为0.01，这个是调整之后的吗？
能大概解释一下吗？

shenweis 发表于 2015-4-19 13:17
自己再根据实际数据情况简单修改就可以了。另，这个用SAS JMP可以很简单实现。

另外，没有显示出精确的P值。只有标注"***"的方式显示是否有显著性差异。
请问在output中显示P值。

精确P值是有的，但当P很小时，软件将以小于0.001的形式给出。
另外，显著性水平你自己根据需要可以定为0.05或0.01的啊。稍改一下就可以了。

精确P值是有的，但当P很小时，软件将以小于0.001的形式给出。
另外，显著性水平你自己根据需要可以定为0.05或0.01的啊。稍改一下就可以了。

shenweis 发表于 2015-4-19 20:30
精确P值是有的，但当P很小时，软件将以小于0.001的形式给出。
另外，显著性水平你自己根据需要可以定为0.0 ...

                                        The GLM Procedure

                                   Dunnett's t Tests for value

   NOTE: This test controls the Type I experimentwise error for comparisons of all treatments
                                        against a control.


                              Alpha                             0.05
                              Error Degrees of Freedom            11
                              Error Mean Square              0.02349
                              Critical Value of Dunnett's t  2.95000


                 Comparisons significant at the 0.05 level are indicated by ***.


                                        Difference      Simultaneous
                       concentration       Between     95% Confidence
                        Comparison           Means         Limits

                           2 - 1           -0.0413    -0.4540  0.3715
                           3 - 1           -0.1125    -0.4817  0.2566
                           4 - 1           -0.2394    -0.6086  0.1298
                           5 - 1           -0.3525    -0.7216  0.0167
                           6 - 1           -0.4512    -0.8203 -0.0820  ***

这是我运行出的结果，结果显示第6组和对照组1有显著性差异，其他没有；但是我想要这5组比较结果的精确P值。SAS只给出了***表示。