[R-sig-teaching] Analyze Nested and Split-Plot Designs with R

[Steven Stoline]

Dear All:

example 14.1, Montgomery, chapter 14. *Supplier* is a *fixed factor*,
*Batches* is a *random factor* nested within the fixed factor Supplier.

How to analyze these data in R in two ways:

1- Nested Design

fit <- aov(purity~Supplier/Batch)


it did not give me the expected output.


2- as a factorial (suppliers Fixed, Batches Random)

fit.out <- aov(Purity~Supplier*Batch, data=have)

it did not give me the expected output.


Here is the data set:
===============

> data
  Supplier Batch  Purity
1     1  1
1     1 -1
1     1  0
1     2 -2
1     2 -3
1     2 -4
1     3 -2
1     3  0
1     3  1
1     4  1
1     4  4
1     4  0
2     1  1
2     1 -2
2     1 -3
2     2  0
2     2  4
2     2  2
2     3 -1
2     3  0
2     3 -2
2     4  0
2     4  3
2     4  2
3     1  2
3     1  4
3     1  0
3     2 -2
3     2  0
3     2  2
3     3  1
3     3 -1
3     3  2
3     4  3
3     4  2
3     4  1


many thanks
Steven

-- 
Steven M. Stoline
1123 Forest Avenue
Portland, ME 04112
sstoline at gmail.com

	[[alternative HTML version deleted]]