[R] lsmeans

John Fox jfox at mcmaster.ca
Fri Jun 6 19:25:40 CEST 2008

Dear Dani,

I intend at some point to extend the effects package to linear and
generalized linear mixed-effects models, probably using lmer() rather
than lme(), but as you discovered, it doesn't handle these models now.

It wouldn't be hard, however, to do the computations yourself, using
the coefficient vector for the fixed effects and a suitably constructed
model-matrix to compute the effects; you could also get standard errors
by using the covariance matrix for the fixed effects.

I hope this helps,

On Fri, 06 Jun 2008 17:05:58 +0200
 Dani Valverde <daniel.valverde at uab.cat> wrote:
> Hello,
> I have the next function call:
> lme(fixed=Error ~ Temperature * Tumour ,random = ~1|ID,
> data=error_DB)
> which returns an lme object. I am interested on carrying out some
> kind of lsmeans on the data returned, but I cannot find any function
> to do this in R. I'have seen the effect() function, but it does not
> work with lme objects. Any idea?
> Best,
> Dani
> -- 
> Daniel Valverde Saubí
> Grup de Biologia Molecular de Llevats
> Facultat de Veterinària de la Universitat Autònoma de Barcelona
> Edifici V, Campus UAB
> 08193 Cerdanyola del Vallès- SPAIN
> Centro de Investigación Biomédica en Red
> en Bioingeniería, Biomateriales y
> Nanomedicina (CIBER-BBN)
> Grup d'Aplicacions Biomèdiques de la RMN
> Facultat de Biociències
> Universitat Autònoma de Barcelona
> Edifici Cs, Campus UAB
> 08193 Cerdanyola del Vallès- SPAIN
> +34 93 5814126
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.

John Fox, Professor
Department of Sociology
McMaster University
Hamilton, Ontario, Canada

More information about the R-help mailing list