I am trying to build a linear mixed effect model assuming both random intercept and random slope. The data is like this: There are 3 groups (1 ,2 ,3) and 15 subjects in group 1, 8 subjects in group 2, 8 subjects in group 3. Each patient has been measured 3 times (time: 0, 1, 2). The outcome is continuous. so I assumed normal distribution and identity link function. I wrote code like this: proc glimmix data=xxx; class group time; model outcome = group time group*time / dist=normal link=identity; random intercept time / subject=id type=un; estimate "group 1 - group 2 at time 0" .................... "group 1 - group 3 at time 0".................... "group 2 - group 3 at time 0".................../ adjust=t; run; After I ran the model, I found the interaction term "group*time" is significant. So I want to estimate the differences between groups at time 0 and adjust p-values by tukey. But I am not sure how to specify the values of fixed effect and random effect in the "estimate" statement. can anyone help with this? In addition, if I specify time as linear instead of categorical variable, would it change the way how we specify the values of fixed effect and random effect? Thanks.
... View more