Rumus Federer
Rumus Federer
Rumus Federer
,/-
~---
,:
by
Walter T. Federer, Cornell University
Abstract
Statistical design has many facets, two of which are experiment design and
treatment design.
treatme~ts
i~
an ex-
;~ed
here because
it has been used in many diffe:r-ent contexts and besides the '':iesign" should :not
be "experimental".
the 1920's.
Since then
connecte~Dess,
resolvability,
Many of the above principles have been used for treatment de-
sign and others have been added, e.g., unbiasedness, rotatability, mean squared
error optimality, alias balance, saturation, and symmetry.
To many writers of
in attaining certain properties for the plans, the last two are not statistical
design.
planning investigations.
by
BU-813-M
Walter T. Federer, Cornell University
May, 1983
Revised July 1983
Introduction
Statistical design for an investigation encompasses many items; these are:
(i)
(s.u. ), the elements making up the population; unit of the experiment (e.u. ),
or experiment unit, the smallest unit to which one treatment is applied; of the
variables to be Eeasured, as well as their 1istributions; the response equations
prior to applying treatments; and the observational
Q~it
(ii)
Survey design.
Sequential design.
previous results determine which treatment(s) is(are) applied to the next e.u. (s)
- 2 -
(vii)
met~e=atical
expression of results is
Selection of an appro~~d
arrange~ent
model design.
(viii)
numoe~,
Statistical
crite~ia,
as well as available
reso~ces,
(ix)
to as experiment technique
that popular in the field
o~
o~
~oo
plot technique).
investigation.
refe~red
Several principles
~ar-domization,
control (blocking) and several others such as orthogonality, confounding, sensitivity, balance, efficiency, connectedness,
and sufficiency.
designs.
~esolvability,
variance
optimali~J,
thogonality, experiment and treatment designs have the property of being orthogonal
designs.
Some basic axioms (rules) to follow in designing an experiment are given
below.
Fail~e
to follow
~~y
expe~ience
experimenters is:
- 3 Axiom I:
complete, precise,
which inferences
fu~d
It is insufficient to say "the results of this experiment apply to the population of which it is a sample".
Textbooks
and literature, statistical or otheTwise, ignore this fact for the most part.
Any reference to the population is more of a "hand-waving" description.
does not
-~iom
~ow
II:
If one
des~gr:.
All too often ar;. experimer:ter and statistician feel their choice
ac.
c:~
th~~
Following
~e
added
II means
JLv.io~
that new designs may have to be constructed to meet the requirements o: the experiment.
usi~g
One can easily construct many types of designs that are ::ot
i~
Wilii"'m"~
tables,
lsr76,
Ttere is no
need to allocate extra degrees of freedom to blocking when they should be allocated to error variance.
seven e.u. 's and has four treatments, the block size should be seven, not four.
~
Each of the treatments but one can be included twice in each block, ~~d a different
treatment is omitted from each block.
- 4 (An experiment
Also, if several
tr~~splanting
blocki~~.
~~d
all be called a block rather than blocking on three separate sources of variation
and using a complex
design.
experime~t
~ill
suffice.
Ease of statistical
a~a~'Ses,
when data are lost, and interpretations are additional reasons for
us~ng
especially
a design
numbers of replicates, will have different e.u. 's, and will have
~
different~
variances.
Consider the two sets of data, 100 measurements on one plant and one measurement from each of 100 plants.
puted using exactly the same formulae for both sets of data.
interested in the
computir~
The indi7idual
plant.
was randomly selected from the population of plants, the sample size is one, no
matter if a million measurements on that one plant are made.
100 randomizations were used and the sample size of plants is 100.
in this example should be obvious to almost everyone, but in
distinction may be rather subtle.
m~~y
The differences
situations the
different treatments are different sizes and involve different blocking patterns.
- 5 -
Axiom V:
differe~ce
must contain all sources of variation in the e.u. 's except that due to the treatments themselves.
This is Sir R. A. Fisher's
(1935)
definitio~
so~e
term designated as
measu~ement
for different
treatments.
Many, many experiment designs and classes of designs have been devised by
statisticians and mathematicians, but they have been constructed without any con-
for~
It should
be noted that the linear model provided may be incorrect and inappropriate, leading
to incorrect statistical analyses and inferences.
the above five axioms should be kept
L~
Fisher
Figure 2
(1972).
(1933, 1937).
section.
- 6 -
Replication
Validity of estimate
of experimental error
Figure 1.
Reduction of
experimental error
- 7 -
randomization
...
valid estimate of
experimental error
'---,--------~
~ ~~
balancing
""'
f?artial
~---------------~Jconfoundinq 1~---
Figure 2 .
4lt
des~gn.
treat~ent
Tee reason
~o
de~ined,
Additional
ar.d then to
this
~stake
va~iance
are not
fo~
is
Some of the
~~ture.
additional ones discussed below are balance, se:r:sitivity, efficiency, coru:ectedness, optimality, resolvability, and
sufficie~cy.
Replication
The replicatio::: principle is a
Experiment
units
necessa~:.- pa~t
of scientific ir:-.e.s-:igation.
t~eatment
variation
The number of e.u. 's, not the number of s.u. 's or o.u. 's,
4lt
of replications in an experiment.
e._;_. to e. u.
:~~om
dete~mine t~e
number
The idea
o~
repeat-
the duplicate plot, give me one and I know where I am" was not much of a scientist.
He wanted to ignore the universal fact that e.u. 's differ and that repeatability
of results of scientific investigations is a necessity.
Randomization
Randomization is the basis for obtaining a valid estimate of a
effect and valid estimate of error variation.
been unclear for
41t
~~y,
The concept of
~~derstand
t~eatment
r~~domization
has
~andom-
permutatio~s
treat~ef-L
~ermutation
results in a
systematically arranged experiment design, which usually results in biased estimates of treatment effects and of error variances.
be avoided.
The randomization principle is difficult for certain scientists because they
feel they know how to place treatments in an experiment.
vestigators can
results.
j~diciously
place treatments in
the experimental
~reat~e~ts
~terial
lli~der
investigation,
as :est as possible.
impl~es
t~e
biasedness.
In order
a~ ex~er~ment
to be fair to the
BlockL~g
experime~t
as discussed above is
techniques.
The
Continuous
variates classified into groups and called blocks, e.g., weight of plants, is an
inefficient method.
the experimental material is controlled and that this variation can be removed
- 10 -
ment
error,
thereby
increasir~
exper~e~tal
Theorems on and
othe~
field.
p~ocedures
p~inciple
designs have the property that the various sets of effects are orthogonal to each
other, that simple statistical analyses result, that valid estimates of effects
and their variances are obtained, and the greatest precision usually results.
Simple arithmetic averages may be used to obtain the differences between estimated
effects, and the variances of these differences are easily computed when a linear
model holds.
experiment designs, there is no loss in information due to this source, and maximum precision results.
As Preece
To illustrate
: r vj
- ll -
are the same for every block j, we say that treatment effects are orthogonal to
~efinition
If differences between arithmetic means for any and all pairs of means, contains
only treatment and random error terms, treatments are orthogonal to all other
effects in the linear model.
If one
ha~
contai~s
o~~y
geo~etric
.
.
. J..l.Tl\( .IIr .A. . . \1/r;
( nr .....A . )1/r
1.gnor2ng
eryor t erms, lS
)
i ...:.1: 2 .
.
=
'J=l J.
J=l J.
1:
I
1 -r 2
Th e comb"2na t or1.c
.
However, the
user must first decide what type of orthogonality is being considered for his
situation.
deter~~e
If two or
more effects are not separable, the effects are said to be completely confounded.
Otherwise, when the effects can be separately estimated and they are not orthogonal,
they are said to be partially confounded.
plans for v
=3
For
ex~ple,
=4
blocks:
- 12 -
Plan I
Plan II
Block
Block
c
c
c
c
c
c
In the first design plan treatment A effects a!:d blocks 1 and 2 are completely
confounded.
In design
pl~
I:~
co~foQ~dei,
block and if treatments B and C had bee!: ir.cluied once in each of the
i~
fo~
every
blocks,
the plan would have been an orthogonal design plar. with blocks and treatment
th~~
can be included in a block, one must use an incomplete block design with partial
confounding of treatment and block effects.
on some effects than others, or when it is
ment, the principle of confounding may be invoked to obtain smaller blocks than
can accommodate all the treatments in each block.
When there is no confounding of various types of effects in a statistical
design, orthogonality results.
be estimated in the experiment.
statistical information.
- 13 -
Balancedness
Balancedness of effects is another useful principle in selecting and constructing design plans for an experiment.
ment effects.
There are many types of balance (see Heiayat and Federer, 1974, and Preece, 1982),
but attention here is confined mostly to variance balance; i.e., when the errors
in the experiment come from a distribution with a single ,variance parameter (homoscedastic), all differences between treatments will have the same variance.
When
There are situations where a balanced design has a lower error variance
ments A, B, and C.
I.
Complete Block
II.
Complete Block
III.
Complete Block
_g_
_3_
c
c
c
c
c
c
c
c
c
c
c
c
c
c
c
competL~
One
=3
as in plan III.
- 14 -
groups.
II would cost
1/3
III would cost 8/3- l = 5/3 more than I and twice as much as II.
If each e. u.
costs the same regardless of which group it is in, then the three plans would
cost the same.
However, there are many situations where the e.u. 's do not cost
the same and some e.u. 's would need to be discaried to obtain plans II and III.
For example, suppose that there were 8 positions in an oven for baking pies.
There would only be 3 bakings for I but 8 for III, resulting in considerably
more expense.
Plan I is variance balanced, and II and III are botq variance balanced and
orthogonal.
pairwise balanced in that every pair of treatments occurs together in the blocks
an equal number of times A; note that A
= 21
(1977) 23 pages to discuss the orthogonality muddle, it took him (1982) 102 pages
to discuss the balance tangle.
terms.
Efficiency
The principle of efficiency is useful in deciding which one or ones of a
class of experiment designs is optimal (best) with respect to a selected criterion.
Making use of this principle, designs can be ranked with respect to their relative
efficiencies using the selected criterion.
relative efficiency with respect to variance (perhaps not if cost and variance
are both considered).
says to utilize the most efficient design at one's disposal and suited for the
experiment.
- 15 -
iesigns.
balance~~ess
and cr2j4.
have
a2/4,
E
in degrees of freedom for the remainder sums of squares (i.e., 19, 18, and 14,
respectively), are 0.989, 1.000, and 0.975.
cost of the additional groups required for designs II and III, a balanced design
is more efficient than a competing orthogonal design, III, and is almost as
efficient, 0.989, as the best competing orthogonal design.
additional cost, this balanced incomplete block design
either of the competing orthogonal designs.
i~
v.
Plan
Partially Balanced
Incomplete Block
Block
Block
_2_
_3_
A
c
c
_]_
_5_
c
c
c
c
c
c
D
D
The average variance of a difference for plan IV is cr2/2 and that for plan V is
E
7cr2/18.
E
is limited to the case where treatment i occurs in block j either zero or once,
i.e., n ..
lJ
0, 1.
qua::era_~,
= 0,
~'
= 0,
~'
2, a binary design.
~ina_~
Plan V is a
~ernayy
iesign, n ..
~J
= 0,
It is not
In plan IV,
1, 2, and has
a smaller variance.
Other criteria are often necessary
i~
experimentation.
For example, in a
nutrition experiment, a variance minimal design was selected for one-half of the
experiment because it was nutritionally
opt~~al
in the
eA~erime~t.
a~d
~ense
of having other
statistical optimality.
Other criteria
Optimality
The optimal principle states that the
utilized in conducting experiments.
~st
eA~nditure
of experimental material.
There are many criteria which could be utilized in setting up optimality measures.
Connectedness
The connectedness principle is related to the blocking and confounding
principles.
desi~s
Unconfounded
experimer.~al
This principle
- 17 -
structures and
con~ectedness.
confounding is not
give~ ~:ier
~ecent
papers on this
princi~le
~~d
Resolvability
The principle of resolvability is to have all treatments in a complete block
~egardless
~.-hen
of how
man~.- inco~plete
~andomized
con:plete blocks
.:iesign analysis for an ex:peri.ment actually laid o..:t as an incomplete block design.
Cccasionally it Tiay be desirable to eli,.,.,inate one
Dr
analysis.
Sensitivity
The sensitivity of a statistical design first reached attention of researchers
~rhen
Fisher (1935) used the example of a lady tasting tea where the t-,.,..o kinds of
tea were milk first and then tea versus tea with milk infused.
on sensitivity since that time include papers by
3~adley
Some literature
Dar ( 1962, 1964), Lashof et al. ( 1956), Mandel and Stiehler ( 1954), and Schumann
and Bradley (1957).
For example,
measurement of time to recurrence of cancer after treatment is a much more sensitive measurement than is measuring whether or not the disease recurred.
the effect of
~
applyir~
Measuring
more effectively in samples from the 8- 10 node ~ea than in the leaves.
node area samples were able to detect a 25 lb.
inc~ease
The 8- 10
in nitrogen fertilizer
- 18 -
with the same precision as leaf samples could detect a 150 lb. increase.
Thus
the 8-10 node samples were six times as sensitive as the leaf samples in detecting additional applications of nitrogen fertilizer.
the sugar cane stalk would have been even more sensitive in measuring changes in
nitrogen levels.
Sufficiency
The sufficiency principle states that an experiment design and a given response
~odel
exper~nt.
For
example, Graybill and Weeks (1959) show that Yates' combined interblock and intrablock estimator is based on a set of minimal sufficient statistics but the set is
not coreplete.
(ii)
(iii)
(iv)
(v)
The principles of treatment design discussed below apply to each of the above
types.
- 19 opt~ality,
For fraction-
In
special cases, a fractional replicate is all that is possible, e.g., diallel cross,
tour~aments,
cult~es
etc.
pri~ciple
The replication
treatment design.
o~
Instead another
orthogonal plans leads to more precise interpretation of effects, i.e., they are
unentangled with other effects.
- 20 -
Con~ounding
When
(aliasing)
~ractions
o~
~ull
e~~ects
will necessarily be
con~ounded
con~ound e~~ects o~
o~
treatment
o~
e~~ects
e~~ects
is unknown,
which are
shoul~
o~
e~~ects.
important
This
When the
e~~ects.
should be avoided.
co~oundir~ o~ e~~ects
con~ounded,
e~~ects.
o~ con~ounding ~or
all important
o~
e~~ects
aliases.
Synmetry
The symmetry principle is to have the same number and type
of the
~actors
~actors
are
o~
o~
levels
equal importance.
~or
each
something about the response function for certain ~actors, one may be interested
in the response
less detailed
t~ue,
~or
one
~actor
in~ormation
on one or more
o~
the
~actors.
If any
o~
these were
Saturation
Saturated designs are parsimonious designs in that the number
taken is equal to the number
o~
parameters to be estimated.
o~
observations
In many situations,
o~
Hence,
o~
o~
~1
- 21 -
Efficiency
The efficiency principle for treatment design applies in various ways.
a treatment design should be selected to meet the goals of the experiment.
First,
Unless
this is done, more than one experiment may be necessary to meet the stated goals,
resulting in wasteful expenditure of resources.
or variables must be selected with care.
may not be achieved or achieved with low precision, resulting in low efficiency.
This principle is closely tied to variance optimality, which follows, when
the goal is to optimize a procedure with respect to a stated criterion.
Balance
As noted for experiment designs and by Preece (1982), there are many types
of balance for treatment designs and many uses of the term.
Others could be added from the literature on linear models (e.g., Searle, 1971).
In utilizing the balance principle, one can say that in general most uses of the
principle have some desirable properties, and it is essential that the writer
precisely defines the context in which the term balance is being used.
Otherwise,
Thus, suppose
one confounds all three-factor interactions once with blocks and all four-factor
interactions three times with blocks; this would be a balanced confounding arrangement.
Unbalanced
The
- 22 -
All
t~ends
of
a specified degree can be confounded equally, assuming that they are of equal
importance to the experimenter.
Raktoe et al. (1981) set up measures of alias optimality and alias balance.
If A is the aliasing matrix, then when the absolute value of the
minimal, this is denoted as alias optimal.
row and j'th column of A.
Let
~j
deteriTi~ant
of s is
below:
variance Optimality
In statistical literature considerable emphasis is placed upon variance
optimality or precision of estimating treatment effects.
criteria have been established.
(ii)
(iii)
(iv)
- 23 -
Unbiasedness
In the absence of knowledge about the size and nature of the bias, the unbiasedness principle would indicate that unbiased treatment designs would be
utilized.
Suppose that one has a full or complete treatment design with parametric vector
~'
!:l
E(!) = E !:2
= X~ =
!:3
xll
x12
xl3
x21
x22
x23
x3l
x32
x33
~3
~
If one takes a fraction !:1 , say, then
1
E(!l)
= [Xll
Xl2
Xl3]
= Xlll
+ Xl22 + Xl33 .
If either, or both,
~2
biased estimate of l
small.
~l'
- 24 -
[A
with respect to 1
A*]
~or
the vector
[::]
~our
situations regarding
Connectedness
As in experiment design, it is desirable that all non-zero
mable.
For the situation in the previous section, one needs to have a !l and a ! 2
that allow
estimation
0~
l and 2
be esti-
e~~ects
scienti~ic
13-2
= o,
e~~ects
'
)-l
then (x11x11
can be estimated
o~
statis-
This is a
waste~
use
o~
resources
principle states that the mean squared error (MSE), which is the variance plus
the square
measurement
o~
o~
Discussion
It seems odd that Sir Ronald A. Fisher did not relate the principles of
statistical design given in his design book (1935), viz. replication, randomization, and local control with the principles of estimation given in his statistical methods book (1925), viz. efficiency, sufficiency, and consistency.
It
Should
the above discussed principles in the real world where these designs are used?
r.'
- co -
An overview of the principles of design and estimation indicates that there are
many unanswered questions and that the philosophical aspects need fUrther study.
Literature Cited
Box, G. E. P. and J. S. Hunter (1957). Multi-factor experimental designs for
exploring response surfaces. Annals of Mathematical Statistics 28, 195-241.
Bradley, R. A. and D. E. W. Schumann (1957). The comparison of sensitivities of
similar experiments: Applications. Biometrics 13, 496-510.
Dar, S. N. (1962). On the comparison of the sensitivities of experiments.
of the Royal Statistical Society ~ 24, 447-453.
Journal
Federer,
w.
T. (1973).
799-805.
Hedayat, A. and W. T. Federer (1974). Pairwise and variance balanced incomplete
block designs. Annals of the L~stitute of Statistical Mathematics 26, 331-
338.
Khare, M. and W. T. Federer (1981). A simple construction procedure for resolvable
incomplete block designs for any number of treatments. Biometrical Journal
23(2), 121-132.
- 27 -
..
Lashof, T. w., J, Mandel, and V. Worthington (1956). Use of the sensitivity criterion for the comparison of the Bekk and Sheffield smoothness testers .
Tappi 39, 532-54 3.
Mandel, J. and R. D. Stiehler (1954). Sensitivity- A criterion for the comparison
of methods of the test. Journal of Research, Na.tio:r:al Bureau of Standards 53,
155-159.
Patterson, H. D. and E. R. Williams (1976).
block designs. Biometrika 63, 83-92.
A terminological muddle.
Another
ter~nological
tangle.
Factorial Designs.
John
Linear Models.
Shafi~
A blocking strategy.
Journal of
Communications in
Welch, W. J. (1983). A mean squared error criterion for the design of experiments.
Biometrika 70, 205-2l3.
Yates, F. (l933).
experiments.
Supplement, Journal of