大家好,请教大家一个问题。我有三组数据,要做方差分析,按不同的visit 输出一些统计量及P值,trtpn是组别,dif是测量值。
1、 要输出组间差异及其95%CI;程序如下:
ods listing close;
proc mixed data=effic;
class visit trtpn;
by visit;
model dif=trtpn;
lsmeans trtpn / cl;
estimate "group1 vs group2" trtpn 1 -1 0 /cl;
estimate "group1 vs group3" trtpn 1 0 -1 /cl;
estimate "group2 vs group3" trtpn 0 1 -1 /cl;
ods output estimates=out1(keep=visit Label Estimate Lower Upper);
run;
ods listing;
该程序是否有问题?
2、 数据同上,要输出三组间比较的统计量与P值;程序如下:
ods listing close;
proc glm data=effic;
class trtpn visit;
by visit;
model dif=trtpn;
ods output overallanova=out2(where=(upcase(source)='MODEL') keep=visit FValue ProbF source);
run;
quit;
ods listing;
该程序是否有问题?
3、 数据同上,要输出两两比较的P值; 程序如下:
ods listing close;
proc glm data=effic;
class visit trtpn;
by visit;
model dif=trtpn;
means trtpn/snk;
estimate "group1 vs group2" trtpn 1 -1 0;
estimate "group1 vs group3" trtpn 1 0 -1;
estimate "group2 vs group3" trtpn 0 1 -1;
ods output estimates=out3(keep=visit Parameter Probt);
run;
quit;
ods listing;
该程序是否有问题?
望大家指教,谢谢!