I'm fitting a simple one-way nested anova and would like to obtain standard errors or confidence limits. I'm using se.contrast to try to get se for the contrast between the two levels of the main effect, but I get an error message (see below). What is going on?
Don
> shp<-factor(rep(c("reserve","strip"),each=96))
> site<-factor(rep(c("1g","1p", "1t","2g","2p","2t","3g","3p","3t","4g","4p","4t"),each=16))
> pit<-factor(rep(1:16,12))
> reptsp<-c(4,5,6,4,6,6,6,7,3,5,2,2,4,8,5,4,2,4,2,2,4,5,2,4,4,4,3,2,3,2,5,3,5,3,4,4,4,3,4,3,4,4,4,3,4,3,6,3,3,5,4,6,4,4,2,4,2,6,5,5,5,7,4,4,5,1,4,5,6,5,5,2,6,3,5,6,4,5,4,8,2,4,2,4,2,4,3,3,4,4,3,2,1,3,4,4,2,2,3,2,4,1,2,2,3,4,5,5,3,5,5,4,1,1,2,1,3,1,4,1,6,1,2,3,2,2,2,1,1,2,2,6,5,3,2,3,5,3,2,3,2,1,3,2,4,4,3,3,3,1,2,4,3,4,5,6,5,2,3,2,2,5,5,5,2,2,5,2,4,4,3,2,2,3,2,2,2,2,5,4,3,3,5,2,5,4,3,2,2,2,1,2)
> ddata<-data.frame(shp,pit,site,reptsp)
> #Fit a Standard Nested Anova Model
> repmod1<-aov(reptsp~shp/site/pit+ Error(shp/site/pit))
> summary(repmod1)
Error: shp
Df Sum Sq Mean Sq
shp 1 53.13 53.13Error: shp:site
Df Sum Sq Mean Sq
shp:site 10 61.885 6.189Error: shp:site:pit
Df Sum Sq Mean Sq
shp:site:pit 180 318.56 1.77
> se.contrast(repmod1, list(shp=="strip", shp=="reserve"),data=ddata)
Error in rep.int(n.object - 1, nrow(c.qr) - length(e.assign)) :
invalid number of copies in "rep"______________________________________________ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
