Как перевести PRO C MIXED из SAS в R с помощью gls? (MMRM, AR1) - PullRequest
Новая
       35

Как перевести PRO C MIXED из SAS в R с помощью gls? (MMRM, AR1)

1 голос
/ 17 июня 2020

Я пытаюсь получить тот же результат, что и sas с R, но, похоже, возникают некоторые трудности.

Например, R: contras: 2-1 AVISITN = 6: оценка -1,81 и SE 1,59 

library(nlme)
library(emmeans) 
gls <- do.call("gls", list(data=data,
                           model=CHG~TRTPN+AVISITN+TRTPN*AVISITN+BASE+COUNTRY,
                           correlation=corAR1(form=~1|SUBJID),
                           weights=varIdent(form=~1|AVISITN), method="REML" )) 
emm <- emmeans(gls, specs=trt.vs.ctrl ~ TRTPN, at=list(AVISITN=c(6,7,8,9,10)),
               by = "AVISITN", level = 0.4) summary(emm)

SAS: оценка -1,2 и SE 1,24 

proc mixed data=data method=reml;
  class subjid avisit trtp country ;
  model chg = trtp avisit trtpavisit country base /CL SOLUTION DDFM=KR; repeated avisit / subject=subjid type=AR(1) r;
  lsmeans trtpavisit / PDIFF CL alpha=0.40;
  ODS output Diffs=diffs01 lsmeans=lsmeans01;
run;

У меня очень мало опыта использования R, есть ли способ получить тот же результат? Большое спасибо!

...