Strip-Plot Designs - Crop and Soil Science

Download Report

Transcript Strip-Plot Designs - Crop and Soil Science

Strip-Plot Designs
 Sometimes called split-block design
 For experiments involving factors that are
difficult to apply to small plots
 Three sizes of plots so there are three
experimental errors
 The interaction is measured with greater
precision than the main effects
For example:
 Three seed-bed preparation methods
 Four nitrogen levels
 Both factors will be applied with large scale machinery
S3
S1
S2
S1
N1
N2
N2
N3
N0
N1
N3
N0
S3
S2
Advantages --- Disadvantages
 Advantages
– Permits efficient application of factors that would be
difficult to apply to small plots
 Disadvantages
– Differential precision in the estimation of interaction
and the main effects
– Complicated statistical analysis
Strip-Plot Analysis of Variance
Source
df
SS
MS
F
Total
rab-1
Block
r-1
SSR
MSR
A
a-1
SSA
MSA
FA
SSEA
MSEA
Factor A error
SSB
MSB
FB
(r-1)(b-1)
SSEB
MSEB
Factor B error
(a-1)(b-1)
SSAB
MSAB
FAB
SSEAB
MSEAB
Subplot error
Error(a)
B
Error(b)
AB
(r-1)(a-1)
b-1
Error(ab) (r-1)(a-1)(b-1)
SSTot
Computations
 There are three error terms - one for each main plot and
interaction plot
SSTot

 i  j  k Yijk  Y

rb   Y
SSR
ab  k Y ..k  Y
SSA

SSEA
i
i..
Y


2
2
2

b  i  k Y i .k  Y

2
 SSA  SSR
SSAB
Y  Y 
a    Y  Y   SSB  SSR
r   Y  Y   SSA  SSB
SSEAB
SSTot-SSR-SSA-SSEA-SSB-SSEB-SSAB
2
SSB
ra 
. j.
j
2
SSEB
j
k
. jk
2
i
j
ij.
F Ratios
 F ratios are computed somewhat differently
because there are three errors
 FA = MSA/MSEA
tests the sig. of the A main effect
 FB = MSB/MSEB
tests the sig. of the B main effect
 FAB = MSAB/MSEAB tests the sig. of the AB
interaction
Standard Errors of Treatment Means
 Factor A Means
M SE A
rb
 Factor B Means
M SEB
ra
 Treatment AB Means
M S E AB
r
SE of Differences for Main Effects
 Differences between 2 A means
2 * M SE A
with (r-1)(a-1) df
rb
 Differences between 2 B means
2 * M SEB
ra
with (r-1)(b-1) df
SE of Differences
 Differences between A means at same level of B
2 *   b  1  M S E A B  M S E A 
rb
 Difference between B means at same level of A
2 *   a  1  M S E A B  M S E B 
ra
 Difference between A and B means at diff. levels
2 *   a b  a  b  M S E A B  a * M S E A  b * M S E B 
ra b
For sed that are calculated from >1 MSE, t tests and df are approximated
Interpretation
Much the same as a two-factor factorial:
 First test the AB interaction
– If it is significant, the main effects have no meaning
even if they test significant
– Summarize in a two-way table of AB means
 If AB interaction is not significant
– Look at the significance of the main effects
– Summarize in one-way tables of means for factors
with significant main effects
Numerical Example
 A pasture specialist wanted to determine the
effect of phosphorus and potash fertilizers on
the dry matter production of barley to be used
as a forage
–
–
–
–
Potash: K1=none, K2=25kg/ha, K3=50kg/ha
Phosphorus: P1=25kg/ha, P2=50kg/ha
Three blocks
Farm scale fertilization equipment
K3
K1
K2
P1
56
32
49
P2
67
54
58
K1
K3
K2
P2
38
62
50
P1
52
72
64
K2
K1
K3
P2
54
44
51
P1
63
54
68
Raw data - dry matter yields
Treatment
I
II
III
P1K1
32
52
54
P1K2
49
64
63
P1K3
56
72
68
P2K1
54
38
44
P2K2
58
50
54
P2K3
67
62
51
Construct two-way tables
Potash x Block
Phosphorus x Block
K
I
II
III
Mean
P
1
43.0
45.0
49.0
45.67
1
45.67 62.67 61.67
56.67
2
53.5
57.0
58.5
56.33
2
59.67 50.00 49.67
53.11
3
61.5
67.0
59.5
62.67
Mean 52.67 56.33 55.67
54.89
Mean 52.67 56.33 55.67
I
II
III
Mean
54.89
SSR=6*devsq(range)
Potash x Phosphorus
P
Main effect of Potash
SSA=6*devsq(range)
SSEA =
2*devsq(range) – SSR – SSA
K1
K2
K3
Mean
1
46.00 58.67 65.33
56.67
2
45.33 54.00 60.00
53.11
Mean 45.67 56.33 62.67
54.89
Construct two-way tables
Potash x Block
Phosphorus x Block
K
I
II
III
Mean
P
1
43.0
45.0
49.0
45.67
1
45.67 62.67 61.67
56.67
2
53.5
57.0
58.5
56.33
2
59.67 50.00 49.67
53.11
3
61.5
67.0
59.5
62.67
Mean 52.67 56.33 55.67
54.89
Mean 52.67 56.33 55.67
I
II
III
Mean
54.89
Main effect of Phosphorus
SSB=9*devsq(range)
Potash x Phosphorus
P
K1
K2
K3
Mean
SSEB =
3*devsq(range) – SSR – SSB
1
46.00 58.67 65.33
56.67
2
45.33 54.00 60.00
53.11
SSAB=
3*devsq(range) – SSA – SSB
Mean 45.67 56.33 62.67
54.89
ANOVA
Source
df
SS
MS
Total
17
1833.78
Block
2
45.78
22.89
Potash (K)
2
885.78
442.89
Error(a)
4
78.22
19.56
Phosphorus (P) 1
56.89
56.89
Error(b)
2
693.78
346.89
KxP
2
19.11
9.56
Error(ab)
4
54.22
13.55
F
22.64**
0.16ns
0.71ns
See Excel worksheet calculations
Interpretation
Potash
None
25 kg/ha
50 kg/ha
SE
Mean Yield
45.67
56.33
62.67
1.80
 Only potash had a significant effect on barley
