From kwright at eskimo.com  Tue May  3 15:17:49 2005
From: kwright at eskimo.com (kwright@eskimo.com)
Date: Tue May  3 15:17:58 2005
Subject: [Rd] Bug #7829 reproducible example (rejected by bug tracker)
Message-ID: <12745.170.54.59.167.1115126269.squirrel@170.54.59.167>


I submitted a bug report (#7829) that is in the bug tracking system, but
bounced off the email list for some reason.  Since then, I've come up with
a reproducible example:

Oats <- Oats[-1,]
Oats$nitro <- as.factor(Oats$nitro)
Oats.aov <- aov(yield~Block+Variety+nitro+Variety:nitro,data=Oats)
model.tables(Oats.aov,"means",cterms=c("nitro","Variety:nitro"))

Gives the following results.  The bug is in the number of replications for
the interaction term.  See the bug tracking system for more details.


Tables of means
Grand mean

103.8732

 nitro
        0  0.2   0.4 0.6
    78.74 98.5 113.8 123
rep 17.00 18.0  18.0  18

 Variety:nitro
             nitro
Variety       0      0.2    0.4    0.6
  Golden Rain  79.53  98.03 114.20 124.37
  rep          79.53  98.03 114.20 124.37
  Marvellous   86.20 108.03 116.70 126.37
  rep          86.20 108.03 116.70 126.37
  Victory      69.77  89.20 110.37 118.03
  rep          69.77  89.20 110.37 118.03


[Rd] Bug #7829 reproducible example (rejected by bug tracker)
kwright@eskimo.com