Re: [R] Re peated Measure ANOVA-Old question

From: <ctu_at_bigred.unl.edu>
Date: Wed, 30 Jul 2008 02:05:21 -0500

Tank you (Anna & Mark) very much for this super information Let me confirm that IF I want to perform the RM-ANOVA I should use "lmer" and perform the post-hoc test by using "glht", right?

Because I am not so familiar with "lmer" so I have two more questions. Let me use that following example again, 1. Is the random effect statement correct? if no what will be? 2. What did I do wrong in glht?

require(lme4);require(nlme)
group <-rep(rep(1:2, c(5,5)), 3);time <-rep(1:3, rep(10,3));subject <- rep(1:10,3)
p.pa <- c(92, 44, 49, 52, 41, 34, 32, 65, 47, 58, 94, 82, 48, 60, 47, 46, 41, 73, 60, 69, 95, 53, 44, 66, 62, 46, 53, 73, 84, 79)

P.PA <- data.frame(subject, group, time, p.pa)
P.PA$group=factor(P.PA$group);P.PA$time=factor(P.PA$time)
P.PA$subject=factor(P.PA$subject)

tlmer<-lmer (p.pa ~ time * group + (time*group | subject), data=P.PA )

> summary(glht(tlmer,linfct=mcp(group="Tukey")))
Error in UseMethod("fixef") : no applicable method for "fixef" In addition: Warning message

Thank you very much
Chunhao



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. Received on Wed 30 Jul 2008 - 07:08:19 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Wed 30 Jul 2008 - 08:33:01 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive