Блог

Главная - Вопросы по программированию - как перевести код смешанной модели SAS в R?

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

#r #model #mixed

#r #Модель #смешанный

Вопрос:

Я попытался применить аналогичную смешанную модель, которую я использовал в 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?
Спасибо

Комментарии:

1. Это может помочь: bbolker.github.io/mixedmodels-misc/notes/corr_braindump.html

2. На этой странице показано множество кодов R и SAS (бок о бок), которые выполняют один и тот же анализ; Я нахожу это очень полезным: sas-and-r.blogspot.com/p/statistics-examples.html .