kwright@eskimo.com
2005-May-03 15:17 UTC
[Rd] Bug #7829 reproducible example (rejected by bug tracker)
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