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.