STAT 5303: Designing Experiments
Spring Semester 2014
Here is the course information handout.
Official University policies on many issues are available
Some good news and some bad news:
The bad news is that the textbook is out of print and hence not readily available from bookstores, although you may find a used copy somewhere.
The good news is that the author (Prof. Oehlert) has generously made the book available as a no-cost download. Go to
for links to the PDF file for the book itself and datafiles for examples and exercises from that book.
You can also buy a hard copy version of that PDF file from Paradigm Copies.
Assignment #8 ( due Friday, April 18):
Assignment #7 (due Wednesday, April 9):
P14.6 Just do parts b, c, and d, but omit part a
These handouts (written by Prof. Oehlert) demonstrate using R to work through examples from the text.
As you work through them, remember you want to see both
(1) how to get the computer to perform certain tasks for you, and
(2) why you want the computer to do those tasks.
In other words, you should learn something about the underlying statistical concepts in addition to learning how to use the computer.
Here is a handout with examples showing how to use R to perform randomization tests.
This R handout uses the resin lifetime data set (Example 3.2 in the text) to demonstrate analysis of a completely randomized design, and also discusses dose-response polynomial modeling.
Our analysis procedures depend on certain assumptions. Some violations of those assumptions are not too serious, but others can cause problems.
This handout discusses how you can detect problems and what you can do about them.
This handout shows how to use R to find power and sample size for the ANOVA F test.
This handout shows how to use R to analyze data from experiments where the treatments are combinations of factors at various levels.
R can also deal with two-series factorials.
Here is a random effects handout. The R command we will use for these models is lmer() (pronounced "elmer") which is part of the "lme4" package. We will also use some functions from the Stat5303 package to perform additional analyses after running lmer().
Here are data files for the calves and resistors examples discussed in the random effects handout.
Here is information about nesting and mixed effects in R and
here is a data file for the particle counts example discussed in this handout.
Here is information about
restricted mixed models in R and here is background information about REML.
Here is a data file for the cheese tasting example discussed in the restricted models handout.
Here is an example of a randomized complete block design in R.
We can analyze Latin square designs in R.
This handout also shows how to deal with an outlier in the data by using a dummy variable to effectively remove that observation from the dataset.
However, doing that also destroys balance, so we need to consider the order in which terms are entered in the model in R. Examples are presented in the handout.
Of course, there is a handout for balanced incomplete block designs.
Naturally, there is a handout for confounding in factorial designs.
Yes, R knows about split-plot and split-split-plot designs.
This handout also deals with repeated-measures designs.
Here are data files for the first (gum) and second (emul2) emulsion examples discussed in the handout.
A dataset for the other example (weeds) is in the oehlert library as emp16.7, or you can find it at
if you have not installed the oehlert package.
R knows about fractional factorial designs.
Basic R information
We will use R.
Here is an introduction to R.
You may download R for Macintosh, Windows, and Linux from the
R package Stat5303
We will use Prof. Oehlert's R package called Stat5303, which
adds some extra commands that we will need.
The Stat5303 package depends on certain other R packages:
The new version of this is actually a pair of packages; one is cfcdae (Companion to a First Course in Design and Analysis of Experiments) and the other is Stat5303libs. Install both packages.
car (version 2.0 or later), effects, FrF2, mvtnorm, perm, RLRsim, and tseries.
You should install these and all the packages they depend on.
Companion (cfcdae) R package (0.8-2) for Mac
Companion (cfcdae) R package (0.8-2) for Windows
Stat5303libs R package (0.7-3) for Mac
Stat5303libs R package (0.7-3) for Windows
Download the packages and save the files into a place where R can find them (e.g., your home directory or the desktop). Start R, set the working directory to that location (e.g., use setwd(), and then use
(The repos=NULL says not to find it online but to look for the package in the local files, and replace XXX by tgz or zip as appropriate.)
install.packages("Stat5303libs_0.7-3.XXX",repos=NULL)to install the other package.
When you start R you then use two library commands: library(cfcdae) and library(Stat5303libs); you do not use library(Stat5303) for the new version.
Because typing in data is tiresome and can
be a source of errors, the data from the homework problems
and examples (from the text) are available as files that R can read.
There are a couple of ways to get data into R.
Individual data files
One approach is a collection of individual data files.
(These are plain text files, so you can also read them with an
editor and cut/paste what is needed into other programs too.)
These files are set up to use the read.table() command in R.
Data sets for examples are named exmpl6.3 for example 3 from chapter 6.
Data sets for exercises are named ex3.1 for chapter 3, exercise 1, and
data sets for problems are named pr3.3 for chapter 3, problem 3.
Here are dataset files for R that you can download.
So, for example---provided your machine is connected to the Internet---from within R you can use
resindata <- read.table("http://www.stat.umn.edu/~corbett/classes/5303/RDataFiles/exmpl3.2",header=TRUE)
to save the data for Example 2 from Chapter 3 as a data frame called "resindata".
You don't even have to type that in; just copy and paste and you're done.
When you are not connected to the Internet, you could use
resindata <- read.table(file="exmpl3.2",header=TRUE)
This would allow you to load data into R without an active Internet connection, but of course you have to already have downloaded the file to your machine.
Another way to get data into R is to load an R package with all the data.
Russ Lenth at the University of Iowa has also provided two R packages that
include the data sets from the book.
For Macintosh/Linux oehlert_1.02.tar.gz
For Windows oehlert_1.02.zip
Download the package and save the file into a place where R can find it
(e.g., your home directory or the desktop).
Start R, set the working
directory to that location (e.g., use setwd(), and then use
(The repos=NULL says not to find it online but to look for the
package in the local files; you need to replace XXX by tar.gz or zip as appropriate.)
Once the package is installed, you can do
from within R to load all of the data. At that point, the command
should give you problem 4 from chapter 17.
Note that the data set names, variable names, and variable codings
in the oehlert data package and the direct-web-accessible data files may not be the same.
This is as much a list of what the project is not as it is a list of what it is.
- There is more information in the text---see chapter 20.
- The project involves a designed experiment, not an observational study or a survey.
- You need to clearly identify experimental units and treatments and, if appropriate, blocking factors and measurement units.
- Describe how you use randomization in your experiment.
- You should analyze the data using methods we have covered in this course, such as ANOVA for a continuous response.
- Remember that we've studied comparative experiments, so you probably want to examine differences between treatments.
- No human subjects, nor anything else that would require review board approval.
- You need to propose a question, plan an experiment to answer that question, perform that experiment, collect and analyze the data, and write a final report covering all of that.
- Remember to include in both the proposal and the final report enough background information so someone who is not trained in your area of expertise can understand what you're trying to do and why someone would want to do it.
- This has to be a new experiment; you cannot just take some old data from work you did before and reanalyze that.
- This is to be an actual experiment, and not just a computer simulation.
- The proposal is not like a contract, but it needs to be specific enough so I can tell tell if what you're planning fits within these guidelines.
Comments? Questions? Send me an e-mail note:
The views and opinions expressed in this page are strictly those
of the page author.
The contents of this page have not been reviewed or approved by the University of Minnesota.