Как добавить групповой код в парное сравнение в SAS Proc mixed

Question

Когда есть много групп для сравнения, опция lsmeans/pdiff дает слишком много парных сравнений. Мне нужен стол типа

, который показывает, что группа 1,2 отличается от группы 5,6, но не отличается от группы 3,4.

Есть ли возможность в смешанной процедуре SAS или других процедурах сделать это?

Answer 1

Вы ищете опцию LINES в операторе LSMEANS, однако она не работает в PROC MIXED, поэтому вам нужно будет использовать PROC PLM.

proc mixed data=sashelp.class;
   class age;
   model weight=age;
   lsmeans age; *for check;
   store out=classmodel;
   run;
   quit;
proc plm restore=classmodel;
   lsmeans age / lines;
   run;
   quit;

Answer 2

Попробуйте использовать пользовательский формат, который применяется во время процесса.

proc format;
  value groupCat
    1 = 'A'
    2 = 'A'
    3 = 'AB'
    4 = 'AB'
    5 = 'B'
    6 = 'B'
  ;
run;

proc mixed …;
  … ;
  format group groupCat.;
  …