dry matter production
 Each increment of added potash resulted in an
increase in the yield of dry matter (~340 g/plot
per kg increase in potash
 The increase took place regardless of the level
of phosphorus
Repeated measurements over time
 We often wish to take repeated measures on experimental units to
observe trends in response over time.
– Repeated cuttings of a pasture
– Multiple harvests of a fruit or vegetable crop during a season
– Annual yield of a perennial crop
– Multiple observations on the same animal (developmental responses)
 Often provides more efficient use of resources than using different
experimental units for each time period.
 May also provide more precise estimation of time trends by reducing
random error among experimental units – effect is similar to blocking
 Problem: observations over time are not assigned at random to
experimental units.
– Observations on the same plot will tend to be positively correlated
– Violates the assumption that errors (residuals) are independent
Analysis of repeated measurements
 The simplest approach is to treat sampling times
as sub-plots in a split-plot experiment.
– Some references recommend use of strip-plot rather
than a split-plot
 Univariate adjustments can be made
 Multivariate procedures can be used to adjust for
the correlations among sampling periods
 Mixed Model approaches can be used to adjust
for the correlations among sampling periods
Split-plot in time
 In a sense, a split-plot is a specific case of repeated
measures, where sub-plots represent repeated
measurements on a common main plot
 Analysis as a split-plot is valid only if all pairs of sub-plots
in each main plot can be assumed to be equally
correlated
– Compound symmetry
– Sphericity
Formal names for required assumptions
 When time is a sub-plot, correlations may be greatest for
samples taken at short time intervals and less for distant
sampling periods, so assumptions may not be valid
– Not a problem when there are only two sampling periods
Univariate adjustments for repeated measures
 Fit a smooth curve to the time trends and analyze a
derived variable
– average
– maximum response
– area under curve
– time to reach the maximum
 Use polynomial contrasts to evaluate trends over time
(linear, quadratic responses) and compare responses for
each treatment
 Reduce df for subplots, interactions, and subplot error
terms to obtain more conservative F tests
Multivariate adjustments for repeated measures
 In PROC GLM, each repeated measure is treated like an
additional variable in a multivariate analysis:
model yield1 yield2 yield3 yield4=variety/nouni;
repeated harvest / printe;
 MANOVA approach is very conservative
– Effectively controls Type I error
– Power may be low
• Many parameters are estimated so df for error may be too low
• Missing values result in an unnecessary loss of available
information
 No real benefit compared to a Mixed Model approach
Covariance Structure for Residuals

 Y  Y   Y   Y  2  Y ,Y
2
2
1
2
sed2
Y  Y
2
2
1
2
2
1
se2
2
se2
 Y ,Y  0
1
2
1
2
covariance
 Y Y 
2
1
2
2M S E
r
Covariance Structure for Residuals
 No correlation (independence)
– 4 measurements per subject
2

– All covariances = 0
0
  
 0

 0
0
0

0
2
0

0
0
2
0 
1


0
0
  2 
0
0 

2
 
0
0
0
1
0
0
1
0
0
0

0

0

1
 Compound symmetry (CS)
– All covariances (off-diagonal elements) are the same
– Often applies for split-plot designs (sub-plots within main plots are equally
correlated)
1


2

  





1


1








1
Covariance Structure for Residuals
 Autoregressive (AR)
– Applies to time series analyses
– For a first-order AR(1) structure, the
within subject correlations drop off
exponentially as the number of time
lags between measurements
increases (assuming time lags are all
the same)
 1
Complex and computer intensive

 12
No particular pattern for the
2
 
covariances is assumed
  13
May have low power due to loss of df

for error
  14
 Unstructured (UN)
–
–
–
1

2  
 
 2
 3



1


1


2
 12
 13
1
 23
 23
1
 24
 34
2
3
 
2




1
 14 

 24

 34 

1 
Mixed Model adjustment for error structure
 Stage one: estimate covariance structure for residuals
1. Determine which covariance structures would make sense for the
experimental design and type of data that is collected
2. Use graphical methods to examine covariance patterns over time
3. Likelihood ratio tests of more complex vs simpler models
4. Information content

2
= (-2 res log likelihood)simple model
minus (-2 res log likelihood)complex model
df = difference in # parameters estimated
AIC, AICC, BIC – information content
adjust for loss in power due to loss of df in
more complex models
Null model - no adjustment for correlated errors
Mixed Model adjustment for error structure
 Stage two:
– include appropriate covariance structure in the model
– use Generalized Least Squares methodology to evaluate
treatment and time effects
 Computer intensive
– use PROC MIXED or GLIMMIX in SAS