[R] (no subject)

[Don Driscoll]

G'day,

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.13

Error: shp:site
          Df Sum Sq Mean Sq
shp:site 10 61.885   6.189

Error: 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"