[prev in list] [next in list] [prev in thread] [next in thread] 

List:       r-sig-mixed-models
Subject:    Re: [R-sig-ME] question about linear mixed model
From:       Lao Meng <laomeng.3 () gmail ! com>
Date:       2011-07-28 2:00:56
Message-ID: CANDXbiq5xoTGN7MC9RkpmmbfO+geHRFjh4Uq+WsuC7VDi1bU+g () mail ! gmail ! com
[Download RAW message or body]

Thanks Yuan.
Here's another question:
What's kind of difference does glht tests?
In other words,
If I get 4 intercepts and 4 slopes for all levels of GROUP,then glht tests
for the difference of the 4 intercepts or 4 slopes,or something else?

Thanks.


2011/7/19 Yuan-Ye Zhang <zhangyuanye0706@gmail.com>

> Hi,
>
> 1. Maybe you could try model<- glmmPQL(log10(CD4) ~ time +
> (0+factor(group))+ time:factor(group), random = ~time|id).
> 0+ some times means that you do NOT fit a intercept, but fit each level a
> slope.
>
> 2. Maybe you could try
> library(multcomp),
> glht(model, linfct=mcp(fgroup="Tukey"))
> fgroup<-factor(group).
> The glht function work for "lmer" or "mer", but I am not sure whether it
> will help in your case.
>
> I hope it may help a bit.
> best,
> Yuanye
>
>
> 2011/7/19 Lao Meng <laomeng.3@gmail.com>
>
>> Hi all:
>> I have a question about linear mixed model.
>>
>> my linear mixed model with randomized slope and intercept with interaction
>> of time and group(g1,g2,g3):
>> model<- glmmPQL(log10(CD4) ~ time + factor(group)+ time:factor(group),
>> random = ~time|id)
>>
>> What I get is only the main and interaction of time and group.My question
>> is:
>>
>> 1. How can I get the g1,g2,g3's slope respectively?In other words,how can
>> I
>> get the slope of g1,the slope of g2,and the slope of g3?
>>
>> 2.How can I pair-wise compare the 3 slopes of g1,g2,g3?
>>
>>
>>
>> Thanks a lot!
>>
>>
>> My best
>>
>>        [[alternative HTML version deleted]]
>>
>> _______________________________________________
>> R-sig-mixed-models@r-project.org mailing list
>> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>>
>
>

	[[alternative HTML version deleted]]

_______________________________________________
R-sig-mixed-models@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
[prev in list] [next in list] [prev in thread] [next in thread]