📄 oats.html
字号:
<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html><head><title>R: Data from an Oats Field Trial</title>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<link rel="stylesheet" type="text/css" href="../../R.css">
</head><body>
<table width="100%" summary="page for oats {MASS}"><tr><td>oats {MASS}</td><td align="right">R Documentation</td></tr></table>
<h2>Data from an Oats Field Trial</h2>
<h3>Description</h3>
<p>
The yield of oats from a split-plot field trial using three varieties and
four levels of manurial treatment. The experiment was laid out in 6 blocks
of 3 main plots, each split into 4 sub-plots. The varieties were applied
to the main plots and the manurial treatments to the sub-plots.
</p>
<h3>Usage</h3>
<pre>
oats
</pre>
<h3>Format</h3>
<p>
This data frame contains the following columns:
<dl>
<dt><code>B</code></dt><dd>Blocks, levels I, II, III, IV, V and VI
</dd>
<dt><code>V</code></dt><dd>Varieties, 3 levels.
</dd>
<dt><code>N</code></dt><dd>Nitrogen (manurial) treatment, levels 0.0cwt, 0.2cwt, 0.4cwt and 0.6cwt,
showing the application in cwt/acre.
</dd>
<dt><code>Y</code></dt><dd>Yields in 1/4lbs per sub-plot, each of area 1/80 acre.
</dd></dl>
<h3>Source</h3>
<p>
Yates, F. (1935) Complex experiments,
<EM>Journal of the Royal Statistical Society Suppl.</EM>
<B>2</B>, 181–247.
</p>
<p>
Also given in
Yates, F. (1970)
<EM>Experimental design: Selected papers of Frank Yates, C.B.E, F.R.S.</EM>
London: Griffin.
</p>
<h3>References</h3>
<p>
Venables, W. N. and Ripley, B. D. (2002)
<EM>Modern Applied Statistics with S.</EM> Fourth edition. Springer.
</p>
<h3>Examples</h3>
<pre>
oats$Nf <- ordered(oats$N, levels = sort(levels(oats$N)))
oats.aov <- aov(Y ~ Nf*V + Error(B/V), data = oats, qr = TRUE)
summary(oats.aov)
summary(oats.aov, split = list(Nf=list(L=1, Dev=2:3)))
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)
</pre>
<hr><div align="center">[Package <em>MASS</em> version 7.2-44 <a href="00Index.html">Index]</a></div>
</body></html>
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -