Displaying 1 result from an estimated 1 matches for "mdrfinal".
Did you mean:
md5final
2012 Feb 24
1
code for mixed model in R?
...lysis in R. Here the SAS code and what I wrote in R. It seems to
work but the results are not the same. I don't know how to specify the class
variable in R or specify the variance matrix. Can you please help me?
Thanks
Jurgen
## SAS:
proc glimmix data=trend method=RSPL;
class pid;
model mdrfinal (event = '1') = time therapy comment / dist=binary
link=logit solution cl;
random intercept / subject=pid type=un;
run;
## R:
model_GLMM<- glmmPQL(mdrfinal ~ time + therapy + comment, data = data,
family = binomial, random = ~1| time)
summary(model_GLMM)
--
View this message in con...