[R] similar limma's contrasts.fit() for lme (mixed effect model)object
shirley0818 at gmail.com
Mon Jul 9 17:49:40 CEST 2007
Thanks for your prompt reply. Actually I have already checked contrast
package, but I still could not figure out how to set the contrast
matrix for a 2x2 factorial design.
I would like to set a contrast exactly similar to the limma's user
guide, 8.7 Factorial Designs (page 45). For example, after a mixed
linear model fit with the default contr.treatment, I want to compute
an estimated coefficient as the sum of the third and fourth
coefficients of the original lme model (2x2 factorial design).
Could you give me an example how to implement this using the contrast package?
On 7/9/07, Kuhn, Max <Max.Kuhn at pfizer.com> wrote:
> The contrast package can do this. The method of specifying the contrast
> conditions/coefficients is different (and I think easier).
> -----Original Message-----
> From: r-help-bounces at stat.math.ethz.ch
> [mailto:r-help-bounces at stat.math.ethz.ch] On Behalf Of shirley zhang
> Sent: Monday, July 09, 2007 11:16 AM
> To: R-help at stat.math.ethz.ch
> Subject: [R] similar limma's contrasts.fit() for lme (mixed effect
> Dear R help,
> In limma package, contrasts.fit() function is very useful. I am
> wondering whether there is a similar function for lme object, which
> means given a mixed linear model fit, compute estimated coefficients
> and standard errors for a given set of contrasts.
> R-help at stat.math.ethz.ch mailing list
> PLEASE do read the posting guide
> and provide commented, minimal, self-contained, reproducible code.
> LEGAL NOTICE
> Unless expressly stated otherwise, this message is confidential and may be privileged. It is intended for the addressee(s) only. Access to this E-mail by anyone else is unauthorized. If you are not an addressee, any disclosure or copying of the contents of this E-mail or any action taken (or not taken) in reliance on it is unauthorized and may be unlawful. If you are not an addressee, please inform the sender immediately.
More information about the R-help