search for: ndrinks

Hi all, In assessing the concordance of two sets of ANOVA results, I would like to extract the p values from the summaries. I had a look at previous answers to this problem, and the suggestions do not seem to work: # toy example test.df<-data.frame(subno=c(1:50,1:50), ndrinks=c(rpois(50,4.5),rpois(50,4)), b4after=c(rep(1,50),rep(-1,50))) test.summary<-summary(aov(ndrinks~b4after+Error(subno),test.df)) test.summary[[2]][,5] Error in test.summary[[2]][, 5] : incorrect number of dimensions and so on. I can get the structure of the object: str(test.summary[[2]]) Lis...