The R Project SVN – R – Diff – /trunk/src/library/stats/man/se.contrast.Rd



summary(fit, split = list(B = 1:2), expand.split = TRUE)
## t.stat^2 is the F value on the A:B: C1 line (with Helmert contrasts)
## Now look at all three interaction contrasts
cont <- c(1, -1)[A] * cbind(c(1, -1, 0), c(1, 0, -1), c(0, 1, -1))[B,]
se.contrast(fit, cont)  # same, due to balance.
rm(A, B, score)
 
 
## multi-stratum example where efficiencies play a role
utils::example(eff.aovlist)
fit <- aov(Yield ~ A + B * C + Error(Block), data = aovdat)

Rev 61160	Rev 61168
Line 94...	Line 94...
94	`summary(fit, split = list(B = 1:2), expand.split = TRUE)`	94	`summary(fit, split = list(B = 1:2), expand.split = TRUE)`
95	`## t.stat^2 is the F value on the A:B: C1 line (with Helmert contrasts)`	95	`## t.stat^2 is the F value on the A:B: C1 line (with Helmert contrasts)`
96	`## Now look at all three interaction contrasts`	96	`## Now look at all three interaction contrasts`
97	`cont <- c(1, -1)[A] * cbind(c(1, -1, 0), c(1, 0, -1), c(0, 1, -1))[B,]`	97	`cont <- c(1, -1)[A] * cbind(c(1, -1, 0), c(1, 0, -1), c(0, 1, -1))[B,]`
98	`se.contrast(fit, cont) # same, due to balance.`	98	`se.contrast(fit, cont) # same, due to balance.`
99	`rm(A,B,score)`	99	`rm(A, B, score)`
100		100
101		101
102	`## multi-stratum example where efficiencies play a role`	102	`## multi-stratum example where efficiencies play a role`
103	`utils::example(eff.aovlist)`	103	`utils::example(eff.aovlist)`
104	`fit <- aov(Yield ~ A + B * C + Error(Block), data = aovdat)`	104	`fit <- aov(Yield ~ A + B * C + Error(Block), data = aovdat)`

The R Project SVN R

(root)/trunk/src/library/stats/man/se.contrast.Rd – Rev 61160 → 61168