On Tue, 24 Oct 2006, A Ezhil wrote: > Hi All, > > I have performed a 3-way ANOVA analysis for my > experimental data using aov function. My simple R > funtion for this is: > > 3aof <- function(x){ > m <- data.frame(R,S,T, x); > anova(aov(x ~ R+S+T+R*S+R*T+S*T+R*S*T, m) ) > } >
you can use R> library("multcomp") ### check for version 0.991-1 for some model R> amod <- aov(x ~ R+S+T+R*S+R*T+S*T+R*S*T, m) with, say, all-pair comparisons of the levels of `T': R> posth <- glht(amod, linfct = mcp(T = "Tukey")) various adjusted p-values can be computed using R> summary(posth) and confidence intervals (adjusted and unadjusted) are available from R> confint(posth) Hope that helps, Torsten > Now, I am getting P values for all the main and > interactions effects. If I want to perform postdoc > test on one of my main effects, say T, what method I > should use (i have unequal sample size)? Is there > anyway I can extend my above '3aof' function to do > this? > > I greatly appreciate your help. Thanks in Advance. > > Regards, > Ezhil > Ph.D Student > School of Biotechnology > AU, India. > > ______________________________________________ > R-help@stat.math.ethz.ch 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. > > ______________________________________________ R-help@stat.math.ethz.ch 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.