An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-sig-mixed-models/attachments/20110727/fac3da57/attachment.pl>
question about the lme result
3 messages · Lao Meng, ONKELINX, Thierry
Dear Lao, Rewrite your model as lme(sqrt(CD4) ~ 0 + group + group:time,random = ~time|id,data=data_ori) Best regards, Thierry
-----Oorspronkelijk bericht----- Van: r-sig-mixed-models-bounces at r-project.org [mailto:r-sig-mixed-models- bounces at r-project.org] Namens Lao Meng Verzonden: woensdag 27 juli 2011 9:26 Aan: r-sig-mixed-models at r-project.org Onderwerp: [R-sig-ME] question about the lme result Hi all. My data is a longitudinal data,which measure everyone(id)'s CD4 level at some specific time point.All persons are grouped into 3 groups(0,1,2,3). y:sqrt(CD4) x:time,GROUP(3levels: 0,1,2,3) data_ori is my original data,which inclues:CD4,group,time,id. I performed a linear mixed model using lme: result <- lme(sqrt(CD4) ~ group*time,random = ~time|id,data=data_ori)
From the summary(result),I can only get the intercept and slope of
GROUP0(which is described as "intercept" and "time"). If I wanna get intercepts and slopes of GROUP1,GROUP2,GROUP3,how can I do? Simply,I can do this manually by the following steps(GROUP1 for instance): "intercept" + GROUP1 = intercept_GROUP1; "time" + time:GROUP1 = slope_GROUP1; But is there any "shortcut" method to choose?Also,I cann't get the corresponding SE and p value for each parameter.So the above "manually"method is not effective to some degree. Any help from you are welcome. Many thanks! [[alternative HTML version deleted]]
_______________________________________________ R-sig-mixed-models at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-sig-mixed-models/attachments/20110727/66c73fcc/attachment.pl>