Lecture 7 Randomized Complete Block Design (RCBD) [ST&D sections (except 9.6) and section 15.8]
|
|
- Rosalind Payne
- 6 years ago
- Views:
Transcription
1 Lecture 7 Randomized Complete Block Design () [ST&D sections (except 9.6) and section 15.8] The Completely Randomized Design () 1. It is assumed that all experimental units (EU's) are uniform.. Treatments are randomly assigned to EUs such that each treatment occurs equally often in the experiment. (1 randomization per experiment) 3. It is advocated to include as much of the native variability of the experiment as possible within each EU. 4. When EU's are not uniform, experimental error () increases, F (MST/) decreases, and the experiment loses sensitivity. If the experiment is replicated in a variety of situations to increase its scope, the variability increases even further. The Randomized Complete Block Design () 1. The population of EU's is divided into a number of relatively homogeneous subpopulations or blocks, and it is assumed that all EU's within a given block are uniform.. Within each block, treatments are randomly assigned to EU's such that each treatment occurs equally often (usually once) in each block. (1 randomization per block) 3. It is advocated to minimize the native variability as much as possible within blocks and to maximize the native variability as much as possible among blocks. 4. Variation among blocks can be partitioned out of the experimental error (), thereby reducing this quantity and increasing the power of the test. Additional variability introduced when increasing the scope of the experiment can also be partitioned out of the. Blocks usually represent levels of naturally-occurring differences or sources of variation that are unrelated to the treatments, and the characterization of these differences is not of interest to the researcher. 1
2 Example: A field trial comparing three cultivars (A, B, and C) of sugar beet with four replications. North end of field 1 3 Hi N South end of field Low N : One randomization per experiment North end of field B C B Hi N C A C A B A C B A South end of field Low N : One randomization per block Block North end of field Hi N South end of field Low N Block North end of field Hi N 1 B A C A B C 3 A C B 4 A C B South end of field Low N
3 The linear model The model underlying each observation in the experiment: Y And the sum of squares: t r Y ij = µ + τ i + β j + ε ij ij = Y.. ij + Y + ( Y i. Y..) + ( Y. j Y..) + ( Y Y i. Y. j..) t ( Yij Y ) = r ( Y i. Y..) + t ( Y. j Y..) +.. ( Y Y i. Y. j+ Y..) i= 1 j= 1 i= 1 j= 1 i= 1 j= 1 r TSS = SST + SSB + SSE t r ij This partitioning of variance is possible because the sums of squares of treatments, blocks, and error are orthogonal to one another. This orthogonality is a direct result of the completeness of the block design. Source df SS MS F Total rt - 1 TSS Treatments t - 1 SST SST/(t-1) MST/ Error t(r - 1) TSS-SST SSE/r(t-1) (one replication per block-treatment combination) Source df SS MS F Total rt - 1 TSS Treatments t - 1 SST SST/(t-1) MST/ Blocks r - 1 SSB SSB/(r-1) Error (t-1)(r-1) TSS-SST-SSB SSE/(t-1)(r-1) 1. has (r - 1) fewer df e than the.. If there are no differences among blocks (SSB = 0), <. 3. If there are large enough differences among blocks (SSB >> 0), >. 3
4 Example: An experiment was conducted to investigate the effect of estrogen on weight gain in sheep. The treatments are combinations of sex of sheep (M, F) and level of estrogen (Est 0, Est 3 ). The sheep are blocked by ranch, with one replication of each treatment level at each ranch. Ranch Trtmt M Est 0 M Est 3 F Est 0 F Est 3 Effect of estrogen on weight gain in sheep (lbs). Ranch (i.e. block) Treatment Treatment I II III IV Total Mean F-S M-S F-S M-S Block Total Block Mean ANOVA (treating blocks as reps) Source df SS MS F Totals Treatment NS Error ANOVA Source df SS MS F Total Treatment ** Blocks ** Error
5 Relative efficiency: When to block? F = MST = SSE F crit = Fα, df df trt, df e e Blocking reduces SSE, which reduces. Blocking reduces df e, which increases and increases F crit. The concept of relative efficiency formalizes the comparison between two experimental methods by quantifying this balance between loss of degrees of freedom and reduction in experimental error. The information per replication in a given design is: I = 1 df df 3 + ε RE 1: = I I 1 df df df df ( df = ( df 1 + 1)( df + 1)( df 1 + 3) + 3) 1 The main complication is how to estimate for the alternative design. If an experiment was conducted as an, can be estimated by the following formula (ST&D ): ˆ df B MSB + ( df df + df B T T + df + df e e ) Assume TSS of the two designs is the same. Rewrite TSS in terms of its components and simplify the expression. 5
6 For the interested: Derivation of the expected 1. Set the total sums of squares of each design equal to each other and rewrite in terms of mean squares and degrees of freedom: TSS SST df T ( MST R ( t 1) MST = TSS + SSB R + df B( + SSE MSB R + ( r 1) MSB R + df = SST R + SSE = df T ( C) + ( t 1)( r 1) R MST C + df C) = ( t 1) MST C C + t( r 1) C. Replace each mean square with the variance components of its expected mean square: ( t 1)( + r T ( ) + ( r 1)( [( t 1) + ( r 1) + ( t 1)( r 1)] t( r 1) C) = B( + ( tr 1) t( r 1) B( + ( tr 1) + t B( + t( r 1) = ( tr 1) C) ) + ( t 1)( r 1) B( + r( t 1) T ( = ( t 1)( C) + r T ( C ) = [( t 1) + t( r 1)] C) ) + t( r 1) + r( t 1) C ) T ( C) 3. Finally, rewrite this expression in terms of mean squares and degrees of freedom: = MSB + t( r 1) t( tr 1) MSB = + ( r 1) ( r 1) tr 1 tr 1 MSB = [( tr 1) ( r 1)] + ( r 1) tr 1 tr 1 r( t 1) + ( r 1) MSB = tr 1 ( dft ( + df ) + df BMSB = df + df + df T ( B( 6
7 Example: From the sheep experiment, = 7.78 and MSB = Therefore: ˆ df BMSB + ( dft + dfe) 3(19) + (3+ 9)7.78 M SE = = 44.6 df + df + df B T e RE ( df ( df + 1)( df + 1)( df + 3) ˆ + 3) (9+ 1)(1+ 3)44.6 = (1+ 1)(9+ 3) : = Interpretation: It takes 5.51 replications in the to produce the same amount of information as one replication in the. Or, the is 5.51 time more efficient than the in this case. 1. When there are no significant differences among blocks distribution of the difference between two means () distribution of the difference between two means () α/ SSE = SSE > 0 Power < Power Power 1 β β 7
8 . When there are significant differences among blocks distribution of the difference between two means () α/ distribution of the difference between two means () 0 df e < df e SSE < SSE Power > Power or Power < Power Power 1 β β Assumptions of the model The model for the with a single replication per block-treatment combination: Y ij = µ + τ i + β j + ε ij 1. The residuals (ε ij ) are independent, homogeneous, and normally distributed.. The variance within each treatment levels is homogeneous across all treatment levels. 3. The main effects are additive. 8
9 Recall that experimental error is defined as the variation among experimental units that are treated alike. Ranch Trtmt M Est 0 M Est 3 F Est 0 F Est 3 There is an expected value for each sheep, given by: Expected Y ij = µ + τ i + β j Observed Y ij = µ + τ i + β j + ε ij With only one replication per cell, the residuals are the combined effects of experimental error and any non-additive treatment*block interactions: ε ij = τ i *β j + error ij So when we use ε ij as estimates of the true experimental error, we are assuming that τ i *β j = 0. This assumption of no interaction is referred to as the assumption of additivity of the main effects. If this assumption is violated, all F-tests will be very inefficient and possibly misleading, particularly if the interaction effect is very large. 9
10 Example: Consider the hypothetical effects of two factors (A, B) on some response variable: Purely additive Factor A Factor B τ 1 = 1 τ = τ 3 = 3 β 1 = β = Purely multiplicative Factor A Factor B τ 1 = 1 τ = τ 3 = 3 β 1 = β = If multiplicative data of this sort are analyzed by a conventional ANOVA, the interaction SS will be large due to the nonadditivity of the treatment effects. In the case of multiplicative effects, there is a simple remedy. Simply transform the variable by taking the log of each mean to restore additivity: Restoring additivity (log of data) Factor A Factor B τ 1 = 0.00 τ = 0.30 τ 3 = 0.48 β 1 = β =
11 Tukey s 1-df test for nonadditivity [ST&D 395] Under our linear model, each observation is characterized as: yij = µ + β + τ + ε i j ij The predicted value of each individual is given by: pred ij = µ + β + τ i j So, if we had no error in our experiment (i.e. if ε = 0 ), the observed data would exactly match its predicted values and a correlation plot of the two would yield a perfect line with slope = 1: ij Observed vs. Predicted Values (, no error) 0 18 Observed Predicted Now let's introduce some error: Observed vs. Predicted (, with error) 0 18 Observed Predicted 11
12 But what happens when you have an interaction (e.g. Block * Treatment) but lack the degrees of freedom necessary to include it in the linear model? ε ε + B*T Interaction Effects ij = RANDOMij Observed vs. Predicted (, with error and B*T) 0 18 Observed Predicted SO, if the observed and predicted values obey a linear relationship, then the nonrandom Interaction Effects buried in the error term are sufficiently small to uphold our assumption of additivity of main effects. This test is easily implemented using SAS: Data Lambs; Input Sex_Est Do Ranch = 1 to 4; Input Output; End; Cards; F M F M ; Proc GLM; Class Sex_Est Ranch; Model Gain = Ranch Sex_Est; Output out = LambsPR p = Pred r = Res; Proc GLM; Class Sex_Est Ranch; Model Gain = Ranch Sex_Est Pred*Pred; Run; Quit; * This is the Tukey 1 df test; NOTE: Pred*Pred is not in the Class statement, and it is the last term in the Model. 1
13 Output from the nd Proc GLM (the Tukey 1 df test): Source DF Type III SS Mean Square F Value Pr > F Ranch Sex_Est Pred*Pred NS This test is necessary ONLY when there is one observation per block-treatment combination. If there are two or more replications per block-treatment combination, the block*treatment interaction can be tested directly in an exploratory model. Example: Yield of penicillin from four different protocols (A D). Blocks are different stocks of an important reagent. The numbers below each observation (O) are the predicted values (P = Grand Mean + Treatment effect + Block effect) and the residuals (. Block Stock 1 Stock Stock 3 Stock 4 Stock 5 Treatment A B C D O: 88 O: 97 O: 94 P: 91 P: 95 P: 9 R: -3 R: R: O: 89 P: 90 R: -1 O: 84 P: 81 R: 3 O: 81 P: 83 R: - O: 87 P: 86 R: 1 O: 79 P: 80 R: -1 O: 77 P: 8 R: -5 O: 87 P: 84 R: 3 O: 9 P: 87 R: 5 O: 81 P: 81 R: 0 O: 9 P: 86 R: 6 O: 87 P: 88 R: -1 O: 89 P: 91 R: - O: 80 P: 85 R: -5 O: 79 P: 83 R: -4 O: 85 P: 85 R: 0 O: 84 P: 88 R: -4 O: 88 P: 8 R: 6 Treatment mean Treatment effect Block Mean Block Effect Mean = 86 13
14 The SAS code for such an analysis: Data Penicillin; Do Block = 1 to 4; Do Trtmt = 1 to 4; Input Output; End; End; Cards; ; Proc GLM Data = Penicillin; Class Block Trtmt; Model Yield = Block Trtmt; Output out = PenPR p = Pred r = Res; Proc Univariate Data = PenPR normal; * Testing for normality of residuals; Var Res; Proc GLM Data = Penicillin; * Testing for variance homogeneity (1 way ANOVA); Class Trtmt; Model Yield = Trtmt; Means Trtmt / hovtest = Levene; Proc Plot Data = PenPR; * Generating a plot of predicted vs. residual values; Plot Res*Pred = Trtmt; Proc GLM Data = PenPR; * Testing for nonadditivity; Class Block Trtmt; Model Yield = Block Trtmt Pred*Pred; Run; Quit; This dataset meets all assumptions: normality, variance homogeneity, and additivity: Tests for Normality Test --Statistic p Value Shapiro-Wilk W Pr < W NS Levene's Test for Homogeneity of Yield Variance ANOVA of Squared Deviations from Group Means Sum of Mean Source DF Squares Square F Value Pr > F Trtmt NS Error Tukey's 1 df test for nonadditivity Source DF Type III SS Mean Square F Value Pr > F Block Trtmt Pred*Pred NS 14
15 No particular pattern presents itself in this graph. 15
16 Nesting within an Data Lambs; Input Sex_Est $ Ranch Animal Cards; f f f f m m m m f f3 1 5 f f m m m m f f f f m m m m f f f f m m m m ; Proc GLM Data = Lambs Order = Data; Class Ranch Sex_Est Animal; Model Gain = Ranch Sex_Est Animal(Ranch*Sex_Est); Random Animal(Ranch*Sex_Est); Test h = Sex_Est e = Animal(Ranch*Sex_Est); Contrast 'sex' Sex_Est / e = Animal(Ranch*Sex_Est); Contrast 'estrogen' Sex_Est / e = Animal(Ranch*Sex_Est); Contrast 'interaction' Sex_Est / e = Animal(Ranch*Sex_Est); Means Sex_Est / Tukey e = Animal(Ranch*Sex_Est); Proc Varcomp Method = Type1; Class Ranch Sex_Est Animal; Model Gain = Ranch Sex_Est Animal(Ranch*Sex_Est); Run; Quit; Ranch Trtmt M Est 0 M Est 3 F Est 0 F Est 3 measurements Type 1 Analysis of Variance Source Ranch Sex_Est Animal(Ranch*Sex_Es) Error Expected Mean Square Var(Error) + Var(Animal(Ranch*Sex_Es)) + 8 Var(Ranch) Var(Error) + Var(Animal(Ranch*Sex_Es)) + 8 Var(Sex_Est) Var(Error) + Var(Animal(Ranch*Sex_Es)) Var(Error) Type 1 Estimates Variance Component Estimate % Var(Ranch) Var(Sex_Est) Var(Animal(Ranch*Sex_Es)) Var(Error) The only reason to analyze this dataset as a nested is to calculate the variance components. If you do not need the variance components, simply average the subsamples for each experimental unit and analyze it as a simple. 16
17 with multiple replications per block-treatment combination Data Lambs; Do Ranch = 1 to 4; Do Sex_Est = 'M0', 'M3', 'F0', 'F3'; Input Output; End; End; Cards; Ranch Trtmt M Est 0 M Est 3 F Est 0 F Est ; Proc GLM Data = Lambs Order = Data; Class Ranch Sex_Est; Model Gain = Ranch Sex_Est Ranch*Sex_Est; e * The exploratory model; Proc GLM Data = Lambs Order = Data; Class Ranch Sex_Est; Model Gain = Ranch Sex_Est; Output out = LambsPR p = Pred r = Res; Contrast 'sex' Sex_Est ; Contrast 'estrogen' Sex_Est ; Contrast 'interaction' Sex_Est ; Means Sex_Est / Tukey; * The ANOVA; Proc Univariate Data = LambsPR normal; Var Res; Proc GLM Data = Lambs Order = Data; Class Sex_Est; Model Gain = Sex_Est; Means Sex_Est / hovtest = Levene; * Testing normality; * Levene's test (1 way ANOVA); Run; Quit; Output from the exploratory model: Sum of Source DF Squares Mean Square F Value Pr > F Model Error Corrected Total Source DF Type III SS Mean Square F Value Pr > F Ranch Sex_Est <.0001 Ranch*Sex_Est NS 17
18 1 rep/cell 1 rep/cell with subsamples T1 T T1 T B1 B1 B B Class Block Trtmt; Model Y = Block Trtmt; Tukey Test Required Class Block Trtmt Pot; Model Y = Block Trtmt Pot(Block*Trtmt); Random Pot(Block*Trtmt); Test h = Trtmt e = Pot(Block*Trtmt); Tukey Test Required >1 rep/cell >1 rep/cell with subsamples T1 T T1 T B1 B1 B B Exploratory model: Class Block Trtmt; Model Y = Block Trtmt Block*Trtmt; Tukey Test not Required ANOVA: Class Block Trtmt; Model Y = Block Trtmt; Exploratory model: Class Bl Trt Pot; Model Y = Bl Trt Bl*Trt Pot(Bl*Trt); Random Pot(Bl*Trt); Test h = Bl*Trt e = Pot(Bl*Trt); Tukey Test not Required ANOVA: Class Bl Trt Pot; Model Y = Bl Trt Pot(Bl*Trt); Random Pot(Bl*Trt); Test h = Trt e = Pot(Bl*Trt); 18
Topic 6. Two-way designs: Randomized Complete Block Design [ST&D Chapter 9 sections 9.1 to 9.7 (except 9.6) and section 15.8]
Topic 6. Two-way designs: Randomized Complete Block Design [ST&D Chapter 9 sections 9.1 to 9.7 (except 9.6) and section 15.8] The completely randomized design Treatments are randomly assigned to e.u. such
More information2. Treatments are randomly assigned to EUs such that each treatment occurs equally often in the experiment. (1 randomization per experiment)
Lecture 9 Topic 6: Randomized Complete Block Designs (RCBD's) The Completely Randomized Design () 1. It is assumed that all experimental units (EU's) are uniform.. Treatments are randomly assigned to EUs
More informationPLS205 Lab 2 January 15, Laboratory Topic 3
PLS205 Lab 2 January 15, 2015 Laboratory Topic 3 General format of ANOVA in SAS Testing the assumption of homogeneity of variances by "/hovtest" by ANOVA of squared residuals Proc Power for ANOVA One-way
More informationPLS205 Lab 6 February 13, Laboratory Topic 9
PLS205 Lab 6 February 13, 2014 Laboratory Topic 9 A word about factorials Specifying interactions among factorial effects in SAS The relationship between factors and treatment Interpreting results of an
More informationPLS205 Winter Homework Topic 8
PLS205 Winter 2015 Homework Topic 8 Due TUESDAY, February 10, at the beginning of discussion. Answer all parts of the questions completely, and clearly document the procedures used in each exercise. To
More informationTopic 13. Analysis of Covariance (ANCOVA) [ST&D chapter 17] 13.1 Introduction Review of regression concepts
Topic 13. Analysis of Covariance (ANCOVA) [ST&D chapter 17] 13.1 Introduction The analysis of covariance (ANCOVA) is a technique that is occasionally useful for improving the precision of an experiment.
More informationPLS205!! Lab 9!! March 6, Topic 13: Covariance Analysis
PLS205!! Lab 9!! March 6, 2014 Topic 13: Covariance Analysis Covariable as a tool for increasing precision Carrying out a full ANCOVA Testing ANOVA assumptions Happiness! Covariable as a Tool for Increasing
More informationPLS205 KEY Winter Homework Topic 3. The following represents one way to program SAS for this question:
PL05 KEY Winter 05 Homework Topic 3 Answers to Question [30 points] The following represents one way to program A for this question: Data Weeds; Input Cover $ Biomass; Cards; 66 634 63 63 633 645 69 63
More informationTopic 13. Analysis of Covariance (ANCOVA) - Part II [ST&D Ch. 17]
Topic 13. Analysis of Covariance (ANCOVA) - Part II [ST&D Ch. 17] 13.5 Assumptions of ANCOVA The assumptions of analysis of covariance are: 1. The X s are fixed, measured without error, and independent
More informationTopic 8. Data Transformations [ST&D section 9.16]
Topic 8. Data Transformations [ST&D section 9.16] 8.1 The assumptions of ANOVA For ANOVA, the linear model for the RCBD is: Y ij = µ + τ i + β j + ε ij There are four key assumptions implicit in this model.
More informationTopic 9: Factorial treatment structures. Introduction. Terminology. Example of a 2x2 factorial
Topic 9: Factorial treatment structures Introduction A common objective in research is to investigate the effect of each of a number of variables, or factors, on some response variable. In earlier times,
More informationThe ε ij (i.e. the errors or residuals) are normally distributed. This assumption has the least influence on the F test.
Lecture 11 Topic 8: Data Transformations Assumptions of the Analysis of Variance 1. Independence of errors The ε ij (i.e. the errors or residuals) are statistically independent from one another. Failure
More informationTopic 12. The Split-plot Design and its Relatives (continued) Repeated Measures
12.1 Topic 12. The Split-plot Design and its Relatives (continued) Repeated Measures 12.9 Repeated measures analysis Sometimes researchers make multiple measurements on the same experimental unit. We have
More informationCOMPLETELY RANDOM DESIGN (CRD) -Design can be used when experimental units are essentially homogeneous.
COMPLETELY RANDOM DESIGN (CRD) Description of the Design -Simplest design to use. -Design can be used when experimental units are essentially homogeneous. -Because of the homogeneity requirement, it may
More informationAssessing Model Adequacy
Assessing Model Adequacy A number of assumptions were made about the model, and these need to be verified in order to use the model for inferences. In cases where some assumptions are violated, there are
More informationIncreasing precision by partitioning the error sum of squares: Blocking: SSE (CRD) à SSB + SSE (RCBD) Contrasts: SST à (t 1) orthogonal contrasts
Lecture 13 Topic 9: Factorial treatment structures (Part II) Increasing precision by partitioning the error sum of squares: s MST F = = MSE 2 among = s 2 within SST df trt SSE df e Blocking: SSE (CRD)
More informationIX. Complete Block Designs (CBD s)
IX. Complete Block Designs (CBD s) A.Background Noise Factors nuisance factors whose values can be controlled within the context of the experiment but not outside the context of the experiment Covariates
More informationOutline Topic 21 - Two Factor ANOVA
Outline Topic 21 - Two Factor ANOVA Data Model Parameter Estimates - Fall 2013 Equal Sample Size One replicate per cell Unequal Sample size Topic 21 2 Overview Now have two factors (A and B) Suppose each
More informationThe legacy of Sir Ronald A. Fisher. Fisher s three fundamental principles: local control, replication, and randomization.
1 Chapter 1: Research Design Principles The legacy of Sir Ronald A. Fisher. Fisher s three fundamental principles: local control, replication, and randomization. 2 Chapter 2: Completely Randomized Design
More informationLecture 3. Experiments with a Single Factor: ANOVA Montgomery 3-1 through 3-3
Lecture 3. Experiments with a Single Factor: ANOVA Montgomery 3-1 through 3-3 Page 1 Tensile Strength Experiment Investigate the tensile strength of a new synthetic fiber. The factor is the weight percent
More informationSTAT 115:Experimental Designs
STAT 115:Experimental Designs Josefina V. Almeda 2013 Multisample inference: Analysis of Variance 1 Learning Objectives 1. Describe Analysis of Variance (ANOVA) 2. Explain the Rationale of ANOVA 3. Compare
More informationLecture 4 Topic 3: General linear models (GLMs), the fundamentals of the analysis of variance (ANOVA), and completely randomized designs (CRDs)
Lecture 4 Topic 3: General linear model (GLM), the fundamental of the analyi of variance (ANOVA), and completely randomized deign (CRD) The general linear model One population: An obervation i explained
More informationAllow the investigation of the effects of a number of variables on some response
Lecture 12 Topic 9: Factorial treatment structures (Part I) Factorial experiments Allow the investigation of the effects of a number of variables on some response in a highly efficient manner, and in a
More informationTopic 17 - Single Factor Analysis of Variance. Outline. One-way ANOVA. The Data / Notation. One way ANOVA Cell means model Factor effects model
Topic 17 - Single Factor Analysis of Variance - Fall 2013 One way ANOVA Cell means model Factor effects model Outline Topic 17 2 One-way ANOVA Response variable Y is continuous Explanatory variable is
More informationLecture 3. Experiments with a Single Factor: ANOVA Montgomery 3.1 through 3.3
Lecture 3. Experiments with a Single Factor: ANOVA Montgomery 3.1 through 3.3 Fall, 2013 Page 1 Tensile Strength Experiment Investigate the tensile strength of a new synthetic fiber. The factor is the
More informationComparison of a Population Means
Analysis of Variance Interested in comparing Several treatments Several levels of one treatment Comparison of a Population Means Could do numerous two-sample t-tests but... ANOVA provides method of joint
More informationTopic 12. The Split-plot Design and its Relatives (Part II) Repeated Measures [ST&D Ch. 16] 12.9 Repeated measures analysis
Topic 12. The Split-plot Design and its Relatives (Part II) Repeated Measures [ST&D Ch. 16] 12.9 Repeated measures analysis Sometimes researchers make multiple measurements on the same experimental unit.
More informationSingle Factor Experiments
Single Factor Experiments Bruce A Craig Department of Statistics Purdue University STAT 514 Topic 4 1 Analysis of Variance Suppose you are interested in comparing either a different treatments a levels
More informationFactorial designs. Experiments
Chapter 5: Factorial designs Petter Mostad mostad@chalmers.se Experiments Actively making changes and observing the result, to find causal relationships. Many types of experimental plans Measuring response
More informationLecture 3: Inference in SLR
Lecture 3: Inference in SLR STAT 51 Spring 011 Background Reading KNNL:.1.6 3-1 Topic Overview This topic will cover: Review of hypothesis testing Inference about 1 Inference about 0 Confidence Intervals
More informationStatistics For Economics & Business
Statistics For Economics & Business Analysis of Variance In this chapter, you learn: Learning Objectives The basic concepts of experimental design How to use one-way analysis of variance to test for differences
More informationChap The McGraw-Hill Companies, Inc. All rights reserved.
11 pter11 Chap Analysis of Variance Overview of ANOVA Multiple Comparisons Tests for Homogeneity of Variances Two-Factor ANOVA Without Replication General Linear Model Experimental Design: An Overview
More informationExample: Four levels of herbicide strength in an experiment on dry weight of treated plants.
The idea of ANOVA Reminders: A factor is a variable that can take one of several levels used to differentiate one group from another. An experiment has a one-way, or completely randomized, design if several
More informationLecture 14 Topic 10: ANOVA models for random and mixed effects. Fixed and Random Models in One-way Classification Experiments
Lecture 14 Topic 10: ANOVA models for random and mixed effects To this point, we have considered only the Fixed Model (Model I) ANOVA; now we will extend the method of ANOVA to other experimental objectives.
More informationCHAPTER 4 Analysis of Variance. One-way ANOVA Two-way ANOVA i) Two way ANOVA without replication ii) Two way ANOVA with replication
CHAPTER 4 Analysis of Variance One-way ANOVA Two-way ANOVA i) Two way ANOVA without replication ii) Two way ANOVA with replication 1 Introduction In this chapter, expand the idea of hypothesis tests. We
More informationLecture 7: Latin Square and Related Design
Lecture 7: Latin Square and Related Design Montgomery: Section 4.2-4.3 Page 1 Automobile Emission Experiment Four cars and four drivers are employed in a study for possible differences between four gasoline
More informationLecture 4. Random Effects in Completely Randomized Design
Lecture 4. Random Effects in Completely Randomized Design Montgomery: 3.9, 13.1 and 13.7 1 Lecture 4 Page 1 Random Effects vs Fixed Effects Consider factor with numerous possible levels Want to draw inference
More informationTwo-factor studies. STAT 525 Chapter 19 and 20. Professor Olga Vitek
Two-factor studies STAT 525 Chapter 19 and 20 Professor Olga Vitek December 2, 2010 19 Overview Now have two factors (A and B) Suppose each factor has two levels Could analyze as one factor with 4 levels
More informationTopic 20: Single Factor Analysis of Variance
Topic 20: Single Factor Analysis of Variance Outline Single factor Analysis of Variance One set of treatments Cell means model Factor effects model Link to linear regression using indicator explanatory
More informationAnalysis of Variance and Design of Experiments-I
Analysis of Variance and Design of Experiments-I MODULE VIII LECTURE - 35 ANALYSIS OF VARIANCE IN RANDOM-EFFECTS MODEL AND MIXED-EFFECTS MODEL Dr. Shalabh Department of Mathematics and Statistics Indian
More information5.3 Three-Stage Nested Design Example
5.3 Three-Stage Nested Design Example A researcher designs an experiment to study the of a metal alloy. A three-stage nested design was conducted that included Two alloy chemistry compositions. Three ovens
More informationdf=degrees of freedom = n - 1
One sample t-test test of the mean Assumptions: Independent, random samples Approximately normal distribution (from intro class: σ is unknown, need to calculate and use s (sample standard deviation)) Hypotheses:
More informationEstimating σ 2. We can do simple prediction of Y and estimation of the mean of Y at any value of X.
Estimating σ 2 We can do simple prediction of Y and estimation of the mean of Y at any value of X. To perform inferences about our regression line, we must estimate σ 2, the variance of the error term.
More informationLecture 4. Checking Model Adequacy
Lecture 4. Checking Model Adequacy Montgomery: 3-4, 15-1.1 Page 1 Model Checking and Diagnostics Model Assumptions 1 Model is correct 2 Independent observations 3 Errors normally distributed 4 Constant
More informationLecture 15 Topic 11: Unbalanced Designs (missing data)
Lecture 15 Topic 11: Unbalanced Designs (missing data) In the real world, things fall apart: plants are destroyed/trampled/eaten animals get sick volunteers quit assistants are sloppy accidents happen
More information4.1. Introduction: Comparing Means
4. Analysis of Variance (ANOVA) 4.1. Introduction: Comparing Means Consider the problem of testing H 0 : µ 1 = µ 2 against H 1 : µ 1 µ 2 in two independent samples of two different populations of possibly
More informationTopic 7: Incomplete, double-blocked designs: Latin Squares [ST&D sections ]
Topic 7: Incomplete, double-blocked designs: Latin Squares [ST&D sections 9.10 9.15] 7.1. Introduction The Randomized Complete Block Design is commonly used to improve the ability of an experiment to detect
More informationThe Random Effects Model Introduction
The Random Effects Model Introduction Sometimes, treatments included in experiment are randomly chosen from set of all possible treatments. Conclusions from such experiment can then be generalized to other
More informationWELCOME! Lecture 13 Thommy Perlinger
Quantitative Methods II WELCOME! Lecture 13 Thommy Perlinger Parametrical tests (tests for the mean) Nature and number of variables One-way vs. two-way ANOVA One-way ANOVA Y X 1 1 One dependent variable
More informationEXST 7015 Fall 2014 Lab 11: Randomized Block Design and Nested Design
EXST 7015 Fall 2014 Lab 11: Randomized Block Design and Nested Design OBJECTIVES: The objective of an experimental design is to provide the maximum amount of reliable information at the minimum cost. In
More informationBALANCED INCOMPLETE BLOCK DESIGNS
BALANCED INCOMPLETE BLOCK DESIGNS V.K. Sharma I.A.S.R.I., Library Avenue, New Delhi -110012. 1. Introduction In Incomplete block designs, as their name implies, the block size is less than the number of
More informationIf we have many sets of populations, we may compare the means of populations in each set with one experiment.
Statistical Methods in Business Lecture 3. Factorial Design: If we have many sets of populations we may compare the means of populations in each set with one experiment. Assume we have two factors with
More informationSTOR 455 STATISTICAL METHODS I
STOR 455 STATISTICAL METHODS I Jan Hannig Mul9variate Regression Y=X β + ε X is a regression matrix, β is a vector of parameters and ε are independent N(0,σ) Es9mated parameters b=(x X) - 1 X Y Predicted
More informationANALYSIS OF VARIANCE OF BALANCED DAIRY SCIENCE DATA USING SAS
ANALYSIS OF VARIANCE OF BALANCED DAIRY SCIENCE DATA USING SAS Ravinder Malhotra and Vipul Sharma National Dairy Research Institute, Karnal-132001 The most common use of statistics in dairy science is testing
More informationInference for Regression Simple Linear Regression
Inference for Regression Simple Linear Regression IPS Chapter 10.1 2009 W.H. Freeman and Company Objectives (IPS Chapter 10.1) Simple linear regression p Statistical model for linear regression p Estimating
More informationLecture 6: Single-classification multivariate ANOVA (k-group( MANOVA)
Lecture 6: Single-classification multivariate ANOVA (k-group( MANOVA) Rationale and MANOVA test statistics underlying principles MANOVA assumptions Univariate ANOVA Planned and unplanned Multivariate ANOVA
More informationLecture 5: Comparing Treatment Means Montgomery: Section 3-5
Lecture 5: Comparing Treatment Means Montgomery: Section 3-5 Page 1 Linear Combination of Means ANOVA: y ij = µ + τ i + ɛ ij = µ i + ɛ ij Linear combination: L = c 1 µ 1 + c 1 µ 2 +...+ c a µ a = a i=1
More informationChapter 11 - Lecture 1 Single Factor ANOVA
Chapter 11 - Lecture 1 Single Factor ANOVA April 7th, 2010 Means Variance Sum of Squares Review In Chapter 9 we have seen how to make hypothesis testing for one population mean. In Chapter 10 we have seen
More informationANOVA (Analysis of Variance) output RLS 11/20/2016
ANOVA (Analysis of Variance) output RLS 11/20/2016 1. Analysis of Variance (ANOVA) The goal of ANOVA is to see if the variation in the data can explain enough to see if there are differences in the means.
More informationY it = µ + τ i + ε it ; ε it ~ N(0, σ 2 ) (1)
DOE (Linear Model) Strategy for checking experimental model assumptions Part 1 When we discuss experiments whose data are described or analyzed by the one-way analysis of variance ANOVA model, we note
More informationRCB - Example. STA305 week 10 1
RCB - Example An accounting firm wants to select training program for its auditors who conduct statistical sampling as part of their job. Three training methods are under consideration: home study, presentations
More informationOutline. Topic 22 - Interaction in Two Factor ANOVA. Interaction Not Significant. General Plan
Topic 22 - Interaction in Two Factor ANOVA - Fall 2013 Outline Strategies for Analysis when interaction not present when interaction present when n ij = 1 when factor(s) quantitative Topic 22 2 General
More informationWITHIN-PARTICIPANT EXPERIMENTAL DESIGNS
1 WITHIN-PARTICIPANT EXPERIMENTAL DESIGNS I. Single-factor designs: the model is: yij i j ij ij where: yij score for person j under treatment level i (i = 1,..., I; j = 1,..., n) overall mean βi treatment
More informationAssignment 6 Answer Keys
ssignment 6 nswer Keys Problem 1 (a) The treatment sum of squares can be calculated by SS Treatment = b a ȳi 2 Nȳ 2 i=1 = 5 (5.40 2 + 5.80 2 + 10 2 + 9.80 2 ) 20 7.75 2 = 92.95 Then the F statistic for
More informationWhat is Experimental Design?
One Factor ANOVA What is Experimental Design? A designed experiment is a test in which purposeful changes are made to the input variables (x) so that we may observe and identify the reasons for change
More informationLec 5: Factorial Experiment
November 21, 2011 Example Study of the battery life vs the factors temperatures and types of material. A: Types of material, 3 levels. B: Temperatures, 3 levels. Example Study of the battery life vs the
More informationChapter 10. Design of Experiments and Analysis of Variance
Chapter 10 Design of Experiments and Analysis of Variance Elements of a Designed Experiment Response variable Also called the dependent variable Factors (quantitative and qualitative) Also called the independent
More informationTwo-Way Analysis of Variance - no interaction
1 Two-Way Analysis of Variance - no interaction Example: Tests were conducted to assess the effects of two factors, engine type, and propellant type, on propellant burn rate in fired missiles. Three engine
More informationSAS Commands. General Plan. Output. Construct scatterplot / interaction plot. Run full model
Topic 23 - Unequal Replication Data Model Outline - Fall 2013 Parameter Estimates Inference Topic 23 2 Example Page 954 Data for Two Factor ANOVA Y is the response variable Factor A has levels i = 1, 2,...,
More information1-Way ANOVA MATH 143. Spring Department of Mathematics and Statistics Calvin College
1-Way ANOVA MATH 143 Department of Mathematics and Statistics Calvin College Spring 2010 The basic ANOVA situation Two variables: 1 Categorical, 1 Quantitative Main Question: Do the (means of) the quantitative
More informationCh 2: Simple Linear Regression
Ch 2: Simple Linear Regression 1. Simple Linear Regression Model A simple regression model with a single regressor x is y = β 0 + β 1 x + ɛ, where we assume that the error ɛ is independent random component
More informationChapter 20 : Two factor studies one case per treatment Chapter 21: Randomized complete block designs
Chapter 20 : Two factor studies one case per treatment Chapter 21: Randomized complete block designs Adapted from Timothy Hanson Department of Statistics, University of South Carolina Stat 705: Data Analysis
More informationAnalysis of Variance. ภาว น ศ ร ประภาน ก ล คณะเศรษฐศาสตร มหาว ทยาล ยธรรมศาสตร
Analysis of Variance ภาว น ศ ร ประภาน ก ล คณะเศรษฐศาสตร มหาว ทยาล ยธรรมศาสตร pawin@econ.tu.ac.th Outline Introduction One Factor Analysis of Variance Two Factor Analysis of Variance ANCOVA MANOVA Introduction
More informationIn social and biological sciences, a useful and frequently used statistical technique is the
Testing the additivity assumption in repeated measures ANOV A Robert I. Zambarano University of Texas Computation Center In social and biological sciences, a useful and frequently used statistical technique
More informationIntroduction to Design and Analysis of Experiments with the SAS System (Stat 7010 Lecture Notes)
Introduction to Design and Analysis of Experiments with the SAS System (Stat 7010 Lecture Notes) Asheber Abebe Discrete and Statistical Sciences Auburn University Contents 1 Completely Randomized Design
More informationMean Comparisons PLANNED F TESTS
Mean Comparisons F-tests provide information on significance of treatment effects, but no information on what the treatment effects are. Comparisons of treatment means provide information on what the treatment
More informationAnswer Keys to Homework#10
Answer Keys to Homework#10 Problem 1 Use either restricted or unrestricted mixed models. Problem 2 (a) First, the respective means for the 8 level combinations are listed in the following table A B C Mean
More informationAssignment 9 Answer Keys
Assignment 9 Answer Keys Problem 1 (a) First, the respective means for the 8 level combinations are listed in the following table A B C Mean 26.00 + 34.67 + 39.67 + + 49.33 + 42.33 + + 37.67 + + 54.67
More information1 Introduction to One-way ANOVA
Review Source: Chapter 10 - Analysis of Variance (ANOVA). Example Data Source: Example problem 10.1 (dataset: exp10-1.mtw) Link to Data: http://www.auburn.edu/~carpedm/courses/stat3610/textbookdata/minitab/
More informationRepeated Measurement ANOVA. Sungho Won, Ph. D. Graduate School of Public Health Seoul National University
1 Repeated Measurement ANOVA Sungho Won, Ph. D. Graduate School of Public Health Seoul National University 2 Analysis of Variance (ANOVA) Main idea Evaluate the effect of treatment by analyzing the amount
More informationWeek 7.1--IES 612-STA STA doc
Week 7.1--IES 612-STA 4-573-STA 4-576.doc IES 612/STA 4-576 Winter 2009 ANOVA MODELS model adequacy aka RESIDUAL ANALYSIS Numeric data samples from t populations obtained Assume Y ij ~ independent N(μ
More informationV. Experiments With Two Crossed Treatment Factors
V. Experiments With Two Crossed Treatment Factors A.The Experimental Design Completely Randomized Design (CRD) Let A be a factor with levels i = 1,,a B be a factor with levels j = 1,,b Y ijt = the response
More informationSleep data, two drugs Ch13.xls
Model Based Statistics in Biology. Part IV. The General Linear Mixed Model.. Chapter 13.3 Fixed*Random Effects (Paired t-test) ReCap. Part I (Chapters 1,2,3,4), Part II (Ch 5, 6, 7) ReCap Part III (Ch
More informationThis is a Randomized Block Design (RBD) with a single factor treatment arrangement (2 levels) which are fixed.
EXST3201 Chapter 13c Geaghan Fall 2005: Page 1 Linear Models Y ij = µ + βi + τ j + βτij + εijk This is a Randomized Block Design (RBD) with a single factor treatment arrangement (2 levels) which are fixed.
More informationOutline. Topic 20 - Diagnostics and Remedies. Residuals. Overview. Diagnostics Plots Residual checks Formal Tests. STAT Fall 2013
Topic 20 - Diagnostics and Remedies - Fall 2013 Diagnostics Plots Residual checks Formal Tests Remedial Measures Outline Topic 20 2 General assumptions Overview Normally distributed error terms Independent
More informationSummary of Chapter 7 (Sections ) and Chapter 8 (Section 8.1)
Summary of Chapter 7 (Sections 7.2-7.5) and Chapter 8 (Section 8.1) Chapter 7. Tests of Statistical Hypotheses 7.2. Tests about One Mean (1) Test about One Mean Case 1: σ is known. Assume that X N(µ, σ
More informationAnalysis of Variance
Statistical Techniques II EXST7015 Analysis of Variance 15a_ANOVA_Introduction 1 Design The simplest model for Analysis of Variance (ANOVA) is the CRD, the Completely Randomized Design This model is also
More informationTopic 28: Unequal Replication in Two-Way ANOVA
Topic 28: Unequal Replication in Two-Way ANOVA Outline Two-way ANOVA with unequal numbers of observations in the cells Data and model Regression approach Parameter estimates Previous analyses with constant
More informationSMA 6304 / MIT / MIT Manufacturing Systems. Lecture 10: Data and Regression Analysis. Lecturer: Prof. Duane S. Boning
SMA 6304 / MIT 2.853 / MIT 2.854 Manufacturing Systems Lecture 10: Data and Regression Analysis Lecturer: Prof. Duane S. Boning 1 Agenda 1. Comparison of Treatments (One Variable) Analysis of Variance
More informationANOVA: Analysis of Variation
ANOVA: Analysis of Variation The basic ANOVA situation Two variables: 1 Categorical, 1 Quantitative Main Question: Do the (means of) the quantitative variables depend on which group (given by categorical
More informationIncomplete Block Designs
Incomplete Block Designs Recall: in randomized complete block design, each of a treatments was used once within each of b blocks. In some situations, it will not be possible to use each of a treatments
More informationGeneral Linear Model (Chapter 4)
General Linear Model (Chapter 4) Outcome variable is considered continuous Simple linear regression Scatterplots OLS is BLUE under basic assumptions MSE estimates residual variance testing regression coefficients
More informationWhat If There Are More Than. Two Factor Levels?
What If There Are More Than Chapter 3 Two Factor Levels? Comparing more that two factor levels the analysis of variance ANOVA decomposition of total variability Statistical testing & analysis Checking
More informationAnalysis of variance
Analysis of variance 1 Method If the null hypothesis is true, then the populations are the same: they are normal, and they have the same mean and the same variance. We will estimate the numerical value
More information2. Two distinct error terms (Error for main plot effects > Error for subplot effects)
Lecture 16 Topic 12: The split-plot design and its relatives (Part I) Definition A split plot design results from a two-stage randomization process of a factorial treatment structure. Because of this two-stage
More informationOne-Way Analysis of Variance. With regression, we related two quantitative, typically continuous variables.
One-Way Analysis of Variance With regression, we related two quantitative, typically continuous variables. Often we wish to relate a quantitative response variable with a qualitative (or simply discrete)
More informationOne-way ANOVA Model Assumptions
One-way ANOVA Model Assumptions STAT:5201 Week 4: Lecture 1 1 / 31 One-way ANOVA: Model Assumptions Consider the single factor model: Y ij = µ + α }{{} i ij iid with ɛ ij N(0, σ 2 ) mean structure random
More informationAnalysing qpcr outcomes. Lecture Analysis of Variance by Dr Maartje Klapwijk
Analysing qpcr outcomes Lecture Analysis of Variance by Dr Maartje Klapwijk 22 October 2014 Personal Background Since 2009 Insect Ecologist at SLU Climate Change and other anthropogenic effects on interaction
More informationBIOL 933!! Lab 10!! Fall Topic 13: Covariance Analysis
BIOL 933!! Lab 10!! Fall 2017 Topic 13: Covariance Analysis Covariable as a tool for increasing precision Carrying out a full ANCOVA Testing ANOVA assumptions Happiness Covariables as Tools for Increasing
More informationAnalysis of Covariance
Analysis of Covariance (ANCOVA) Bruce A Craig Department of Statistics Purdue University STAT 514 Topic 10 1 When to Use ANCOVA In experiment, there is a nuisance factor x that is 1 Correlated with y 2
More information