On this Picostat.com statistics page, you will find information about the oats data set which pertains to Data from an Oats Field Trial. The oats data set is found in the MASS R package. You can load the oats data set in R by issuing the following command at the console data("oats"). This will load the data into a variable called oats. If R says the oats data set is not found, you can try installing the package by issuing this command install.packages("MASS") and then attempt to reload the data. If you need to download R, you can go to the R project website. You can download a CSV (comma separated values) version of the oats R data set. The size of this file is about 2,155 bytes.
Data from an Oats Field Trial
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.
This data frame contains the following columns:
Blocks, levels I, II, III, IV, V and VI.
Varieties, 3 levels.
Nitrogen (manurial) treatment, levels 0.0cwt, 0.2cwt, 0.4cwt and 0.6cwt,
showing the application in cwt/acre.
Yields in 1/4lbs per sub-plot, each of area 1/80 acre.
Yates, F. (1935) Complex experiments,
Journal of the Royal Statistical Society Suppl.
Also given in
Yates, F. (1970)
Experimental design: Selected papers of Frank Yates, C.B.E, F.R.S.
Venables, W. N. and Ripley, B. D. (2002)
Modern Applied Statistics with S. Fourth edition. Springer.
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, split = list(Nf=list(L=1, Dev=2:3)))
par(mfrow = c(1,2), pty = "s")
abline(h = 0, lty = 2)
oats.pr <- proj(oats.aov)
qqnorm(oats.pr[][,"Residuals"], ylab = "Stratum 4 residuals")
qqline(oats.pr[][,"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)
Dataset imported from https://www.r-project.org.