大家好，请教大家一个问题。我有三组数据，要做方差分析，按不同的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;
该程序是否有问题？

望大家指教，谢谢！