Hi Paul, Perhaps
# Data set set.seed(123) x=rnorm(100) y=x+2*rnorm(100) # ANOVA AOV=anova(lm (y ~ x)) p1=AOV$"Pr(>F)"[1] p1 or # ANOVA 2 AOV1=aov(y ~ x) p2=unlist(summary(AOV1))["Pr(>F)1"] p2 HTH, Jorge On Wed, May 7, 2008 at 8:47 PM, Paul Benton <[EMAIL PROTECTED]> wrote: > hello all, > > Quick question, how do I get the p value out of the anova? > > Thanks, > > Paul > > > pb<-aov(as.numeric(diff[5,16:33]) ~ grF) > > summary(pb) > Df Sum Sq Mean Sq F value Pr(>F) > grF 3 2.7860e+10 9.2867e+09 4.2236 0.02534 * > Residuals 14 3.0783e+10 2.1988e+09 > --- > Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 > > str(summary(pb)) > List of 1 > $ :Classes 'anova' and 'data.frame': 2 obs. of 5 variables: > ..$ Df : num [1:2] 3 14 > ..$ Sum Sq : num [1:2] 2.79e+10 3.08e+10 > ..$ Mean Sq: num [1:2] 9.29e+09 2.20e+09 > ..$ F value: num [1:2] 4.22 NA > ..$ Pr(>F) : num [1:2] 0.0253 NA > - attr(*, "class")= chr [1:2] "summary.aov" "listof" > > attr(summary(pb), "Pr(>F)") > NULL > > ______________________________________________ > R-help@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. [[alternative HTML version deleted]] ______________________________________________ R-help@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.
