(Anderson) Design of Experiments Realistic Approach
(Anderson) Design of Experiments Realistic Approach
A Realistic Approach
STATISTICS
A SERIES EDITED BY
D. B. O W E N , Coordinating Editor
Department of Statistics
Southern Methodist University
Dallas , Texas
PETER LEWIS P A U L D. M I N T O N
Naval Postgraduate School Virginia Commonwealth University
Monterey ,California Richmond ,Virginia
J O H N W. P R A T T
Harvard University
Boston, Massachusetts
O t h e r Vo l u m e s in P r e p a r a t io n
DESIGN OF EXPERIMENTS
A Realistic Approach
IT * __*1 V A 1 _ 1 n L I I II V
^
Department of Statistics Department of Statistics
Purdue University University of Tennessee
Lafayette, Indiana Knoxville, Tennessee
MARCEL
WL M a r c e l D e k k e r, In c. N e w Y o rk • B a s e l • H o n g K o n g
COPYRIGHT © 1974 by MARCEL DEKKER, INC. ALL RIGHTS RESERVED.
ISBN 978-0-8247-6131-8
20 19 18 17 16
PROBLEMS xi
PREFACE xiii
v
CONTENTS
APPENDICES 388
Problem
Number Applied Area Page
1.1.1 Civil Engineering 6
1.3.1 Pharmacy 21
1.3.2 Pharmacy 24
1.3.3 Pharmacy 27
1.4.1 Engineering 37
1.4.2 General 37
1.4.3 General 37
4.1.1 Engineering 97
4.2.1 Engineering 98
4.2.2 Animal Nutrition 110
4.2.3 General 110
4.2.4 Agronomy 110
4.2.5 General 111
4.2.6 General 111
4.2.7 Engineering 116
4.2.8 Engineering 116
4.2.9 Engineering 119
4.2.10 Genetics 120
xi
7.1.2 Medical 198
7.2.1 Alloy (Engineering) 202
7.2.2 Battery (Engineering) 203
7.2.3 Steel (Engineering) 204
7.3.1 Highway (Engineering) 206
xiii
xiv PREFACE
The book is arranged so that the reader may go from the simple
to the complex designs and grasp the appropriate analyses from the
resulting data. In Chapter 1 we review analyses for one way
classified data and in Chapter 2 we introduce two way analyses.
Random sampling concepts are used in these two chapters, assuming
PREFACE xv
1
2 1. BASIC STATISTICAL CONCEPTS
confounds the two factors. ’’Confounds’1 merely means that the effects
of the two factors are not separable.
TABLE 1.1
4% cement 8% cement
and and
21 days curing time 7 days curing time
(psi) (psi)
70 89
81 79
82 85
68 76
74 74
76 81
65 70
74 79
60 65
70 72
10 10
I y ,• = 720 = 770
i = i 11
yx = 72 = 77
^2
10 - 9 10 9
I (yu - y p - 422 = 460
l& n - y2J
i=l 11 1
_ yi ' X2 _ 72 - 77 I 5
tdf = S - - -,v 3.13 =
yl ' y 2 /49(°-2)
10 9 10
I (y li - yP + ^
i=l 1=1
n^ + n^ - 2
10 9 10 9
I Cyu - vx) + I (y2i -
i=l 11 1 i=l _ 882 _
n, + n - 2 18 ” 4y
y = overall mean
TABLE 1.2.1
A 6, 2, 5, 4, 6, 7 30 5.0
B 10, 8, 11, 7, 7, 9 52 8.7
C 3, 7, 6, 4, 8, 6 34 5.7
D 10, 4, 6, 6, 7, 8 41 6.8
TABLE 1.2.2
(0.05)
Source df SS MS F
crit
Total 23 114.0
Conclusions:
find this out. One test that allows investigation of all possible
pairs of means in a sequential manner, has very good power (1-3)
and keeps the a level constant for investigation of all pairs of
means is the Newman-Keuls test. The two basic references for this
test (See the preface for a philosophical discussion of multiple
comparison tests) are Newman (1939) and Keuls (1952). In addition
to the Newman-Keuls test the researcher is apt to encounter other
multiple comparison procedures such as: the least significant
difference, Scheffe's method, Tukeyfs honestly significant difference,
and Duncan's multiple range test. Two excellent references which
compare these various methods are Bancroft (1968, Chapter 8), and
Carmer and Swanson (1973).
1 2 3 4
B D C A
8.7 6.8 5.7 5.0
where column headings indicate the number of means spanned for each
test on that diagonal.
’ -y ' ■ / ¥ ■
*k_i - qa0c-Mf) s-
R2 = %(2,df) s-
(0.025) _
20 df “ 2*086
or
(7.1 < y < 10.3)
Observation
The symbol T^# represents the total of all observations for the ith
treatment and the "dot" indicates a summation over the subscript that
has been replaced by the dot, that is T2# = y^y In addition,
y\# = T^#/n, the sample mean of the i ^ treatments. Note that the
grand total of all observations is T#> and the grand mean is
y
• • • = T • •/tn.
1 n 7
SS total =I I (y. . - y . J
i=l j=l 13
2
By expanding ][ J (y.^ - y ##) this expression can be seen to be equal
t n ?
SS total *I I yf. - CT
i=l j=l 13
In addition, by writing
t n
+ I I2(y.. -yi#)(y. - y.J
i=l j=l ij i i
Note that in the last term of the above equation that 2 (y\# - y ##) is
a constant with respect to the inner summation over j and that
l (y ^ - yj.) = o
j= i 13 1
consequently
Observing that
= SS treatments
we see that
t n 2
SS total = I I ( Y ii - Yim) + SS treatments
i=l j=l 1J 1
so that
t n 2
SS error = I I (y.. - y. )
i=l j=l 1J 1
It should be noted that the above data array and the algebraic
relation among the SS total, SS treatments, and SS error was presented
assuming an equal number of observations per cell. This assumption
was only made for ease of notation and is not a mathematical
requirement. This statement remains valid for the following
discussion on the F-test utilized to test the hypothesis of no
significant differences among treatment means.
The F-test is well known for its use in testing for equality of
two variances. This test is carried out by computing the ratio of
two independent variance estimates. If this ratio is greater than
16 1. •BASIC STATISTICAL CONCEPTS
, t 2 1
rr r I (Xi. - y.J = f7ft~ T T css treatments)
i=l v
= i (MS treatments)
2
and hence must be an estimate of a /n so that MS treatments is an
2
estimate of a under the null hypothesis. Intuitively, it has been
demonstrated that the ratio
MS treatments
MS error
*
Robust for a in this situation means that even if the variances
are heterogeneous the a-level for the F-test on treatment means is not
changed materially.
18 1. BASIC STATISTICAL CONCEPTS
TABLE 1.3.1
Types of tablets
Mean 28 13 82 154
Variance 59.55 26.22 430.89 436.22
1.3 UNEQUAL VARIANCES 19
2
but through the log s technique one can also study the relationship
between the factors and variability of the process.
TABLE 1.3.2
1 A 59.55 1.77488
2 B 26.22 1.41863
3 C 430.89 2.63437
4 D 436.22 2.63971
Total 952.88 8.46759
T Zs2 952.88
s = 238.22
k 4
log s2 = 2.37698
2 M 21.56 *
*3df = C “ 0 4 6 " 20-6
and
xcrit 3df(°-001) = 16-3
TABLE 1.3.3
2
Log s Analysis of Homogeneity
Type of tablet
A B C D
2 . s2 2 . 2 2 2 2 2
s log s log s s log s s log s
59 1.7709 10 1.0000 449.5 2.6527 493.0 2.6928
52.5 1.7202 39 1.5911 510.0 2.7076 478.5 2.6799
2
ANOVA of log s
Source df MS F (0.01)
^calc crit(3,4)
Types of tablet 3 0.96876 21.8 16.7
Within error _4_ 0.04439
Total 7
2
For unequal sample sizes we specify that each sample variance, s^,
be calculated by dividing by the degrees of freedom, v^, rather than
by the sample size, n^ (where = n^ - 1, for i = 1, ..., p ) . Let
v denote the arithmetic average of the degrees of freedom. In this
case we have:
2 4
Types_______ s___________ s
A 59.55 3,546.20
B 26.22 687.49
C 430.89 185,666.19
D 436.22 190,287.89
Total 952.88 380,187.77
= 380,187,77 „ Q>41fl
q =
■ (952.88)
(°;0D a 0 398
ncrit
q <°:?01) - 0 . 4 M
ncrit
1.3.4 Transformation of y
TABLE 1.3.4
Type log [ s ^ ]
A 0.538 -0.2692
B 0.441 -0.3566
C 1.279 +0.1069
D 0.734 -0.1343
Total 2.992 -0.6522
Problem 1.3.2. (a) Carry out the Q-test on the data furnished
in Table 1.3.4. (b) Run the ANOVA on Jy in Table 1.3.1, p. 18
and explain your results, (c) Determine the 95% confidence
interval for the mean of the Type D population.
2. Compute E(yi - y) .
k
b = V i + i (V i +i - y p
k
b = Vi+i (V i +i -
1. Ordered observations: 8, 9, 10, 10, 10, 12, 12, 16, 19, 24.
Problem 1.3.5. Use the W test to test the data utilized in the
above example for normality after carrying out the square root
transformation.
k k
where the Zfs are functions of X and the afs are the corresponding
coefficients.
Z = Z
1 linear
2
= X2 [(X - X)2 -
In order to obtain the estimate of the response surface one must also
determine estimates of the coefficients of the polynomials. For
example,
a.
1
When one recalls that we are going to assume that the X values
take on the values 1, 2, ..., n, it is readily seen that one is able
to create a table of Z values in order to simplify calculations.
Appendix 11 gives an abbreviated table of Z's along with the
2
corresponding Xfs and E(Z) . The case where n = 6 is reproduced in
Table 1.4.1 for illustrative purposes.
TABLE 1.4.1
Notice for example that these values are orthogonal, e.g., EZ^ Z^^
® (-5)(-5) + (-3)(7) + ... + (5)(5) = 0. In analyzing a set of
data for trends of various degrees it becomes necessary to be able
to compute the sums of squares which is explained by each of the
terms included in the model of Eq. (1.4.2). For example, the
formula required for the sum of squares (SS) for the quadratic term
is
SSquadratic = = “l l ^ l P
The same procedure is used for any degree terms given in Eq. (1.4.2).
The following example concerns the observation of hardness as
measured every 2 centimeters across a block of synthetic material.
The data and the method of calculating the afs are shown in
Table 1.4.2.
1.4 CURVE FITTING 31
TABLE 1.4.2
x» y X
Z1 ^1 Z11 yZll Zlll ^Zlll
2 16 1 -5 -80 5 80 -5 -80
4 14 2 -3 -42 -1 -14 7 98
o
6 10 3 -1 -10 -4 4 40
00
i
8 12 4 1 12 -4 -48 -4
10 13 5 3 39 -1 -13 -7 -91
12 15 6 5 75 5 75 5 75
80 -6 40 -6
O
olvo
O1
13.33 a -6
ii
ao = -0.0857
1 70 ~
volo
40
1l
0.4762 -0.0333
II
ii
|rH
®ir 84 ~ alll
oo
= 13.93 - 0.1714 X
= 13.93 - 0.1714
TABLE 1.4.3
Source df SS MS R2
,514)
Linear regression 1 0.514 0.514 0.022 = | 2 ^;
,333)
Error 4 22.819 5.705
Total 5 23.333
y = y + «i z i + ° n zn
2
= 13.33 - 0.0857 [X (X-Y)] + 0.4762 X2[(X-X)2 -
TABLE 1.4.4
Source df SS MS F R2
Linear 1 0.514
Quadratic _1 19.048
2
The figures in Table 1.4.4 still do not indicate a significant R
2
when one uses the F-test; however, the relatively large jump in R
from 0.022 for the linear model to 0.838 for the quadratic model would
probably lead one to consider the quadratic model as the appropriate
model. If one continued this example and calculated the cubic sum of
squares one would not find much of an increase in the regression sum
of squares. In order to obtain any further information on the fit of
this model we need additional observations at the same X values which
is the topic of the next section.
m.
■i ( I Xj)2
y v 2 j=l
SSpooled error “ x yj m7 ^
n
SS due to the quadratic term = a.,(m J y. Z.J
i=l 1 L
The ANOVA for this example appears in Table 1.4.6. Note that
the lack of fit is carried out twice and in general it would be
carried out on successive terms until the hypothesis of an adequate
model is accepted or there are no degrees of freedom remaining for
lack of fit. One must recognize that this testing procedure involves
36 1. BASIC STATISTICAL CONCEPTS
TABLE 1.4.5
79
13.17
ao 55 6 ~
"4h-*
11
o|o
=
1!
-0.14
al
= 36.5 = 0.43
all 84 *
2 X (-10)2
SSCap SSL = 2.86
70
2 X (36.5)2
SS(on ) SSQ = = 31.72
84
TABLE 1.4.6
Source df SS MS F R2
y = °F
Positions
1 2 3 4 5 6 7
1750 1700 1730 1710 1630 1665 1775
1775 1750 1700 1665 1660 1690 1710
1730 1720 1690 1705 1645 1710 1735
1725 1725 1700 1680 1600 1705 1740
1770 1705 1680 1690 1615 1690 1740
1.5 REFERENCES
Anderson, R. L. and Bancroft, T. A. Statistical Theory in Research
McGraw-Hill, New York, 1952.
Bancroft, T. A. Topics in Intermediate Statistical Methods, Iowa
State University Press, Ames, Iowa, 1968.
Bartlett, M. S. Proc. Roy. Soc.A, 160:268(1937).
Bartlett, M. S. Biometrics 3:39 (1947).
Bartlett, M. S. and Kendall, D. G. JRSS Suppl.8:128(1946).
Box, G. E. P. Ann.Math. Stat. 25:290 (1954).
Bozivich, H., Bancroft, T. A., and Hartley, H. 0. Ann. Math. Stat.
27:1017 (1956).
Bratcher, T. L., Moran, M. A., and Zimmer, W. J. J. Quality Tech.
2:156 (1970).
Burr, I. W. and Foster, L. A. A Test for Equality ofVariances,
Department of Statistics Mimeo Series No.282, Purdue University,
Lafayette, Indiana, 1972.
Carmer, S. G. and Swanson, M. R. JASA 68:66 (1973).
Cochran, W. G. Biometrics 3:22 (1947).
Dixon, W. J. and Massey, F. J., Jr. Introduction to Statistical
Analysis, 3rd ed., McGraw-Hill, New York, 1969.
Draper, N. R. and Smith, H. Applied Regression Analysis, Wiley,
New York, 1966.
Eisenhart, C. Biometrics 3:1 (1947).
1.5 REFERENCES 39
y. . = y + M. + T. + MT. . + e r . (2.1.1)
ij i 3 13 (13)
i = 1, 2, 3, 4 j = 1, 2, 3
where
y = overall mean
42 2. INTERMEDIATE DATA ANALYSIS
In this model the subscripts for MT and e are the same, thus
indicating that the effects due to the interaction and the error
effects will be confounded. The brackets on the subscripts of e
indicate there is only one observation per cell or only one
experimental unit (piece of metal) for each material-temperature
combination. We could write the error symbol indicating one
observation per cell but the e , . i s sufficient if understood.
The data for this part of the experiment are given in Table
2. 1. 1.
TABLE 2.1.1
Temperatures
Materials 1 2 3 Total
Using a = 0.05 the ANOVA for these data is given in Table 2.1.2.
2.1 ONE OBSERVATION PER CELL 43
TABLE 2.1.2
F(0.05)
Source df SS MS F
crit
= 354.67
= 654.67
SS MT and/or
a SS total - SS materials - SS temperatures
error
= 654.67 - 8.67 - 354.67
= 291.33
TABLE 2.1.3
Source EMS
F - MS temperatures
2,6 ” MS MT and/or error
2.1 ONE OBSERVATION PER CELL 45
TABLE 2.1.4
Temperature
:erial i 2 3 <y±. - y,
*i.
a [1 yii • y.J]2
aSS (for 1 df) i,j 13 1____________I _________
Nonadditivity ‘ [E (- . - _)2][I (- . -_,2,
i j J
73.01
TABLE 2.1.5
F P(0.0£
Source df SS MS EMS
crit
Total 11 654.67
It still turns out that all F-tests are not significant at the
0.05 level, including the interaction of MT, an indication that the
MT term could be deleted from Eq. (2.1.1).
2.1 ONE OBSERVATION PER CELL 47
TABLE 2.1.6
Temperatures
orthogonal polynomials
(coefficients)
Materials -1 0 1
a
SS of materials by linear temperatures:
The ANOVA for the data from Table 2.1.6 is given in Table 2.1.7.
TABLE 2.1.7
(Materials) by
(linear temp.) 3 207.00 69.00 a2 + <|>(TM) 2.45 9.28
Total 11 654.67
Utilizing the same data as was used in Sections 1.2.1 and 1.2.2
we use the orthogonal polynomials in Appendix 11 and refer to Section
1.4 as a reminder of the technique. Here, however, we must consider
the orthogonal polynomials for both variables and cross multiply them
to get the appropriate coefficients for each observation. The
procedure is shown in Table 2.1.8.
TABLE 2.1.8
Temperatures
Orthogonal
coefficients -1 0 1
= 202.5
TABLE 2.1.9
F(0.0S)
Source df SS MS EMS F
crit
Total 11 654.67
Temperatures
Materials 1 2 3
a 169 182 188
b 174 183 190
c 181 188 194
d 184 191 195
2.2 TWO-FACTOR EXPERIMENTS WITH MORE THAN ONE OBSERVATION PER CELL
Example 2.1; Fixed Model. Consider the case where both factors
are fixed and the model is y.M = p + a. + B. + aB.- + i = 1,
ijk i j lj (ijjk
...,a j = 1, ..., b k = 1, ..., n where
til
y . a variable to be analyzed from the k experimental unit
th th
associated with the i level of factor a and j level
of factor B
y = overall mean
I*l.
ou * effect ofthe i level of factor a (fixed)
a b n
F F R
Source
0
0
0 0
1 1 1
£(ij)k
where
a b n
F F R
Source
0 b n
a 0 n
0 0 n
1 1 1
E(ij)k
The product turns out to be 1. The EMS for this source is then
2 2
l«o » a . Next to determine the EMS for aB interaction, we only
consider those rows which contain at least the same subscripts as
a3..> in this case only z.,. ... and aB... Then cover all columns with
13 (i3)k iJ
headings the same as those contained in the aB row heading, i.e.,
columns i and j. For each rcw under consideration, determine the
product of the remaining coefficients. In this case it is the number
1 for the e r . ... row and n for the aB.• row. These values then
become the coefficient for the corresponding component. The sum of
these terms is then the EMS for aB.
a b n
F F R
Source i j k EMS
a. 0 b n
l
a 0 n
63 2
aB 0 0 n a +
1 1 1 a2
6(ij)k
We will use the notation given by Cochran (1951) for fixed components
as <p which is a non-negative function dependent upon the effect
contained within the parentheses. Hence, following the same rules
of allowing only those components that have at least the same subscripts
as the source being investigated to enter the EMS and using the same
column covering rule with multiplication of coefficients we finally
complete the EMS as:
Source i j k EMS
2
a. 0 b n a + bn<J>(a)
l 2
a 0 n a + an<J>(B)
63
0 0 n + n<|>(aB)
2
1 1 1 a
e(ij)k
56 2. INTERMEDIATE DATA ANALYSIS
Example 2.2: Random Model. For this case all factors are
random and we place R!s in all column headings. Using the same rules
as before we obtain the following:
a b n
R R R
Source i P k EMS
2 2
a. 1 b n +
a2 na 0 + bna
l a$ a
a 1 n a2 + naa$ + ana2
6i
a$. . 1 1 n a2 + na2
ij aB
2
1 1 1 a
2
Notice that since all effects are random the a notation
(rather than any <J> for fixed) is used for all components. This
indicates sampling variance and the experimenter's interest is in
estimation of these components of variance as well as testing the
2 2 2
hypotheses that the components, a^, a^ and a ^ are zero.
2
It can easily be seen that the test for 0 requires the e
or within error mean square as the denominator and both main effects
are tested using the interaction (aB^j) mean square. The latter is
quite different from the previous tests for fixed models (the means
are of interest) in which the within error mean square is used as
the denominator for all tests.
2
In this case <J> is used only for fixed and a is used for both
2
random and mixed (a^ is mixed here) components.
Hence all tests made within the inference space are based on the
within error.
I 2 .... j
Experim ent
Inferen ce-----
w here:a « A
Space and b « B
0
1 2 ... j b=B
1
2
i E xp erim en t
9 ------------ In fe re n c e ------------------------
Space
A» a
2.2 MORE THAN ONE OBSERVATION PER CELL 59
TABLE 2.2.1
4 3 1
F F R
df_____________ Source______________ i j k________ EMS
3 Ma 0 3 1 a2 + 3*(M)
2 T\ 4 0 1 a2 +
6 MT. . 0 0 1 a2 + <HMT)
0 i l l
We could use this result for the EMS, but since there are no
degrees of freedom for the within error, many authors prefer
dropping this term in the ANOVA. From Table 2.2.1 above it can be
seen that Table 2.1.3 combines the components of both MI\^ and
in the same source, MI\^. This consolidation of effects is necessary
since there are no degrees of freedom for error. Notice, however,
in Tables 2.1.5, 2.1.7 and 2.1.9 that there are separate estimates
of MT and e mean squares each of which has nonzero degrees of freedom.
Example 2.5: Fixed Model (More than one observation per cell).
If in the same example given in Section 2.1, there were two
experimental units per treatment combination, then the algorithm
method of obtaining the EMS would be that presented in Table 2.2.2.
Notice that the two main effects and interaction are tested
using the within error.
60 2. INTERMEDIATE DATA ANALYSIS
TABLE 2.2.2
4 3 2
F F R
df Source i j k EMS
0 3 2 a + 6<J>(M)
4 0 2 a2 + 8<j>(T)
0 0 2 a2 + 2<KMT)
2
■
1 1 a
—i
Example 2.6: Random Model. Using the same example but assuming
both factors are random we have Table 2.2.3.
TABLE 2.2.3
4 3 2
R R R
df Source i ) k EMS
3 M. 1 3 2 a2 +
l 2° m * “ m
2 4 1 2 a2 +
tj 2ow * *>?
a2 +
£
6 1 1 2
2°S r
•H
•r-»
12 1 1 1 a2
e(ij)k
TABLE 2.2.4
4 3 2
F R R
df Source_____ i j k___________ EMS______
3 0 3 2 a2 + 2a^, + 6<|>(M)
2 T\ 4 1 2 a2 + 8a2
6 0 1 2 a2 + 20 ^,
12 .(..)k 1 1 1 «2
/
/Mean square used to test the fixed treatment_______W1 2
INumber of observations in each mean of the treatment#
TABLE 2.2.5
1 2 3 Total
1 10,12
© 15,19 © 10, 8 © 74
i = 1, 2, .., 5 j = 1, 2, 3 k « 1, 2 where
th th
y. = coded value of strength of the k weld at the i
1*J time of welding and jth gage bar setting
y = overall mean
All assumptions for the ANOVA are assumed met, and the ANOVA
data are given in Table 2.2.6
64 2. • INTERMEDIATE DATA ANALYSIS
TABLE 2.2.6
5 3 2
F F R
df Source i j k EMS MS F
Within
2
15 error (e,. ...) 1 1 1 a 10.9
(ij)k
29 Total
96.4
^within
SS 9 9
within _ error = 1 2 2_ 2+ 2_ J S N j
error “ df = 15 ~ 15LJ‘w T
Jerror 9
2 2 J
= 10.9
TABLE 2.2.7
a /in Q ~
Standard errors: cell (body) = s- =J — j- = 2.33
Xij*
Overall = s- = = 0.60
X... y 30
Figure 2.2.1
A - 1, 2, ..., 15 (2 . 2 . 2)
m = 1, 2
2.2 MORE THAN ONE OBSERVATION PER CELL 67
where
Using Eq. (2.2.2) we can rank the cell means and use the
Newman-Keuls test as we did in Section 1.2 (See Fig. 2.2.2).
(i.J) (3^2) (3,1) (5^1) (1,2) (4,1) (2.1) (5.3) (4.3) (2.2) (4 .2 ) (5 .2 ) 0,1) (2 .3 ) (1,3) (3.3)
Mean 34 25.5 19 17 17 15 15 14.5 13 13 13 il 10.5 9 7.5
Figure 2.2.2
2.2 MORE THAN ONE OBSERVATION PER CELL 69
R6 = q0.05
(6,15) s- = (4.60) (2.33) = 10.7
yij-
AJ
(5,15) s- = (4.37) (2.33) = 10.2
RS = q0.05
yy.± •
r
R4 = q0.05
(4,15) s- = (4.08)(2.33) = 9.5
y- •
(3,15) s- (3.67) (2.33) = 8 . 6
R3 = q0.05 =
yir^j
(2,15) s- = (3.01)(2.33) = 7.0
R2 = q0.05 yn.
th
The basic model would now be only appropriate for the 5
level of time of welding:
y = y + G + e, x P - 1, 2, 3 (2,2.3)
pq p (p )q q - 1, 2
where
th
y = variable to analyze the q experimental unit in
the p gage bar setting
Means 19 13 15
(3) Differences:
Levels
Conclusion:
Time
i 2 3 4 5
Z1 -2 -1 0 1 2
2 -1 -2 -1 2
Z11
-1 2 0 -2 1
zm
y 34 26 68 26 26
SS ,
-K»l 1 . 12.8
[(-2r * (-1) + r + 2 ] • (2)
2 2 2 2 2 2
34 + 26 + 6 8 + 26 + 26 (180)
= ----------- -2-------------------- = 664
TABLE 2.2.8
Source df SS MS
Linear 1 12.8
Quadratic 1 : 165.1
Lack of fit 2 jI 486.1 243.0*
V
Cubic lll) 3.2
Lack of fit 1 1 482.9 482.9*
Conclusion:
By this we mean for the example on time and gage that the 8
degrees of freedom for interaction may be broken into the following
parts each with 1 df:
1. Time linear x gage linear
2. Time linear x gage quadratic
3. Time quadratic x gage linear
and so on to
8. Time quartic x gage quadratic
the within error mean square of Table 2.2.6 will be used for any
tests of hypothesis. Reference should also be made to Section 1.4
for methods of calculation.
TABLE 2.2.9
The ANOVA for the data of Table 2.2.9 is given in Table 2.2.10
for reference purposes in the derivation of the response surface
approximation.
TABLE 2.2.10
Source df SS MS
Within error
frdnuTable 2.2.6 15 164.0 10.9
TABLE 2.2.11
SS for 1 df T. x G
I q
G
Orthogonal q
coefficients 1 -2 1
-1 30 (-1) 26 ( 2) 21 (-1)
0 51 ( 0) 68 (0) 15 ( 0)
1 34 ( 1) 26 (-2) 29 ( 1)
= 6.0
TABLE 2.2.12
12.0 1.00
T* zi al
289.0 -2.83
Tq zn all
301.0
208.3 -4.17
ga Z2 a2
I CM
Nl
G 100.0 -1.67
CM
q a 22
308.3
2.0 0.50
t ji x G* Z1 Z2 a12
6.0 0.50
T* x G q Z1 Z22 a 122
140.2 2.42
Tq x G * Z11 Z2 a112
T x G 312.5 2.08
q q Z11 Z22 a 1122
460.7
2.2 MORE THAN ONE OBSERVATION PER CELL 77
y = ao ♦ a1Z1 ♦ Z2 ♦ an 2 2 Zn Z22
zn z 2 = x2 t(T ■ T) 2 ■ ^ n ] xi ( G " G)
= 3 t(T - 2)2 - |] (G - 2)
Temperatures
Supplements 40 80 120
2.3 REFERENCES
79
80 3. A SCIENTIFIC APPROACH TO EXPERIMENTATION
The change cost the company $200,000, and after two years in the
field there was as much trouble with the product made from the new
alloy as there had been with the old product. The vice president was
disgusted and called one of the authors of this book to say he would
never allow his company to use designed experiments again. After
some discussion, the vice president allowed the author to talk with
the experimenter to find out how the experiment was conducted.
TABLE 3.1.1
Source df EMS
Hence, rather than testing that the alloy effect was zero
[<f>(A) = 0], he was really testing that the total effect for ingots,
2 2
heats, and alloy was zero [300^ + 30cfh + 30<J>(A) = 0]. Since the
long-run production of the new alloy did not produce the improvement
seen in the experiment, <KA) must equal zero. Thus, it must have
2 2
been that a j. and/or a were/was not zero.
When the results were explained, the vice president was willing
to use design of experiments again on this type of problem but
insisted upon having more than one heat for each alloy.
3.1 RECOGNITION THAT A PROBLEM EXISTS 81
The authors feel that through committee action of this type one
is almost always able to find realistic levels for all major factors.
Occasionally considerable effort is needed to find out just how
shoddy one's manufacturing operation really is, and this example
shows that actual production operations will allow their process to
operate at various levels as long as it works. The determination
of the optimal levels is, of course, the desired end of the
experimental investigation.
84 3. A SCIENTIFIC APPROACH TO EXPERIMENTATION
this point, the statistician must write down the mathematical model
that has evolved as a result of the committee activity in the
preceding sections. This mathematical model will give rise to the
ANOVA table. The ANOVA table will now consist of the degrees of
freedom and the expected mean squares for each of the specific
factors selected in Section 3.3 above. These expected mean squares
will indicate to the statistician and experimenter all of the various
factors which will have tests available. In addition, if the test
is not directly available, it will be clear as to what assumptions
must be made in order to obtain some type of conservative test. Once
the statistician has reviewed this ANOVA table he must then present
his findings to the committee to see if these assumptions and
conclusions are realistic from a practical point of view. Should
the committee decide that some of the necessary assumptions (which
must be made in order to test the effects of various factors) are not
realistic, the design will have to be changed before proceeding
further in the experimentation process.
3.11 CONCLUSIONS
3.12 IMPLEMENTATION
3.13 SUMMARY
1. Inference Space
2. Randomization
3. Replication
3.14 REFERENCES
94
4.1 ONE FACTOR 95
Packing Pressure
7 14 21
and draw a random number between 1 and 21 again. Say the number is
12 , then.the number 12 explosive switch (2 0 , 0 0 0 psi pressure) is
fired first. Continue in this manner until all 21 switches have an
order number to be fired. All of this randomization must be
completed before the firing is begun so the experiment can be carried
off smoothly. One possible firing order is
2 2 12,000
3 11 20,000
4 20 28,000
5 15 28,000
6 3 12,000
7 9 20,000
8 17 28,000
21 8 2 0 ,0 0 0
4.1 ONE FACTOR 97
where
\i = overall mean
th
P^ = effect of the i packing pressure
TABLE 4.1.1
Problem 4.1.1. Show how you would examine the trend of firing
time over packing pressures.
98 4. COMPLETELY RANDOMIZED DESIGN
Draw a random number between 1 and 18, say 3. Then the first
combination or explosive switch to be fired will be number 3
100 4. COMPLETELY RANDOMIZED DESIGN
y i - U ♦ Ti * ; i = 1, 2, ..., 18 (4.2.1)
where
y^ = firing time of the switch with the i treatment
(combination) considered the variable to be analyzed
y = overall mean
TABLE 4.2.1
Source df EMS
Treatments 17 a2 + « (T)
Experimental Error 0 o2
where y. ., = the firing time required for the switch made up of the
th J ^ ^
i metal, j initiator and k pressure, as indicated in the
treatment combination table, and all components on the right-hand
side of the equation are the main effects and interactions of metals,
initiators, and pressures plus the mean and error.
TABLE 4.2.2
Source df EMS
2
Metals (M) 1 0 + 9<j>(M)
Initiators (I) 2 a 2 + 6 <J>(I)
M x I 2 a2 + 3<|>(MI)
Pressures (P) 2 a2 + 6 *(P)
M x P 2 a2 + 3<j>(MP)
Ix P 4 a2 + 2*(IP)
M x Ix P 4 a2 + <#>(MIP)
2
Error __0 a
Total 17
the MIP mean square the tests are conservative and those effects and
interactions that show significance should be significant.
2
If the experimenter had no prior knowledge of a and could not
assume <|>(MIP) = 0, he may consider repeating treatment combinations
on randomly chosen experimental units or repeating a few of the
treatment combinations in the experiment. The usual practice in
analyzing the experiment with only a few of the treatment combinations
repeated on other experimental units is to use the additional
observations to estimate the experimental error, but not to use these
additional observations in estimating the effects of the treatments.
If a program for a digital computer were available to use the
appropriate weights, it may be desirable to utilize all the
information; but if the investigator does not have the program and
is in a hurry for the results, he usually does not use these
additional observations in the analysis of effects.
where the components are the same as in Eq. (4.2.2) except that
I = 1 , 2 for only six of the treatment combinations, and I = 1 for
the remaining twelve treatment combinations. The analysis of
variance is the same as that given in Table 4.2.2 except that we use
only the first observation for the factorial analysis, but use the
two observations per treatment combination to calculate the error
sums of squares (SS). For each treatment combination with two
observations, calculate the SS with 1 df by squaring the difference
between the observations and dividing by 2 [see Example 7.2 for
proof]. Then merely add the six SS to obtain the SS for error. The
ANOVA is given in Table 4.2.3.
4.2 MORE THAN ONE FACTOR 103
TABLE 4.2.3
Source df EMS
2
Error 6 a
where
y = firing time (assumed normally distributed)
TABLE 4.2.4
ai
Metals (M) 1 a2 + <|>(M)
Linear initiators (1^) 1 *2 + *(1 *)
a2
Quadratic initiators (I ) 1 a2 + *(l )
a 22
2 "
ai2 M x I, 1 o + ♦(MI1)
Linear pressure (P^) 1 a2 + ♦(P1)
a3
Quadratic pressure (P ) 1 a2 + <t>CP )
a33
2 q
MxP, 1 + <KMPa)
ai3
a23 1 °2 * ♦(!,?,)
df
M x I
q
Residual (higher order M x P
q
terms from the P0 x I
interactions) _9 I q > a2
P x I0
17 q *
P X I
q q
MIP
4. 2 MORE THAN ONE FACTOR 105
)
and similarly for other means.
Hence, with the higher order terms in the model being zero,
there is an estimate of the error variance from the data with nine
degrees of freedom. If the assumptions are correct, the design of
this experiment is very "efficient,” that is, few observations have
provided the estimates desired for the model; therefore there is a
reasonably good error estimate.
degrees of freedom, the residual from regression (or the lack of fit)
could be tested using the estimate of the experimental error from
the six repeated treatment combinations as the error. It would be
advisable to use the probability of a Type I error of 0.25 before
pooling the residual error with the experimental error because under
these experimental conditions and reasonable alternative hypotheses
the probability of the Type II error, 3, is usually quite large for
small values of a. Pooling is accomplished by adding together the
SS of each source and the corresponding degrees of freedom.
regression model for a completely randomized design, and the data for
this two factor design with coded levels of temperature and calcium
content are given below. The measured response, y, is time-measured
in hours (cell totals are circled).
i = 1, 2, 3 j = 1, 2, 3, 4 k = 1, 2
Temperature
1 2 3 4 Total
26 ^
c
1
: © ; © © 27 © T."*• = 117
1
o
-M
e
o
u 2
; © : © :: © © T 02 *• = 113
G)
©21© W
u
3
: © : 23
;; © T.3* • = 143
Total T.i. = 43 T 0 = 51 T •4*
A =169 T>-# = 373
•2 * T*3* " 110
CT
L = 2x3x4
-*■*— = = 5797 04
24 S/S/.U4
Source df SS MS EMS F
Total 23 1885.96
1591.41 3.64
Z2 a2
T 108.38 2.12
q Z 22 a 22
T 21.67 -0.42
c Z 222 a 222
1721.46
Cl x Tl 33.80 0.65
Z1 Z 2 a 12
0.00 0.00
C* X T q Z1 Z 2 2 a 122
0.45 -0.08
C* X T C Z1 Z 2 2 2 a 1222
35.27 0.38
Cq X T * Z1 1 Z 2 a 112
C x T 0.00 0.00
q q Z1 1 Z 2 2 all2 2
C x T Z1 1 Z 2 2 2 0.15 -0 . 0 2
q c a 11222
69.67
y = 19.38 - 6 .8 8 C + 2.12C2
Source______df_____ SS_______ MS
Regression 3 1721.46 573.82
Residual 164.50 8.22
Total 23 1885.96
As seen in the above ANOVA table, the sum of square due to regression
is 1773, not 69.67. The reason for this result is that the cross
product sum of squares includes some of the main effects of both
calcium content and temperature as well as their interaction effect
for this nonorthogonal regression case.
The sum of squares due to regression is now the same as the sum of
squares due to the calcium content, temperature main effects and
the interaction combined.
For mixed models see Anderson and Bancroft (1952, Chapter 23),
Hicks (1973, Chapter 10), Winer (1971, Chapter 7), Edwards (1962,
Chapter 17) and John (1971, Section 4.7).
112 4. COMPLETELY RANDOMIZED DESIGN
TABLE 4.2.5
ANOVA
Source df EMS
2
Metals (M) 1 a +
9 °M
2
Initiators (I) 2 a + 3 ojjj * 64.(13
M x I 2 a2 +
" mi
Pressures (P) 2 a2 + 3 o9p • 6 *(P)
M x P 2 a2 +
3 4
Ix P 4 a2 +
°MIP * 2*<IP>
2
M x I x P 4 a2 +
°MIP
Error 0 a2
2
to take a partial replication to obtain an estimate of a directly
(not assuming the higher order terms are zero for the residual).
Metals (random)
1 2
1
44
39
© 27©
r \
20
35
30
© 13©
10
22
I S ©
r \
00
25
6©
©
Pressures ~ 12 /~\
(fixed) 1 I ? © 27 ©
O
21 3 2 © 11 ^ 22
3 43 /^\ 28
22
r \
31 ©
38 W
7 ©
12
“ ©
13
21©
19
i = 1,2 j = 1, 2, 3 k = 1, 2, 3 t = 1, 2
114 4. COMPLETELY RANDOMIZED DESIGN
Totals:
Metals Initiators Pressure Grand
T 78 T2.2. - 95
T1 2 .. * 143 22 *• t i .2. = 197
T = 126
T13*• = 197 23* • t i .3. = 203 T2.3. - 83
Initiators by Pressure
Z Z Z Z y2 . = 25196
i j k it
T2
CT = — = 20258.78
n
Sums of squares:
= 969.50
- C T - S S Metals - SS Pressure
= 14.00
- C T - S S Initiators - SS Pressures
= 20788 - CT - SS I - SS P
= 4.61
S S M x I x P = i s £ ST?., - CT - SS M - SS I - SS P - SS M x I
i j k 1Jk*
-SSMxP-SSIxP
= 24946 - 24907.50
= 38.50
SS Error = SS Total - SS M - SS I - SS P - SS M x I
- SS M x P - SS I x P - SS M x I x P
= 4937.22 - 4687.22
= 250.00
116 4. COMPLETELY RANDOMIZED DESIGN
Source df SS MS EMS F
2
Metals 1 3136.00 3136.00 a + 225.8*
18 4
2
Initiators 2 513.72 256.86 0 + 6 12 *(I) <1
’mi *
2
M x I 2 969.50 484.75 a + 6 34.9*
•2 ,
2
Pressure 10.89 5.44 a2 +
QMP + 12*(P)
2 6 <1
2
M x P 2 14.00 7.00 a2 + 6 <1
aMP
I x P 4 4.61 1.15 a2 +
aMIP + 4<j>(IP)
2 <1
Total 35 4937.22
TABLE 4.2.6
ANOVA
Source df EMS
2 2 2 2
Metals (M) 1 a2 + + 3 3 + 9
°MIP aMI + aMP aM
Initiators (I) 2 a2 +
2
°MIP
+ 2
«?p *
3
■s, + 6
2
al
M x I 2 a2 + a2
MIP
+ 3
•2,
Pressures (P) 2 a2 + 2
MIP
2
+ 2
2
3
2
aMP
+ 6
<•1
M x P 2 a2 + + 3
°MIP MP
2
Ix P 4 a2 + + 2
aMIP alP
2
M x I x P 4 a2 +
°MIP
2
Error __0 a
Total 17
The two-factor interactions are all tested against the three factor
interaction, and all main effects require approximate F tests which
are sometimes called F 1 tests, Cochran (1951), Anderson and Bancroft
(1952, p. 350), Ostle (1963, p. 302). For example, to test
2
H: aw = 0 use
M
r, _ MS(metals) + MS(M x I x P)
“ MS(M x I) + MS(M x P)
2
since the EMS (metals) + EMS(M x I x P) contains only the term 9 aXi
M
in addition to the terms in EMS(M x I) + EMS(M x P).
where
k
M = 7 a.M.
i=i 11
th
M^ = i mean square in the numerator if the degrees of freedom
are being calculated for the numerator or similarly if the
denominator degrees of freedom are being calculated
Of all the models fixed, mixed, and random; the random model
allows a design with the fewest assumptions after the factors are
declared random. Thus it is absolutely necessary for the experimenter
to know whether each factor will be random or fixed before an
appropriate design is prescribed for factorial experiments.
ANOVA
Source df MS
Metals 1 530
Initiators 2 225
M x I 2 46
Pressures 2 35
M x P 2 28
I x P 4 10
M x Ix P 4 5
2 2 2 2
(a) Test the various hypotheses, = aj = = ap =
aMP = aIP = using the F f test and pooling
techniques where appropriate.
(b) Explain the results.
(c) Comment on a possible new model.
2 a2 - 2
maximum information because a, = — z and I(X.-X) is
1 I(Xi-X) 1
The experimenter may select doses of 100, 200, 400, and 800
roentgens for a fixed time and examine the trend by testing a
departure from the exponential. Of course it can be seen that if
the logarithm of dosages are taken and the departure from a straight
line is then investigated, the experimenter has really examined the
departure from the exponential. If the trend is not found to depart
significantly from the exponential, the coefficients a and b can be
estimated using the transformed X variable.
4.3 REFERENCES
123
124 5. RANDOMIZED COMPLETE BLOCK DESIGN
where
y ^ = the response from the experimental unit treated with the
i**1 level of treatment B and the j**1 level of treatment T
Vi = overall mean
An estimate of the error mean square may come from the interaction
source. Assume all the analysis of variance assumptions hold, so
2
that the is NID (0, a ), i.e., normally and independently
distributed with mean zero and variance a^. For this case the
appropriate analysis of variance (ANOVA) is given in Table 5.1.1.
TABLE 5.1.1
ANOVA for the Data from the CRD Using Equation (5.1.1)
Source df EMS
y. .. n = y + B. + + T. +
'ljk l l (i)k j (ijk)Z
(5.1.2)
i = l , . . . ,b j = 1 ,..., t k = 1 1 = 1
where
t»h
y - i v o = t*ie response from the % experimental unit in block
l j KXr th th
i, the k randomization, and given the j treatment
y = overall mean
th
B. = effect of the i block. In this case B. may be
i l 7
random or fixed, but we will assume it is fixed for
the ANOVA
128 5. RANDOMIZED COMPLETE BLOCK DESIGN
The analysis of variance for this model, along with the EMS
derivation is given in Table 5.1.2.
The corresponding model for Eq. (5.1.2) not using the subscripts
k and I is:
5.1 BLOCKS FIXED 129
TABLE 5.1.2
b 1 t 1
F R F R
df Source i k j Z EMS
b - 1 B. 0 1 t 1 a2 + ta2 + t<KB)
l
0 6 f ^i 1 1 t 1 a2 + ta2
(i)k
t - 1 T. b 1 0 1 a2 + bKT)
3
2
(b - 1 )(t - 1 ) 1 1 1 1 a
e(ijk)£
i = 1 , 2, .
,.,. + T.
, 2 , .;;; bt
y. . = U + B. + 6 cs.i.3)
0) J + £(ij) j = 1
l e r e the terms are defined as they were for Eq. (5.,1 .2 ) except that
k and Z are deleted.
It can be seen from Table 5.1.3 that the correct error for
blocks is the restriction error with zero df not the within error
with (b - 1)(t - 1) df as suggested in Table 5.1.1.
TABLE 5.1.3
Source df EMS
*
For example, Ostle (1963, p. 368).
5.1 BLOCKS FIXED 131
and is estimable.
The purpose of the experiment was to select the best valve type
out of the four for all pulse rates and/or best valve type for
particular pulse rates. To a statistician this means that the main
effect of valve types and the interaction of valve types by pulse
rates should be examined carefully.
132 5. RANDOMIZED COMPLETE BLOCK DESIGN
48. 1 6
TABLE 5.1.4
Pulse rates
1 3 2 1 4 2 3 5 4 6 6 5
I
Run c i) ( 2 ) ( 3) ( 4) ( 5) ( 6 ) ( 7) ( 8 ) ( 9) (1 0 ) (ID (1 2 )
5 2 3 4 2 6 1 4 5 3 6 1
(/> 4
<D
Run (13) (14) (IS) (16) (17) (18) (19) (2 0 ) (2 1 ) (2 2 ) (23) (24)
4&->
<D
>*4
r- 4 2 1 6 4 3 5 2 1 6 5 3
a
> 1
Run (25) (26) (27) (28) (29) (30) (31) (32) (33) (34) (35) (36)
6 1 5 4 1 2 3 3 5 6 2 4
3
Run (37) (38) (39) (40) (41) (42) (43) (44) (45) (46) (47) (48)
i = 1, 2, 3, 4 j = 1, 2, ..., 6 k = 1, 2
134 5. RANDOMIZED COMPLETE BLOCK DESIGN
where
th th
y . = k maximum flow gradient from the j pulse rate of
1 JK th
the i valve type
y = overall mean
TABLE 5.1.5
Source df EMS
Total 47
1
;® 7 © j© ,i© s® 104
2
3 ® 15® ‘?© :® 3 © a© 170
3
3 ® 3 ® ;® ;® s ® s® 174
4
!S@ ;s© !® 5© IS© n® 141
Source df SS MS EMS F
Total 47 1709.48
= 6.54
Thus the standard error to use for the multiple range test is
^6754
= /0782 = 0.91
8
T• •• = 589 CT = - = 7227.52
S.1 BLOCKS FIXED 137
- CT
= 7489.42 - CT = 261.90
= 8419.87 - CT = 1192.35
SS VP s i £ T?. - CT - SS V - SS P
2 lj*
- SS V - SSP
= 55.73
SS error = SS total - SS V - SS P - SS VP
= 199.50
5.1.3 RCBD When Interaction Is Present and One Observation Per Cell
i = 1 , 2 , ..., b j = 1 , 2 , ..., t
TABLE 5.1.6
Source df EMS
Total (bt - 1 )
TABLE 5.2.1
Source df EMS
Treatment (T) (t - 1 ) a2 +
°BT + b<KT)
2 2
B x T (b - 1 ) (t - 1 ) a +
aBT
Within error 0 a2
Total (bt - 1 )
2
where: cr^ is the ’’mixed” interaction component and cannot be
estimated with only one experimental unit per treatment, unless one
is willing to use a test for nonadditivity.
Let us now extend the mixed model case discussed in Example 2.7
to factorials. Consider the model
y. „ = y + B. + 6 ...
'ljkfc l (l)
+ p* + BPu + MP
where i = 1, 2, 3; j = 1, 2; k = 1, 2, 3, and Jl = 1, 2, 3
df: 3 4 = 2 + 4 + 4 + 4 + 4 + 8 + 8
y. . = y + G. + + H. + GH. . + z f . .v (5.2.3)
'ij i (i) j ij (ij) v J
i = 1, 2, 3, 4 j = 1, 2, 8
TABLE 5.2.2
Source df EMS
2
Blocks (genome sample) 3 a + 8 a2 ♦
8aG
Restriction error 0 a2 + 8 a2
a2 +
+
4*(H)
Total 31
142 5. RANDOMIZED COMPLETE BLOCK DESIGN
Example 5.2: Consider the case of Eq. (5.2.3) when blocks are
genome samples and treatments are heterozygous states. The coded
value of Drosophila body weight is as follows:
1 16 13 8 7 15 3 14 7 83
2 15 15 10 6 13 6 15 10 90
3 20 18 14 10 22 9 20 14 127
4 12 14 9 5 13 1 12 5 71
Totals 63 60 41 28 63 19 61 36 371
Source df SS MS EMS F
2
Genome sample 3 218.59 72.86 a + ♦ 8 a2 None
Total 31 797.72
s- = A . 67/4 = C 4 2 .
4
4
SS blocks = 1 V T? - CT = 1 (832 + 902 + 1272 + 712)
8 .I
4-1 i* 8
i=l
- CT = 218.59
g
SS treatments = -j I T2^ - CT = i (632 + ... + 362)
- CT = 543.97
For all RCBD there may be more than one observation for each
treatment in each block if there is enough time (and/or money) and
a reason from an efficient design point of view to handle the extra
observations. An extra observation may just be another measurement
of the same experimental unit, possibly by another technician. This
type of observation may not have much value for the inference because
the variation measured may be only an operator and measurement error
variation. Under some circumstances it is of value to have this
144 5. RANDOMIZED COMPLETE BLOCK DESIGN
y + B. +5
yijkJl l
where
i = 1 , 2 , ..., b j = 1 , 2 , ..., t k = 1 £ = 1 , 2
th
y. 0 = cellulose content of the JI observation of the
ij t!h tlr
plants in the k pot from the j treatment in the
i location in the greenhouse
1.L
B^ = effect of the i location
pr.... = random error due to the pot in the i**1 block and
.th
j treatment
th th
nr* •V'w = random error due to the I observation in the k
th th
pot in the j treatment in the i location.
TABLE 5.3.1
Source df EMS
2
= variance component between samples within experimental
units (pots) (a type of sampling error)
2
= variance component between experimental units (pots)
treated alike
2
cjg = variance component due to extraneous variables on the
blocks
146 5. RANDOMIZED COMPLETE BLOCK DESIGN
and the components are defined the same as for Eq. (5.3.1) except
that k = 1 , 2 here and I = 1 only here.
TABLE 5.3.2
Source df EMS
2
Blocks (B) (b - 1 ) a + a2 + 2ta2 + 2ti(>(B)
n P
2 2
Restriction error 0 a + a2 + 2 ta^
n P 6
Total 2 bt - 1
where
the variable to be analyzed from the setup of
yijk
the j method in the i laboratory
u = overall mean
TABLE 5.3.3
Source df EMS
2
Laboratories (L) 19 0 + 50o2 + 50a2
o L
Restriction error 0 a2 + 50a2
2
Methods (M) 4 a + 10a + 200<)>(M)
L x M 76 a2 +
Total 999
TABLE 5.3.4
Source df MS EMS
Laboratories (L) 4 70 a2 + 12 a2 ♦
12aL
Restriction error 0 a2 + 12 o2
L x M 8 40 a2 +
Within error 45 25 a2
Total 59
F = = 3 25
2,8 40
(n1 + l)(n2 + 3) s2
(n2 + l)fn1 + 3) s2
5.5 REFERENCES
154
6.1 NESTED (HIERARCHICAL) 155
where
y = the overall mean
Notice that the example shows the impossibility of the level of any
factor such as countiesbeing the samefor the levels of the factor,
state. Since there are never the same counties within the various
states in this investigation, it is impossible to obtain an
interaction between counties and states. This concept of never being
able to obtain an interaction is always present in a nested design.
TABLE 6.1.1
Source df EMS
, 2
°l + hoT + htoc
2
States s - 1 itcas
, 2
Counties in states s(c - i) haT +
•2 * htac
Towns in counties in states sc(t - i)
y. = u + S. + T r.. . + H r....
ijk H l (i);j (ij)k
i = 1, 2, 3 j = 1, 2 k = 1, 2, 3, 4
where the terms in the model are as defined in Eq. (6.1.1). Two
methods of calculation are presented: one gives the sum of squares
for nested factors directly, whereas the other utilizes the method
of calculation for a factorial design and then combines appropriate
sums of squares to obtain the correct values. The data for this
example are given in the following tabulation.
State
1 2 3
1 2 1 2 1 2
Household
1 10 7 6 6 15 12
2 13 12 5 12 18 15
3 16 11 9 7 20 18
4 12 9 3 10 19 16
Totals 51 39 23 35 72 61
T. = 90 T0 = 58 T-
1 *• 2 *• 3* • = 133
T• • • = 281
6.1 NESTED (HIERARCHICAL) 157
Source df SS MS EMS
2 2
States 2 354.08 177.04
4 + 4 + *4
2
Towns/states 3 51.13 17.04
4 + °T
Households/towns/states 18 107.75 5.99
•2
Total 23 512.96
where
2
i
1 3
n 0
/ T•••
SS states = j ) T. - «--- =-
2x4 i** 2 x 4 x 3
1=1
1 3 2 , 2
SS towns/states = j I I T. - £ T...
i=l j=l J i
= (i T2 . At2 )
>.*■4 J=1
1=1 i ij-
J
8 i-*J
1 2 2 2 1 2 2
= j (51Z + 39 + ... + 6 1Z) - j (90Z + 58
+ 1332)
51.13
2 1 2
SS households/towns/states = Z EZ y...- 7 Z Z T,.
i .k ^ j k 4 . . ir
1 2 2 2
- ~ (51 + 39 + ... + 61 )
= 3803 - 3695.25
= 107.75
158 6. NESTED AND NESTED FACTORIAL DESIGNS
2 •••
SS total * E n y...--
i j k 1JK
= SS states + SS towns/states + SS households/
towns/states
= 3803 - 3290.04
= 512.96
y. ., = \i + S. + T. + ST. . +
ijk l j ij (ij)k
T2
SS "towns" = [(51 + 23 + 72)2 + (39 + 35 + 61)2] -
= 46.09
6.1 NESTED (HIERARCHICAL) 159
Note that 46.09 + 5.04 =» 51.13, which is the same as that obtained
for the towns within states sum of squares. Calculation of sums of
squares for any nested type of sums of squares in the nested
factorial models may be done in a similar fashion.
' 5'99
i = 1, 2, 3, 4 j = 1, 2, 3, 4 k = 1, 2, 3 I = 1, 2 m = 1,
One proposal was to make four ingots, one with each chemistry,
and compare the strengths of the finished product. If this were
done the inference could be made to this one ingot for each
chemistry made from only one heat (batch of molten aluminum). If
each ingot could possibly represent an adequate random sample of
ingots for each heat, the inference could be made to all ingots from
that one heat with that chemistry. To infer indefinitely over heats
one must replicate the experiment by having at least one more heat
run per chemistry or assume the variation between heats was no
greater than the variation between ingots within heats. This latter
assumption, however, is usually very poor for problems of this type.
The cost of heats was quite large and the research worker
decided to take three ingots per chemistry chosen at random from
each heat-chemistry combination. Further he took four random
samples of the material in each ingot to measure the variation
6.1 NESTED (HIERARCHICAL) 161
between pieces within ingots. The appropriate model and ANOVA for
this experiment is:
i = 1, 2, 3, 4 j = 1 k = 1, 2, 3 I = 1, 2, 3, 4
where
y. ., 0 = strength of the ft**1 piece of metal in the k ^ ingot
13 th th
from the j heat with the i chemistry
y = overall mean
TABLE 6.1.2
Source df EMS
Problem 6.1.3. Using the EMS results from Table 6.1.2 and the
information in the following tabulation
162 6. NESTED AND NESTED FACTORIAL DESIGNS
Source df MS
Chemistries 3 160
Ingots 8 20
Pieces 36 10
Valve types
1 2 3 4
4 5 2 7 6 8 1 3
TABLE 6.2.1
Valve types
1 2 3 4
1 2 3 4 5 6 7 8
Pulse rates
1
3
4
5
6
The model for the analysis of the data from this experiment,
assuming the errors are not correlated, is:
where
\i = overall mean
TABLE 6.2.2
Source df EMS
2
Valve types (V^) 3 a2 + 6 a2 + 6 a2 + 12$(V)
Occurrences in valve 9 9 9
types (O^jj) 4 a2 +
<*« + 6%
Restriction error (6 ^ ^ ) 0 a2 +
6»s
MS(Pk)
(1) Pulse rates using F = j
(i)3 kJ
MS (VP..)
and (2) (VP.k) using F3 = MS(OP
(i)jk)
(see Morrison, 1967, pp. 194-196; Winer, 1971, pp. 522-524; and
Greenhouse and Geisser, 1959).
6.2 NESTED FACTORIAL 167
(2) Use the results of the direct ANOVA tests if they show no
significance
Valve types
1 2 3 4
6 6 6 6 7 7 8 8 9 57
Totals 27 32 31 24 40 43 45 44 286
59 55 83 89
i = 1 , 2 , 3, 4 j = 1 , 2 k = 1 , 2 , •••, 6
ANOVA
SS MS p(0.05)
Source df EMS F
crit
where
2 2 2 2 2
no ir _ 59Z + 55 + 83Z + 89 286
i " 12 “ 48
= 1776.33 - 1704.08
= 72.25
2 2 2 2
oo ^ _ 34 + 35 + ... + 57z 286
Fk “ 8 “ 48
= 1809.50 - 1704.08
= 105.42
2 2 2 2 2 2 2
„ _ 21l + 32* + ... + 44* 59 + 55* + 83* + 89*
bb U (i)j 6 “ 12
= 1783.33 - 1776.33
= 7.00
SS vp.k - 52 • 62 ~ ....... 1 7 2 - - SS V. - SS p k
= 38.25
where
2 2 2
SS OP(V1) = 2 + 3 + ... + 6 -
2
2 2 2 55
SS OP(V2) = 4* + 2* + ... + 7* -
SS OP(V3) = 62 + 52 + ... + 82 -
SS OP(V4) = 72 + S2 + ... + 92 -
= 4.41
Therefore
= 237.92
y.
'ljk
= y + M. + Sr.v. +
H i (i)j (1 3 )
+ T. + MT., + ST,.*., +
k lk (1 )jk
.n
(ijk)
(6.2.2)
^ '
i = 1, 2, 3 j = 1, 2, 10 k = 1, 2, 3, 4
172 6. NESTED AND NESTED FACTORIAL DESIGNS
where
th
variable to be analyzed, time to perform the k task
yijk
th th
by the j student in the i training method
y = overall mean
within error
e(ijk)
Note that the model and analysis for this example is the same as that
given in Section 6.2.1 (also see Winer, 1971).
The model, assuming the correlations among all errors are zero,
(Mandel, 1957), for the gain variable to be analyzed follows:
6.2 NESTED FACTORIAL 173
i = 1, 2, 3 j = 1 , 2, ..., 8 k = 1, 2, ..., 5
where
y. = gain in weight of the j***1 hen on the -ration
th
during the k week
y = overall mean
TABLE 6.2.3
Source df EMS
2
Rations (R) 2 a2 + So2 + 5a„ + 40<j)(R)
n
2
Hens in Rations (H) 21 a2 + 5a'
♦ *>2
c 2
Restriction error 0 a2 + 6
2 2
Weeks (W) 4 a + °HW + 24<t>(W)
2 2
R x W 8 a + + 81)1 (RW)
°HW
2 2
H x W 84 a +
°HW
2
Within error 0 a
174 6. NESTED AND NESTED FACTORIAL DESIGNS
s-
y*.
i = 1, 2, 3 k = 1, 2, ..., S
y =overall mean
TABLE 6.2.4
Source df EMS
2
Weeks (W) 4 a + 3*00
R x W 8 a2 + ♦ (MO
2
Within error 0 a
The research workers for this problem should have been aware of
the difficulty in the analysis of lfeggs laid” and have designed the
experiment differently. One major contribution at the design stage
to a good analysis would have been to 11trap nest” the hens and
record individual henfs performance. The analysis, then, would have
been similar to the analysis for the weight varaible, if the
assumption of no correlated errors holds.
Discs (thickness)
Powder lots
1 2
Container Container
A B C D
0.2 0.4 0.8 0.2 0.4 0.8 0.2 0.4 0.8 0.2 0.4 0.8
1 6
2 7
3 8
4 9
5 10
Discs (thickness)
Powder lots
1 2
Container Container
A B C D
0.2 0.4 0.8 0.2 0.4 0.8 0.2 0.4 0.8 0.2 0.4 0.8
116
119
120
i = 1, 2 j = 1, 2 k = 1, 2 Z = 1, 2, 3 m = 1, 2, 5
where
^ijk£m = var:*-a^le t0 analyzed if each igniter performed
properly and the length of time is used
y = overall mean
x.L
D^ = i thickness of disc effect
th th
C(j)k = e^ ect ^
container in the j powder lot
th
DC.,.,*, = effect of the interaction of the i thickness of
1(;,)k th th
disc with the k container in the j powder lot
and so on
th th
= effect of the m igniter in the ijkfc cell assumed
"(ijkJl)m
NID (0, a 2)
Powder lots
Container Container
A B C D
Moisture ^ ^
1, 2, 3, 4, 5
content
1 0.4
0.8
Discs
0.2
2 0.4
0.8 116, ..., 120
+ e(ijkil)m
COMMENT
It has been our experience that to understand the restriction
on randomization (if one exists) in nested factorial experiments, the
experimenter should place the factor that has the nesting at the top
of the design layout (horizontally). This allows one to see the
repeated measures (if they exist) over the other factor, down the
layout (vertically). As a result of all the above one can write
180 6. NESTED AND NESTED FACTORIAL DESIGNS
down the model in a manner that provides an ANOVA (with the EMS)
that is easy to interpret.
6.3 REFERENCES
181
182 7. SPLIT PLOT TYPE DESIGNS
A B C
1. 1 1 1
2. 1 1 2
3. 1 1 3
60. 3 4 5
+ B C + ABC. ., + z f . .|x
}k ljk (ljk)
i = 1, 2, 3 j = 1, 2, 3, 4 k = 1, 2, 3, 4, 5
TABLE 7.1.1
Source df EMS
A 2 a 2 + 20<KA)
B 3 a2 + 154i(B)
AB 6 a2 + 54(AB)
C 4 a 2 + 124(C)
AC 8 a2 + 4<j>(AC)
BC 12 a2 + 34(BC)
ABC 24 a2 + 4 (ABC)
Error 0 a2
Total 59
where
th th
= variable to be analyzed from the i level of A, j
level of B and k**1 level of C
y = overall mean
The tests for all main effects and interactions are not
straightforward in Table 7.1.2. Some assumptions on interactions
are needed or conservative tests must be made when this fixed type
design is used. If, however, it is known by the experimenter that
all interactions with A in this experiment are zero, the tests are
all obvious.
TABLE 7.1.2
Source df EMS
First restriction o 9
error 0 a2 + 5a + 20a*
0) 0
B 3 a2 + Sa2 + 15*(B)
AB 6 a2 + 5a2 + 5* (AB)
Second restriction 9
error 0 a2 + 5a2
w
C 4 a2 + 12*(C)
AC 8 a2 + 4* (AC)
BC 12 a2 + 3<f>(BC)
ABC 24 a2 + *(ABC)
2
Within error _0 a
Total 59
TABLE 7.1.3
Source df MS
Blocks 1
Valve types
Whole plot error
Pulse rates
Valve types x pulse rates
Split plot error
Harter (1961), Chew (1958, p. 48), and others have doubts that
one should always obtain the split plot error by pooling. In
biological studies using the "sometimes pooling" technique described
by Bozivich et al. (1956) and Bancroft (1968, p. 11), there seems to
be good reason to pool "almost always." Some authors have used the
sometimes pooling criterion when the split plot treatments are
random.
(7.1.3)
i = 1, 2 j * 1, 2, 3, 4 k = 1, 2, ..., 6
7.1 SPLIT PLOT 187
where
U ® overall mean
TABLE 7.1.4
Source df EMS
2 2 2 2
Blocks (Bi) a + 6cT + 24a* + 24a^
a) 6 B
First restriction 2 2 2
error C5^ ) (T + 6<T + 24a*
(0 0
Pulse rates 5 a2 +
" b p * 8<KP)
Interaction (BPik) 5 a2 +
<
Interaction a2 +
<vv
15 °BVP + 2+(VP)
2
Interaction (Bw ijk) 15 a2 +
BVP
Within error 0 a2
(e(ijk)}
Example 7.1: Using the actual data from Beeson (1965), the
analysis turned out to be as that given in Table 7.1.5. When the
correct tests are made, only one source, pulse rates, is significant
at the a = 0.05 level. However, the mean squares of the
interactions with blocks at the various stages are about the same
size, indicating that the actual variance components with blocks
are possibly zero [from Eq. (7.1.3)]. Using a certain criterion
for pooling, the more liberal ’’sometime poolers” may even obtain
the analysis of Table 7.1.6. This shows that all effects and the
interaction are significant; however, this last analysis is valid
only if the errors are poolable and pooling is correct. In general,
for a correct analysis, all restrictions on randomization must be
considered before pooling procedures are undertaken.
190 7. SPLIT PLOT TYPE DESIGNS
TABLE 7.1.5
Source df MS
♦Significance at a = 0.05.
TABLE 7.1.6
Source df MS
♦♦Significance at a = 0.01.
♦Significance at a = 0.05.
7.1 SPLIT PLOT 191
TABLE 7.1.7
Source df EMS
It was shown that temperatures had a major effect and the best
temperature from this analysis was used in production. Some time
later in the plant production, it was shown that the new temperature
used in this experiment had no significant effect over the old one
on strength of the alloys and the experimenter could not understand
why his experiment had given him the wrong results.
192 7. SPLIT PLOT TYPE DESIGNS
where
i = 1, 2, 3, 4 j = 1, 2, 3 k = 1, 2 £ = 1, 2
and all effects and interactions are defined in the usual manner.
The corresponding analysis of variance is given in Table 7.1.8.
TABLE 7.1.8
Source df EMS
Alloys 2 a2 + 164(A)
T x A 6 a2 ♦ 44 (TA)
Orientation 1 a2 + 244(0)
2 +
T x 0 3 a 64(TO)
A x 0 2 a2 + 84(AO)
T x A x 0 6 a2 + 24 (TAO)
Total 47
7.1 SPLIT PLOT 193
Temperature
1 2 3 1 2 3 1 2 3 1 2 3
Orientation 2
15 20 27 35 48 47 55 62 23 33 25
i
i
19 28 26 40 39 55 55 63 58 34 39 17
15 25 25 48 48 55 48 63 68 28 33 14
2
23 31 33 55 62 64 60 68 62 37 38 12
Source df SS MS EMS F
Total 47 14508.67
S H I y* = 88142
i j k 4 1JJa
- CT
= 84,499.16 - CT
= 10865.83
= 74067.62 - CT
= 434.29
= 74,102.08 - CT
= 468.75
Now,
S S T x A -ilS T ? . -CT-SSA-SST
4 . . li••
i 3 J
- SS A - SS T
= 86,174 - CT - SS A - SS T
* 1240.55
SS T x 0 4 Z I T? , - CT - SS T - SSO
6 ., i*k*
ik
* 85368.33 - CT - SS T - SS 0
* 400.42
s s a x o 4 z z t 2 ., - c t - s s a - sso
» j k ' 3K*
* 74,611.00 - CT - SS A - SS 0
= 74.63
S S T x A x O = i z Z Z T?..- CT - SS T - SS A - SS 0 - SS TxA
2 . . . ilk*
1 3 k J
-SSTxO-SSAxO
1 2 2 2 2
=j (21z ♦ 38^ + 43^ ♦ ...♦ 26 ) - CT
- SS T ..... S S A x O
« 87,193 - CT - SS T .......SS A x 0
* 75.20
= 88,142.00 - 87193.00
= 949.00
Since there areonly two observations per cell for thisproblem, the
error sum of squares may be calculated by summing the squareof the
differences of the two observations in each cell and dividing by
two, i.e.,
=j (1898.00)
= 949.00
- =-
y yl
- + y2 9 the cell mean. Note that
- 2 - 2 yl + y2 2 yl + y2 2
(yx - y) + (y2 - /) s Cxx - 3 + (y2 " " S " 3
2yr yr y2 2 2y2_yr y2 2
= +( 2 )
( y r y 2) 2 (y2- y p 2
4
cy r y 2)2
2
Summing over all cells gives the sum of squares for error for all
cells.
198 7. SPLIT PLOT TYPE DESIGNS
+ HWcr.. .. „ + e r .
(i)jk£ (ijk£)m
i = 1, 2, 3 j = 1, 2 k = 1, 2, 3 Z = 1, 2 m = 1, 2, 3
where
^ijkfcm = var^-a^^e t0 analyzed
y s overall mean
TABLE 7.2.1
Source df EMS
Whole Plot
Chemistry (C) 2 a2 + 6a2 + 18a2 +18a2 + 36<|>(TV
C)
a) <5 H J
Split Plot
Wheel heights (W) 2 a2 + 6a2 + 6a ^ + 36<KW)
Split-Split Plot
2
Conditions (c) 1 a + 9a2c + 54*(c)
2
C x c 2 a + 9a2c + 18* (Cc)
H x c 3 a2 +
2
W x c 2 a + ' i * »«♦<*>
C x W x c 4 a2 +
3<THWc + 6^ CWc^
H x W x c 6 a2 +
3«mc
Within error 72 a2
Total 107
202 7. SPLIT PLOT TYPE DESIGNS
Columns
1 ooooo
2 ooooo
Rows 3 ooooo
4
ooooo
5 ooooo
The active ingredients are most important to the
experimenter and the effects of rows and columns in describing
the response surface of the sheets are less important. Each
sheet was selected at random for a run through the equipment
and was made up of one combination of three active ingredients
(A at 4 levels, B at 3 levels and C at 2 levels). This made a
total of 24 sheets ( 4 x 3 x 2 ) for one experiment.
(a) Show the design of the experiment, emphasizing the
randomization and inference space.
(b) Show the model and analysis of variance for your
experiment.
204 7. SPLIT PLOT TYPE DESIGNS
h — ii
Cut
Source MS
(a) Mix types: (1) sand with cement; (2) 3/8 inch
rock with cement; (3) 3/4 inchrock with cement.
(b) Percent volumeof fibers: 0, 1, 1.5, 2, 2.5.
(c) Size of cylinders: 3” x 6”, 4”x 8", 6” x 12”.
The experiment was conducted in the following manner:
Day 1: One batch of mix (sand with cement) was made up and
split into 5 parts.
The first part had 0% fibers added,
the second part had 1% fibers added,
the third part had 1.5% fibers added,
the fourth part had 2% fibers added,
and the fifth part had 2.5% fibers added.
Each part was made into 6 cylinders, two each for the
three sizes: 3" x 6”, 4" x 8", 6” x 12”.
Day 2: Repeat of Day 1.
Day 3: The batch of mixwas made of 3/8” rock with cement and
the same procedure was followed for the other two
factors as in Days 1 and 2.
Day 4: Repeat of Day 3.
Day 5: The batch of mixwas made of 3/4” rock with cement and
the same procedure as above was followed for the other
two factors.
Day 6: Repeat of Day 5.
Source df
Mixes (M) 2
Days (D) 1
M x D 2
Fibers (F) 4
MF 8
DF 4
MDF 8
S 2
MS 4
DS 2
MDS 4
FS 8
MFS 16
DFS 8
MDFS 16
within 90
(f) If it turned out that the main effects were all significant
and no interactions were significant show what cautions you
would make in the interpretations and what further analyses you
would do on each.
(g) What further analyses would be appropriate if the
interaction of mixes x fibers turned out to be significant?
Consider regression analyses.
(h) Consider the whole experimental procedure used and
recommend a better design explaining at each stage why you
think your design is better. (Keep the cost of running the
experiment in mind, so you do not run a completely randomized
design each day.)
7.4 REFERENCES
210
8.1 LATIN SQUARE 211
Columns (Blocks)
Fertility Trend ------------
1 2
Rows (Blocks)
y. ., = y + R. + 6... + C. + r),.^ + T, + £r .n
•'ljk i (i) j Q) k (ljk)
j = k = 1, 2, ..., t
where
y.., = variable to be analyzed from the i row j column
1*J th
and k treatment
y = overall mean
.th
R. = effect of the i row
l
6... = row error, due to restriction on randomization
(1) (Anderson, 1970) of treatments in the i
th row
•f”L
Cj = effect of the j column
TABLE 8.1.1
Source df EMS
Experimental error
2
(residual) (t - 1) (t - 2) a
Total t2 - 1
8.1 LATIN SQUARE 213
TABLE 8.1.2
Men (m)
Machines
1 2 3 4 5
(M)
where
The ANOVA for the part of the factorial experiment actually run
from Table 8.1.2 and using a CRD as given for Eq. (8.1.2), is given
in Table 8.1.3.
8.1 LATIN SQUARE 215
TABLE 8.1.3
Source df EMS
Within error 0 a2
The Latin square design is not too useful for fewer than four
treatments nor more than eight (Kempthorne, 1952).
Men
Machines Total
1 C 25 B 17 A 13 D 32 E 42 129
2 B 2 D 26 E 35 C 8 A 8 79
3 A 6 E 43 C 16 B 14 D 25 104
4 D 27 A 9 B 17 E 31 C 19 103
5 E 31 C 22 D 21 A __2 B 13 89
>tal 91 117 102 87 107 504
A B C D E
13 17 25 32 42
8 2 8 26 35
6 14 16 25 43
9 17 19 27 31
_2 13 22 21 31
38 63 90 131 182
Source df SS MS F
Columns (Men) 4 117.76 29.44
Rows (Machines) 4 284.96 71.24
Pressures 4 2598.96 649.74 42.44
Residual 12 183.68 15.31
Total 24 3185.36
where
Z I y. ., = 504
i j
CT="inr=10»160-64
Z Z y ? = 13,346
i j J
218 8. LATIN SQUARE TYPE DESIGNS
= 10278.4 - CT
= 117.76
= 284.96
= 2598.96
= 183.68
y. = y + S. + <$... + 0. + + F, + e,..,.
ijk i (i ) j '(j) k (ljk)
i = 1, 2, 18 j = 1, 2 k = 1, 2
where
th
^ijk = CUTnu^ative amount of drug in the blood of the i
subject after 8 hours given order j and formula k.
Subjects
Order 1 2 3 4 5 6 7 8 9
2 (Crossover) 43.5 45.4 35.4 43.7 49.8 39.9 41.4 37.5 39.5
Total 95.4 80.5 74.0 79.8 84.4 79.6 79.2 76.3 78.6
Subj ects
Order 10 11 12 13 14 15 16 17 18
2 (Crossover) 44.2 33.4 32.7 33.1 33.4 27.1 33.1 39.5 40.1
Total 95.0 74.5 71.8 68.8 67.1 58.3 67.4 71.3 74.4
, ( 6 6 8 V . (7 .8 .1 )2 . 52i62136 , ^
- 52,624.36 = 717.79
8.2 ASSOCIATED DESIGNS 221
- 52,624.36 = 1,093.98
Source df SS MS EMS
To find the smallest error from these data we may test for
order
H : a + (j)(0) = 0
o n
2.25
< 1
*1,16 ” 20.62
Hence accept Hq at a > 0.25 and pool mean squares, which actually
results in a RCBD where subjects act as blocks and formulas are
treatments. The error mean square with 17 df is:
F - = 2 25
1,17 19.54
ToKi^ c.
Tabled F a = 0.10
n = 3.03
1,17 a = 0.25 = 1.42
Order 1 2 3 4 5 6 7 8 9
Formula Two
1 Initial 51.9 35.1 38.6 36.1 34.6 39.7 37.8 38.8 39.1
Formula One
2 Crossover 58.5 60.4 50.4 58.7 64.8 54.9 56.4 52.5 54.5
Subjects
Order 10 11 12 13 14 15 16 17 18
Formula One
1 Initial 50.8 41.1 39.1 35.7 33.7 31.2 34.3 31.8 34.3
Formula Two
2 Crossover 62.2 51.4 50.7 51.1 51.4 45.1 51.1 57.5 58.1
8.3 REFERENCES
In the past it was felt that two leveled factorials were good
only for exploratory work by some research workers. Recently it
seems that the philosophy of experimenters has changed in this
regard. Many investigators believe it is better to introduce many
factors each with a high and low level into the experiment rather
225
226 9. 2n FACTORIAL EXPERIMENTS
than choosing arbitrarily a few factors and run many levels on each.
The reasoning for this whole approach seems to be that experimenters
who encounter many factors (independent variables) in their
investigations would prefer to have experimental, along with their
theoretical, evidence to discard independent variables rather than
using the theoretical approach only. The reason for this approach
is that the actual theory may be unknown but assumed based on past
experience. Thus, with the availability of rather efficient designs
of experiments involving two-leveled factorials this ’’theoretical
only” approach does not seem very scientific. In any case the
inference space must be carefully considered.
Notations
1 2 3
( 1) 00
a 10
b 01
ab 11
__________ b__________
a __________ 0____________ 1_____________ Total____
0 00 or (1) 01 or b (1) + b
1 10 or a 11 or ab a + ab
TABLE 9.1.1
Source df SS
A 1 [a + ab - (1) - b]2/4
B 1 [b + ab - (1) - a]2/4
AB 1 [ab + (1) - a - b]2/4
TABLE 9.1.2
2
Treatment Combinations and Effects in 2
Treatment combination
Average effect
Effect (i) a b ab divisor
A - + + 2
B - + + 2
AB + - + 2
Mean + + + + 4
TABLE 9.1.3
3
Treatment Combinations and Effects in 2
Treatment combination
Effect d)_ a b ab c ac be abc Divisor
A + - + - + - + 4
B - - + + - - + + 4
AB + - - + + - - + 4
C - - - - + + + + 4
AC + - + - - + - + 4
BC + + - - - - + + 4
ABC - + + - + - - + 4
Mean + + + + + + + + 8
As before all main effects have a plus (+) sign when the small
letter is. present in the treatment combination and a minus (-) when
absent. To get the signs for the interactions multiply the
corresponding main effect rows. Multiplication of like signs is
defined to give a positive sign and multiplication of unlike signs
is defined to give a negative sign. The analysis of variance for a
3
2 factorial experiment is given in Table 9.1.4.
TABLE 9.1.4
Source df
A 1
B 1
AB 1
C 1
AC 1
BC 1
ABC __1
Total 7
9.1 COMPLETE BLOCKS 231
A = [(a-l)(b+1)(c+1)(d+1)(e+1) ....]
TABLE 9.1.5
Source df
Main effects n
n factor interaction 1
Total 2n - 1
232 9. 2n FACTORIAL EXPERIMENTS
Factor Levels
Solution temperature (T) Low High
Solution concentration CQ Low High
Roll size (R) 1 2
Roll tension (F) Low High
Replication (blocks) (B) I II
If all the assumptions for an ANOVA are met, the ANOVA in this
case is similar to that given in Table 9.1.5 except that one must
account for the restriction on randomization in each replicate. In
addition one will have a measurement of error from the block by
treatment interactions.
(y = smoothness score).
9.1 COMPLETE BLOCKS 233
and
SS T = (64)2/32 = 128.0
C
L (-1) H (+1)
(+) 105 (-) 74
T C L(-l) H(+l)
Note that the sum of the values of the cells prior to the sign
alteration checks out to a total of 422. The TCF contrast is given
by
[TCF] = -56 + 49 + 34 - 40 + 58 - 77 - 40 + 68
= -4
and
SS TCF = (-4)2/32 = 0.5
„ .Based on.
Source df MS F(verror J
J
Replicates (B) 1 0.5
Restriction error 0
Treatments 15
T 1 128.000 15.4*
C 1 105.125 12.6*
TC 1 0.500 < 1
R 1 24.500 2.9
TR 1 36.125 4.3
CR 1 32.000 3.8
TRC 1 136.125 16.4*
F 1 66.125 8.0*
TF 1 72.000 8.7*
CF 1 15.125 1.8
TCF 1 0.500 < 1
RF 1 40.500 4.9*
TRF 1 3.125 < 1
CRF 1 2.000 < 1
TCRF 1 6.125 < 1
B x treatments 15 8.833
Replication I Replication II
CL b C d £ Strength (coded) a 6 c d £ Strength 1
1. 1 7 3H 1 62 3. 2 7 4 L 1 42
2. 2 7 4 L 3 45 6. 2 7 4H 3 39
3. 2 7 4 L 1 50 10. 2 5 4 L 3 62
4. 2 5 3 L 1 49 16. 1 7 4 L 3 48
5. 1 5 3H 3 60 21. 1 5 4H 1 47
6. 2 7 4H 3 61 29. 1 7 3H 3 61
7. 2 7 4 H 1 60 26. 2 5 3H 3 57
8. 2 5 4H 3 56 19.--1 5 4 L 3 61
9. 1 5 3H 1 59 13. 2 5 3 L 3 49
10. 2 5 4 L 3 61 2. 2 7 4 L 3 60
11. 1 7 4H 1 52 8. 2 5 4 H 3 43
12. 1 7 3 L 3 60 15. 1 7 3 L 1 61
13. 2 5 3 L 3 61 18. 2 5 4 H 1 41
14. 2 5 4 L 1 60 1. 1 7 3H 1 62
15. 1 7 3 L 1 62 20. 2 5 3H 1 61
16. 1 7 4 L 3 47 28. 1 7 4H 3 45
17. 2 7 3 L 3 58 31. 2 7 3 L 1 56
18. 2 5 4H 1 60 25. 1 5 3 L 1 59
19. 1 5 4 L 3 57 17. 2 7 3 L 3 61
20. 2 5 3H 1 61 7. 2 7 4H 1 61
21. 1 5 4H 1 61 4. 2 5 3 L 1 49
22. 2 7 3H 1 62 9. 1 5 3H 1 60
23. 2 7 3H 3 57 22. 2 7 3 H 1 56
24. 1 7 4 L 1 44 27. 1 5 4 L 1 60
25. 1 5 3 L 1 60 12. 1 7 3 L 3 61
26. 2 5 3H 3 61 14. 2 5 4 L 1 46
27. 1 5 4 L 1 58 23. 2 7 3 H 3 61
28. 1 7 4H 3 56 30. 1 5 4H 3 61
29. 1 7 3H 3 60 32. 1 5 3 L 3 60
30. 1 5 4H 3 61 24. 1 7 4 L 1 30
31. 2 7 3 L 1 53 5. 1 5 3H 3 60
32. 1 5 3 L 3 61 11. 1 7 4 H 1 50
9.2 INCOMPLETE BLOCKS 241
2
9.2.1 The 2 Factorial
Day 1 (1) ab
Day 2 a b
[(1) + ab - a - b]
TABLE 9.2.1
2
ANOVA Using 2 Blocks of 2 in a 2 Factorial
Source df
AB and/or days 1
Restriction error 0
A 1
B J.
Total 3
Note: Once we have block one we can take any remaining treatment
combination, say a, and multiply the treatment combinations
in block 1 by a to generate block 2;
a • (1) = a
a • ab = a^b = b
2
a • ac = a c = c
a • be = abc
TABLE 9.2.2
3
ANOVA Using 2 Factorial in 2 Blocks of 4 Each
Source df
ABC and/or Blocks 1
Restriction error 0
A 1
B 1
AB 1
C 1
AC 1
BC _1
Total 7
_________ AB__________
BC - +
Consider the interaction of this table, i.e., [the block with (b)
and (ac) plus the block with (1) and (abc)] minus [the block with
(a) and (be) plus the block with (ab) and (c)] is the same as
9.2 INCOMPLETE BLOCKS 245
ABD : x, + x0 + x. = 0, mod 2
1 2 4
ACE : + x^ + X(- = 0, mod 2,
TABLE 9.2.3
A Design of a 25 Factorial
Experiment in 4 Blocks of 8
ABD: 0 1 1 0
ACE: 0 1 0 1
a b c d e Block 1 Block 2 Block 3 Block 4
0 0 0 0 0 = (1) a b c
0 0 1 0 1 = ce ace bee e
0 1 0 1 0 = bd abd d bed
0 1 1 1 1 = bcde abcde cde bde
1 0 0 1 1 = ade de abde acde
X 0 1 1 0 = acd cd abed ad
1 1 0 •o 1 = abe be ae abce
1 1 1 0 0 = abc be ac ab
TABLE 9.2.4
Source df
Main effects 5
2 factor interactions 10
Error 13
(assume 3 factor and higher
interactions are zero)
Total 31
9.2 INCOMPLETE BLOCKS 247
TABLE 9.2.5
Replication
I II III
Block 1 Block 3 Block 5
a ab b ab a b
TABLE 9.2,.6
df
Replications 2
Blocks in Replications 3
Restriction error 0
A 1 Two-thirds information on
each. The sum of squares
B 1
are computed from those
AB 1 replicates in which the
effect or interaction is
Error _3_
not confounded.
Total 11
Replication
II III
Block 1 Block 3 Block 5
a =8 ab = 17 b = 13 ab = 17 a =9 b = 10
I y = 123
£ y2 = 1437
CT = ^12* = 1260-75
2 2 2 2 2 2
SS ui ir m• reps = --------------
cc blocks 18 + 25 + 13 + 30z + 18 + 19
*--------------- 10„
1266.75
= 1351.50 - 1266.75
= 84.75
250 9. 2n FACTORIAL EXPERIMENTS
= 8.00
Source df SS MS F
Replication 2 6.00
Blocks in reps 3 84.75
Restriction error 0
A 1 12.50 12.50 4.68
B 1 60.50 60.50 22.66*
AB 1 4.50 4.50 1.69
Error _3 8.00 2.67
Total 11 176.25
♦Significant at a = 0.05
9.3 REFERENCES
This chapter will cover two-level factors only, but one of the
methods of fractional replication given here can be applied to any
prime number equal leveled factorial such as 3, 5, 7, and so on.
Besides, if a nonprime number such as 4, 8, 9, and so on occur then
2 3 2
the experimenter can convert the factorial into a 2 , 2 , 3 , and
so on, which can be fractionated easily (see John, 1961 and 1966).
252
10.1 THE 1/2 REPLICATION 253
TABLE 10.1.1
Treatment combination
Effect a) a b ab c ac be abc
V + + + + + + + +
A - + - + - + - +
B - - + + - - + +
AB + - - + + - - +
C - - - - + + + +
AC + - + - - + - +
BC + + - - - - + +
ABC - + + - + - - +
254 10. FRACTIONALS FOR TWO-LEVELED FACTORS
(-) (+)
Block 1 Block 2
(1) a
ab b
ac c
be abc
TABLE 10.1.2
Treatment combination
Effect a b c abc
A + - - +
BC + - - +
B - + - +
AC - + - +
C - - + +
AB - - + +
U + + + +
ABC + + + +
Hence, the effect for A will be calculated the same as the BC effect;
similarly for B and AC; C and AB; and y and ABC. This feature is
the basis for saying A is an alias of BC, and so on to y is an alias
of ABC, because in this particular fractional replication whenever we
measure the effect of A we also measure the effect of BC and so on.
10.1 THE 1/2 REPLICATION 255
y = I = ABC
Where I means the identity and all effects and interactions may be
multiplied by it to show the results of the design. Thus we would
have the following display of confounded effects:
A = A2BC = BC
B = AB2C = AC
C = ABC2 = AB
TABLE 10.1.3
Source df
A and/or BC 1
B and/or AC 1
C and/or AB J.
Total 3
256 10. FRACTIONAL^ FOR TWO-LEVELED FACTORS
had been used, the relationship between A and BC would have been
negative as given in Table 10.1.4.
TABLE 10.1.4
A + +
BC + +
I = -ABC B = -AC
A = -BC C = -AB
Level
Factor 1 2
a. Aggregate gradation Fine Coarse
b. Compaction temperature 250°F 300°F
c. Asphalt content 5% 7%
d. Curing condition Unwrapped Wrapped
e. Curing temperature 45°F 72°F
AB:
B
199 = 54 + 14 + 73 + 58 49 + 85 + 17 + 29 = 180
+
ae + ac + ad + acde ab + abce + abde + abed
AC:
C +
Of course one can see that the SS for the main effects are
easily found from the marginal totals of the two-way tables. For
A we obtain
SS A . (1« • >80 - M 9 - i s i ) 2 , i 9 0 1
lo
10.1 THE 1/2 REPLICATION 259
for B
s s i s - t 1 5 1 * . 18.p . - .6 14 » .- . A 9a . i , 18 ,1
and for C
TABLE 10.1.5
Source df MS F
A 1 390.1 1.95
B 1 18.1 < 1
C 1 232.6 1.16
D 1 612.6 3.06
E 1 76.6 < 1
AB 1 27.6 < 1
AC 1 138.1 < 1
AD 1 1387.6 6.94^
AE 1 1105.6 5. 53^
BC 1 68.1 < 1
BD 1 3052.6 15.26^
BE 1 0.6 < 1
CD 1 0.1 < 1
CE 1 410.1 2.05
DE 1 3570.1 17.85^
Error 00 200
TABLE 10.2.1
9
ANOVA for 1/4 Fractional Replicated 2 Factorial
Source df
Main effects 9
Two factor interactions 36
(Remainder) error 82
Total 127
A rule that may be used for determining the signs of the effects in
the identity relation is that a minus sign is placed before all
effects which have an odd number of letters provided the intrablock
subgroup is used for the experiment. If we are interested in all
main effects we can obtain the following aliases for the 7 degrees
of freedom from the experiment:
(1) C1) 6 16 24 45
a b ab 10 8 21 -9 A = BC = BDE*
c be cbe 8 9 -4 -5 C = AB = DE
d e de 4 4 -8 -3 D = CE = ABE
10.
*
Fourth and higher order interactions are omitted.
FRACTIONALS FOR TWO-LEVELED
FACTORS
10.2 THE 1/4 REPLICATION 265
we obtain
[A] = -9
[BC] = 9
[ACDE] = 9
[BDE] = -9
Next we need to show the sums of squares for the main effects
and two factor interactions, AD and AE. These sums of squares can
now be obtained as follows:
266 10. FRACTIONALS FOR TWO-LEVELED FACTORS
SS - (Total effect)2
?L 2
c-
i=l 1
Sums of squares:
A - — - 10.1 AD = - 0.1
45.1 E.iiiii.is.i
C.i^.5,1 AE - - 3.1
4
Problem 10.2.3. Given a 1/2 replicate of a 2 factorial,
(a) Discuss the quality of the best design.
(b) Set up the identity and aliases.
(c) Show your intrablock subgroup.
(d) If the yields for the corresponding treatment
combinations are (1) =6; ab = 10; ac = 4, be = 5, ad = 9,
bd = 3; cd = 8 and abed = 0, use the Yates method and the
Berger procedure to find the best ANOVA assuming the error
mean square is known to be 5.0.
268 10. FRACTIONALS FOR TWO-LEVELED FACTORS
To write out the blocks we may use the four defining equations
(see Kempthorne, 1952, Chapter 20).
and using levels such that the equations are all set equal to zero
gives the intrablock subgroup. Since we will deal only with the
intrablock subgroup of the 1/4 replicate, the levels of ABCDE and
ABFGH will always be zero, but the levels of the ACF and BDG will
be 0 and 1 to make up the four blocks of 16 in the 1/4 replicate.
In other words, the four blocks may be represented as:
10.3 BLOCKING IN FRACTIONAL 269
ABCDE : 0 ACF :: 0
ABFGH : 0 BDG : 0
ABCDE : 0 ACF : 0
ABFGH : 0 BDG : 1
ABCDE :: 0 ACF :: 1
ABFGH :: 0 BDG : 0
ABCDE : 0 ACF :: 1
ABFGH : 0 BDG : 1
TABLE 10.3.1
o
ANOVA for 1/4 Fractional Replicated 2 in Blocks of 16
Source df
Total 63
I = ABC
A = BC
B = AC
C = AB
TABLE 10.4.1
ANOVA of y
Source df
A 1
B 1
AB 1
TABLE 10.4.2
ANOVA of y2
Source df
B 1
C 1
BC 1
I = ABCD
I II III IV
(1) a b d
ac c abc acd
bed abed cd be
abd bd ad ab
TABLE 10.4.3
A = ABD in I I and IV
B = ABCD in I and IV
C = BCD in I I and IV
D = ABCD in I and III
AB = BCD in I and IV
AC = ABCD in I I and IV
AD = BCD in I and III
BC = ABD in I and IV
BD = ACD in I and III Use the mean of these
or ABC in I and IV two estimates for BD
CD = ABD in I and III
ADF BEF
ADG BEG
ADH BEH
AEF CDF
AEG CDG
AEH CDH
BDF CEF
BDG CEG
BDH CEH
10.6 REFERENCES
280
11.1 CONFOUNDING IN 3n 281
b
0 1 2
II
0 00 01 02-- 0, mod 3
H-X
1
a l 10 11 12-- mod 3
— Xl = 1,
2 20 21 22-- mod 3
1 | 1 — Xl = 2,
I
X2 = fo x2 = 1
X2 = 2
mod 3 mod 3 mod 3
2, mod 3
TABLE 11.1.1
Source df
A 2
B 2 This breakdown of the
AB 4 I 2 Hinteraction has no )
*AB (2 df)* *experimental meaning.*
where
th
A^ = sumof observations in i level of a (i = 0, 1, 2)
th
Bj = sumof observations in j level of b (j= 0, 1, 2)
2 2
(AB C D)2 = total of the observations from the defining
equation X^ + 2X2 + 2X^ + X^ = 2, mod 3
2
(AC )1 = total of the observations from the defining
equation X^ + 2X^ = 1, mod 3
11.1 CONFOUNDING IN 3n 285
4
Problem 11.1.1. For a 3 factorial experiment, show the sum
2 2
of squares for AB CD . (Write out the equation, establish the
correct treatment combinations for the various totals and
indicate the way to obtain the sum of squares.)
For all designs with prime and equal number of levels, the
first letter of any part of an interaction may have an exponent of
unity. This condition is easily seen for the three-leveled
2 2
factorials. For example: AB = A B since the treatment combinations
for each interaction are combined in the same manner as shown below
(AB2)0 = 00 + 11 + 22 (a2b)0 = 00 + 11 + 22
(AB2)1 = 02 + 1 0 + 2 1 (A2B)2 = 0 2 + 1 0 + 21
(AB2)2 = 0 1 + 1 2 + 2 0 (A2B)1 = 01 + 12 + 20
TABLE 11.1.2
3
ANOVA of a 3 Factorial Experiment
Pieces that can
Source df
be confounded
A A 2
B B 2
AB 2
AB
AB 2
C C 2
AC 2
AC
AC 2
BC 2
BC
BC 2
ABC 2
ABC2 2
ABC 2
AB C 2
a b 2c 2 2
and so on
11.1 CONFOUNDING IN 3n 287
X Y XY XY
X Y XY XY
2 2
(a) ABC a b 2c A C = AC B
2 2 5 2
(b) ABC AB2 A C = AC B^C = BC
Results:
ABC : Xx + X2 + X3 = 0, 1, 2 mod 3
levels of ABC
levels of AB2
Blocks
0 0 0 0 0 1 0 0 2 0 2 1 2 0 0
1 1 1 1 1 2 1 1 0 10 2 0 1 1
2 2 2 2 2 0 2 2 1 2 10 1 2 2
TABLE 11.1.3
3
ANOVA for 3 Factorial in Blocks of 3
Source df
A 2
B 2
AB 2
C 2
AC 2
BC 2
Error and/or AB2C, ABC2, AB2C2 _6
Total 26
2
11.2.1 A3 Factorial Experiment in Blocks of 3
0 0 1 0 2 0
1 2 2 2 0 2
2 1 0 1 1 1
TABLE 11.2.1
Source df
Blocks and/or AB 2
Restriction error 0
A 2
B 2
AB (actually AB2) 4-2 = 2
a 2
AB (actually AB ) is the estimate of the whole interaction
2
AB coming from the part, AB .
variance for the data from this design is given in Table 11.2.2.
Refer to Example 9.3 for the method of calculation of SS.
TABLE 11.2.2
Source df Information
Replications (R) 3
Blocks in R 8
Restriction error 0
A 2 3/4 Balanced
B 2 3/4 information
. AB 3/4
AB 4 9
3/4 on all effects
AB 3/4
Error 16
B1 1 00 = 2 01 = 4 02 = 3 B1 4 00 = 3 10 = 5 20 = 6
B1 2 10 = 3 11 = 6 12 = 5 B1 5 01 = 5 11 = 7 21 = 8
B1 3 20 = 7 21 = 6 22 = 5 B1 6 02 = 2 12 = 7 22 = 3
III IV
B1 7 00 = 1 12 = 5 21 = 9 B1 10 00 = 4 11 = 6 22 = 1
B1 8 02 = 3 11 = 6 20 = 6 B1 11 01 = 6 12 = 5 20 = 4
B1 9 01 = 4 10 = 7 22 = 4 B1 12 02 = 7 10 = 8 21 = 7
(c) Show how you calculated the SS for all sources and
give the amount of information on each of A, B, and AB.
3
11.2.2 Two 3 Factorial Experiments
_X _Y XY XY2
AB BC AB2C AC2
4
11.2.3 Two 3 Factorial Experiments
X Y XY XY2
2 2 2 2 2
ABC AC D AB D BCV
TABLE 11.2.3
Source df
2
Blocks and/or ABC, AC D,
2 2 2 2
AB D , and BC D 8
Restriction error
Main effects 8
2-factor interactions 24
3-factor interactions 24
4-factor interactions 16
I = ABCDE = (ABCDE)2
2 2 2 2
A = AB C D E = BCDE
B = AB2CDE = ACDE
AB = ABC2D2E2 = CDE
and so on.
TABLE 11.3.1
Source df
This analysis shows that one obtains full information on the main
effects and two factor interactions from the one block (this is the
1/3 fraction).
243 = total df
-__9 = df confounded on line I
234 = df for main effects and interactions remaining.
234 df
= 13, the number of main effects
18 df for each set of aliases
TABLE 11.4.1
Source df
Blocks and/or
(a b 2c , bd 2e 2, ac d 2e 2) 2
Restriction error 0
Main effects 10
Two factor interactions 40
Residual (error) 28 May be broken into components
if the mean square is much
larger than within, then look
at the larger and smaller
pieces of these interactions
separately
Within duplicate error _3
Total 83
It follows that
2 2222
A = ABC = A B C D E
2
= AC DE = ABD E
22 2
= BC = BCDE = CD^E
22 2
= AB DE
and similarly for the other main effects. In general, ifthe two
factor interactions are considered smaller than the main effects,
a conclusion may be made to discard some factors from the next two
blocks if the main effects are quite small.
Columns
1 2 3
1 A B C
Rows 2 B C A
3 C A B
abc
where:
0 0 0
Xx + X2 + 2X3 = 0, mod 3
10 1
2 0 2 satisfies the block. Hence
0 11
I = ABC2 = (ABC2)2
1 1 2
A = AB2C = BC2
2 1 0
2 2 2
0 2 2 B = AB C = KC
12 0 C = AB = ABC
2 2 1
AB = AC = BC
and it can be seen that all main effects are confounded with parts
of two factor interactions. If interactions are zero, this
fractional factorial is of value to investigate main effects. In
large investigations, the fractional setup as a Latin square has
advantages as a screening device if the lower order terms express
the larger source of variation of the variable being analyzed. In
other words, even if there is confounding of the two factor
interactions with main effects, the main effects ought to be
larger and indicate the major source of variation.
= a b 3c 3 = a b 4c 4 = a b 2c 2 = BC
11.7 REFERENCES
Many other incomplete block designs exist and special plans are
available for some of them in Cochran and Cox (1957). The treatment
of all incomplete block designs can follow a different procedure
from the one presented here but the results are essentially the
same. In all instances in this chapter it is assumed that there is
complete randomization in the blocks considered.
302
12.1 SAME PRIME POWER 303
Treatment Pseudofactor
tf
CM
a ai *3
1 0 0 0
2 0 0 1
3 0 1 0
4 0 1 1
5 1 0 0
6 1 0 1
7 1 1 0
8 1 1 1
Factor 1 Factor 2
a b
*1 *2 bi fa2 b3
1 0 0 0 1 0 0 0
2 0 0 1 2 0 0 1
3 0 1 0 3 0 1 0
4 0 1 1 4 0 1 1
5 1 0 0 5 1 0 0
6 1 0 1 6 1 0 1
7 1 1 0 7 1 1 0
8 1 1 1 8 1 1 1
TABLE 12.1.1
Original Pseudofactorial
Source df Source df
A 7 1
Ai
1
A2
1
A3
1
A1A2
1
A1A3
1
A2A3
1
A1A2A3
B 7 1
B1
1
*2
1
B3
1
B1B2
1
B1B3
1
B2B3
B,B„B„ 1
12 3
AB 49 1
AiBi
A1B2 1
1
A1A2A3B1B2B3
Total 63 63
306 12. MIXED FACTORIALS
X Y XY
TABLE 12.1.2
Source df Information
Blocks 3
Restriction n
error u
A 6 6/7
B 6 6/7
AB 48 48/49
:2) is:
X Y XY
A1B1 A
2B2 A1A2B1B2
TABLE 12.1.3
Source df Information
Blocks 3
Restriction o
error
A 7 Full
B 7 Full
AB 46 46/49
x2 + x5 = °, 1 mod 2
Pseudofactor Factor
1 X2 X3 X4 X5 X6
0 0 0 0 0 0 1 1
1 0 0 1 0 0 5 5
0 1 0 0 1 1 3 4
1 1 1 1 1 0 8 7
TABLE 12.2.1
A 1 A
B 1 B
AB 1 AB
C 2 C
AC 2 AC
BC 2 BC
ABC 2 ABC
D 2 D
AD 2 AD
BD 2 BD
ABD 2 ABD
CD 4 CD and CD2
ACD 4 ACD and ACD2
BCD 4 BCD and BCD2
ABCD __4 ABCD and ABCD2
Total 35
X Y XY XY
c D CD CD'
df = 2 2 2 2
This is a very poor design because all main effects of C and D are
lost and so is their interaction.
X Y XY
AB CD ABCD
df = 1 2 2
This is a fairly good design because all main effects and 1/2 of CD
are estimable. The informationon AB,however, is lost.
X Y XY
A B AB
df = 1 1 1
X
CD or CD2
df = 2
This is not as bad a design as the one given in Section 12.2.3 since
one obtains half information on CD, and full information on all main
effects and the other interactions.
310 12. MIXED FACTORIALS
X
AB
df = 1
df — 1 2 2
JC Y XY
Replication 1 AB C ABC
Replication 2 A CD ACD
Replication 3 B CD2 BCD-
Replication 4 AB D ABD
a k ci c2 ^1 ^2
2 2 2 3 2 3
c
C1 C2
1 0 0
2 0 1
3 0 2
4 1 0
5 1 1
6 1 2
Factors Levels
a
f
t
CM I
1
0° 0 0
45° 0 1
90° 1 0
105° 1 1
if — ► 2 1 2 Number of
Condition X Y XY subjects
I. BC A1DE A1BCDE 2
II. BC A2DE A2BCDE 4
III. BC2 A^DE AjBC2DE 2
IV. BC2 A^DE A1A2BC2DE 4
12.3 LEVELS, PRODUCT OF PRIMES 313
The reader must understand that for each subject blocks are
chosen at random to run first, second and so on. Then, given a
block the treatment combinations within that block are chosen at
random to find out which will be run first, second and so on for
each sitting.
BC : X3 + X4 = 0 mod 3
A.DE : X. + X. + X. = 0 mod 2
1 1 5 6
13. I 0 0 0 1 0 13. 30 20 I F
14. 1 0 0 0 0 1 14. 90° 30 20 R B
15. 1 0 1 2 1 0 15. 90° 40 60 I F
16. 1 0 1 2 0 1 16. 90° 40 60 R B
to
oo
17. 1 0 2 1 1 0 17. 50 40 I F
o
o
18. 1 0 2 1 0 1 18. 50 40 R B
19. 1 1 0 0 1 0 19. 105° 30 20 I F
20. 1 1 0 0 0 1 20. 105° 30 20 R B
21. 1 1 1 2 1 0 21. 105° 40 60 I F
22. 1 1 1 2 0 1 22. 105° 40 60 R B
23. 1 1 2 1 1 0 23. 105° 50 40 I F
24. 1 1 2 1 0 1 24. 105° 50 40 R B
Assuming that all 24 different blocks have been set up for the
whole experiment we next outline the ANOVA showing the source, df
and amount of information for each source.
.3 LEVELS, PRODUCT OF PRIMES 315
TABLE 12.3.1
Subj ects 11
Blocks and/or interactions 60
Restriction error 0
A 3 Full
B 2 Full
AB 6 Full
C 2 Full
AC 6 Full
BC 4 Half
ABC 12 Full
D 1 Full
AD 3 Full
BD 2 Full
ABD 6 Full
CD 2 Full
ACD 6 Full
BCD 4 Full
ABCD 12 Full
E 1 Full
AE 3 Full
BE 2 Full
ABE 6 Full
CE 2 Full
ACE 6 Full
BCE 4 Full
ABCE 12 Full
DE 1 Full
ADE 3 2/3
BDE 2 Full
ABDE 6 Full
CDE 2 Full
ACDE 6 Full
BCDE 4 Full
ABCDE 12 5/6
Error (Residual) 1513
Total 1727
316 12. MIXED FACTORIALS
An example is 9 x 8 or
ai *2 61 b2 b3
3 3 2 2 2
12.6 FRACTIONALS 317
2 3
which results in a pseudofactorial experiment of 3 2 if incomplete
block designs are required.
ai a2 bl b2 C1 C2 dl d 2 ei £2 h <2 91 92 fel k 2
2 x 3 x 2 x 3 x 2x3x 2x3x 2x5x 2x5x 2x5x 2x5
8 4 4
or a pseudofactorial experiment 2 3 5 .
3
Problem 12.5.1. In a 6 experiment, show a good system of
confounding blocks of 36.
ai *2 *3 *1 *2 *3
0 0 0 1 1 1
1 1 0 1 0 0
1 0 1 0 1 0
0 1 1 0 0 1
2
Similarly using B2 for confounding the 3 part, l
Wx + W2 = 0, 1,
2, mod 3 (using Wfs insteadI of Xfs to identity the equation for
mod 3), the three resulting groups of treatment combinations are
r2
Ri R3
b 6, b, 6 bn
1 2 1 2 1 2
1 0
0 1
CM
2
12.7 LATTICE DESIGNS 319
TABLE 12.7.1
Source df
Replications 1
Restriction error 0
Treatments k2 - 1
R x T (error) k2 - 1
322 12. MIXED FACTORIALS
Similarly for the double simple lattice the analysis for RCBD is
given in Table 12.7.2
TABLE 12.7.2
Source df
Replications 3
Restriction error 0
Treatments k2 - 1
R x T (error) 3(k2 - 1)
(nx + l)(n2 + 3) s2
_ 2
(nx + 3)(n2 + 1) s j
where
n^ = degrees of freedom for the error in.the lattice
2
s^ « error mean square for the lattice
Effect Replication
confounded
Wlth______ 1 2 3 ... (k + lj/2
Rows A AB AB3 ... ABk"2
Columns B AB2 AB4 ... ABk_1
Effect Replication
confounded ,
with 1 ... — ^— (— 2 — + ^ ... (k + 1)
Designs that have only one factor with more than one level and
all the levels of the factor cannot be handled in a single block have
12.9 REFERENCES 325
12.9 REFERENCES
326
13.1 FIXED DESIGNS 327
Levels 3 2 7 4 2 3 6 5 2 3 5 6 3 4 5 2
328 13. RESPONSE SURFACE EXPLORATION
Then the first treatment combination for (X., X0, X,, ..., X,.) to
1 Z O ID
be used in the experiment is (3, 1, 6, ..., 2). For the next
treatment combination take a random number between 1 and 3, say 2
for Xj. Continue as above to obtain the second treatment combination.
This method is continued for all 25 treatment combinations with the
restriction that the number of treatment combinations containing
each level of each factor be as equal as possible. For example,
consider factor X ^ each of the three levels will be represented
at least eight times, but one level must be represented nine times;
whereas for factor Xg, each of the five levels must be present in
the experiment five times. This completes the design.
Factors (q)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
1 + + + + + + + + + + + - - - - -
2 + - + + +
3 - + .+ + - - - + - - + + + - + +
4 + + + - - - + - - + - - + + + +
g
5 + + - - - + - - + - + + + - + -
cd
C
•H + + + + + +
6 + - - - - - + - + -
1
CJ
+ + + +
7 - - - + - - + - + + + -
P
c
<D
a
p
8 - - + - - + - + + + - + - + - +
cd
a) + + +
U 9 - + - - + - + + + - - + -
H
10 + - - + - + + + - - - - + + - -
11 - - + - + + + - - - + - - - + +
12 _ + + + + _ _ _ + _ _ _ _ -
i i* Number of pairs
|s|
Designs q n
of columns = (^)
0 2 4 6 8
I (Plackett-Burman) 16 12 67 53 120
II 20 12 75 115 190
III 24 12 101 175 276
IV 24 18 193 83 276
V 30 18 281 154 435
VI 36 18 385 245 630
VII (Plackett-Burman) 30 24 295 81 59 435
Booth and Cox (1962) indicate that since the number of -4!s and +4*s
are equally likely over the set of designs from which this particular
design was evolved, the E(s) = 0 and the variance of s is
E(s^) = (53/120)(4^) + (67/120)(0^). Hence the variance for this
design is 7.07.
Systematic
Design q n Random balance supersatured
I 16 12 13.1 7.07
II 20 12 13.1 9.68
III 24 12 13.1 10.1
IV 24 18 19.1 13.6
V 30 18 19.1 15. 3
VI 36 18 19.1 17.4
VII 30 24 25.0 11.4
TABLE 13.1.1
0.100 16 22 29 44
0.050 32 45 59 90
0.025 64 91 119 182
0.010 161 230 299 459
0.005 322 460 598 919
In general,
, _ (m + q - 1) !
K " m! Cq - 1) I
k = + q ( q - 1) (g - 2)
2 6
y2 y23 y3
338 13. RESPONSE SURFACE EXPLORATION
+ Y 123X 1X 2X 3
. * W s < V X/ * 52 3 W X 2 - V
+ e 22x 2 + e33x 3
where y is the dependent variable and the X^, X2, X3 are the
independent variables. It is required that X1 + X2 + X3 = 1 in the
simplex lattice design. We can form the following four identities:
13.1 FIXED DESIGNS 339
2
Substituting these identities into the general equation for 30, X-,
2 2
X2, and X^, we will get, after simplification,
+ ^B12 " B11 ' B22^X1X2 + ^B13 ' B11 ' B33^X1X3
B0 + Bl + B11
I—1
>-
II
CM
B0 + B2 + B22
Y3 = B0 + B3 B33
Y12 = B12* Bll‘ B22
Yi13
z = B13 ■M
7 ' i 1 - ^33
B11 B,
and
-2L xi X2 X3
i 0 0
yi
0 1 0
y2
0 0 1
y3
1/2 1/2 0
y12
1/2 0 1/2
y13
y23 0 1/2 1/2
Y13 = 4y13 - 2yl - 2y3 and Y23 = 4y23 ’ 2y2 ' 2y3
2 r i [xi(2xi-i)]2 2 t4xixi)2i
v(y)=a2 I * I — ^ — \
Li=l i (i<j)=l ij J
2
where a = error variance from a previous experiment or estimated
from the within treatment combinations: r. or r..
l ij
(i < j = 1, 2, 3) = the number of observations at the
designated treatment combination (i or ij).
of the factors
i - i 1
q
I X. = 1,
i=l
and
0 < a. < X. < b. < 1
— i — l— l—
the a. and b. levels for all but one factor which is left blank,
1 l
e«g»> (a!> b2, a3* * as* b^), for a six factor experiment. This
Treatment Treatment
>ination combination
X1 X2 X3 X4 X1 X2 X3 X4
and the final centroid of the polyhedron, of course,, comes; from the
average of all eight treatment combinations and is
2
The squared multiple correlation coefficient (R ) for this
model is 0.9833, with only five degrees of freedom for residual.
2
If only Xj, x y x i x 2> xix3 > xlx4* X2X3 were u s e &> t*ie R would
be' 0.9829, with eight degrees of freedom. Since all four variables
still appeared in the latter model, the authors decided to retain
the full model. The reader should recognize that, as in any model
13.1 FIXED DESIGNS 347
where points (1), (2), and (3) make up the design of the experiment,
p = distance from the center point (0, 0) in standardized units
which the experimenter believes will give information on
the response surface; 8 = the angle from one point,
treatment combination to another. (In this case 0 = 360°/3
= 120°), and the resulting design is an equilaterial triangle
13.1 FIXED DESIGNS 349
y = e0 + eiXl + 62X2
where the usual assumptions apply to the nonorthogonal regression
analysis except that there is no estimate of error. To obtain an
error estimate, at least one of the three treatment combinations
must be repeated in the experiment. Usually forthis type of
experiment previous knowledge ofthe error is available.
y = 60 ♦ ♦ e2X2 ♦ 3n x2 + ^ ♦ B^X,
with no error estimate from the data unless there are repeated
observations for at least one of the treatment combinations. Box
and Wilson (1951, pp. 8 and 9), compare pentagonal and factorial.
There is little to choose between the two. Quadratic effects are
estimated more accurately with pentagonal and interactions less
accurately.
Treatment
combination X^ X2
1 p COS 6 p sin 0
2 p COS 20 p sin 20
•
.
•
CD
p COS sin n 2 0
c
p
CM
n2
n2+l 0 0
0 0
V ni
y = I = Z lZ 2 Z 3
Z1 = Z2Z3
Z2 = Z1Z3
Z3 " Z1Z2
y = a0 * OjZj ♦ a2 Z2 ♦
Z1 Z2 Z3
-1 -1 -1
Same as I = -ABC
1 1 -1
from fractional replicated
1 -1 1 factorial experiments
- 1 1 1
^2 ^3
2 15 5
14 23 5
14 15 27
2 23 27
352 13. RESPONSE SURFACE EXPLORATION
z
h fs h h
-1 -l -1 2 15 5
-1 -l 1 2 15 27
l -1 2 23 5
-1 l 1 2 23 27
or
1 -l -1 14 15 5
1 -l 1 14 15 27
1 l -1 14 23 5
1 l 1 14 23 27
The composite design Box and Wilson, (1951); Myers, (1971, Ch. 7;
and Davies, (1971, Ch. 11) has three parts: a basic two-leveled
factorial or fractional factorial, an extra point at the center of
the entire design and 2k (where k = number of factors) extra points,
one at either extreme of each factor and at the center of all other
factors. Hence in a composite design with a complete factorial
experiment in it there are 2 + 2k + 1 treatment combinations.
Correspondingly, if there was a fractional factorial instead of a
k
complete factorial experiment in the design, the 2 would be reduced
as required. The particular type of composite design depends on the
location of the extreme points. If the extreme points are located
at the same standardized distance from the center point as the
factorial points, the design is called a rotatable composite design
(sometimes the word “central" is included in the title of these
designs to indicate that there is a center point). If the extremes
are located at a distance that makes the squared terms in the model
orthogonal to each other, the design is called an orthogonal composite
design. Any other location of the extreme points may be used, but
the analysis is just a nonorthogonal regression analysis.
Number of Three-leveled
v
factors k factorial, 3 Composite, 2 + 2k + 1
2 9 9
3 27 15
4 81 25
5 243 43
5 (1/3 fractional) 81 (1/2 fractional) 27
6 729 77
6 (1/3 fractional) 243 (1/2 fractional) 45
y ■ e„ . CjXj . e 2* 2 . 6 jX 3 . 6 3 3 X 3 X 3 . 6 3 3 X 3X 3 . 6 3 3 X 3X 3
♦ 633 X 3 . 6 3 3 X 3 * 6 3 3 X 3 . k
1 2 3 W 3 + 61123X1X2X3
1233X1X2X3 +
612233X1X2X3
This design requires that the distance from the center point
to each of the star points be such that all effects and interactions
estimated in the second order model are orthogonal to one another.
This distance, d, is
1/2 1/4
[(F + T) - F1/2]2
4n
1/4
[(8 ♦ 7)
1/2
- 81 / 2 ] 2 1.215
■i
z. Treatment comt
Z2 Z3
1 -1 -1 -1 (1)
2 1 -1 -1 zi
3 -1 1 -1 Z2
4 1 1 -1 Z1Z2
5 -1 -1 1
Z3
6 1 1
Z1Z3
7 -1 1 1
Z2Z3
8 1 1 1
Z1Z2Z3
9 0 0 0
10 1.215 0 0
dl
11 -1.215 0 0
-dl
12 0 1.215 0
d2
13 0 -1.215 0
-d2
14 0 0 1.215
d3
15 0 0 -1.215
*d3
Treatment Combination
Effect 1 2 3 4 5 6 7 8 9 10 r 11 12 13 14 15
z1z2z3
V(a)
2
where a is the error variance, and Z's are as described in the
layout. As examples consider the variance of a^:
a2 n2
v(« ) = - 2 j = ----- 2 ------
ZZX 8 + 2(1.215)
2 2
XTf > _ a________________a
^alP 2 2 2 2
ZZn 8(. 27) + 5(-. 73) + 2(.746)Z
where only the three extreme points are used in the design.
13.2 SEQUENTIAL DESIGNS 363
1 1
2/3 2/6 2/l0 /2n(n+l)
1 1 __1_
2/3 2 /6 2/10 /2n(n+l)
2 1
2/3 2/6 2/10 /2n(n+l)
3
2/6 2/fi ) /2n(n+l)
ij
(n + 1)
/2n(n+l)
Independent variable
Point
364 13. RESPONSE SURFACE EXPLORATION
Independent variable
Point
1 . JL
2 2/3
1 . JL
2 2/3
The next step is to add the levels of each factor and multiply each
by 2/n (where n = the number of X's) or 2/2 = 1 for this case, and
subtract the level of the factor deleted. For variable 1 we have
| c- i + = °
| (. J —) = - —
2/3 2/3 2/3
has a yield of 58. For X^, X2, and X^ the experimenter felt that
increments of 10, 2, and 0.5, respectively, were equally important.
Using these increments to code about the centroid, we obtain the
four treatment combinations
1 _ JL_
’ 2/6
1 _ JL_
2/J " 2/6
2 _ _ l_
2/3 2/6
3
o,
2/6
The yields for the four treatment combinations were 52, 62, 61, 57,
respectively. Hence delete treatment combination 1 and replace it
by the three levels:
This new treatment combination (118.3, 44.0, 86.3), along with the
remaining three treatment combinations, forms the new simplex.
After running the new treatment combination, the sequential
procedure continues as discussed previously.
procedure if there is no error and somewhat poorer (but not too bad)
when error is present.
To obtain the next best (in the sense of giving most information
using the distance criterion) point from the remaining (n - 2)
points, one must find the point that is farthest from the first two.
This point may not be unique either (as for the first two points),
and the same criterion of taking the point with the smallest index
is utilized.
-2 -1 0 1 2
+2 1 2 3 4 5
+1 6 7 8 9 10
0 11 12 13 14 15
-1 16 17 18 19 20
-2 21 22 23 24 25
13.2 SEQUENTIAL DESIGNS 369
The pair of points (1, 25) and the pair (5, 21) are tied for
being farthest apart. Since the index on the pair (1, 25) is
smaller, namely 1, this pair is used in the design.
The same type of calculations are made for all 23 points and
it is found that 5 and 21 are tied for being farthest from design
points 1 and 25. Choosing the smaller index point, 5, the third
design point is 5.
One can continue this until all 25 points are in the design.
This concludes the demonstration of the design.
A + B + C — ► D + other products
Factor level
-1 +1
assuming a = 1.
xi X2 X3 X4 xs
In this calculation line (1) shows the base level for the
factors, line (2) shows the unit, i.e., the change in level for the
X-variable corresponding to the change from 0 to 1 for the
Z-variable, and line (3) shows the estimated slopes calculated on
the assumption that three and higher factor interactions are
negligible. The levels of the factors are then changed in proportion
to these slopes. For every 7.9 units increase in Z^, should be
decreased 2.2 units, Z^ increased 6.0 units, and so on. Hence, for
each 25 x 7.9 cc. that is changed, X2 should be changed by (.25)
x (-2.2) mole. These changes are given on line (4) and define the
374 13. RESPONSE SURFACE EXPLORATION
y = 6 + 3 Xj + 5 X2 - 4 Xj - 3 Xj - 2 Xj X2 (13.2.1)
one may find the contour for y = 8 by setting Eq. (13.2.1) equal to
8 and substituting a value for X^. This makes the equation a
quadratic in X2 which allows one to solve for two points. Next,
change the value of X^ and solve for two new X2 points. Continue
this until one obtains the contour of 8 and the picture shown in
Fig. 13.2.1.
f£ - . 3 - 8 h - 2 X2 ■ 0
2 17
Solve and obtain the point of maximum, (X^ = — , X2 = -jj ), and
refer to it as M. Substituting these solutions into Eq. (13.2.1),
we obtain the maximum y^ = 8.2.
bu - » 1 /2 b,2 -4-A -1
r<
-1 -3-A
1
CM
CM
1/2 b12
X2 + 7X + 11 = 0.
where the and axes are labeled by observing that the contours
drop off faster in the direction than in the SC^ direction.
13.2 SEQUENTIAL DESIGNS 377
y - yM =
An example:
If
y - yM = B11& 1 * B2 & 2
50-
Q>
E
40- •2
_L_ _L
1220 1260
Temperature
Operating condition 1 2 3 4 5
The form for Cycle 1 has very little filled in since there was
no previous information and no estimate of the standard deviation
can be made with only one set of observations. The calculations
shown for Cycle 2 in Table 13.2.2 shows the form much more complete.
All data shown under the heading of "Previous” have been transferred
from the form used for the calculations of Cycle 1. The value of
the range in line (iv) is calculated by finding the range of the
five differences shown on that same line. In order to find the
New s use the value of f^ n given in Table 13.2.4. Note that the
two sigma error limits For New Effects is 0.89 which is smaller than
the temperature and time effects indicating that there is
13.2 SEQUENTIAL DESIGNS 381
* EX.y re i
^ = ~2 (temperature effect).
Using the data of Cycle 3 then, the calculations are the following:
X1 X2
Base level 1240 45
Unit 20 5
Estimated slope -1.45 -0.85
Unit x slope -29.0 -4.25
Amount of change per 20° -20.0 -2.93
TABLE 13.2.1 382
Cycle: n = 1 Phase: 1
Time 1
9 . Response: Tensile strength (1000 lb) Date: 5/20 - 5/22
Temp 4
Calculation of averages Calculation of standard deviation s
Operating Conditions (1) (2) (3) (4) (5)
(v) New sums [(i)+(iii)] 27.0 28.6 25.1 26.0 26.8 New sum s =
(vi) New average y^[(v)/n] 27.0 28.6 25.1 26.0 26.8 New average s = (New sum s)/(n-l) =
Cycle: n = 2 Phase: 1
Time l
9 . Response: Tensile strength (1000 lb) Date: 5/23 - 5/25
Temp 4
Calculation of averages Calculation of standard deviation s
Operating Conditions CD (2) (3) (4) (5)
SEQUENTIAL DESIGNS
(i) Previous cycle sum 27.0 28.6 25.1 26.0 26.8 Previous sum s =
(ii) Previous cycle average 27.0 28.6 25.1 26.0 26.8 Previous average s =
(iii) New observations 26.2 29.1 24.0 26.4 27.8 New s = range x f^ n = 2.1 x .30 = 0.63
5*,2
(iv) Differences [(ii)-(iii)] 0.8 -0.5 1.1 -0.4 -1.0 Range of (iv) * 1.1 -■ (-1.0) = 2.1
(v) New sums [(i)+(iii)] 53.2 57.7 49.1 52.4 54.6 New sum s = 0.63
(vi) New averages y^[(v)/n] 26.6 28.8 24.6 26.2 27.3 New average s = (New sum s)/(n-l)= 0.63
Cycle: n = 3 Phase: 1
Time 1
^ A Response: Tensile strength (1000 lb) Date: 5/27 - 5/29
1 Temp 4
Calculation of averages Calculation of standard deviation s
Operating conditions (1) (2) (3) (4) (5)
(i) Previous cycle sum 53.2 57.7 49.1 52.4 54.6 Previous sum s =0.63
(ii) Previous cycle average 26.6 28.8 24.6 26.2 27.3 Previous average s =0.63
(iii) New observations 26.5 29.2 24.2 25.6 27.0 New s = range x f^ = 0.35
(iv) Differences [(ii)-(iii)] 0.1 -0.4 0.4 0.6 0.3 Range of (iv) = 1.0
13.
(v) New sums [(i)+(iii)] 79.7 86.9 73.3 78.0 81.6 New sum s =0.98
(vi) New averages y^[(v)/n] 26.6 29.0 24.4 26.0 27.2 New average s = (New sum s)/(n-l) = 0.49
TABLE 13.2.4
Values of f,
k,n
n
2 3 4 5 6 7 8 9 10
k
4 0.34 0.40 0.42 0.43 0.44 0.45 0.45 0.46 0.46
5 0.30 0.35 0.37 0.38 0.39 0.40 0.40 0.40 0.41
6 0. 28 0.32 0.34 0.35 0.36 0.37 0.37 0.37 0.37
7 0.26 0.30 0.32 0.33 0.34 0.34 0.35 0.35 0.35
8 0.25 0.29 0.30 0.31 0.32 0.32 0.33 0.33 0.33
9 0.24 0.27 0.29 0.30 0.31 0.31 0.31 0.32 0.32
10 0.23 0.26 0.28 0.29 0.30 0.30 0.30 0.31 0.31
386 13. RESPONSE SURFACE EXPLORATION
13.3 REFERENCES
Random Numbers
0 0 -0 4 0 5 -0 9 1 0 -1 4 1 5 -1 9 2 0 -2 4 2 5 -2 9 3 0-3 4 3 5 -3 9 4 0 -4 4 4 5 -4 9
00 01826 72696 67261 13748 57834 27748 47472 43428 85524 19311
01 70731 12890 90395 45245 71282 15960 02749 86763 80564 02651
02 46616 84522 17249 78172 14197 84272 53226 96719 83462 05628
03 46384 26607 53444 68780 00458 27488 04949 30717 14304 60222
04 58645 90453 83413 15983 42345 27118 29425 42172 58412 31268
05 41874 83916 07454 42647 09410 12882 80001 98932 63277 94255
06 26032 68402 68176 87347 52572 79056 22703 83175 95807 60134
07 56142 74950 98878 36160 28616 77908 04908 63013 30555 51831
08 57893 57918 56290 00612 53884 76761 47934 20050 84977 11053
09 46729 47281 30473 94911 69061 87486 96782 90636 91405 01219
10 14107 46068 15859 99140 29872 70750 35757 25344 84845 95688
11 27713 13273 09015 42262 24580 30925 96900 56246 10613 17230
12 49886 62560 23023 09812 91948 92265 77407 93047 45352 61689
13 86651 47164 79270 45746 73141 67388 08454 68246 39046 57933
14 96133 52026 25837 64698 06911 21730 86390 19749 25859 80005
15 00926 46212 00204 58304 73907 97914 06786 47324 81225 50754
16 89036 39245 74371 87236 92131 38908 22146 15710 36858 35904
17 79592 66982 87984 30503 86953 56490 20414 07662 11122 95517
18 25084 83491 79984 64583 33924 36146 54692 44989 71603 31274
19 90261 79382 24483 75975 14765 49897 75468 66940 82187 68875
20 94937 36303 85514 76537 86864 88600 51843 18255 49660 87326
21 10178 16103 33283 38239 32402 02950 90428 24405 38580 61198
22 05983 63553 93546 66079 18389 57705 86746 62498 71642 60219
23 80832 36191 32042 76015 51331 40592 04887 50771 23810 64650
24 85889 08924 41542 04758 98753 52695 28165 87031 00479 95974
25 77032 15760 83026 92461 51806 65495 32148 40714 49107 05758
26 84080 67254 89239 13272 50218 17737 00242 15203 69060 39047
27 41022 70272 73710 17625 51325 69525 63464 89526 76738 80210
28 72398 03426 36476 62922 43624 98779 63289 34550 94270 01263
29 52224 04363 92979 26520 91076 13849 17740 25964 22169 23110
30 26223 68525 82483 15232 93800 78903 18831 19546 50469 74328
31 57872 36109 88383 83512 58763 57230 35952 75716 57094 53951
32 44529 90455 96666 31804 24979 51863 42983 13367 14111 38541
33 13586 36649 80041 54602 17281 22392 63074 53831 85782 94644
34 49084 08133 45510 21472 14644 38592 04490 52187 33945 78491
35 42309 57412 33314 07820 31345 45074 93547 30023 47632 93222
36 08174 49231 02588 36639 53978 40177 45572 44117 97946 44214
37 58679 57695 94704 16260 41928 27300 14053 93050 87103 34434
38 46847 17986 19347 74125 64945 97496 00465 65830 56564 46884
39 91053 64137 26104 29911 01242 58960 60582 83119 35911 68859
40 98796 24540 75578 79529 63199 78959 60734 73433 91995 81783
41 66824 16251 63616 91787 77519 14592 15333 63525 88097 79340
42 35492 57116 78056 72780 44502 97058 25342 53590 29196 06552
43 53794 57593 27187 99355 12652 77611 95229 41602 99946 94801
44 86332 05856 07786 78241 13584 29243 81609 55078 25937 90051
15 56530 97998 92028 18839 05879 43900 91769 87557 50750 05667
46 22966 67453 46806 57010 49278 99735 84996 74079 77917 48474
47 17068 92733 58021 81102 19871 17442 45375 27328 58055 20323
48 21570 52311 25131 76662 18619 91141 89382 41329 91702 39975
49 09699 95392 78594 40396 83605 01008 48694 85571 64427 23547
388
A p pen d ix 2
z 0.00 0 .0 1 0 .0 2 0 .0 3 0 .0 4 0 .0 5 0 .0 6 0 .0 7 0 .0 8 0 .0 9
0 .0 0 .0 0 0 0 0 .0 0 4 0 0 .0 0 8 0 0 .0 1 2 0 0 .0 1 6 0 0 .0 1 9 9 0 .0 2 3 9 0 .0 2 7 9 0 .0 3 1 9 0 .0 3 5 9
0.1 .0 3 9 8 .0438 .0 4 7 8 .05 1 7 .0 5 5 7 .0 5 9 6 .0 6 3 6 .06 7 5 0 .7 1 4 .0 7 5 3
0 .2 .0793 .0 8 3 2 .0871 .0 9 1 0 .0 9 4 8 .0 9 8 7 .1 0 2 6 .1 0 6 4 .1 1 0 3 .1141
0 .3 .1 1 7 9 .1 2 1 7 .1 2 5 5 .1 2 9 3 .1331 .1 3 6 8 .14 0 6 .1 4 4 3 .1 4 8 0 .1 5 1 7
0 .4 .1 5 5 4 .1591 .1628 .16 6 4 .1 7 0 0 .1 7 3 6 .1 7 7 2 .1 8 0 8 .1 8 4 4 .1 8 7 9
2 .0 .4 7 7 2 .47 7 8 .4 7 8 3 .4 7 8 8 .4 7 9 3 .4 7 9 8 .4 8 0 3 .4 8 0 8 .4 8 1 2 .4 8 1 7
2.1 .4821 .4 8 2 6 .4 8 3 0 .4 8 3 4 .4838 .4 8 4 2 .4 8 4 6 .4 8 5 0 .4 8 5 4 .4 8 5 7
2.2 .48 6 1 .4 8 6 4 .4 8 6 8 .4871 .48 7 5 .4 8 7 8 .48 8 1 .4 8 8 4 .4 8 8 7 .4 8 9 0
2 .3 .4 8 9 3 .4 8 9 6 .4 8 9 8 .4901 .4904 .4 9 0 6 .4 9 0 9 .4911 .49 1 5 .4 9 1 6
2 .4 .4 9 1 8 .4 9 2 0 .4 9 2 2 .49 2 5 .4 9 2 7 .4 9 2 9 .4931 .49 3 2 .49 3 4 .4 9 3 6
2 .5 .4 9 3 8 .4 9 4 0 .4 9 4 1 .4 9 4 3 .4 9 4 5 .4 9 4 6 .4 9 4 8 .4 9 4 9 .4951 .4 9 5 2
2 .6 .4 9 5 3 .4 9 5 5 .4 9 5 6 .49 5 7 .4 9 5 9 .4 9 6 0 .4961 .4 9 6 2 .4 9 6 3 .4 9 6 4
2 .7 .4 9 6 5 .4 9 6 6 .4 9 6 7 .49 6 8 .4 9 6 9 .4 9 7 0 .4971 .4 9 7 2 .4 9 7 3 .4 9 7 4
2.8 .4 9 7 4 .49 7 5 .4 9 7 6 .4 9 7 7 .4 9 7 7 .4 9 7 8 .4 9 7 9 .4 9 7 9 .4 9 8 0 .4981
2 .9 .4981 .49 8 2 .49 8 2 .4 9 8 3 .4 9 8 4 .4 9 8 4 .4 9 8 5 .4 9 8 5 .4 9 8 6 .4 9 8 6
3 .6 .4998 .4 9 9 8 .4 9 9 9 .4 9 9 9 .4 9 9 9 .4 9 9 9 .49 9 9 .4 9 9 9 .4 9 9 9 .4 9 9 9
3 .9 .5 0 0 0
389
A p pen d ix 3
P e r c e n ta g e p o i n t s o f th e t - d i s t r i b u t i o n
Q=0. 4 0 .2 5 0. 1 0 .0 5 0 .0 2 5 0.01 0 .0 0 5 0 .0 0 2 5 0 .0 0 1 0 .0 0 0 5
V
2Q=0. 8 0 .5 0.2 0.1 0 .0 5 0.02 0 .0 1 0 .0 0 5 0.002 0 .0 0 1
1 0 .3 2 5 1 .0 0 0 3 .0 7 8 6 .3 1 4 1 2 .7 0 6 3 1 .8 2 1 6 3 .6 5 7 1 2 7 .3 2 3 1 8 .3 1 6 3 6 .6 2
2 .2 8 9 0 .8 1 6 1 .8 8 6 2 .9 2 0 4 .3 0 3 6 .9 6 5 9 .9 2 5 1 4 .0 8 9 2 2 .3 2 7 3 1 .5 9 8
3 .2 7 7 .7 6 5 1 .6 3 8 2 .3 5 3 3 .1 8 2 4 .5 4 1 5 .8 4 1 7 .4 5 3 1 0 .2 1 4 1 2 .9 2 4
4 .2 7 1 .7 4 1 1 .5 3 3 2 .1 3 2 2 .7 7 6 3 .7 4 7 4 .6 0 4 5 .5 9 8 7 .1 7 3 8 .6 1 0
5 0 .2 6 7 0 .7 2 7 1 .4 7 6 2 .0 1 5 2 .5 7 1 3 .3 6 5 4 .0 3 2 4 .7 7 3 5 .8 9 3 6 .8 6 9
6 .2 6 5 .7 1 8 1 .4 4 0 1 .9 4 3 2 .4 4 7 3 .1 4 3 3 .7 0 7 4 .3 1 7 5 .2 0 8 5 .9 5 9
7 .2 6 3 .711 1 .4 1 5 1 .8 9 5 2 .3 6 5 2 .9 9 8 3 .4 9 9 4 .0 2 9 4 .7 8 5 5 .4 0 8
8 .2 6 2 .7 0 6 1 .3 9 7 1 .8 6 0 2 .3 0 6 2 .8 9 6 3 .3 5 5 3 .8 3 3 4 .5 0 1 5 .0 4 1
9 .2 6 1 .7 0 3 1 .3 8 3 1 .8 3 3 2 .2 6 2 2 .8 2 1 3 .2 5 0 3 .6 9 0 4 .2 9 7 4 .7 8 1
10 0 .2 6 0 0 .7 0 0 1 .3 7 2 1 .8 1 2 2 .2 2 8 2 .7 6 4 3 .1 6 9 3 .5 8 1 4 .1 4 4 4 .5 8 7
11 .2 6 0 .6 9 7 1 .3 6 3 1 .7 9 6 2 .2 0 1 2 .7 1 8 3 .1 0 6 3 .4 9 7 4 .0 2 5 4 .4 3 7
12 .2 5 9 .6 9 5 1 .3 5 6 1 .7 8 2 2 .1 7 9 2 .6 8 1 3 .0 5 5 3 .4 2 8 3 .9 3 0 -4 .3 1 8
13 .2 5 9 .694 1 .3 5 0 1 .7 7 1 2 .1 6 0 2 .6 5 0 3 .0 1 2 3 .3 7 2 3 .8 5 2 4 .2 2 1
14 .2 5 8 .6 9 2 1 .3 4 5 1 .761 2 .1 4 5 2 .6 2 4 2 .9 7 7 3 .3 2 6 3 .7 8 7 4 .1 4 0
15 0 .2 5 8 0 .6 9 1 1 .3 4 1 1 .7 5 3 2 .1 3 1 2 .6 0 2 2 .9 4 7 3 .2 8 6 3 .7 3 3 4 .0 7 3
16 .2 5 8 .6 9 0 1 .3 3 7 1 .7 4 6 2.120 2 .5 8 3 2 .9 2 1 3 .2 5 2 3 .6 8 6 4 .0 1 5
17 .2 5 7 .6 8 9 1 .3 3 3 1 .7 4 0 2.110 2 .5 6 7 2 .8 9 8 3 .2 2 2 3 .6 4 6 3 .9 6 5
18 .2 5 7 .68 8 1 .3 3 0 1 .7 3 4 2.1 0 1 2 .5 5 2 2 .8 7 8 3 .1 9 7 3 .6 1 0 3 .9 2 2
19 .2 5 7 .68 8 1 .3 2 8 1 .7 2 9 2 .0 9 3 2 .5 3 9 2 .8 6 1 3 .1 7 4 3 .5 7 9 3 .8 8 3
20 0 .2 5 7 0 .6 8 7 1 .3 2 5 1 .7 2 5 2 .0 8 6 2 .5 2 8 2 .8 4 5 3 .1 5 3 3 .5 5 2 3 .8 5 0
21 .2 5 7 .68 6 1 .3 2 3 1 .7 2 1 2 .0 8 0 2 .5 1 8 2 .8 3 1 3 .1 3 5 3 .5 2 7 3 .8 1 9
22 .2 5 6 .6 8 6 1 .3 2 1 1 .7 1 7 2 .0 7 4 2 .5 0 8 2 .8 1 9 3 .1 1 9 3 .5 0 5 3 .7 9 2
23 .2 5 6 .6 8 5 1 .3 1 9 1 .7 1 4 2 .0 6 9 2 .5 0 0 2 .8 0 7 3 .1 0 4 3 .4 8 5 3 .7 6 7
24 .2 5 6 .6 8 5 1 .3 1 8 1.7 1 1 2 .0 6 4 2 .4 9 2 2 .7 9 7 3 .0 9 1 3 .4 6 7 3 .7 4 5
25 0 .2 5 6 0 .6 8 4 1 .3 1 6 1 .7 0 8 2 .0 6 0 2 .4 8 5 2 .7 8 7 3 .0 7 8 3 .4 5 0 3 .7 2 5
26 .2 5 6 .6 8 4 1 .3 1 5 1 .7 0 6 2 .0 5 6 2 .4 7 9 2 .7 7 9 3 .0 6 7 3 .4 3 5 3 .7 0 7
27 .2 5 6 .6 8 4 1 .3 1 4 1 .7 0 3 2 .0 5 2 2 .4 7 3 2 .7 7 1 3 .0 5 7 3 .4 2 1 3 .6 9 0
28 .2 5 6 .6 8 3 1 .3 1 3 1 .7 0 1 2 .0 4 8 2 .4 6 7 2 .7 6 3 3 .0 4 7 3 .4 0 8 3 .6 7 4
29 .2 5 6 .6 8 3 1 .3 1 1 1 .6 9 9 2 .0 4 5 2 .4 6 2 2 .7 5 6 3 .0 3 8 3 .3 9 6 3 .6 5 9
30 0 .2 5 6 0 .6 8 3 1 .3 1 0 1 .6 9 7 2 .0 4 2 2 .4 5 7 2 .7 5 0 3 .0 3 0 3 .3 8 5 3 .6 4 6
40 .2 5 5 .6 8 1 1 .3 0 3 1 .6 8 4 2 .0 2 1 2 .4 2 3 2 .7 0 4 2 .9 7 1 3 .3 0 7 3 .5 5 1
60 .2 5 4 .6 7 9 1 .2 9 6 1 .6 7 1 2.000 2 .3 9 0 2 .6 6 0 2 .9 1 5 3 .2 3 2 3 .4 6 0
120 .2 5 4 .6 7 7 1 .2 8 9 1 .6 5 8 1 .9 8 0 2 .3 5 8 2 .6 1 7 2 .8 6 0 3 .1 6 0 3 .3 7 3
OO .2 5 3 .6 7 4 1 .2 8 2 1 .6 4 5 1 .9 6 0 2 .3 2 6 2 .5 7 6 2 .8 0 7 3 .0 9 0 3 .2 9 1
i s t h e u p p e r - t a i l a r e a o f th e d i s t r i b u t i o n f o r v d e g r e e s o f fre e d o m ,
a p p r o p r i a t e f o r u se in a s i n g l e - t a i l te st. F or a t w o - t a i l t e s t , 2Q must be u s e d .
390
A p p e n d ix 4
Number o f O b s e r v a t io n s p e r S am ple:
U sing t f o r d i f f e r e n c e o f means
S in g le -S id e d T e st a = 0 .0 0 5 a * 0.01 a = 0 .0 2 5 0 .0 5
D o u b le -S id e d T e s t a = 0 .0 1 a = 0 .02 a = 0 .0 5 0. 1
0 .0 1 0 .0 5 0 .1 0 .2 0 .5 0 .0 1 0 .0 5 0 .1 0 .2 0 .5 0 .0 1 0 .0 5 0 .1 0 .2 0 .5 0 .0 1 0 .0 5 0 .1 0 .2 0 .5
0 .3 0 123 87 61 0 .3 0
0 .4 0 85 70 100 50 108 78 35 0 .4 0
0 .5 0 96 55 106 82 45 106 86 64 32 88 70 51 23 0 .5 0
0 .6 0 101 85 67 39 90 74 58 32 104 74 60 45 23 89 61 49 36 16 0 .6 0
0 .7 0 100 75 63 50 29 90 66 55 43 24 76 55 44 34 17 66 45 36 26 12 0 .7 0
0 .7 5 88 66 55 44 26 79 58 48 38 21 67 48 39 29 15 57 40 32 23 11 0 .7 5
0 .8 0 77 58 49 39 23 70 51 43 33 19 59 42 34 26 14 50 35 28 21 10 0 .8 0
0 .8 5 69 51 43 35 21 62 46 38 30 17 52 37 31 23 12 45 31 25 18 9 0 .8 5
0 .9 0 62 46 39 31 19 55 41 34 27 15 47 34 27 21 11 40 28 22 16 8 0 .9 0
0 .9 5 55 42 35 28 17 50 37 31 24 14 42 30 25 19 10 36 25 20 15 7 0 .9 5
1.0 0 50 38 32 26 15 45 33 28 22 13 38 27 23 17 9 33 23 18 14 7 1.00
V alu e o f 1.1 42 32 27 22 13 38 28 23 19 11 32 23 19 14 8 27 19 15 12 6 1 .1
1 .2 36 27 23 18 11 32 24 20 16 9 27 20 16 12 7 23 16 13 10 5 1 .2
D - ia 1 .3 31 23 20 16 10 28 21 17 14 8 23 17 14 11 6 20 14 11 9 5 1 .3
1 .4 27 20 17 14 9 24 18 15 12 8 20 15 12 10 6 17 12 10 8 4 1 .4
1 .5 24 18 15 13 8 21 16 14 11 7 18 13 11 9 5 15 11 9 7 4 1 .5
1.6 21 16 14 11 7 19 14 12 10 6 16 12 10 8 5 14 10 8 6 4 1.6
1 .7 19 15 13 10 7 17 13 11 9 6 14 11 9 7 4 12 9 7 6 3 1 .7
1 .8 17 13 11 10 6 15 12 10 8 5 13 10 8 6 4 11 8 7 5 1.8
1 .9 16 12 11 9 6 14 11 9 8 5 12 9 7 6 4 10 7 6 5 1 .9
2.0 14 11 10 8 6 13 10 9 7 5 11 8 7 6 4 9 7 6 4 2 .0
2 .5 10 8 7 6 4 9 7 6 5 4 8 6 5 4 6 5 4 3 2 .5
3 .0 8 6 6 5 4 7 6 5 4 3 6 5 4 4 5 4 3 3 .0
3 .5 6 5 5 4 3 6 5 4 4 5 4 4 3 4 3 3 .5
4 .0 6 5 4 4 5 4 4 3 4 4 3 4 4 .0
Appendix 5
F-Distribution
Upper 25%
4 5 6 7 8 9 10 12 15 20 24 30
V2
1 5,,83 7..50 8,.20 8.58 8.82 8.98 9. 10 9.19 9. 26 9.32 9.41 9. 49 9. 58 9. 63 9. 67
2 2.,57 3,.00 3,.15 3.23 3.28 3.31 3.34 3.35 3. 37 3.38 3.39 3.41 3.43 3. 43 3. 44
3 2,.02 2,.28 2,.36 2.39 2.41 2.42 2.43 2.44 2.44 2.44 2.45 2.46 2. 46 2. 46 2. 47
4 1.,81 2,.00 2,.05 2. 06 2.07 2.08 2.08 2.08 2.08 2.08 2.08 2.08 2.08 2.08 2. 08
5 1.,69 1..85 1,,88 1,.89 1,.89 1,.89 1.,89 1.,89 1..89 1,.89 1.,89 1,.89 1,.88 1..88 1,.88
6 1.,62 1.,76 1,.78 1..79 1,.79 1,.78 1.,78 1,,78 1.,77 1,.77 1.,77 1,.76 1,.76 1,.75 1..75
7 1.,57 1.,70 1.,72 1..72 1,.71 1,.71 1,.70 1.,70 1.,69 1,.69 1.,68 1,.68 1,.67 1..67 1,.66
8 1..54 1..66 1..67 1,.66 1,.66 1,.65 1..64 1.,64 1.,63 1,.63 1.,62 1,.62 1,.61 1,.60 1,.60
9 1.,51 1.,62 1,.63 1..63 1,.62 1,.61 1..60 1,,60 1.,59 1,.59 1.,58 1,.57 1,.56 1..56 1,.55
10 1,.49 1.,60 1,.60 1,.59 1,.59 1 .58 1..57 1,,56 1..56 1,.55 1,.54 1,.53 1,.52 1,.52 1,.51
11 1..47 1.,58 1..58 1,.57 1,.56 1 .55 1.,54 1.,53 1.,53 1,.52 1,.51 1,.50 1,.49 1,.49 1,.48
12 1..46 1,.56 1,.56 1,.55 1 .54 1 .53 1..52 1,.51 1,.51 1,.50 1,.49 1,.48 1,.47 1,.46 1,.45
13 1..45 1..55 1,.55 1,.53 1,.52 1,.51 1,.50 1,.49 1.,49 1 .48 1,.47 1 .46 1,.45 1,.44 1 .43
14 1..44 1..53 1,.53 1,.52 1,.51 1 .50 1.,49 1,.48 1,.47 1 .46 1,.45 1 .44 1 .43 1 .42 1 .41
15 1,.43 1,.52 1,,52 1..51 1,.49 1,.48 1..47 1..46 1..46 1,.45 1.,44 1,.43 1,.41 1..41 1,.40
16 1..42 1..51 1,.51 1,.50 1,.48 1,.47 1..46 1,.45 1.,44 1,.44 1,.43 1,.41 1,.40 1,.39 1,.38
17 1.,42 1.,51 1,.50 1.,49 1..47 1,.46 1..45 1,.44 1.,43 1,.43 1..41 1,.40 1,.39 1,.38 1,.37
18 1..41 1..50 1,.49 1,.48 1,.46 1 .45 1.,44 1,.43 1.,42 1 .42 1,,40 1,.39 1,.38 1,.37 1,.36
19 1..41 1,.49 1,.49 1,.47 1,.46 1,.44 1.,43 1,.42 1.,41 1,.41 1.,40 1,.38 1,.37 1,.36 1,.35
20 1,.40 1..49 1,.48 1,.47 1 .45 1,.44 1..43 1,.42 1..41 1,.40 1..39 1 .37 1,.36 1 .35 1,.34
21 1.,40 1,.48 1,.48 1,.46 1,,44 1,.43 1.,42 1,.41 1..40 1,.39 1.,38 1,.37 1,.35 1,.34 1,.33
22 1..40 1..48 1 .47 1 .45 1,.44 1 .42 1,.41 1,.40 1,.39 1,.39 1..37 1 .36 1,.34 1,.33 1,.32
23 1,.39 1,.47 1,.47 1,.45 1,.43 1,.42 1,.41 1,.40 1,.39 1 .38 1,.37 1,.35 1 .34 1 .33 1,.32
24 1,.39 1..47 1,.46 1,.44 1,.43 1,.41 1..40 1,.39 1..38 1,.38 1,.36 1,.35 1,,33 1,.32 1,.31
25 1.,39 1..47 1.,46 1..44 1,.42 1,.41 1..40 1,.39 1.,38 1,.37 1.,36 1,.34 1,.33 1..32 1,.31
26 1..38 1..46 1,.45 1,.44 1,.42 1 .41 1..39 1,.38 1.,37 1,.37 1..35 1,.34 1,.32 1,.31 1,.30
27 1..38 1,.46 1,.45 1,.43 1,.42 1,.40 1.,39 1,.38 1,.37 1,.36 1,.35 1,.33 1,.32 1,.31 1,.30
28 1,.38 1,.46 1,.45 1,.43 1 .41 1,.40 1..39 1,.38 1,.37 1 .36 1,.34 1 .33 1,.31 1,.30 1 .29
29 1..38 1,.45 1 .45 1,.43 1,.41 1 .40 1,.38 1,.37 1.,36 1,.35 1,.34 1,.32 1,.31 1,.30 1,.29
30 1..38 1,.45 1 .44 1 .42 1 .41 1,.39 1,.38 1,.37 1.,36 1,.35 1,.34 1 .32 1 .30 1,.29 1,.28
40 1,.36 1,.44 1,.42 1 .40 1,.39 1,.37 1,.36 1,.35 1.,34 1,.33 1.,31 1 .30 1,.28 1..26 1,.25
60 1..35 1,.42 1,.41 1,.38 1,.37 1,.35 1,.33 1,.32 1.,31 1,.30 1,,29 1,.27 1,.25 1,.24 1,.22
120 1,.34 1..40 1,.39 1,.37 1,.35 1,.33 1.,31 1,.30 1.,29 1,.28 1.. f 6 1,.24 1,.22 1..21 1,.19
00 1..32 1..39 1,.37 1,.35 1,.33 1,.31 1.,29 1,.28 1.,27 1..25 1.,24 1,.22 1,.19 1..18 1,.16
Appendix 5 (Continued)
F-Distribution
Upper 10%
V1 1 2 3 4 5 6 7 8 9 10 12 15 20 24 30 40 60 120
V2 ____________
1 3 9 .8 6 4 9 .5 0 5 3 .5 9 5 5 .8 3 5 7 .2 4 5 8 .2 0 58..91 59.,44 59.,86 60.,19 60.,71 61.,22 61. 74 62.,00 62.,26 62.,53 62. 79 63.,06 63.,33
2 8 .5 3 9 .0 0 9 .1 6 9 .2 4 9 .2 9 9 .3 3 9.,35 9.,37 9.,38 9.,39 9.,41 9.,42 9. 44 9. 45 9.,46 9.,47 9. 47 9.,48 9.,49
3 5 .5 4 5 .4 6 5 .3 9 5 .3 4 5 .3 1 5 .2 8 5.,27 5.,25 5.,24 5.,23 5.,22 5.,20 5. 18 5. 18 5., 17 5.,16 5. 15 5.,14 5. 13
4 4 .5 4 4 .3 2 4 .1 9 4 .1 1 4 .0 5 4 .0 1 3.,98 3.,95 3.,94 3.,92 3..90 3.,87 3. 84 3. 83 3.,82 3.,80 3.,79 3.,78 3. 76
5 4 .0 6 3 .7 8 3 .6 2 3 .5 2 3 .4 5 3 .4 0 3.,37 3.,34 3..32 3.,30 3.,27 3,,24 3.,21 3.,19 3.,17 3.,16 3.,14 3. 12 3. 10
6 3 .7 8 3 .4 6 3 .2 9 3 .1 8 3 .1 1 3 .0 5 3.,01 2 .,98 2 .,96 2 .,94 2 ..90 2 .,87 2 .,84 2 .,82 2 .,80 2 .,78 2 ..76 2 .,74 2 .,72
7 3 .5 9 3 .2 6 3 .0 7 2 .9 6 2.88 2 .8 3 2 .,78 2 .,75 2 .,72 2 .,70 2 .,67 2 . 63 2 .,59 2 .,58 2 ..56 2 .,54 2 .,51 2 .,49 2 .,47
8 3 .4 6 3 .1 1 2 .9 2 2 .8 1 2 .7 3 2 .6 7 2 .,62 2 . 59 2 .,56 2 .,54 2 . 50 2 , 46 2 .,42 2 . 40 2 .,38 2 .,36 2 .,34 2 . 32 2 .,29
9 3 .3 6 3 .0 1 2 .8 1 2 .6 9 2 .6 1 2 .5 5 2 ..51 2 .,47 2 ..44 2 .,42 2 .,38 2 , 34 2 .,30 2 .,28 2 .,25 2 , 23 2 .,21 2 .,18 2 .,16
10 3 .2 9 2 .9 2 2 .7 3 2 .6 1 2 .5 2 2 .4 6 2 .,41 2 .,38 2 .,35 2 ..32 2 .,28 2 ..24 2 ..20 2 ,,18 2 ., 16 2 . 13 2 .,11 2 .,08 2 .,06
11 3 .2 3 2.86 2.66 2 .5 4 2 .4 5 2 .3 9 2 .,34 2 ..30 2 .,27 2 .,25 2 . 21 2 ..17 2 .,12 2 .,10 2 ,,08 2 .,05 2 . 03 2 .,00 1.,97
12 3 .1 8 2 .8 1 2 .6 1 2 .4 8 2 .3 9 2 .3 3 2 .,28 2 .,24 2 ,. 21 2 .,19 2 .,15 2 ,. 10 2 .,06 2 .,04 2 .,01 1..99 1.,96 1.,93 1.,90
13 3 .1 4 2 .7 6 2 .5 6 2 .4 3 2 .3 5 2 .2 8 2 ..23 2 .,20 2 .,16 2 .,14 2 .,10 2 ., 05 2 . 01 1.,98 1.,9 6 1.,93 1.,90 1.,88 1 .,85
14 3 .1 0 2 .7 3 2 .5 2 2 .3 9 2 .3 1 2 .2 4 2 .,19 2 .,15 2 .,12 2 .,10 2 ., 05 2 ,. 01 1.,96 1.,94 1 ,.91 1.,89 1.,86 1.,83 1 .,80
15 3 .0 7 2 .7 0 2 .4 9 2 .3 6 2 .2 7 2 .2 1 2 .,16 2 .,12 2 ,,09 2 ..0 6 2 ., 02 1,.97 1.,92 1.,90 1.,87 1,.85 1..82 1,.79 1..76
16 3 .0 5 2 .6 7 2 .4 6 2 .3 3 2 .2 4 2 .1 8 2 .,13 2 .,09 2 ,. 06 2 ,. 03 1..99 1,.94 1.,89 1.,87 1.,84 1,.81 1,,78 1.,75 1,.72
17 3 .0 3 2 .6 4 2 .4 4 2 .3 1 2 .2 2 2 .1 5 2 ..10 2 .,06 2 ,. 03 2 ,. 00 1..96 1.,91 1.,86 1,.84 1,.81 1,.78 1,.75 1.,72 1..69
18 3 .0 1 2 .6 2 2 .4 2 2 .2 9 2.20 2 .1 3 2 .,08 2 ..04 2 ,. 00 1.,98 1..93 1.,89 1,,84 1.,81 1,.78 1,.75 1,.72 1.,69 1,.66
19 2 .9 9 2 .6 1 2 .4 0 2 .2 7 2 .1 8 2 .1 1 2 ..06 2 .,02 1.,98 1 ,. 9 6 1..91 1,.86 1.,81 1.,79 1,.7 6 1,.73 1,,70 1..67 1,.63
20 2 .9 7 2 .5 9 2 .3 8 2 .2 5 2 .1 6 2 .0 9 2 .,04 2 .,00 1.,96 1.,94 1 ,,89 1,.84 1.,79 1.,77 1,.7 4 1,.71 1,.68 1,.64 1,.61
21 2 .9 6 2 .5 7 2 .3 6 2 .2 3 2 .1 4 2 .0 8 2 .,02 1.,98 1..95 1.,92 1,.87 1,.8 3 1.,78 1,.75 1,.72 1,.69 1..66 1,.62 1,.59
22 2 .9 5 2 .5 6 2 .3 5 2 .2 2 2 .1 3 2 .0 6 2 .,01 1.,97 1..93 1,.9 0 1,.86 1,.81 1.,76 1..73 1,.7 0 1,.67 1,.6 4 1,.60 1,.57
23 2 .9 4 2 .5 5 2 .3 4 2 .2 1 2 .1 1 2 .0 5 1.,9 9 1.,95 1.,92 1,.8 9 1..84 1,.8 0 1..74 1.,72 1,.6 9 1,.66 1..62 1..59 1,.55
24 2 .9 3 2 .5 4 2 .3 3 2 .1 9 2.10 2 .0 4 1.,98 1.,94 1.,91 1.,88 1,.83 1,.78 1,,73 1,.70 1,.6 7 1,.64 1,.61 1.,57 1,.53
25 2 .9 2 2 .5 3 2 .3 2 2 .1 8 2 .0 9 2 .0 2 1.,97 1 ,,93 1..89 1.,8 7 1,.82 1,.7 7 1,.72 1..69 1,.66 1.6 3 1,.5 9 1.56 1 .52
26 2 .9 1 2 .5 2 2 .3 1 2 .1 7 2 .0 8 2 .0 1 1.,9 6 1.,92 1.,88 1.,86 1.,81 1,.7 6 1,.71 1,.68 1 ,. 65 1,.61 1,.58 1,.54 1,.5 0
27 2 .9 0 2 .5 1 2 .3 0 2 .1 7 2 .0 7 2.00 1..95 1.,91 1..87 1.,85 1,.8 0 1 ,. 75 1,.70 1,,67 1,.64 1,.6 0 1,.5 7 1,.53 1,.49
28 2 .8 9 2 .5 0 2 .2 9 2 .1 6 2 .0 6 2.00 1..94 1.,90 1..87 1.,8 4 1,.79 1 ,. 7 4 1,.69 1..66 1..6 3 1,.5 9 1,.5 6 1,.52 1,.4 8
29 2 .8 9 2 .5 0 2 .2 8 2 .1 5 2 .0 6 1 .9 9 1..9 3 1.,89 1,.86 1.,8 3 1..78 1,.7 3 1,,68 1.,65 1,.62 1,.5 8 1,,55 1,.51 1,.4 7
30 2.88 2 .4 9 2 .2 8 2 .1 4 2 .0 5 1 .9 8 1..93 1.,88 1..85 1.,82 1..77 1,.72 1..67 1,.64 1,.61 1,.5 7 1,.5 4 1,.50 1,.4 6
40 2 .8 4 2 .4 4 2 .2 3 2 .0 9 2.00 1 .9 3 1.,87 1.,83 1,.79 1.,7 6 1,.71 1,.66 1,.61 1,.57 1.5 4 1,.51 1,.4 7 1.42 1.38
60 2 .7 9 2 .3 9 2 .1 8 2 .0 4 1 .9 5 1 .8 7 1.,82 1.,77 1..74 1..71 1,.66 1,.6 0 1,.54 1,.51 1,.4 8 1,.4 4 1,.4 0 1,.35 1,.2 9
120 2 .7 5 2 .3 5 2 .1 3 1 .9 9 1 .9 0 1 .8 2 1.,77 1.,72 1..68 1..65 1,.60 1,.55 1..48 1,.45 1,.41 1.3 7 1,.32 1.26 1.1 9
00 2 .7 1 2 .3 0 2 .0 8 1 .9 4 1 .8 5 1 .7 7 1..72 1.,67 1..63 1.,6 0 1,.55 1,.4 9 1.,42 1.,38 1,.3 4 1,.3 0 1,.2 4 1.17 1.00
Appendix 5 (Continued)
F-Distribution
Upper 5%
V2
l 8 10 12 15 20 24 30 40 60 120
1 1 61. 200. 216. 225. 230 . 234. 237 . 239 . 241. 242. 2 44 . 246. 248. 2 49 . 2 50 . 2 51 . 252. 2 53 . 254 .
2 1 8 .5 1 1 9 .0 0 1 9 .1 6 1 9 .2 5 1 9 .3 0 1 9 .3 3 1 9 .3 5 1 9 .3 7 1 9 .3 8 1 9 .4 0 19.,41 19.,43 1 9 .4 5 19., 45 19.,46 19,.47 19..48 19.,49 19..50
3 1 0 .1 3 9 .5 5 9 .2 8 9 .1 2 9 .0 1 8 .9 4 8 .8 9 8 .8 5 8 .8 1 8 .7 9 8 .,74 8 .,70 8 .6 6 8 ..64 8 ..62 8 ,. 5 9 8 ..57 8 .,55 8 ,. 53
4 7 .7 1 6 .9 4 6 .5 9 6 .3 9 6 .2 6 6 .1 6 6 .0 9 6 .0 4 6 .0 0 5 .9 6 5.,91 5..86 5 .8 0 5,.77 5,.75 5,.72 5,.69 5..66 5,.6 3
5 6 .6 1 5 .7 9 5 .4 1 5 .1 9 5 .0 5 4 .9 5 4 .8 8 4 .8 2 4 .7 7 4 .7 4 4.,68 4..62 4 .5 6 4..53 4.,50 4..4 6 4..4 3 4..40 4,.3 6
6 5 .9 9 5 .1 4 4 .7 6 4 .5 3 4 .3 9 4 .2 8 4 .2 1 4 .1 5 4 .1 0 4 .0 6 4..00 3..94 3 .8 7 3..84 3.,81 3..77 3,.74 3.,70 3,.67
7 5 .5 9 4 .7 4 4 .3 5 4 .1 2 3 .9 7 3 .8 7 3 .7 9 3 .7 3 3 .6 8 3 .6 4 3.,57 3..51 3 .4 4 3..41 3..38 3,.34 3,.3 0 3..27 3..2 3
8 5 .3 2 4 .4 6 4 .0 7 3 .8 4 3 .6 9 3 .5 8 3 .5 0 3 .4 4 3 .3 9 3 .3 5 3..28 3..22 3 .1 5 3,.12 3,.08 3 .04 3,.01 2 .,97 2 ,. 9 3
9 5 .1 2 4 .2 6 3 .8 6 3 .6 3 3 .4 8 3 .3 7 3 .2 9 3 .2 3 3 .1 8 3 .1 4 3.,07 3,.01 2 .9 4 2 ..90 2 . 86 2 ,. 8 3 2 ,. 79 2 ..75 2 .71
10 4 .9 6 4 .1 0 3 .7 1 3 .4 8 3 .3 3 3 .2 2 3 .1 4 3 .0 7 3 .0 2 2 .9 8 2 ..91 2 ,.85 2 .7 7 2 ,. 74 2 ,. 70 2 ,. 66 2 ,. 62 2 ..58 2 ,.54
11 4 .8 4 3 .9 8 3 .5 9 3 .3 6 3 .2 0 3 .0 9 3 .0 1 2 .9 5 2 .9 0 2 .8 5 2 ..79 2 .,72 2 .6 5 2 ,. 61 2 ,. 57 2 ,. 53 2 ,. 49 2 ..45 2 .40
12 4 .7 5 3 .8 9 3 .4 9 3 .2 6 3 .1 1 3 .0 0 2 .9 1 2 .8 5 2 .8 0 2 .7 5 2 ..69 2 ,. 62 2 .5 4 2 ,.51 2 . 47 2 ,. 4 3 2 ,. 38 2 ..34 2 .30
13 4 .6 7 3 .8 1 3 .4 1 3 .1 8 3 .0 3 2 .9 2 2 .8 3 2 .7 7 2 .7 1 2 .6 7 2, .60 2 ,. 53 2 .4 6 2 ., 42 2 ..38 2 .34 2 ,. 3 0 2 ,.25 2 ,. 21
14 4 .6 0 3 .7 4 3 .3 4 3 .1 1 2 .9 6 2 .8 5 2 .7 6 2 .7 0 2 .6 5 2 .6 0 2 .,53 2 ..4 6 2 .3 9 2 ..35 2 ..31 2 ,. 2 7 2 ,. 22 2 .,18 2 ,. 1 3
15 4 .5 4 3 .6 8 3 .2 9 3 .0 6 2 .9 0 2 .7 9 2 .7 1 2 .6 4 2 .5 9 2 .5 4 2 .,48 2 .,4 0 2 .3 3 2 .,29 2 ,. 25 2 .,20 2 ..1 6 2 .,11 2 . 07
16 4 .4 9 3 .6 3 3 .2 4 3 .0 1 2 .8 5 2 .7 4 2.66 2 .5 9 2 .5 4 2 .4 9 2 .,42 2 ,. 35 2 .2 8 2 ..24 2 ..19 2 ,. 15 2 ,. 11 2 .,06 2 . 01
17 4 .4 5 3 .5 9 3 .2 0 2 .9 6 2 .8 1 2 .7 0 2 .6 1 2 .5 5 2 .4 9 2 .4 5 2 . 38 2 ..31 2 .2 3 2 .,19 2 .,15 2 ,. 10 2 ,. 0 6 2 ..01 1 ,. 9 6
18 4 .4 1 3 .5 5 3 .1 6 2 .9 3 2 .7 7 2.66 2 .5 8 2 .5 1 2 .4 6 2 .4 1 2 .,34 2 ..27 2 .1 9 2 ,. 15 2 , 11 2 ,. 06 2 ,. 02 1 ..97 1 .92
19 4 .3 8 3 .5 2 3 .1 3 2 .9 0 2 .7 4 2 .6 3 2 .5 4 2 .4 8 2 .4 2 2 .3 8 2 .,31 2 .,23 2 .1 6 2 ..11 2 .,07 2 ,. 0 3 1 ..98 1 .,93 1 ..88
20 4 .3 5 3 .4 9 3 .1 0 2 .8 7 2 .7 1 2 .6 0 2 .5 1 2 .4 5 2 .3 9 2 .3 5 2 . 28 2 ,. 20 2.12 2 . 08 2 ,. 04 1 ,. 9 9 1 .95 1 ..90 1 .84
21 4 .3 2 3 .4 7 3 .0 7 2 .8 4 2.68 2 .5 7 2 .4 9 2 .4 2 2 .3 7 2 .3 2 2 . 25 2 ..18 2.10 2 ,. 05 2 .,01 1 .9 6 1 .92 1 ..87 1 ,.81
22 4 .3 0 3 .4 4 3 .0 5 2 .8 2 2.66 2 .5 5 2 .4 6 2 .4 0 2 .3 4 2 .3 0 2 ,. 23 2 ,. 15 2 .0 7 2 ,. 03 1 ,. 98 1 .9 4 1 .8 9 1 ,.84 1 .7 8
23 4 .2 8 3 .4 2 3 .0 3 2 .8 0 2 .6 4 2 .5 3 2 .4 4 2 .3 7 2 .3 2 2 .2 7 2 ,. 20 2 ,. 1 3 2 .0 5 2 ,. 01 1 ,. 96 1 .91 1 ,. 86 1 ..81 1 .7 6
24 4 .2 6 3 .4 0 3 .0 1 2 .7 8 2 .6 2 2 .5 1 2 .4 2 2 .3 6 2 .3 0 2 .2 5 2 ;.18 2 ,. 11 2 .0 3 1 ,.98 1 ,. 94 1 .8 9 1 .84 1 ,. 79 1 . 73
25 4 .2 4 3 .3 9 2 .9 9 2 .7 6 2 .6 0 2 .4 9 2 .4 0 2 .3 4 2 .2 8 2 .2 4 2 ,. 16 2 ..0 9 2 .0 1 1 ,. 9 6 1 ..92 1 .87 1 .82 1 ,.77 1 .71
26 4 .2 3 3 .3 7 2 .9 8 2 .7 4 2 .5 9 2 .4 7 2 .3 9 2 .3 2 2 .2 7 2.22 2 ..15 2 ,.0 7 1 .9 9 1 ,. 95 1 .,90 1 ,. 85 1 ,. 8 0 1 .,75 1 ,. 6 9
27 4 .2 1 3 .3 5 2 .9 6 2 .7 3 2 .5 7 2 .4 6 2 .3 7 2 .3 1 2 .2 5 2.20 2 ,. 13 2 .,0 6 1 .9 7 1 ,. 93 1 . 88 1 ,. 84 1 ,. 7 9 1 .,73 1 ,. 6 7
28 4 .2 0 3 .3 4 2 .9 5 2 .7 1 2 .5 6 2 .4 5 2 .3 6 2 .2 9 2 .2 4 2 .1 9 2 ..12 2 ,.04 1 .9 6 1 . 91 1 ..87 1 .82 1 ,. 77 1 ..71 1 .65
29 4 .1 8 3 .3 3 2 .9 3 2 .7 0 2 .5 5 2 .4 3 2 .3 5 2 .2 8 2 .2 2 2 .1 8 2, .10 2 ,. 0 3 1 .9 4 1 ,. 9 0 1 ,. 85 1 .81 1 .7 5 1 ,. 7 0 1 .64
30 4 .1 7 3 .3 2 2 .9 2 2 .6 9 2 .5 3 2 .4 2 2 .3 3 2 .2 7 2 .2 1 2 .1 6 2. ,09 2 ..01 1 .9 3 1 ,. 8 9 1 .,84 1 ,. 7 9 1 ,. 7 4 1 ..68 1 .62
40 4 .0 8 3 .2 3 2 .8 4 2 .6 1 2 .4 5 2 .3 4 2 .2 5 2 .1 8 2 .1 2 2 .0 8 2, .00 1 ,.92 1 .8 4 1 ,. 79 1 ,. 74 1 ,. 6 9 1 ,.64 1 ..5 8 1 ,.51
60 4 .0 0 3 .1 5 2 .7 6 2 .5 3 2 .3 7 2 .2 5 2 .1 7 2.10 2 .0 4 1 .9 9 1 .,92 1 ..8 4 1 .7 5 1 .,7 0 1 ..65 1 ,. 5 9 1 ,. 5 3 1 .,47 1 ,. 3 9
120 3 .9 2 3 .0 7 2.68 2 .4 5 2 .2 9 2 .1 7 2 .0 9 2.02 1 .9 6 1 .9 1 1 .,83 1 .,75 1.66 1 ..61 1 .,55 1 ,. 5 0 1 ,. 4 3 1 ..35 1 ,.25
3 .8 4 3 .0 0 2 .6 0 2 .3 7 2.21 2.10 2 .0 1 1 .9 4 1 .8 8 1 .8 3 1 .,75 1 .,67 1 .5 7 1 .,52 1 .,46 1 ,. 3 9 1 ,.32 1 .,22 1 ..00
Appendix 5 (Continued)
F-Distribution
Upper 1%
V1 r
1 2 3 4 5 6 8 9 10 12 15 20 24 30 40 60 120 00
V2
1 4 052. 5 00 0 . 5 40 3 . 5625. 5 764. 5859. 5 928. 5 982. 6 02 2 . 6056. 6106. 615 7 . 6209. 6235. 6261. 6287. 631 3 . 633 9 . 6366.
2 9 8 .5 0 9 9 .0 0 9 9 .1 7 9 9 .2 5 9 9 .3 0 9 9 .3 3 99. 36 9 9 .3 7 9 9 .3 9 9 9 .4 0 9 9 .4 2 9 9 .4 3 9 9 .4 5 9 9 .4 6 9 9 .4 7 9 9 .4 7 9 9 .4 8 9 9 .4 9 99.,50
3 3 4 .1 2 3 0 .8 2 2 9 .4 6 2 8 .7 1 2 8 .2 4 2 7 .9 1 27. 67 2 7 .4 9 2 7 .3 5 2 7 .2 3 2 7 .0 5 2 6 .8 7 2 6 .6 9 2 6 .6 0 2 6 .5 0 2 6 .4 1 2 6 .3 2 2 6 .2 2 26.,13
4 21.20 1 8 .0 0 1 6 .6 9 1 5 .9 8 1 5 .5 2 1 5.2 1 14. 98 1 4 .8 0 1 4 .6 6 1 4 .5 5 1 4 .3 7 1 4 .2 0 1 4 .0 2 1 3 .9 3 1 3 .8 4 1 3 .7 5 1 3 .6 5 1 3 .5 6 13..46
5 1 6 .2 6 1 3 .2 7 1 2 .0 6 1 1 .3 9 1 0 .9 7 1 0 .6 7 10 .,46 1 0 .2 9 1 0 .1 6 1 0 .0 5 9 .8 9 9 .7 2 9 .5 5 9 .4 7 9 .3 8 9 .2 9 9 .2 0 9 .1 1 9.,02
6 1 3 .7 5 1 0 .9 2 9 .7 8 9 .1 5 8 .7 5 8 .4 7 8 .,26 8 . 1 0 7 .9 8 7 .8 7 7 .7 2 7 .5 6 7 .4 0 7 .3 1 7 .2 3 7 .1 4 7 .0 6 6 .9 7 6 .,88
7 1 2 .2 5 9 .5 5 8 .4 5 7 .8 5 7 .4 6 7 .1 9 6 . 99 6 .8 4 6 .7 2 6 .6 2 6 .4 7 6 .3 1 6 .1 6 6 .0 7 5 .9 9 5 .9 1 5 .8 2 5 .7 4 5..65
8 1 1 .2 6 8 .6 5 7 .5 9 7 .0 1 6 .6 3 6 .3 7 6 .,18 6 .0 3 5 .9 1 5 .8 1 5 .6 7 5 .5 2 5 .3 6 5 .2 8 5 .2 0 5 .1 2 5 .0 3 4 .9 5 4.,86
9 1 0 .5 6 8 .0 2 6 .9 9 6 .4 2 6 .0 6 5 .8 0 5. 61 5 .4 7 5 .3 5 5 .2 6 5 .1 1 4 .9 6 4 .8 1 4 .7 3 4 .6 5 4 .5 7 4 .4 8 4 .4 0 4..31
10 1 0 .0 4 7 .5 6 6 .5 5 5 .9 9 5 .6 4 5 .3 9 5. 20 5 .0 6 4 .9 4 4 .8 5 4 .7 1 4 .5 6 4 .4 1 4 .3 3 4 .2 5 4 .1 7 4 .0 8 4 .0 0 3.,91
11 9 .6 5 7 .2 1 6 .22 5 .6 7 5 .3 2 5 .0 7 4.,89 4 .7 4 4 .6 3 4 .5 4 4 .4 0 4 .2 5 4 .1 0 4 .0 2 3 .9 4 3 .8 6 3 .7 8 3 .6 9 3,.60
12 9 .3 3 6 .9 3 5 .9 5 5 .4 1 5 .0 6 4 .8 2 4.,64 4 .5 0 4 .3 9 4 .3 0 4 .1 6 4 .0 1 3 .8 6 3 .7 8 3 .7 0 3 .6 2 3 .5 4 3 .4 5 3.,36
13 9 .0 7 6 .7 0 5 .7 4 5 .2 1 4 .8 6 4 .6 2 4.,44 4 .3 0 4 .1 9 4 .1 0 3 .9 6 3 .8 2 3 .6 6 3 .5 9 3 .5 1 3 .4 3 3 .3 4 3 .2 5 3.,17
14 8 .86 6 .5 1 5 .5 6 5 .0 4 4 .6 9 4 .4 6 4.,28 4 .1 4 4 .0 3 3 .9 4 3 .8 0 3 .6 6 3 .5 1 3 .4 3 3 .3 5 3 .2 7 3 .1 8 3 .0 9 3,.00
15 8 .68 6 .3 6 5 .4 2 4 .8 9 4 .5 6 4 .3 2 4.,14 4 .0 0 3 .8 9 3 .8 0 3 .6 7 3 .5 2 3 .3 7 3 .2 9 3 .2 1 3 .1 3 3 .0 5
395
2 .9 6 2 ,. 87
16 8 .5 3 6 .2 3 5 .2 9 4 .7 7 4 .4 4 4 .2 0 4.,03 3 .8 9 3 .7 8 3 .6 9 3 .5 5 3 .4 1 3 .2 6 3 .1 8 3 .1 0 3 .0 2 2 .9 3 2 .8 4 2 ,. 75
17 8 .4 0 6.11 5 .1 8 4 .6 7 4 .3 4 4 .1 0 3. 93 3 .7 9 3 .6 8 3 .5 9 3 .4 6 3 .3 1 3 .1 6 3 .0 8 3 .0 0 2 .9 2 2 .8 3 2 .7 5 2 ..65
18 8 .2 9 6.01 5 .0 9 4 .5 8 4 .2 5 4 .0 1 3.,84 3 .7 1 3 .6 0 3 .5 1 3 .3 7 3 .2 3 3 .0 8 3 .0 0 2 .9 2 2 .8 4 2 .7 5 2.66 2 ,. 57
19 8 .1 8 5 .9 3 5 .0 1 4 .5 0 4 .1 7 3 .9 4 3.,77 3 .6 3 3 .5 2 3 .4 3 3 .3 0 3 .1 5 3 .0 0 2 .9 2 2 .8 4 2 .7 6 2 .6 7 2 .5 8 2 , 49
20 8.10 5 .8 5 4 .9 4 4 .4 3 4 .1 0 3 .8 7 3.,70 3 .5 6 3 .4 6 3 .3 7 3 .2 3 3 .0 9 2 .9 4 2.86 2 .7 8 2 .6 9 2 .6 1 2 .5 2 2 ,. 42
21 8.02 5 .7 8 4 .8 7 4 .3 7 4 .0 4 3 .8 1 3.,64 3 .51 3 .4 0 3 .3 1 3 .1 7 3 .0 3 2.88 2 .8 0 2 .7 2 2 .6 4 2 .5 5 2 .4 6 2 ,,36
22 7 .9 5 5 .7 2 4 .8 2 4 .3 1 3 .9 9 3 .7 6 3.,59 3 .4 5 3 .3 5 3 .2 6 3 .1 2 2 .9 8 2 .8 3 2 .7 5 2 .6 7 2 .5 8 2 .5 0 2 .4 0 2 ..31
23 7 .8 8 5 .6 6 4 .7 6 4 .2 6 3 .9 4 3 .7 1 3.,54 3 .4 1 3 .3 0 3 .2 1 3 .0 7 2 .9 3 2 .7 8 2 .7 0 2 .6 2 2 .5 4 2 .4 5 2 .3 5 2 ,. 26
24 7 .8 2 5 .6 1 4 .7 2 4 .2 2 3 .9 0 3 .6 7 3.,50 3 .3 6 3 .2 6 3 .1 7 3 .0 3 2 .8 9 2 .7 4 2.66 2 .5 8 2 .4 9 2 .4 0 2 .3 1 2 .,21
25 7 .7 7 5 .5 7 4 .6 8 4 .1 8 3 .8 5 3 .6 3 3. 46 3 .3 2 3 .2 2 3 .1 3 2 .9 9 2 .8 5 2 .7 0 2 .6 2 2 .5 4 2 .4 5 2 .3 6 2 .2 7 2 .,17
26 7 .7 2 5 .5 3 4 .6 4 4 .1 4 3 .8 2 3 .5 9 3.,42 3 .2 9 2 .1 8 2 .0 9 2 .9 6 2 .8 1 2 .66 2 .5 8 2 .5 0 2 .4 2 2 .3 3 2 .2 3 2 ,. 13
27 7 .6 8 5 .4 9 4 .6 0 4 .1 1 3 .7 8 3 .5 6 3.,39 3 .2 6 3 .1 5 3 .0 6 2 .9 3 2 .7 8 2 .6 3 2 .5 5 2 .4 7 2 .3 8 2 .2 9 2 .2 0 2 ., 10
28 7 .6 4 5 .4 5 4 .5 7 4 .0 7 3 .7 5 3 .5 3 3. 36 3 .2 3 3 .1 2 3 .0 3 2 .9 0 2 .7 5 2 .6 0 2 .5 2 2 .4 4 2 .3 5 2 .2 6 2 .1 7 2 .,06
29 7 .6 0 5 .4 2 4 .5 4 4 .0 4 3 .7 3 3 .5 0 3. 33 3 .2 0 3 .0 9 3 .0 0 2 .8 7 2 .7 3 2 .5 7 2 .4 9 2 .4 1 2 .3 3 2 .2 3 2 .1 4 2 ,. 03
30 7 .5 6 5 .3 9 4 .5 1 4 .0 2 3 .7 0 3 .4 7 3. 30 3 .1 7 3 .0 7 2 .9 8 2 .8 4 2 .7 0 2 .5 5 2 .4 7 2 .3 9 2 .3 0 2 .2 1 2.1 1 2 ,. 01
40 7 .3 1 5 .1 8 4 .3 1 3 .8 3 3 .5 1 3 .2 9 3.,12 2 .9 9 2 .8 9 2 .8 0 2.66 2 .5 2 2 .3 7 2 .2 9 2.20 2 .1 1 2.02 1 .9 2 1 ,. 80
60 7 .0 8 4 .9 8 4 .1 3 3 .6 5 3 .3 4 3 .1 2 2 . 95 2 .8 2 2 .7 2 2 .6 3 2 .5 0 2 .3 5 2.20 2.12 2 .0 3 1 .9 4 1 .8 4 1 .7 3 1 ,. 60
120 6 .8 5 4 .7 9 3 .9 5 3 .4 8 3 .1 7 2 .9 6 2 .,79 2.66 2 .5 6 2 .4 7 2 .3 4 2 .1 9 2 .0 3 1 .9 5 1.86 1 .7 6 1.66 1 .5 3 1 ..38
00 6 .6 3 4 .6 1 3 .7 8 3 .3 2 3 .0 2 2 .8 0 2 .,64 2 .5 1 2 .4 1 2 .3 2 2 .1 8 2 .0 4 1 .88 1 .7 9 1 .7 0 1 .5 9 1 .4 7 1 .3 2 1 ..00
Appendix 6
Upper Percentage Points of the Studentized Range
a = .10
1 3 4 5 6 7 8 9 10 12 L5
1 8 .9 2 9 1 3 .4 4 1 6 .3 6 1 8 .4 9 2 0 .1 5 2 1 .5 1 2 2 .6 4 2 3 .6 2 2 4 .4 8 2 5 .9 2 2 7 .6 2
2 4 .1 3 0 5 .7 3 3 6 .7 7 3 7 .5 3 8 8 .1 3 9 8 .6 3 3 9 .0 4 9 9 .4 0 9 9 .7 2 5 1 0 .2 6 1 0 .8 9
3 3 .3 2 8 4 .4 6 7 5 .1 9 9 5 .7 3 8 6 .1 6 2 6 .5 1 1 6 .8 0 6 7 .0 6 2 7 .2 8 7 7 .6 6 7 8.1 2 0
4 3 .0 1 5 3 .9 7 6 4 .5 8 6 5 .0 3 5 5 .3 8 8 5 .6 7 9 5 .9 2 6 6 .1 3 9 6 .3 2 7 6 .6 4 5 7 .0 2 5
5 2 .8 5 0 3 .7 1 7 4 .2 6 4 4 .6 6 4 4 .9 7 9 5 .2 3 8 5 .4 5 8 5 .6 4 8 5 .8 1 6 6.101 6 .4 4 0
6 2 .7 4 8 3 .5 5 9 4 .0 6 5 4 .4 3 5 4 .7 2 6 4 .9 6 6 5 .1 6 8 5 .3 4 4 5 .4 9 9 5 .7 6 2 6 .0 7 5
7 2 .6 8 0 3 .4 5 1 3 .9 3 1 4 .2 8 0 4 .5 5 5 4 .7 8 0 4 .9 7 2 5 .1 3 7 5 .2 8 3 5 .5 3 0 5 .8 2 6
8 2 .6 3 0 3 .3 7 4 3 .8 3 4 4 .1 6 9 4 .4 3 1 4 .6 4 6 4 .8 2 9 4 .9 8 7 5 .1 2 6 5 .3 6 2 5 .6 4 4
9 2 .5 9 2 3 .3 1 6 3 .7 6 1 4 .0 8 4 4 .3 3 7 4 .5 4 5 4 .7 2 1 4 .8 7 3 5 .0 0 7 5 .2 3 4 5 .5 0 6
10 2 .5 6 3 3 .2 7 0 3 .7 0 4 4 .0 1 8 4 .2 6 4 4 .4 6 5 4 .6 3 6 4 .7 8 3 4 .9 1 3 5 .1 3 4 5 .3 9 7
11 2 .5 4 0 3 .2 3 4 3 .6 5 8 3 .9 6 5 4 .2 0 5 4 .4 0 1 4 .5 6 8 4 .7 1 1 4 .8 3 8 5 .0 5 3 5 .3 0 9
12 2 .5 2 1 3 .2 0 4 3 .6 2 1 3 .9 2 2 4 .1 5 6 4 .3 4 9 4 .5 1 1 4 .6 5 2 4 .7 7 6 4 .9 8 6 5 .2 3 6
13 2 .5 0 5 3 .1 7 9 3 .5 8 9 3 .8 8 5 4 .1 1 6 4 .3 0 5 4 .4 6 4 4 .6 0 2 4 .7 2 4 4 .9 3 0 5 .1 7 6
14 2 .4 9 1 3 .1 5 8 3 .5 6 3 3 .8 5 4 4 .0 8 1 4 .2 6 7 4 .4 2 4 4 .5 6 0 4 .6 8 0 4 .8 8 2 5 .1 2 4
15 2 .4 7 9 3 .1 4 0 3 .5 4 0 3 .8 2 8 4 .0 5 2 4 .2 3 5 4 .3 9 0 4 .5 2 4 4 .6 4 1 4 .8 4 1 5 .0 7 9
16 2 .4 6 9 3 .1 2 4 3 .5 2 0 3 .8 0 4 4 .0 2 6 4 .2 0 7 4 .3 6 0 4 .4 9 2 4 .6 0 8 4 .8 0 5 5 .0 4 0
17 2 .4 6 0 3 .1 1 0 3 .5 0 3 3 .7 8 4 4 .0 0 4 4 .1 8 3 4 .3 3 4 4 .4 6 4 4 .5 7 9 4 .7 7 4 5 .0 0 5
18 2 .4 5 2 3 .0 9 8 3 .4 8 8 3 .7 6 7 3 .9 8 4 4 .1 6 1 4 .3 1 1 4 .4 4 0 4 .5 5 4 4 .7 4 6 4 .9 7 5
19 2 .4 4 5 3 .0 8 7 3 .4 7 4 3 .7 5 1 3 .9 6 6 4 .1 4 2 4 .2 9 0 4 .4 1 8 4 .5 3 1 4 .7 2 1 4 .9 4 8
20 2 .4 3 9 3 .0 7 8 3 .4 6 2 3 .7 3 6 3 .9 5 0 4 .1 2 4 4 .2 7 1 4 .3 9 8 4 .5 1 0 4 .6 9 9 4 .9 2 4
24 2 .4 2 0 3 .0 4 7 3 .4 2 3 3 .6 9 2 3 .9 0 0 4 .0 7 0 4 .2 1 3 4 .3 3 6 4 .4 4 5 4 .6 2 8 4 .8 4 7
30 2 .4 0 0 3 .0 1 7 3 .3 8 6 3 .6 4 8 3 .8 5 1 4 .0 1 6 4 .1 5 5 4 .2 7 5 4 .3 8 1 4 .5 5 9 4 .7 7 0
40 2 .3 8 1 2 .9 8 8 3 .3 4 9 3 .6 0 5 3 .8 0 3 3 .9 6 3 4 .0 9 9 4 .2 1 5 4 .3 1 7 4 .4 9 0 4 .6 9 5
60 2 .3 6 3 2 .9 5 9 3 .3 1 2 3 .5 6 2 3 .7 5 5 3 .9 1 1 4 .0 4 2 4 .1 5 5 4 .2 5 4 4 .4 2 1 4 .6 1 9
120 2 .3 4 4 2 .9 3 0 3 .2 7 6 3 .5 2 0 3 .7 0 7 3 .8 5 9 3 .9 8 7 4 .0 9 6 4 .1 9 1 4 .3 5 3 4 .5 4 3
2 .3 2 6 2 .9 0 2 3 .2 4 0 3 .4 7 8 3 .6 6 1 3 .8 0 8 3 .9 3 1 4 .0 3 7 4 .1 2 9 4 .2 8 5 4 .4 6 8
Appendix 6 (Cont'd .)
Upper Percentage Points of the Studentized Range
a = .05
2 3 4 5 6 7 8 9 10 12 15 30 60
1 1 7 .9 7 2 6 .9 8 3 2 .8 2 3 7 .0 8 4 0 .4 1 4 3 .1 2 4 5 .4 0 4 7 .3 6 4 9 .0 7 5 1 .9 6 5 5 .3 6 6 5 .1 5 7 3 .9 7
2 6 .0 8 5 8 .3 3 1 9 .7 9 8 10.88 1 1 .7 4 1 2 .4 4 1 3 .0 3 1 3 .5 4 1 3 .9 9 1 4 .7 5 1 5 .6 5 1 8 .2 7 20.66
3 4 .5 0 1 5 .9 1 0 6 .8 2 5 7 .5 0 2 8 .0 3 7 8 .4 7 8 8 .8 5 3 9 .1 7 7 9 .4 6 2 9 .9 4 6 1 0 .5 3 12.21 1 3 .7 6
4 3 .9 2 7 5 .0 4 0 5 .7 5 7 6 .2 8 7 6 .7 0 7 7 .0 5 3 7 .3 4 7 7 .6 0 2 7 .8 2 6 8 .2 0 8 8 .6 6 4 10.00 1 1 .2 4
5 3 .6 3 5 4 .6 0 2 5 .2 1 8 5 .6 7 3 6 .0 3 3 6 .3 3 0 6 .5 8 2 6 .8 0 2 6 .9 9 5 7 .3 2 4 7 .7 1 7 8 .8 7 5 9 .9 4 9
6 3 .4 6 1 4 .3 3 9 4 .8 9 6 5 .3 0 5 5 .6 2 8 5 .8 9 5 6.122 6 .3 1 9 6 .4 9 3 6 .7 8 9 7 .1 4 3 8 .1 8 9 9 .1 6 3
7 3 .3 4 4 4 .1 6 5 4 .6 8 1 5 .0 6 0 5 .3 5 9 5 .6 0 6 5 .8 1 5 5 .9 9 8 6 .1 5 8 6 .4 3 1 6 .7 5 9 7 .7 2 8 8 .6 3 2
8 3 .2 6 1 4 .0 4 1 4 .5 2 9 4 .8 8 6 5 .1 6 7 5 .3 9 9 5 .5 9 7 5 .7 6 7 5 .9 1 8 6 ,1 7 5 6 .4 8 3 7 .3 9 5 8 .2 4 8
9 3 .1 9 9 3 .9 4 9 4 .4 1 5 4 .7 5 6 5 .0 2 4 5 .2 4 4 5 .4 3 2 5 .5 9 5 5 .7 3 9 5 ,9 8 3 6 .2 7 6 7 .1 4 5 7 .9 5 8
10 3 .1 5 1 3 .8 7 7 4 .3 2 7 4 .6 5 4 4 .9 1 2 5 .1 2 4 5 .3 0 5 5 .4 6 1 5 .5 9 9 5 .8 3 3 6 .1 1 4 6 .9 4 8 7 .7 3 0
11 3 .1 1 3 3 .8 2 0 4 .2 5 6 4 .5 7 4 4 .8 2 3 5 .0 2 8 5 .2 0 2 5 .3 5 3 5 .4 8 7 5 .7 1 3 5 .9 8 4 6 .7 9 0 7 .5 4 6
12 3 .0 8 2 3 .7 7 3 4 .1 9 9 4 .5 0 8 4 .7 5 1 4 .9 5 0 5 .1 1 9 5 .2 6 5 5 .3 9 5 5 .6 1 5 5 .8 7 8 6 .6 6 0 7 .3 9 4
13 3 .0 5 5 3 .7 3 5 4 .1 5 1 4 .4 5 3 4 .6 9 0 4 .8 8 5 5 .0 4 9 5 .1 9 2 5 .3 1 8 5 .5 3 3 5 .7 8 9 6 .5 5 1 7 .2 6 7
14 3 .0 3 3 3 .7 0 2 4 .1 1 1 4 .4 0 7 4 .6 3 9 4 .8 2 9 4 .9 9 0 5 .1 3 1 5 .2 5 4 5 .4 6 3 5 .7 1 4 6 .4 5 9 7 .1 5 9
15 3 .0 1 4 3 .6 7 4 4 .0 7 6 4 .3 6 7 4 .5 9 5 4 .7 8 2 4 .9 4 0 5 .0 7 7 5 .1 9 8 5 .4 0 4 5 .6 4 9 6 .3 7 9 7 .0 6 5
16 2 .9 9 8 3 .6 4 9 4 .0 4 6 4 .3 3 3 4 .5 5 7 4 .7 4 1 4 .8 9 7 5 .0 3 1 5 .1 5 0 5 .3 5 2 5 .5 9 3 6 .3 1 0 6 .9 8 4
17 2 .9 8 4 3 .6 2 8 4 .0 2 0 4 .3 0 3 4 .5 2 4 4 .7 0 5 4 .8 5 8 4 .9 9 1 5 .1 0 8 5 .3 0 7 5 .5 4 4 6 .2 4 9 6 .9 1 2
18 2 .9 7 1 3 .6 0 9 3 .9 9 7 4 .2 7 7 4 .4 9 5 4 .6 7 3 4 .8 2 4 4 .9 5 6 5 .0 7 1 5 .2 6 7 5 .5 0 1 6 .1 9 5 6 .8 4 8
19 2 .9 6 0 3 .5 9 3 3 .9 7 7 4 .2 5 3 4 .4 6 9 4 .6 4 5 4 .7 9 4 4 .9 2 4 5 .0 3 8 5 .2 3 1 5 .4 6 2 6 .1 4 7 6 .7 9 2
20 2 .9 5 0 3 .5 7 8 3 .9 5 8 4 .2 3 2 4 .4 4 5 4 .6 2 0 4 .7 6 8 4 .8 9 6 5 .0 0 8 5 .1 9 9 5 .4 2 7 6 .1 0 4 6 .7 4 0
24 2 .9 1 9 3 .5 3 2 3 .9 0 1 4 .1 6 6 4 .3 7 3 4 .5 4 1 4 .6 8 4 4 .8 0 7 4 .9 1 5 5 .0 9 9 5 .3 1 9 5 .9 6 8 6 .5 7 9
30 2.888 3 .4 8 6 3 .8 4 5 4 .1 0 2 4 .3 0 2 4 .4 6 4 4 .6 0 2 4 .7 2 0 4 .8 2 4 5 .0 0 1 5 .2 1 1 5 .8 3 3 6 .4 1 7
40 2 .8 5 8 3 .4 4 2 3 .7 9 1 4 .0 3 9 4 .2 3 2 4 .3 8 9 4 .5 2 1 4 .6 3 5 4 .7 3 5 4 .9 0 4 5 .1 0 6 5 .7 0 0 6 .2 5 5
60 2 .8 2 9 3 .3 9 9 3 .7 3 7 3 .9 7 7 4 .1 6 3 4 .3 1 4 4 .4 4 1 4 .5 5 0 4 .6 4 6 4 .8 0 8 5 .0 0 1 5 .5 6 6 6 .0 9 3
120 2 .8 0 0 3 .3 5 6 3 .6 8 5 3 .9 1 7 4 .0 9 6 4 .2 4 1 4 .3 6 3 4 .4 6 8 4 .5 6 0 4 .7 1 4 4 .8 9 8 5 .4 3 4 5 .9 2 9
CO 2 .7 7 2 3 .3 1 4 3 .6 3 3 3 .8 5 8 4 .0 3 0 4 .1 7 0 4 .2 8 6 4 .8 3 7 4 .4 7 4 4 .6 2 2 4 .7 9 6 5 .3 0 1 5 .7 6 4
Appendix 6 (Cont1d .)
Upper Percentage Points of the Studentized Range
ols .01
2 3 4 5 6~ 7 8 9 10 12 15 30 60 100
1 90,.03 135..0 1 6 4 .3 1 8 5 .6 2 0 2 ..2 215.,8 2 27.,2 237..0 245..6 260..0 277..0 326 ..0 3 70 .,1 400..1
2 14,.04 19,.02 2 2 .2 9 2 4 .7 2 26,.63 28..20 29..53 30..68 31..69 3 3 ..4 0 3 5 ..43 41..32 4 6 .,70 50..38
3 8 ,. 261 10, . 6 2 1 2 .1 7 1 3 .3 3 14,.24 15,.00 15..64 16..20 16,.69 17,.53 18,.52 2 1 ,. 44 24..13 25,.99
4 6 ,. 512 8 . 120 9 .1 7 3 9 .9 5 8 1 0 ,. 58 1 1 ..10 1 1 ..55 1 1 ,. 93 1 2 ,. 27 1 2 ,.84 13,.53 15,.57 17..46 18,.77
5 5,.702 6 ,. 9 7 6 7 .8 0 4 8 .4 2 1 8 ,.913 9..321 9.,669 9,.972 1 0 ,. 24 1 0 ,. 7 0 1 1 ,. 24 12 .87 14..39 15 .45
6 5,.243 6 ..331 7 .0 3 3 7 .5 5 6 7,.973 8 ..318 8 ..613 8 ..8 6 9 9,.097 9,.485 9,.951 1 1 ..34 1 2 .,65 13,.55
7 4,.949 5..9 1 9 6 .5 4 3 7 .0 0 5 7..373 7..6 7 9 7..939 8 ..166 8 ,. 368 8 ,. 711 9,.1 2 4 1 0 ,. 3 6 1 1 ..52 1 2 ,.34
8 4,.7 4 6 5..6 3 5 6 .2 0 4 6 .6 2 5 6 ..960 7,.237 7..474 7,.681 7,.863 8 ,. 1 7 6 8 ,. 552 9,.678 1 0 ,.75 1 1 ,. 49
9 4,.5 9 6 5,.428 5 .9 5 7 6 .3 4 8 6 ,.658 6 ..915 7..134 7..325 7,.495 7,.784 8 ,. 132 9 .177 1 0 .,17 1 0 ,. 87
10 4,.482 5..2 7 0 5 .7 6 9 6 .1 3 6 6 ,.428 6 ..6 6 9 6 .,875 7,.055 7,.213 7,.485 7..812 8 .794 9..726 1 0 ,. 39
11 4..392 5..1 4 6 5 .6 2 1 5 .9 7 0 6 .,247 6 .,4 7 6 6 ..672 6 ..842 6 ..992 7,.2 5 0 7,.5 6 0 8 ,.491 9..377 1 0 ,. 00
12 4,.3 2 0 5,.0 4 6 5 .5 0 2 5 .8 3 6 6 ..101 6 ..321 6 ..507 6 ..6 7 0 6 ,. 814 7..0 6 0 7..3 5 6 8 ,. 246 9..094 9,.693
13 4,.2 6 0 4..964 5 .4 0 4 5 .7 2 7 5,.981 6 ..192 6 ..372 6 ,. 528 6 . 667 6 ,.9 0 3 7..188 8 ,. 043 8 ..859 9 .436
14 4..21 0 4..8 9 5 5 .3 2 2 5 .6 3 4 5..881 6 ..085 6 ..258 6 ., 4 0 9 6 ,. 543 6 ,. 772 7,.047 7,.873 8 ..661 9 .219
15 4..168 4..836 5 .2 5 2 5 .5 5 6 5..796 5.,994 6 ..162 6 ..3 0 9 6 ,. 4 3 9 6 ..6 6 0 6 ,. 927 7,.728 8 ..492 9,.035
16 4..131 4 ..7 8 6 5 .1 9 2 5 .4 8 9 5.,722 5..915 6 .,079 6 . 222 6 ,.3 4 9 6 . 564 6 ..823 7,.602 8 ..347 8 ,. 874
17 4,.099 4 ..742 5 .1 4 0 5 .4 3 0 5..659 5..847 6 ..007 6 ,. 147 6 ,. 2 7 0 6 . 480 6 ..734 7,.493 8 ..219 8 ,. 735
18 4..071 4..703 5 .0 9 4 5 .3 7 9 5,.603 5,.788 5..944 6 ,. 081 6 ,. 201 6 ,. 407 6 ..655 7,.398 8 ,. 107 8 .611
19 4..046 4..6 7 0 5 .0 5 4 5 .3 3 4 5,.554 5..735 5.,889 6 . 022 6 ..141 6 ,. 342 6 ..585 7,.313 8 ,. 008 8 .502
20 4,.024 4.,6 3 9 5 .0 1 8 5 .2 9 4 5,.510 5..688 5..8 3 9 5..9 7 0 6 ,. 087 6 ,. 285 6 ,. 523 7,.237 7,.919 8 .404
24 3..956 4..5 4 6 4 .9 0 7 5 .1 6 8 5..374 5..542 5.,685 5..8 0 9 5,.9 1 9 6 ..1 0 6 6 ..3 3 0 7,.001 7..642 8 .097
30 3..889 4,.455 4 .7 9 9 5 .0 4 8 5,.242 5.,401 5..536 5..653 5,.756 5,.932 6 ..143 6 ,. 772 7,.370 7,.796
40 3..825 4,.367 4 .6 9 6 4 .9 3 1 5,.114 5..265 5..392 5..502 5,.5 9 9 5..764 5,.961 6 .547 7,.104 7 .5 0 0
60 3..762 4 ,.2 8 2 4 .5 9 5 4 .8 1 8 4..991 5..133 5..253 5..3 5 6 5,.447 5..601 5,.7 8 5 6 ,. 3 3 0 6 . 843 7 .207
120 3,,702 4..2 00 4 .4 9 7 4 .7 0 9 4,.872 5..005 5..118 5..214 5,.2 9 9 5..443 5,,614 6 ,. 117 6 ,.588 6 .9 1 9
00 3..643 4..1 20 4 .4 0 3 4 .6 0 3 4,.757 4..882 4..987 5..078 5,.157 5..2 9 0 5..4 4 8 5 .911 6 ,.338 6 .636
A p p e n d ix 6 ( C o n t 'd . )
U pper P e r c e n ta g e P o i n t s o f t h e S t u d e n t iz e d Range
a = .0 01
d f/k 2 3 4 5 6 7 8 9 10 12 15 30 60 100
1 9 0 0 .3 1351. 164 3 . 1856. 2022. 2158. 2272. 2370. 2455. 2600. 2770. 3260. 3701. 4002.
2 4 4 .6 9 6 0 .4 2 7 0 .7 7 7 8 .4 3 8 4 .4 9 8 9 .4 6 9 3 .6 7 9 7 .3 0 1 0 0 .5 1 0 5 .9 1 1 2 .3 1 3 1 .0 1 4 8 .0 1 5 9 .7
3 1 8 .2 8 2 3 .3 2 2 6 .6 5 2 9 .1 3 3 1 .1 1 3 2 .7 4 3 4 .1 2 3 5 .3 3 3 6 .3 9 3 8 .2 0 4 0 .3 5 4 6 .6 8 5 2 .5 1 5 6 .5 3
4 1 2 .1 8 1 4 .9 9 1 6 .8 4 1 8 .2 3 1 9 .3 4 2 0 .2 6 2 1 .0 4 2 1 .7 3 2 2 .3 3 2 3 .3 6 2 4 .5 9 2 8 .2 4 3 1 .6 5 3 4 .0 0
5 9 .7 1 4 1 1 .6 7 1 2 .9 6 1 3 .9 3 1 4 .7 1 1 5 .3 5 1 5 .9 0 1 6 .3 8 1 6 .8 1 1 7 .5 3 1 8 .4 1 21.01 2 3 .4 5 2 5 .1 5
6 8 .4 2 7 9 .9 6 0 1 0 .9 7 1 1 .7 2 1 2 .3 2 1 2 .8 3 1 3 .2 6 1 3 .6 3 1 3 .9 7 1 4 .5 4 1 5 .2 2 1 7 .2 8 1 9 .2 2 2 0 .5 8
7 7 .6 4 8 8 .9 3 0 9 .7 6 8 1 0 .4 0 1 0 .9 0 1 1 .3 2 11.68 1 1 .9 9 1 2 .2 7 1 2 .7 4 1 3 .3 2 1 5 .0 5 1 6 .6 9 1 7 .8 5
8 7 .1 3 0 8 .2 5 0 8 .9 7 8 9 .5 2 2 9 .9 5 8 1 0 .3 2 1 0 .6 4 1 0 .9 1 1 1 .1 5 1 1 .5 6 1 2 .0 6 1 3 .5 7 1 5 .0 1 1 6 .0 2
9 6 .7 6 2 7 .7 6 8 8 .4 1 9 8 .9 0 6 9 .2 9 5 9 .6 1 9 9 .8 9 7 1 0 .1 4 1 0 .3 6 1 0 .7 3 1 1 .1 8 1 2 .5 3 1 3 .8 2 1 4 .7 4
10 6 .4 8 7 7 .4 1 1 8 .0 0 6 8 .4 5 0 8 .8 0 4 9 .0 9 9 9 .3 5 2 9 .5 7 3 9 .7 6 9 10.11 1 0 .5 2 1 1 .7 5 1 2 .9 4 1 3 .7 8
11 6 .2 7 5 7 .1 3 6 7 .6 8 7 8 .0 9 8 8 .4 2 6 8 .6 9 9 8 .9 3 3 9 .1 3 8 9 .3 1 9 9 .6 3 0 10.01 1 1 .1 6 1 2 .2 5 1 3 .0 4
12 6 .1 0 6 6 .9 1 7 7 .4 3 6 7 .8 2 1 8 .1 2 7 8 .3 8 3 8 .6 0 1 8 .7 9 3 8 .9 6 2 9 .2 5 4 9 .6 0 6 10.68 1 1 .7 1 1 2 .4 5
399
13 5 .9 7 0 6 .7 4 0 7 .2 3 1 7 .5 9 5 7 .8 8 5 8 .1 2 6 8 .3 3 3 8 .5 1 3 8 .6 7 3 8 .9 4 8 9 .2 8 1 1 0 .2 9 1 1 .2 7 1 1 .9 7
14 5 .8 5 6 6 .5 9 4 7 .0 6 2 7 .4 0 9 7 .6 8 5 7 .9 1 5 8.1 1 0 8 .2 8 2 8 .4 3 4 9 .6 9 6 9 .0 1 2 9 .9 7 2 1 0 .9 1 1 1 .5 7
15 5 .7 6 0 6 .4 7 0 6 .9 2 0 7 .2 5 2 7 .5 1 7 7 .7 3 6 7 .9 2 5 8 .0 8 8 8 .2 3 4 8 .4 8 3 8 .7 8 6 9 .7 0 3 1 0 .5 9 1 1 .2 3
16 5 .6 7 8 6 .3 6 5 6 .7 9 9 7 .1 1 9 7 .3 7 4 7 .5 8 5 7 .7 6 6 7 .9 2 3 8 .0 6 3 8 .3 0 3 8 .5 9 3 9 .4 7 5 1 0 .3 4 1 0 .9 5
17 5 .6 0 8 6 .2 7 5 6 .6 9 5 7 .0 0 5 7 .2 5 0 7 .4 5 4 7 .6 2 9 7 .7 8 1 7 .9 1 6 8 .1 4 8 8 .4 2 7 9 .2 7 7 10.10 1 0 .7 0
18 5 .5 4 6 6 .1 9 6 6 .6 0 4 6 .9 0 5 7 .1 4 3 7 .3 4 1 7 .5 1 0 7 .6 5 7 7 .7 8 8 8.0 1 2 8 .2 8 3 9 .1 0 6 9 .9 0 4 1 0 .4 8
19 5 .4 9 2 6 .1 2 7 6 .5 2 5 6 .8 1 7 7 .0 4 9 7 .2 4 2 7 .4 0 5 7 .5 4 9 7 .6 7 6 7 .8 9 3 8 .1 5 6 8 .9 5 5 9 .7 3 0 1 0 .2 9
20 5 .4 4 4 6 .0 6 5 6 .4 5 4 6 .7 4 0 6 .9 6 6 7 .1 5 4 7 .3 1 3 7 .4 5 3 7 .5 7 7 7 .7 8 8 8 .0 4 4 8 .8 2 1 9 .5 7 5 10.12
24 5 .2 9 7 5 .8 7 7 6 .2 3 8 6 .5 0 3 6 .7 1 2 6 .8 8 4 7 .0 3 1 7 .1 5 9 7 .2 7 2 7 .4 6 7 7 .7 0 1 8 .4 1 1 9 .1 0 0 9 .5 9 6
30 5 .1 5 6 5 .6 9 8 6 .0 3 3 6 .2 7 8 6 .4 7 0 6 .6 2 8 6 .7 6 3 6 .8 8 0 6 .9 8 4 7 .1 6 2 7 .3 7 5 8.021 8 .6 4 7 9 .0 9 6
40 5 .0 2 2 5 .5 2 8 5 .8 3 8 6 .0 6 3 6 .2 4 0 6 .3 8 6 6 .5 0 9 6 .6 1 6 6 .7 1 1 6 .8 7 2 7 .0 6 7 7 .6 5 1 7 .2 1 4 8 .6 1 8
60 4 .8 9 4 5 .3 6 5 5 .6 5 3 5 .8 6 0 6.022 6 .1 5 5 6 .2 6 8 6 .3 6 6 6 .4 5 1 6 .5 9 8 6 .7 7 4 7 .2 9 9 7 .8 0 2 8 .1 6 1
120 4 .7 7 1 5 .2 1 1 5 .4 7 6 5 .6 6 7 5 .8 1 5 5 .9 3 7 6 .0 3 9 6 .1 2 8 6 .2 0 6 6 .3 3 9 6 .4 9 6 6 .9 6 6 7 .4 1 1 7 .7 2 6
00 4 .6 5 4 5 .0 6 3 5 .3 0 9 5 .4 8 4 5 .6 1 9 5 .7 3 0 5 .8 2 3 5 .9 0 3 5 .9 7 3 6 .0 9 2 6 .2 3 4 6 .6 5 1 7 .0 4 1 7 .3 1 4
A p pen d ix 7
C u m ulative D i s t r i b u t i o n o f C h i-S q u a re
D e g re e s o f P r o b a b i l i t y o f a G r e a t e r V alu e
Freedom
0 .9 9 5 0 .9 9 0 0 .9 7 5 0 .9 5 0 0 .9 0 0 0 .7 5 0 0 .5 0 0 0 .2 5 0 0 .100 0 .0 5 0 0 .0 2 5 0.010 0 .0 0 5 0 .0 0 1
1 0.02 0.10 0 .4 5 1 .3 2 2 .7 1 3 .8 4 5 .0 2 6 .6 3 7 .8 8 1 0 .8 3
2 0 .01 0 .02 0 .0 5 0.10 0.21 0 .5 8 1 .3 9 2 .7 7 4 .6 1 5 .9 9 7 .3 8 9 .2 1 1 0 .6 0 1 3 .8 2
3 0 .0 7 0.11 0.22 0 .3 5 0 .5 8 1 .2 1 2 .3 7 4 .1 1 6 .2 5 7 .8 1 9 .3 5 1 1 .3 4 1 2 .8 4 1 6 .2 7
4 0.21 0 .3 0 0 .4 8 0 .7 1 1 .0 6 1 .9 2 3 .3 6 5 .3 9 7 .7 8 9 .4 9 1 1 .1 4 1 3 .2 8 1 4 .8 6 1 8 .4 7
5 0 .4 1 0 .5 5 0 .8 3 1 .1 5 1 .6 1 2 .6 7 4 .3 5 6 .6 3 9 .2 4 1 1 .0 7 1 2 .8 3 1 5 .0 9 1 6 .7 5 2 0 .5 2
6 0.6 8 0 .8 7 1 .2 4 1 .6 4 2 .20 3 .4 5 5 .3 5 7 .8 4 1 0 .6 4 1 2 .5 9 1 4 .4 5 1 6 .8 1 1 8 .5 5 2 2 .4 6
7 0 .9 9 1 .2 4 1 .6 9 2 .1 7 2 .8 3 4 .2 5 6 .3 5 9 .0 4 12.02 1 4 .0 7 1 6 .0 1 1 8 .4 8 2 0 .2 8 2 4 .3 2
8 1 .3 4 1 .6 5 2 .1 8 2 .7 3 3 .4 9 5 .0 7 7 .3 4 10.22 1 3 .3 6 1 5 .5 1 1 7 .5 3 2 0 .0 9 2 1 .9 6 2 6 .1 2
9 1 .7 3 2 .0 9 2 .7 0 3 .3 3 4 .1 7 5 .9 0 8 .3 4 1 1 .3 9 1 4 .6 8 1 6 .9 2 1 9 .0 2 2 1 .6 7 2 3 .5 9 2 7 .8 8
10 2 .1 6 2 .5 6 3 .2 5 3 .9 4 4 .8 7 6 .7 4 9 .3 4 1 2 .5 5 1 5 .9 9 1 8 .3 1 2 0 .4 8 2 3 .2 1 2 5 .1 9 2 9 .5 9
11 2 .6 0 3 .0 5 3 .8 2 4 .5 7 5 .5 8 7 .5 8 1 0 .3 4 1 3 .7 0 1 7 .2 8 1 9 .6 8 2 1 .9 2 2 4 .7 2 2 6 .7 6 3 1 .2 6
12 3 .0 7 3 .5 7 4 .4 0 5 .2 3 6 .3 0 8 .4 4 1 1 .3 4 1 4 .8 5 1 8 .5 5 2 1 .0 3 2 3 .3 4 2 6 .2 2 2 8 .3 0 3 2 .9 1
13 3 .5 7 4 .1 1 5 .0 1 5 .8 9 7 .0 4 9 .3 0 1 2 .3 4 1 5 .9 8 1 9 .8 1 2 2 .3 6 2 4 .7 4 2 7 .6 9 2 9 .8 2 3 4 .5 3
400
14 4 .0 7 4 .6 6 5 .6 3 6 .5 7 7 .7 9 1 0 .1 7 1 3 .3 4 1 7 .1 2 2 1 .0 6 2 3 .6 8 2 6 .1 2 2 9 .1 4 3 1 .3 2 3 6 .1 2
15 4 .6 0 5 .2 3 6 .2 6 7 .2 6 8 .5 5 1 1 .0 4 1 4 .3 4 1 8 .2 5 2 2 .3 1 2 5 .0 0 2 7 .4 9 3 0 .5 8 3 2 .8 0 3 7 .7 0
16 5 .1 4 5 .8 1 6 .9 1 7 .9 6 9 .3 1 1 1.9 1 1 5 .3 4 1 9 .3 7 2 3 .5 4 2 6 .3 0 2 8 .8 5 3 2 .0 0 3 4 .2 7 3 9 .2 5
17 5 .7 0 6 .4 1 7 .5 6 8 .6 7 1 0 .0 9 1 2 .7 9 1 6 .3 4 2 0 .4 9 2 4 .7 7 2 7 .5 9 3 0 .1 9 3 3 .4 1 3 5 .7 3 4 0 .7 9
18 6 .2 6 7 .0 1 8 .2 3 9 .3 9 10.86 1 3 .6 8 1 7 .3 4 2 1 .6 0 2 5 .9 9 2 8 .8 7 3 1 .5 3 3 4 .8 1 3 7 .1 6 4 2 .3 1
19 6 .8 4 7 .6 3 8 .9 1 10. 1 2 1 1 .6 5 1 4 .5 6 1 8 .3 4 2 2 .7 2 2 7 .2 0 3 0 .1 4 3 2 .8 5 3 6 .1 9 3 8 .5 8 4 3 .8 2
20 7 .4 3 8 .2 6 9 .5 9 1 0 .8 5 1 2 .4 4 1 5 .4 5 1 9 .3 4 2 3 .8 3 2 8 .4 1 3 1 .4 1 3 4 .1 7 3 7 .5 7 4 0 .0 0 4 5 .3 2
21 8 .0 3 8 .9 0 1 0 .2 8 1 1 .5 9 1 3 .2 4 1 6 .3 4 2 0 .3 4 2 4 .9 3 2 9 .6 2 3 2 .6 7 3 5 .4 8 3 8 .9 3 4 1 .4 0 4 6 .8 0
22 8 .6 4 9 .5 4 1 0 .9 8 1 2 .3 4 1 4 .0 4 1 7 .2 4 2 1 .3 4 2 6 .0 4 3 0 .8 1 3 3 .9 2 3 6 .7 8 4 0 .2 9 4 2 .8 0 4 8 .2 7
23 9 .2 6 10.20 1 1 .6 9 1 3 .0 9 1 4 .8 5 1 8 .1 4 2 2 .3 4 2 7 .1 4 3 2 .0 1 3 5 .1 7 3 8 .0 8 4 1 .6 4 4 4 .1 8 4 9 .7 3
24 9 .8 9 10.86 1 2 .4 0 1 3 .8 5 1 5 .6 6 1 9 .0 4 2 3 .3 4 2 8 .2 4 3 3 .2 0 3 6 .4 2 3 9 .3 6 4 2 .9 8 4 5 .5 6 5 1 .1 8
25 1 0 .5 2 1 1 .5 2 1 3 .1 2 1 4.6 1 1 6 .4 7 1 9 .9 4 2 4 .3 4 2 9 .3 4 3 4 .3 8 3 7 .6 5 4 0 .6 5 4 4 .3 1 4 6 .9 3 5 2 .6 2
30 1 3 .7 9 1 4 .9 5 1 6 .7 9 1 8 .4 9 2 0 .6 0 2 4 .4 8 2 9 .3 4 3 4 .8 0 4 0 .2 6 4 3 .7 7 4 6 .9 8 5 0 .8 9 5 3 .6 7 5 9 .7 0
40 2 0 .7 1 2 2 .1 6 2 4 .4 3 2 6 .5 1 2 9 .0 5 3 3 .6 6 3 9 .3 4 4 5 .6 2 5 1 .8 0 5 5 .7 6 5 9 .3 4 6 3 .6 9 6 6 .7 7 7 3 .4 0
50 2 7 .9 9 2 9 .7 1 3 2 .3 6 3 4 .7 6 3 7 .6 9 4 2 .9 4 4 9 .3 3 5 6 .3 3 6 3 .1 7 6 7 .5 0 7 1 .4 2 7 6 .1 5 7 9 .4 9 86.66
60 3 5 .5 3 3 7 .4 8 4 0 .4 8 4 3 .1 9 4 6 .4 6 5 2 .2 9 5 9 .3 3 6 6 .9 8 7 4 .4 0 7 9 .0 8 8 3 .3 0 8 8 .3 8 9 1 .9 5 9 9 .6 1
70 4 3 .2 8 4 5 .4 4 4 8 .7 6 5 1 .7 4 5 5 .3 3 6 1 .7 0 6 9 .3 3 7 7 .5 8 8 5 .5 3 9 0 .5 3 9 5 .0 2 1 0 0 .4 2 1 0 4 .2 2 1 1 2 .3 2
80 5 1 .1 7 5 3 .5 4 5 7 .1 5 6 0 .3 9 6 4 .2 8 7 1 .1 4 7 9 .3 3 8 8 .1 3 9 6 .5 8 101.88 1 0 6 .6 3 1 1 2 .3 3 1 1 6 .3 2 1 2 4 .8 4
90 5 9 .2 0 6 1 .7 5 6 5 .6 5 6 9 .1 3 7 3 .2 9 8 0 .6 2 8 9 .3 3 9 8 .6 4 1 0 7 .5 6 1 1 3 .1 4 1 1 8 .1 4 1 2 4 .1 2 1 2 8 .3 0 1 3 7 .2 1
100 6 7 .3 3 7 0 .0 6 7 4 .2 2 7 7 .9 3 8 2 .3 6 9 0 .1 3 9 9 .3 3 1 0 9 .1 4 1 1 8 .5 0 1 2 4 .3 4 1 2 9 .5 6 1 3 5 .8 1 1 4 0 .1 7 1 4 9 .4 5
A p pen d ix 8
PERCENTILE POINTS FOR Q-TEST, FOR EQUAL DEGREES OF FREEDOM v , AND FOR p SAMPLES
v = 1 v = 2 v = 3 v = 4
p .9 9 .9 9 9 .9 9 .9 9 9 .9 9 .9 9 9 .9 9 .9 9 9
3 * * .8 6 3 * .7 5 7 .9 1 9 .6 8 4 .8 2 8
4 .9 2 0 * .7 2 0 .8 9 8 .6 0 5 .7 5 4 .5 4 9 .675
5 .8 2 8 * .6 0 8 .7 7 3 .5 1 2 .6 4 4 .4 4 3 .5 5 2
6 .7 4 4 .9 4 9 .5 3 9 .6 9 0 .4 3 0 .5 4 6 .3 6 9 .461
9 .5 7 6 .7 5 0 .371 .481 .2 8 7 .3 6 3 .2 4 4 .3 0 0
10 .5 2 8 .6 9 4 .3 3 3 .4 3 3 .2 5 7 .3 2 4 .2 1 8 .2 6 7
12 .4 4 8 .5 9 8 .2 7 6 .358 .211 .2 6 5 .1 7 9 .2 1 7
15 .3 6 5 .4 9 0 .2 1 7 .2 8 0 .1 6 5 .205 .1 4 0 .167
16 .3 4 3 .4 6 0 .202 .261 .1 5 4 .1 9 0 .1 3 0 .1 5 5
20 .2 7 3 .3 6 7 .1 5 8 .2 02 .12 0 .1 4 6 .101 .1 1 9
22 .2 4 6 .332 .1 4 2 .1 8 0 .1 0 8 .1 3 0 .0 9 0 .1 0 6
24 .2 2 4 .3 0 2 .1 2 9 .1 6 2 .0 9 8 .1 1 7 .0 8 2 .0 9 6
26 .2 0 6 .2 7 6 .1 1 8 .148 .0 9 0 .1 0 7 .0 7 5 .0 8 7
28 .1 9 0 .2 5 4 .1 0 8 .135 .0 8 2 .0 9 8 .0 6 9 .0 8 0
30 .1 7 6 .234 .10 0 .1 2 4 .0 7 5 .0 9 0 .0 6 4 .0 7 4
32 .1 6 3 .2 1 8 .0 9 3 .1 1 5 .0 7 0 .0 8 3 .0 6 0 .0 6 8
36 .1 4 3 .1 8 9 .0 8 2 .1 00 .0 6 2 .0 7 2 .0 5 2 .0 6 0
40 .1 2 7 .1 6 7 .0 7 2 .0 8 8 .0 5 5 .0 6 4 .0 4 7 .0 5 3
45 .111 .1 4 5 .0 6 3 .0 7 6 .0 4 8 .0 5 5 .041 .0 4 6
50 .0 9 8 .1 2 7 .0 5 6 .0 6 7 .0 4 3 .0 4 9 .0 3 7 .041
60 .0 8 0 .10 2 .045 .0 5 3 .0 3 5 .0 3 9 .0 3 0 .0 3 3
64 .0 7 4 .094 .0 4 2 .0 4 9 .0 3 3 .0 3 7 .0 2 8 .031
♦These e n t r i e s e x c e e d e d 1 u s in g t h e a p p ro x im a tin g d i s t r i b u t i o n . S in c e
Q ^ 1, th e y a r e o m it t e d .
401
A p pen d ix 8 (C o n tin u e d )
v = 5 v = 6 v = 8 v = 10
p .9 9 .9 9 9 .9 9 .9 9 9 .9 9 .9 9 9 .9 9 .9 9 9
3 .631 .7 6 0 .5 9 3 .7 0 8 .5 3 9 .6 3 3 .5 1 2 .5 9 6
4 .4 9 8 .6 0 8 .4 6 1 .5 5 8 .4 1 3 .4 9 0 .3 8 3 .4 4 6
5 .3 9 9 .4 9 0 .3 6 8 .4 4 6 .3 2 8 .3 8 8 .3 0 3 .3 5 1
6 .3 3 4 .4 0 7 .3 0 7 .3 6 8 .271 .318 .2 5 0 .2 8 8
7 .2 8 4 .3 4 5 .261 .3 1 1 .2 3 0 .2 6 8 .21 2 .2 4 2
8 .2 4 6 .2 9 8 .2 2 6 .2 6 8 .1 9 9 .231 .1 8 4 .2 0 9
9 .2 1 7 .261 .1 9 9 .2 3 5 .1 7 6 .20 2 .1 6 2 .1 8 3
10 .1 9 4 .2 3 2 .1 7 8 .2 0 8 .1 5 7 .1 7 9 .1 4 5 .1 6 3
15 .1 2 3 .1 4 5 .1 1 3 .1 3 1 .101 .1 1 3 .0 9 4 .1 0 3
20 .0 9 0 .1 0 4 .0 8 3 .0 9 4 .0 7 4 .0 8 2 .0 6 9 .0 7 5
30 .0 5 8 .0 6 5 .0 5 3 .0 5 9 .0 4 8 .0 5 2 .0 4 5 .0 4 8
40 .0 4 2 .0 4 7 .0 3 9 .0 4 3 .0 3 5 .0 3 8 .0 3 3 .0 3 5
50 .0 3 3 .0 3 6 .031 .0 3 3 .0 2 8 .0 3 0 .0 2 6 .0 2 8
60 .0 2 7 .0 2 9 .0 2 5 .0 2 7 .0 2 3 .0 2 4 .02 2 .0 2 3
v = 12 v = 14 v = 16 v = 20
P .9 9 .9 9 9 .9 9 .9 9 9 .9 9 .9 9 9 .9 9 .9 9 9
3 .4 8 6 .5 5 8 .4 6 6 .5 3 0 .4 5 1 .5 0 8 .4 2 9 .4 7 6
4 .3 6 2 .4 1 5 .3 4 7 .3 9 3 .3 3 5 .3 7 5 .3 1 9 .351
5 .2 8 7 .3 2 6 .2 7 5 .3 0 8 .2 6 5 .2 9 5 .2 5 2 .2 7 6
6 .2 3 6 .2 6 7 .2 2 7 .2 5 3 .2 1 9 .2 4 2 .2 0 9 .2 2 6
7 .201 .2 2 5 .1 9 2 .2 1 3 .1 8 6 .2 0 4 .1 7 8 .191
8 .1 7 4 .1 9 4 .1 6 7 .1 8 4 .1 6 2 .1 7 6 .1 5 4 .1 6 6
9 .1 5 4 .1 7 0 .1 4 8 .1 6 2 .1 4 3 .1 5 5 .1 3 6 .1 4 6
10 .1 3 7 .1 5 2 .1 3 2 .1 4 4 .1 2 8 .1 3 8 .12 2 .1 3 0
15 .0 8 9 .0 9 7 .0 8 6 .0 9 2 .0 8 3 .0 8 9 .0 8 0 .0 8 4
20 .0 6 6 .0 7 0 .0 6 3 .0 6 7 .0 6 2 .0 6 5 .0 5 9 .0 6 2
30 .0 4 3 .0 4 5 .0 4 2 .0 4 3 .0 4 0 .0 4 2 .0 3 9 .0 4 0
40 .0 3 2 .0 3 3 .031 .0 3 2 .0 3 0 .0 3 1 .0 2 9 .0 3 0
50 .0 2 5 .0 2 6 .0 2 4 .0 2 5 .0 2 4 .025 .0 2 3 .0 2 4
60 .021 .022 .02 0 .021 0.20 .0 2 0 .0 1 9 .02 0
2
For v > 60, c a l c u l a t e p v ( p q - l ) and com pare w ith x w ith (p -1 ) degrees o f
freed om in A p pen d ix 7.
402
Appendix 9
Coefficients {a^ for the W test for normality,
f o r n = 2 ( 1 )5 0
X 2 3 4 5 6 7 8 9 10
1 0 .7 0 7 1 0 .7 0 7 1 0 .6 8 7 2 0 .6 6 4 6 0 .6 4 3 1 0 .6 2 3 3 0 .6 0 5 2 0 .5 8 8 8 0 .5 7 3 9
2 .0 0 0 0 .1 6 7 7 .2 4 1 3 .2 8 0 6 .30 3 1 .31 6 4 .3 2 4 4 .3291
3 .0 00 0 .0 8 7 5 .1401 .1 7 4 3 .1 9 7 6 .2141
4 .0 0 0 0 .0561 .0 9 4 7 .1 2 2 4
5 .00 0 0 .0 3 9 9
X 11 12 13 14 15 16 17 18 19 20
1 0 .5 6 0 1 0 .5 4 7 5 0 .5 3 5 9 0 .5 2 5 1 0 .5 1 5 0 0 .5 0 5 6 0 .4 9 6 8 0 .4 8 8 6 0 .4 8 0 8 0 .4 7 3 4
2 .3 3 1 5 .3 3 2 5 .3 3 2 5 .3 3 1 8 .3 3 0 6 .3 2 9 0 .3 2 7 3 .3 2 5 3 .3 2 3 2 .3 2 1 1
3 .2 2 6 0 .2 3 4 7 .2 4 1 2 .2 4 6 0 .2 4 9 5 .2521 .2 5 4 0 .2 5 5 3 .25 6 1 .2 5 6 5
4 .1 4 2 9 .1 5 8 6 .1 7 0 7 .1 8 0 2 .1 8 7 8 .1 9 3 9 .1 9 8 8 .2 0 2 7 .2 0 5 9 .2 0 8 5
5 .0 6 9 5 .0 9 2 2 .1 0 9 9 .1 2 4 0 .1 3 5 3 .1 4 4 7 .1 5 2 4 .1 5 8 7 .16 4 1 .1 6 8 6
6 0 .0 0 0 0 0 .0 3 0 3 0 .0 5 3 9 0 .0 7 2 7 0 .0 8 8 0 0 .1 0 0 5 0 .1 1 0 9 0 .1 1 9 7 0 .1 2 7 1 0 .1 3 3 4
7 .0 0 0 0 .0 2 4 0 .0 4 3 3 .0 5 9 3 .0 7 2 5 .0 8 3 7 .0 9 3 2 .1 0 1 3
8 .00 0 0 .0 1 9 6 .0 3 5 9 .0 4 9 6 .0 6 1 2 .0 7 1 1
9 .00 00 .0 1 6 3 .0 3 0 3 .04 2 2
10 .00 0 0 .0 1 4 0
X" 21 22 23 24*' 25 26 27 28 29 30
1 0 .4 6 4 3 0 .4 5 9 0 0 .4 5 4 2 0 .4 4 9 3 0 .4 4 5 0 0 .4 4 0 7 0 .4 3 6 6 0 .4 3 2 8 0 .4 2 9 1 0 .4 2 5 4
2 .3 1 8 5 .3 1 5 6 .3 1 2 6 .3 0 9 8 .3 0 6 9 .3 0 4 3 .3 0 1 8 .2 9 9 2 .2 9 6 8 .2 9 4 4
3 .2 5 7 8 .2571 .2 5 6 3 .2 5 5 4 .2 5 4 3 .2 5 3 3 .2 5 2 2 .2 5 1 0 .2 4 9 9 .2 4 8 7
4 .2 1 1 9 .2131 .2 1 3 9 .2145 .2 1 4 8 .21 5 1 .2 1 5 2 .2 1 5 1 .2 1 5 0 .2 1 4 8
5 .1 7 3 6 .17 6 4 .1 7 8 7 .1 8 0 7 .1 8 2 2 .1 8 3 6 .18 4 8 .1 8 5 7 .1 8 6 4 .1 8 7 0
6 0 .1 3 9 9 0 .1 4 4 3 0 .1 4 8 0 0 .1 5 1 2 0 .1 5 3 9 0 .1 5 6 3 0 .1 5 8 4 0 .1 6 0 1 0 .1 6 1 6 0 .1 6 3 0
7 .1 0 9 2 .1 1 5 0 .1201 .1245 .1 2 8 3 .1 3 1 6 .1 3 4 6 .13 7 2 .13 9 5 .1 4 1 5
8 .0 8 0 4 .0 8 7 8 .0941 .0997 .1 0 4 6 .1 0 8 9 .1 1 2 8 .11 6 2 .11 9 2 .1 2 1 9
9 .0 5 3 0 .0 6 1 8 .0 6 9 6 .0 7 6 4 .0 8 2 3 .0 8 7 6 .0 9 2 3 .0 9 6 5 .10 02 .1 0 3 6
10 .0 2 6 3 .0 3 6 8 .0 4 5 9 .0 5 3 9 .0 6 1 0 .0 6 7 2 .0 7 2 8 .0 7 7 8 .08 2 2 .0 8 6 2
11 0.0000 0.0 1 2 2 0 .0 2 2 8 0 .0 3 2 1 0 .0 4 0 3 0 .0 4 7 6 0 .0 5 4 0 0 .0 5 9 8 0 .0 6 5 0 0 .0 6 9 7
12 .0 0 0 0 .0 1 0 7 .020 0 .0 2 8 4 .0 3 5 8 .0424 .0 4 8 3 .05 3 7
13 .00 00 .0 0 9 4 .0 1 7 8 .0 2 5 3 .0 3 2 0 .03 8 1
14 .00 0 0 .0084 .0 1 5 9 .02 2 7
15 .000 0 .0 0 7 6
403
Appendix 9 Continued
Coefficients (a^ for the W test for normality,
f o r n « 2 (1 )5 0 (c o n t.)
31 32 33 34 35 36 37 38 39
1 0 .4 2 2 0 0 .4 1 8 8 0 .4 1 5 6 0 .4 1 2 7 0 .4 0 9 6 0 .4 0 6 8 0 .4 0 4 0 0 .4 0 1 5 0 .3 9 8 9
2 .2921 .2898 .2 8 7 6 .2 8 5 4 .2 8 3 4 - .2 8 1 3 .27 9 4 .2 7 7 4 .2 7 5 5
3 .2 4 7 5 .2 4 6 3 .2451 .2 4 3 9 .2 4 2 7 .2 4 1 5 .2 4 0 3 .2391 .2 3 8 0
4 .2 1 4 5 .2141 .2 1 3 7 .2 1 3 2 .2 1 2 7 .2121 .2 1 1 6 .2110 .2104
5 .1 8 7 4 .18 7 8 .1 8 8 0 .1 8 8 2 .1 8 8 3 .1 8 8 3 .1 8 8 3 .18 8 1 .1 8 8 0
6 0 .1 6 4 1 0 .1 6 5 1 0 .1 6 6 0 0 .1 6 6 7 0 .1 6 7 3 0 .1 6 7 8 0 .1 6 8 3 0 .1 6 8 6 0 .1 6 8 9
7 .1 4 3 3 .1 4 4 9 .1 4 6 3 .1 4 7 5 .1 4 8 7 .1 4 9 6 .1505 .1 5 1 3 .1 5 2 0
8 .1 2 4 3 .1 2 6 5 .1 2 8 4 .1301 .1 3 1 7 .13 3 1 .1 3 4 4 .1 3 5 6 .1 3 6 6
9 .1 0 6 6 .1 0 9 3 .1 1 1 8 .1 1 4 0 .1 1 6 0 .1 1 7 9 .1 1 9 6 .1 211 .1 2 2 5
10 .0 8 9 9 .0931 .0961 .09 8 8 .1 0 1 3 .1 0 3 6 .1 0 5 6 .1 0 7 5 .1 0 9 2
11 0 .0 7 3 9 0 .0 7 7 7 0 .0 8 1 2 0 .0 8 4 4 0 .0 8 7 3 0 .0 9 0 0 0 .0 9 2 4 0 .0 9 4 7 0 .0 9 6 7
12 .05 8 5 .0 6 2 9 .0 6 6 9 .0 7 0 6 .0 7 3 9 .0 7 7 0 .07 9 8 .0 8 2 4 .0 8 4 8
13 .0 4 3 5 .0 4 8 5 .0 5 3 0 .0 5 7 2 .0 6 1 0 .0 6 4 5 .06 7 7 .0 7 0 6 .0 7 3 3
14 .0 2 8 9 .03 4 4 .0 3 9 5 .0441 .0 4 8 4 .0 5 2 3 .0 5 5 9 .0 5 9 2 .0 6 2 2
15 .0 1 4 4 .0 2 0 6 .0 2 6 2 .0 3 1 4 .0 3 6 1 .0 4 0 4 .0 4 4 4 .0481 .0 5 1 5
16 0.0000 0 .0 0 6 8 0 .0 1 3 1 0 .0 1 8 7 0 .0 2 3 9 0 .0 2 8 7 0 .0 3 3 1 0 .0 3 7 2 0 .0 4 0 9
17 .0 0 0 0 .00 6 2 .0 1 1 9 .0 1 7 2 .02 2 0 .0 2 6 4 .03 0 5
18 .0 00 0 .0 0 5 7 .0 11 0 .0 1 5 8 .0 2 0 3
19 .00 00 .0 0 5 3 .0101
20 .0 0 0 0
41 42 43 44 45 46 47 48 49
i 0 .3 9 4 0 0 .3 9 1 7 0 .3 8 9 4 0 .3 8 7 2 0 .3 8 5 0 0 .3 8 3 0 0 .3 8 0 8 0 .3 7 8 9 0 .3 7 7 0
2 .2 7 1 9 .2701 .2 6 8 4 .2 6 6 7 .2651 .26 3 5 .2 6 2 0 .2 6 0 4 .2 5 8 9
3 .23 5 7 .23 4 5 .2 3 3 4 .2 3 2 3 .2 3 1 3 .2 3 0 2 .2291 .2 2 8 1 .2271
4 .20 9 1 .2085 .20 7 8 .20 7 2 .20 6 5 .2 0 5 8 .2 0 5 2 .2 0 4 5 .2 0 3 8
5 .1 8 7 6 .1 8 7 4 .1871 .1 8 6 8 .18 6 5 .1 8 6 2 .1 8 5 9 .1 8 5 5 .1851
6 0 .1 6 9 3 0 .1 6 9 4 0 .1 6 9 5 0 .1 6 9 5 0 .1 6 9 5 0 .1 6 9 5 0 .1 6 9 5 0 .1 6 9 3 0 .1 6 9 2
7 .1531 .15 3 5 .1 5 3 9 .1 5 4 2 .1 5 4 5 .1 5 4 8 .1 5 5 0 .1551 .15 5 3
8 .13 8 4 .1392 .1 3 9 8 .14 0 5 .1 4 1 0 .1 4 1 5 .1 4 2 0 .1 4 2 3 .1 4 2 7
9 .1 2 4 9 .1 2 5 9 .1 2 6 9 .1 2 7 8 .1 2 8 6 .1 2 9 3 .1 3 0 0 .1 3 0 6 .1 3 1 2
10 . 1123 .1 1 3 6 .1 1 4 9 .1 1 6 0 .1 1 7 0 .1 1 8 0 .1 1 8 9 .1 1 9 7 .1 2 0 5
11 0 .1 0 0 4 0.1020 0 .1 0 3 5 0 .1 0 4 9 0 .1 0 6 2 0 .1 0 7 3 0 .1 0 8 5 0 .1 0 9 5 0 .1 1 0 5
12 .0891 .0 9 0 9 .0 9 2 7 .0 9 4 3 .0 9 5 9 .0 9 7 2 .0 9 8 6 .0 9 9 8 .10 1 0
13 .0 7 8 2 .0 8 0 4 .0 8 2 4 .08 4 2 .0 8 6 0 .0 8 7 6 .0892 .0 9 0 6 .0 9 1 9
14 .0 6 7 7 .0701 .0 7 2 4 .07 4 5 .07 6 5 .0 7 8 3 .0801 .0 8 1 7 .0 8 3 2
15 .0 5 7 5 .0602 .0 6 2 8 .0651 .0 6 7 3 .0 6 9 4 .0 7 1 3 .07 3 1 .0 7 4 8
16 0 .0 4 7 6 0 .0 5 0 6 0 .0 5 3 4 0 .0 5 6 0 0 .0 5 8 4 0 .0 6 0 7 0 .0 6 2 8 0 .0 6 4 8 0 .0 6 6 7
17 .0 3 7 9 .0411 .0 4 4 2 .0471 .0 4 9 7 .0 5 2 2 .0 5 4 6 .0 5 6 8 .0 5 8 8
18 .0 2 8 3 .0 3 1 8 .03 5 2 .0 3 8 3 .04 1 2 .0 4 3 9 .0 4 6 5 .0 4 8 9 .0511
19 .0 1 8 8 .0227 .0 2 6 3 .0 2 9 6 .03 2 8 .0 3 5 7 .03 8 5 .0 4 1 1 .0 4 3 6
20 .0 0 9 4 .0 1 3 6 .0 1 7 5 .02 11 .0 2 4 5 .0 2 7 7 .0 3 0 7 .0 3 3 5 .0361
21 0.0000 0 .0 0 4 5 0 .0 0 8 7 0 .0 1 2 6 0 .0 1 6 3 0 .0 1 9 7 0 .0 2 2 9 0 .0 2 5 9 0 .0 2 8 8
22 .0 0 0 0 .0 0 4 2 .0081 .0 1 1 8 .0 1 5 3 .0 1 8 5 .02 1 5
23 .00 00 .0 0 3 9 .0 0 7 6 .0111 .0 1 4 3
24 .00 0 0 .0 0 3 7 .0071
25 .00 0 0
404
A p pen d ix 10
_______________ P e rc e n ta g e p o i n t s o f th e W t e s t f o r n = 3 ( 1 ) 5 0 ____________
L evel
n 0.01 0.02 0.05 0.10 0.50 0.90 0.95 0.98
3 0 .7 5 3 0 .7 5 6 0 .7 6 7 0 .7 8 9 0 .9 5 9 0 .9 9 8 0 .9 9 9 1 .0 0 0
4 .687 .7 0 7 .7 4 8 .792 .9 3 5 .9 8 7 .9 9 2 .9 9 6
5 .6 8 6 .7 1 5 .7 6 2 .8 0 6 .9 2 7 .9 7 9 .9 8 6 .991
6 0 .7 1 3 0 .7 4 3 0 .7 8 8 0 .8 2 6 0 .9 2 7 0 .9 7 4 0 .9 8 1 0 .9 8 6
7 .7 3 0 .7 6 0 .8 0 3 .8 3 8 .9 2 8 .972 .9 7 9 .9 8 5
8 .7 4 9 .7 7 8 .8 1 8 .851 .932 .972 .9 7 8 .9 8 4
9 .764 .791 .8 2 9 .8 5 9 .9 3 5 .9 7 2 .9 7 8 .9 8 4
10 .781 .8 0 6 .8 4 2 .8 6 9 .9 3 8 .9 7 2 .9 7 8 .9 8 3
11 0 .7 9 2 0 .8 1 7 0 .8 5 0 0 .8 7 6 0 .9 4 0 0 .9 7 3 0 .9 7 9 0 .9 8 4
12 .8 0 5 .8 2 8 .8 5 9 .8 8 3 .9 4 3 .9 7 3 .9 7 9 .9 8 4
13 .8 1 4 .8 3 7 .8 6 6 .8 8 9 .9 4 5 .9 7 4 .9 7 9 .9 8 4
14 .825 .8 4 6 .8 7 4 .895 .9 4 7 .9 7 5 .9 8 0 .9 8 4
15 .835 .8 5 5 .8 8 1 .901 .9 5 0 .975 .9 8 0 .9 8 4
16 0 .8 4 4 0 .8 6 3 0 .8 8 7 0 .9 0 6 0 .9 5 2 0 .9 7 6 0 .9 8 1 0 .9 8 5
17 .851 .8 6 9 .8 9 2 .910 .9 5 4 .9 7 7 .981 .985
18 .8 5 8 .874 .8 9 7 .9 1 4 .9 5 6 .9 7 8 .9 8 2 .9 8 6
19 .8 6 3 .8 7 9 .901 .9 1 7 .9 5 7 .9 7 8 .9 8 2 .9 8 6
20 .8 6 8 .8 8 4 .9 0 5 .9 2 0 .9 5 9 .9 7 9 .9 8 3 .9 8 6
21 0 .8 7 3 0 .8 8 8 0 .9 0 8 0 .9 2 3 0 .9 6 0 0 .9 8 0 0 .9 8 3 0 .9 8 7
22 .878 .8 9 2 .911 .9 2 6 .961 .9 8 0 .9 8 4 .9 8 7
23 .881 .8 9 5 .9 1 4 .9 2 8 .9 6 2 .981 .9 8 4 .9 8 7
24 .884 .8 9 8 .9 1 6 .9 3 0 .9 6 3 .981 .9 8 4 .9 8 7
25 .8 8 8 .901 .9 1 8 .931 .9 6 4 .981 .9 8 5 .9 8 8
26 0 .8 9 1 0 .9 0 4 0 .9 2 0 0 .9 3 3 0 .9 6 5 0 .9 8 2 0 .9 8 5 0 .9 8 8
27 .8 9 4 .9 0 6 .9 2 3 .935 .9 6 5 .9 8 2 .9 8 5 .9 8 8
28 .8 9 6 .9 0 8 .9 2 4 .9 3 6 .9 6 6 .9 8 2 .985 .9 8 8
29 .8 9 8 .9 1 0 .9 2 6 .9 3 7 .9 6 6 .982 .9 8 5 .9 8 8
30 .9 0 0 .9 1 2 .9 2 7 .9 3 9 .9 6 7 .9 8 3 .9 8 5 .9 8 8
31 0 .9 0 2 0 .9 1 4 0 .9 2 9 0 .9 4 0 0 .9 6 7 0 .9 8 3 0 .9 8 6 0 .9 8 8
32 .9 0 4 .9 1 5 .9 3 0 .941 .9 6 8 .9 8 3 .9 8 6 .988
33 .9 0 6 .9 1 7 .931 .942 .9 6 8 .9 8 3 .9 8 6 .9 8 9
34 .9 0 8 .9 1 9 .9 3 3 .9 4 3 .9 6 9 .9 8 3 .9 8 6 .9 8 9
35 .9 1 0 .9 2 0 .9 3 4 .9 4 4 .9 6 9 .984 .9 8 6 .9 8 9
36 0 .9 1 2 0 .9 2 2 0 .9 3 5 0 .9 4 5 0 .9 7 0 0 .9 8 4 0 .9 8 6 0 .9 8 9
37 .9 1 4 .9 2 4 .9 3 6 .9 4 6 .9 7 0 .9 8 4 .9 8 7 .9 8 9
38 .9 1 6 .9 2 5 .9 3 8 .9 4 7 .971 .9 8 4 .9 8 7 .9 8 9
39 .9 1 7 .9 2 7 .9 3 9 .9 4 8 .971 .9 8 4 .9 8 7 .989
40 .9 1 9 .9 2 8 .9 4 0 .9 4 9 .9 7 2 .9 8 5 .9 8 7 .989
41 0 .9 2 0 0 .9 2 9 0 .9 4 1 0 .9 5 0 0 .9 7 2 0 .9 8 5 0 .9 8 7 0 .9 8 9
42 .922 .9 3 0 .9 4 2 .951 .972 .985 .9 8 7 .9 8 9
43 .9 2 3 .9 3 2 .9 4 3 .951 .9 7 3 .9 8 5 .9 8 7 .9 9 0
44 .9 2 4 .9 3 3 .9 4 4 .9 5 2 .9 7 3 .9 8 5 .9 8 7 .9 9 0
45 .9 2 6 .9 3 4 .945 .9 5 3 .9 7 3 .9 8 5 .9 8 8 .990
46 0 .9 2 7 0 .9 3 5 0 .9 4 5 0 .9 5 3 0 .9 7 4 0 .9 8 5 0 .9 8 8 0 .9 9 0
47 .9 2 8 .9 3 6 .9 4 6 .9 5 4 .9 7 4 .9 8 5 .9 8 8 .9 9 0
48 .9 2 9 .9 3 7 ,9 4 7 .9 5 4 .9 7 4 .9 8 5 .9 8 8 .9 9 0
49 .9 2 9 .937 .9 4 7 .955 .9 7 4 .985 .9 8 8 .9 9 0
50 .9 3 0 .9 3 8 .9 4 7 .955 .9 7 4 .9 8 5 .9 8 8 .9 9 0
405
A p pen d ix 11
C o e f f i c i e n t s o f O r th o g o n a l P o ly n o m ia ls
n P o ly n o m ia l X=1 2 3 4 5 6 7 8 9 10 IZ2 X
2 L in e a r -1 1 2 2
3 L in e a r -1 0 1 2 1
Q u a d r a tic 1 -2 1 6 3
4 L in e a r -3 -1 1 3 20 2
Q u a d r a tic 1 -1 -1 1 4 1
C u b ic -1 3 -3 1 20 1 0 /3
5 L in e a r -2 -1 0 1 2 10 1
Q u a d r a tic 2 -1 -2 -1 2 14 1
C u b ic -1 2 0 -2 1 10 5 /6
Q u a r t ic 1 -4 6 -4 1 70 3 5 /1 2
6 L in e a r -5 -3 -1 1 3 5 70 2
Q u a d r a tic 5 -1 -4 -4 -1 5 84 3 /2
C u b ic -5 7 4 -4 -7 5 180 5 /3
Q u a r t ic 1 -3 2 2 -3 1 28 7 /1 2
7 L in e a r -3 -2 -1 0 1 2 3 28 1
Q u a d r a tic 5 0 -3 -4 -3 0 5 84 1
C u b ic -1 1 1 0 -1 -1 1 6 1 /6
Q u a r tic 3 -7 1 6 1 -7 3 154 7 /1 2
8 L in e a r -7 -5 -3 -1 1 3 5 7 168 2
Q u a d r a tic 7 1 -3 -5 -5 -3 1 7 168 1
C u b ic -7 5 7 3 -3 -7 -5 7 264 2 /3
Q u a r tic 7 -1 3 -3 9 9 -3 -1 3 7 616 7 /1 2
Q u in t ic -7 23 -1 7 -1 5 15 17 -2 3 7 2184 7 /1 0
9 L in e a r -4 -3 -2 -1 0 1 2 3 4 60 1
Q u a d r a tic 28 7 -8 -1 7 -20 -1 7 -8 7 28 2772 3
C u b ic -1 4 7 13 9 0 -9 -1 3 -7 14 990 5 /6
Q u a r tic 14 -2 1 -11 9 18 9 -11 -21 14 2002 7 /1 2
Q u in t ic -4 11 -4 -9 0 9 4 -1 1 4 468 3 /2 0
10 L in e a r -9 -7 -5 -3 -1 1 3 5 7 9 330 2
Q u a d r a tic 6 2 -1 -3 -4 -4 -3 -1 2 6 132 1/2
C u b ic -4 2 14 35 31 12 -12 -3 1 -3 5 -1 4 42 8580 5 /3
Q u a r t ic 18 -22 -1 7 3 18 18 3 -1 7 -2 2 18 2860 5 /1 2
Q u in t ic -6 14 -1 -1 1 -6 6 11 1 -1 4 6 780 1/10
406
A p pen d ix 12
V a r ia n c e o f Means and C o n t r a s t s
fo r
M ixed M odels
f o l l o w i n g d e r i v a t i o n s a r e a p p r o p r ia t e [A nd erson and B a n c r o ft (1 9 5 2 ), p . 3 4 0 -3 4 4 ]:
C o n s id e r th e la y o u t
1 2 j . . . . c
n
yi..
yi..
Rows (p )
hi-
1
yr ..
Y. = y + p. + y . + (p y ) . . + e . .. ,
13k l *3 v" 1 3 13k ’
le t pa be ( 0 , a R) ,
2 (p y )^
2
b e ( 0 , a RC) and e . j k b e ( 0 ,a ) .
2
is :
r n
I X Y. n ][ p .
Y . = i=l k=l
■ ljk = p + — i------1 + r n Yi
L
•3 • rn m rn
n \ (p y ) • • I \ e.
i= l l> .i = l k=l
S in c e E (Y j) * Yy E (y ) = y , and E ( a l l random v a r i a b l e s in th e m o d e l) = 0
E(CP) * 0
E (Y -j.) = M + Yj
407
The variance of the estimate of the column mean is:
Z p. Z (p y ) . . Z Z £.
- 2 (i * 1A i . . ljk 2
'J ik
V (Y -j-) = E[Y
• r - E(v,j,)1 =E| T - (i- 7 - * V |
bu t E(p?) = a2, E (p y )?j = o ^ c > and E (e i ;j k) 2 = ^ •
Hence // 7
2 . 02 (C PI)
CCP1) v
+ ... + p + C ro ss p r o d u c t s 1\
V(Y
■i ■ -?- - - - - - - )
(PY)j j ♦ . . . + (PY)^ ♦ CP2
2
r
( i i + ... + e2 . + CP A x
*( — 7 ^ — 7 * cp<l
^r a R r 0 RC r n a2 aR aRC a2
and V(Y . ) » - V 1 + — ~ + ■■ ■ = -
*J • I 2 2 22 r r rn
J yr r r n
2. V a r ia n c e o f th e mean i n s i d e th e t a b l e .
E (Y ..J = v + Yj
and
V (Y . . ) = E [Y . . - E (Y . . ) ]
- 2 (
= E I p . + ( p y ) . . + — -------- J
k=lEi:ik {
1 1i j3 • 1 ij • v ( i VKWi j n j
E {pi * (p Y )i j * ~ i± 1— * CP6 )
2 2 no2
°R + °RC + I "
n
F in a lly
2
2 2 .a
408
3. Variance of £ contrast.
Y «q * , then th e d i f f e r e n c e i s r e a l l y th e c o n t r a s t :
E p. Z (p y )• Z Z z . .
( i 1 i ip i k ip k
(Y•p. - Y
-q.) = lu + + y + +-— ----------
v ( r p r rn
Z p. Z ( py) . Z Z z. .
i 1 i it 1Rki
r “ Yq r ” rn j
2 (P Y )ip E (P Y )i q
(Y
•p.
- Y
»q»
) = f CY
r Tp
- Y ) + -
q
Z Z e. , Z Z z. .
i k ip!c ■k iqk
rn ~ rn j
The v a r ia n c e o f t h i s c o n t r a s t i s
V(Y - Y ) = E[ (Y - Y ) - E(Y - Y ) ] 2.
1 *p* - q* 1 -p * -q- - p* #q* J
Now
E (V .p . - V .q . ) - (Yp - Yq )
and
/Z ,r (p y ) • y (p y ) • v
E \ //v
z r
Z c
z. , Z r
r Z c
z-
. , \\ ^
V(Y - Y
•p* -q
„fv
V(V -v
V i
)-' r°R
—C * r°RC * ~T~2
2~ rn°2 . ~2~2
rn<y2
r n r r r n r n
_ 2gRC + 2o^
r m
2 2
Knowing t h a t E.M .S. (R .C .) = naRC + a , th e e s tim a t e o f t h i s v a r ia n c e is
409
be q u i t e d i f f e r e n t from th e v a r ia n c e o f th e c o n t r a s t . To b e s a f e f o r m ixed m o d e ls ,
th e e x p e r im e n te r must d e r i v e t h e a p p r o p r ia t e s ta n d a rd e r r o r and v a r ia n c e o f th e
c o n t r a s t from t h e m od el.
O f c o u r s e t h e r e i s no i n t e r e s t in means f o r th e random m o d e l; c o n s e q u e n t ly
e s p e c ia lly a p p lic a b le . C o n s id e r th e r e p e a t e d t - t e s t s
(/2 s -) i s r e a l l y th e s ta n d a rd e r r o r o f th e c o n t r a s t and we m e r e ly r e p l a c e t h i s
com ponent by
f o r t h e m ixed m od el. O f c o u r s e th e e x p e r im e n te r w i l l t h o r o u g h ly i n v e s t i g a t e th e
tre a tm e n t means.
Newman-Keuls t e s t u ses
Rk = <*f) ( s - ) •
’ ( f o r a f i x e d m o d e l).
But th e sta n d a rd e r r o r f o r t h e c o n t r a s t i s
Hence
qa (k, d f)
(N ote t h a t when k = 2 , = t ) and, f o r th e m ixed m o d e l, s u b s t itu tin g
./T a
(mean sq u a re RC)
fo r / 2 s - , we g e t
qa ( k , d f )
sq u a re RC)
7T~
or
410
RC
Rk • rn
N ested F a c t o r i a l f o r S u b s e c t io n 6 . 2 . 3 .
E q u a tion (6 .2 .3 ):
jh HW i , j^(U)
Yi.. 8 8
5 55 8 5
I \k I RW.k
v D“l k I I HW I I e 1
k=l k k=l j =l k=l j =l k=l UJKJ I
+ 5 * 40 + 40 J
U 7~0
S in c e weeks (W^) a r e f i x e d and hens [H ^ ^ ] a re random , th e summation
o f t h e m ixed i n t e r a c t i o n o v e r t h e s u b s c r ip t k i s d e fi n e d as z e r o .
THu* (j i Ji^cnjkV 0 •
E (y . ) = y + R. 8
u „
8 ~
8 5 ,
411
5. E x t e n s io n s .
412
AUTHOR INDEX
Addelman, S., 126, 153, 272, Conner, W. S., 295, 301, 317,
274, 278 319, 325
Anderson, R. L., 27, 28, 29, 38, Cox, D. R., 85, 93, 94, 121, 330,
52, 78, 111, 117, 118, 121, 331, 332, 386
156, 160, 180, 213, 223 Cox, G. M., 94, 121, 151, 153,
Anderson, V. L., 40, 78, 87, 93, 209, 216, 219, 223, 302, 320,
124, 153, 165, 180, 184, 212, 325, 348, 349, 386
223, 275, 276, 278, 342, 386
Daniel, C., 206, 209, 261, 270,
Bancroft, T. A., 10, 27, 28, 29, 271, 278
36, 38, 52, 78, 105, 111, 117, Davies, 0. L. (Ed.), 104, 121,
118, 121, 141, 153, 186, 208, 237, 245, 251, 252, 255, 256,
213, 223 278, 281, 301, 353, 370, 373,
Bartlett, M. S., 19, 20, 25, 27, 386
38 Dixon, W. J., 7, 23, 38
Beeson, J., 131, 153, 184, 189, Draper, N. R., 27, 28, 38, 374,
208 379, 386
Behnken, D. W., 367, 386 Dykstra, 0., 370, 386
Bennett, C. A., 52, 78, 125,
153, 208 Edwards, A. L., 104, 111, 121,
Berger, P. D., 263, 278 180, 213, 215, 223
Bliss, C. I., 208 Eisenhart, C., 7, 38, 159, 180
Booth, K. H. V., 330, 331, 332,
386 Federer, W. T., 185, 209, 302,
Bowman, K. 0., 7, 39 320, 321, 322, 323, 324, 325
Box, G. E. P., 17, 38, 272, 278, Finney, D. J., 87, 93, 120, 121
348, 349, 350, 353, 367, 370, Fisher, R. A., 87, 93, 151, 153,
374, 379, 385, 386 209, 211, 218, 223, 224
Bozivich, H., 36, 38, 105, 121, Foster, L. A., 20, 38
141, 153, 186, 208 Franklin, N. L., 52, 78, 125,
Bratcher, T. L., 7, 38 153
Brooks, S. H., 333, 361, 362,
386 Geisser, S., 166, 180
Burman, J. P., 330, 386 Gorman, J. W., 335, 340, 386
Burr, I. W., 20, 38 Greenhouse, S. W., 166, 180
413
414 AUTHOR INDEX
415
416 SUBJECT INDEX