[R] random error with lme for repeated measures anova

Andrew Robinson A.Robinson at ms.unimelb.edu.au
Wed Aug 27 12:32:28 CEST 2008


Hi JP,

I suggest that you read the literature cited in the help file, most
particularly 

     Pinheiro, J.C., and Bates, D.M. (2000) "Mixed-Effects Models in S
     and S-PLUS", Springer.  

MASS (Modern Applied Statistics in S) also has some useful things in
it.

Cheers

Andrew


On Wed, Aug 27, 2008 at 03:23:58AM -0700, Jean-Pierre Bresciani wrote:
> 
> Hi,
> 
> what is the appropriate syntax to get the random error correct when
> performing repeated measures anova with 'lme'.
> 
> let's say i have 3 independent variables, with 'aov', i would write
> something like: aov(dep_var~(indep_var1*indep_var2*indep_var3) +
> Error(subject/(indep_var1*indep_var2*indep_var3)).
> 
> With 'lme' however, i can't find the right formula. i tried things like:
> lme(dep_var~(indep_var1*indep_var2*indep_var3), random = ~1|subject) or
> nesting my independent variables in 'subject', but those are obviously wrong
> with my design.
> 
> i'm quite clueless (and i haven't found any convincing piece of information
> about how to correctly use 'lme' or 'lmer'). So, any advice along that line
> is more than welcome.
> 
> JP 
> -- 
> View this message in context: http://www.nabble.com/random-error-with-lme-for-repeated-measures-anova-tp19178239p19178239.html
> Sent from the R help mailing list archive at Nabble.com.
> 
> ______________________________________________
> 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.

-- 
Andrew Robinson  
Department of Mathematics and Statistics            Tel: +61-3-8344-6410
University of Melbourne, VIC 3010 Australia         Fax: +61-3-8344-4599
http://www.ms.unimelb.edu.au/~andrewpr
http://blogs.mbs.edu/fishing-in-the-bay/



More information about the R-help mailing list