EVOLUTION-MANAGER
Edit File: oats.html
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Strict//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd"><html xmlns="http://www.w3.org/1999/xhtml"><head><title>R: Data from an Oats Field Trial</title> <meta http-equiv="Content-Type" content="text/html; charset=utf-8" /> <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 style="text-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: </p> <dl> <dt><code>B</code></dt><dd> <p>Blocks, levels I, II, III, IV, V and VI. </p> </dd> <dt><code>V</code></dt><dd> <p>Varieties, 3 levels. </p> </dd> <dt><code>N</code></dt><dd> <p>Nitrogen (manurial) treatment, levels 0.0cwt, 0.2cwt, 0.4cwt and 0.6cwt, showing the application in cwt/acre. </p> </dd> <dt><code>Y</code></dt><dd> <p>Yields in 1/4lbs per sub-plot, each of area 1/80 acre. </p> </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 style="text-align: center;">[Package <em>MASS</em> version 7.3-51.4 <a href="00Index.html">Index</a>]</div> </body></html>