search for: tolerance1

Hi, Try either: tolerance <- read.csv("http://www.ats.ucla.edu/stat/r/examples/alda/data/tolerance1.txt") ?aggregate(exposure~male,data=tolerance,mean) ?# male exposure #1??? 0 1.246667 #2??? 1 1.120000 #or ?library(plyr) ?ddply(tolerance,.(male),summarize,exposure=mean(exposure)) #? male exposure #1??? 0 1.246667 #2??? 1 1.120000 #or with(tolerance,tapply(exposure,list(male),FUN=mean))...

I am having a hard time understanding how to perform a "repeated measures" type of ANOVA with R. When reading the document found here: http://cran.r-project.org/doc/contrib/Lemon-kickstart/kr_repms.html I find that there is a reference to a function make.rm () that is supposed to rearrange a "one row per person" type of frame to a "one row per observation" type

...havior. For example, it works for the first interaction.plot below, with the example from the ALDA book, but not with the other plots, from the NPK dataset: # from http://www.ats.ucla.edu/stat/r/examples/alda/ch2.htm tolerance <- read.csv("http://www.ats.ucla.edu/stat/r/examples/alda/data/tolerance1.txt") tolerance.pp <- read.csv("http://www.ats.ucla.edu/stat/r/examples/alda/data/tolerance1_pp.txt") attach(tolerance.pp) # fitting the linear model by id fit <- by(tolerance.pp, id, function(bydata) fitted.values(lm(tolerance ~ time, data=bydata))) fit <- unlist(fit) # pl...