Как перевести код смешанной модели SAS в R? - PullRequest
Новая
       100

Как перевести код смешанной модели SAS в R?

0 голосов
/ 23 апреля 2019

Я пытался применить аналогичную смешанную модель, которую я использовал в SAS для другого анализа в R. Мне нужно понять, как будет написан исходный код в R. Это будет смешанная модель с повторяющимся измерением и случайным множителем.

Вот что я использовал в SAS: 

proc sort data=comp1; by age order subj light clock;
  PROC MIXED DATA = comp1;by age;
  CLASS  subj order light clock;

  MODEL mel = light clock light*clock/DDFM = KR;
  RANDOM subj order/ subject=subj ;
  Repeated clock /TYPE = ar(1)  ;


  LSMEANS light * clock/PDIFF CL;
ods output  diffs=diffcm1 ;
run;

Вы знаете, как написать это в R?Спасибо

...