Skip to content
Prev 670 / 20628 Next
Original Raw

Split-plot Design

Kevin Wright 
Your question is not very clear, but if you are trying to match the
results in Kuehl, you need a fixed-effects model:

dat <- read.table("expl14-1.txt", header=TRUE)
dat$block <- factor(dat$block)
dat$nitro <- factor(dat$nitro)
dat$thatch <- factor(dat$thatch)

colnames(dat) <- c("block","nitro","thatch","chlor")
m1 <- aov(chlor~nitro*thatch+Error(block/nitro), data=dat)
summary(m1)

Mixed-effects models and degrees of freedom have been discussed many
times on this list....search the archives.

K Wright
On Thu, Mar 20, 2008 at 12:39 PM, <marcioestat at pop.com.br> wrote:

Thread (12 messages)

marcioestat at pop.com.br Split-plot Design Mar 20 Kevin Wright Split-plot Design Mar 20 John Maindonald Split-plot Design Mar 20 Andy Fugard Split-plot Design Mar 21 Robert Kushler Split-plot Design Mar 21 Douglas Bates Split-plot Design Mar 21 John Maindonald Split-plot Design Mar 21 John Maindonald Split-plot Design Mar 22 Andy Fugard Split-plot Design Mar 22 Douglas Bates Split-plot Design Mar 22 marcioestat at pop.com.br Split-plot Design Mar 23 Andrew Robinson Split-plot Design Mar 23