R npk 經典 N、P、K 階乘實驗

說明

在 6 個區塊上進行豌豆生長的經典 N、P、K(氮、磷酸鹽、鉀)析因實驗。其中 3 個圖使用了混雜 NPK 相互作用的部分因子設計的每一半。

用法

npk

格式

npk DataFrame 有 24 行和 5 列：

block

哪個塊(標簽 1 到 6)。

N

施氮指示符 (0/1)。

P

磷酸鹽應用的指示劑(0/1)。

K

鉀施用指標(0/1)。

yield

豌豆產量，以磅/地塊為單位(地塊麵積為 (1/70) 英畝)。

例子

options(contrasts = c("contr.sum", "contr.poly"))
npk.aov <- aov(yield ~ block + N*P*K, npk)
npk.aov
summary(npk.aov)
coef(npk.aov)
options(contrasts = c("contr.treatment", "contr.poly"))
npk.aov1 <- aov(yield ~ block + N + K, data = npk)
summary.lm(npk.aov1)
se.contrast(npk.aov1, list(N=="0", N=="1"), data = npk)
model.tables(npk.aov1, type = "means", se = TRUE)

來源

倫敦帝國理工學院，理學碩士練習表。

參考

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