R help - May 2003 - variance components

Dear All,

I need to calculate the variance components in a mixed effect model (one
fixed and one random effect) with REML (maximizing the proportion of
the likelihood that does not depend on the fixed effects). In S+ there is
the varcomp function, but I would like to do it in R. Is there a way to do
that?

Thanks!
Katalin


___
Katalin Csillery
Division of Biological Sciences
University of Montana, Missoula MT 59801
Phone: 406 243 6106, E-mail: csillery at selway.umt.edu

Have you considered "lme" (documented especially in Pinhiero and
Bates,
Mixed Effects Models in S and S-Pus)?

hth.  spencer graves

Katalin Csillery wrote:
> Dear All,
> 
> I need to calculate the variance components in a mixed effect model (one
> fixed and one random effect) with REML (maximizing the proportion of
> the likelihood that does not depend on the fixed effects). In S+ there is
> the varcomp function, but I would like to do it in R. Is there a way to do
> that?
> 
> Thanks!
> Katalin
> 
> 
> ___
> Katalin Csillery
> Division of Biological Sciences
> University of Montana, Missoula MT 59801
> Phone: 406 243 6106, E-mail: csillery at selway.umt.edu
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://www.stat.math.ethz.ch/mailman/listinfo/r-help