R oats 来自燕麦田间试验的数据

说明

split-plot 田间试验使用三个品种和四个粪肥处理水平得出的燕麦产量。该实验分为 6 个区块，每块 3 个主要地块，每个地块分为 4 个sub-plots。这些品种被应用于主要地块，肥料处理被应用于sub-plots。

用法

oats

格式

该 DataFrame 包含以下列：

B

区块，I、II、III、IV、V 和 VI 层。

V

品种，3个级别。

N

氮(肥料)处理，水平为 0.0cwt、0.2cwt、0.4cwt 和 0.6cwt，显示以 cwt/英亩为单位的应用情况。

Y

每 sub-plot 的产量为 1/4 磅，每个面积 1/80 英亩。

例子

oats$Nf <- ordered(oats$N, levels = sort(levels(oats$N)))
oats.aov <- aov(Y ~ Nf*V + Error(B/V), data = oats, qr = TRUE)
## IGNORE_RDIFF_BEGIN
summary(oats.aov)
summary(oats.aov, split = list(Nf=list(L=1, Dev=2:3)))
## IGNORE_RDIFF_END
par(mfrow = c(1,2), pty = "s")
plot(fitted(oats.aov[[4]]), studres(oats.aov[[4]]))
abline(h = 0, lty = 2)
oats.pr <- proj(oats.aov)
qqnorm(oats.pr[[4]][,"Residuals"], ylab = "Stratum 4 residuals")
qqline(oats.pr[[4]][,"Residuals"])

par(mfrow = c(1,1), pty = "m")
oats.aov2 <- aov(Y ~ N + V + Error(B/V), data = oats, qr = TRUE)
model.tables(oats.aov2, type = "means", se = TRUE)

来源

Yates, F. (1935) 复杂实验，皇家统计学会杂志增刊。 2、181-247。

Yates, F. (1970) 实验设计：Frank Yates, C.B.E, F.R.S. 的论文精选伦敦：格里芬。

参考

Venables, W. N. and Ripley, B. D. (2002) Modern Applied Statistics with S. Fourth edition. Springer.