Lecture 5: Comparing Treatment Means Montgomery: Section 3-5
|
|
- Tamsin James
- 5 years ago
- Views:
Transcription
1 Lecture 5: Comparing Treatment Means Montgomery: Section 3-5 Page 1
2 Linear Combination of Means ANOVA: y ij = µ + τ i + ɛ ij = µ i + ɛ ij Linear combination: L = c 1 µ 1 + c 1 µ c a µ a = a i=1 c iµ i Want to test: H 0 : L = c i µ i = L 0 pairwise comparison: µ i µ j =0 treatment vs. control: µ i µ 1 =0if treatment 1 is control comparing combinations of treatment effects: µ 1 2µ 2 + µ 3 =0. high order effects (orthogonal polynomials): linear, quadratic, etc. Page 2
3 Estimate of L: ˆL = c i ȳ i. Var(ˆL) = c 2 i Var(ȳ i.) =σ 2 c 2 i n i (= σ2 n c 2 i ) Standard Error for ˆL S.E. ˆL = MSE c 2 i n i Test statistic t 0 = (ˆL L 0 ) S.E. ˆL t(n a) under H 0 Page 3
4 Example: Lambs Diet Experiment Recall there are three diets and their treatment means are denoted by µ 1, µ 2 and µ 3. Suppose one wants to consider and test H 0 : L =60. L = µ 1 +2µ 2 +3µ 3 =6µ + τ 1 +2τ 2 +3τ 3 data lambs; input diet wtgain@@; cards; ; proc glm; class diet; model wtgain=diet; means diet; estimate l1 intercept 6 diet 1 2 3; run; Page 4
5 Example: Lambs Diet Experiment SAS output Level of wtgain diet N Mean Std Dev Dependent Variable: wtgain Standard Parameter Estimate Error t Value Pr > t l <.0001 t 0 = ( )/8.89 = 1.91 P value = P (t 1.91 or t 1.91 t(12 3)) =.088 Accept H 0 : µ 1 +2µ 2 +3µ 3 =60at α =5%. Page 5
6 Contrasts Γ= a i=1 c iµ i is a contrast if a i=1 c i =0. Equivalently, Γ= a i=1 c iτ i. Examples 1. Γ 1 = µ 1 µ 2 = µ 1 µ 2 +0µ 3 +0µ 4, c 1 =1,c 2 = 1,c 3 =0,c 4 =0 Comparing µ 1 and µ Γ 2 = µ 1 0.5µ 2 0.5µ 3 = µ 1 0.5µ 2 0.5µ 3 +0µ 4 c 1 =1,c 2 = 0.5,c 3 = 0.5,c 4 =0 Comparing µ 1 and the average of µ 2 and µ 3. Estimate for Γ: C = a i=1 c iȳ i. Page 6
7 Test: H 0 :Γ=0 t 0 = C S.E. C t(n a) t 2 0 = ( c i ȳ i. ) 2 MSE c 2 i Under H 0, t 2 0 F 1,N a. Contrast Sum of Squares n i = ( c i ȳ i. ) 2 / c 2 i /n i MSE = SS C/1 MSE SS C = ( ci y i. ) 2 / (c 2 i /n i ) SS C is the amount of variation caused by Γ. Page 7
8 SAS Code (cont.sas) Tensile Strength Example options ls=80; title1 Contrast Comparisons ; data one; infile c:\saswork\data\tensile.dat ; input percent strength time; proc glm data=one; class percent; model strength=percent; contrast C1 percent ; contrast C2 percent ; contrast C3 percent ; contrast C4 percent ; Page 8
9 Dependent Variable: STRENGTH Sum of Mean Source DF Squares Square F Value Pr > F Model Error Corrected Total Source DF Type I SS Mean Square F Value Pr > F PERCENT Contrast DF Contrast SS Mean Square F Value Pr > F C C C C Page 9
10 Orthogonal Contrasts Two contrasts {c i } and {d i } are Orthogonal if a i=1 c i d i n i =0 Example Γ 1 = µ 1 + µ 2 µ 3 µ 4,Soc 1 =1,c 2 =1,c 3 = 1,c 4 = 1. Γ 2 = µ 1 µ 2 + µ 3 µ 4.Sod 1 =1,d 2 = 1,d 3 =1,d 4 = 1 It is easy to verify that both Γ 1 and Γ 2 are contrasts. Furthermore, c 1 d 1 + c 2 d 2 + c 3 d 3 + c 4 d 4 = ( 1) + ( 1) 1+( 1) ( 1)=0. Hence, Γ 1 and Γ 2 are orthogonal to each other. A complete set of orthogonal contrasts C = {Γ 1, Γ 2,...,Γ a 1 } if contrasts are mutually orthogonal and there does not exist a contrast orthogonal to all contrasts in C. Page 10
11 If there are a treatments, C must contain a 1 contrasts. Complete set is not unique. For example, in the tensile strength example C 1 : includes : Γ 1 = (0, 0, 0, 1, 1) Γ 2 = (1, 0, 1, 1, 1) Γ 3 = (1, 0, 1, 0, 0) Γ 4 = (1, 4, 1, 1, 1) C 2 : includes : Γ 1 = ( 2, 1, 0, 1, 2) Γ 2 = (2, 1, 2, 1, 2) Γ 3 = ( 1, 2, 0, 2, 1) Γ 4 = (1, 4, 6, 4, 1) Page 11
12 Orthogonal Contrasts Orthoogonal contrasts (estimates) are independent with each other. Suppose C 1,C 2,...,C a 1 are the estimates of the contrasts in a complete set of comtrasts {Γ 1, Γ 2,...,Γ a 1 }, then SS Treatment = SS C1 + SS C2 + + SS Ca 1 Recall in ANOVA, F 0 = MS Treatment MSE, F 0 = SS C 1 /MSE + + SS Ca 1 /MSE a 1 = F 10 + F F (a 1)0 a 1 where F i0 is the test statistic used to test contrast Γ i. Example on Slide 9 Page 12
13 Tensile Example Try to model mean response as a function of treatments Page 13
14 Orthogonal contrasts and orthogonal polynomial model Treatments are quantitative (assume a =4) One can use general polynomial model to fit the trend (t: level or treatment). f(t) =a 0 + a 1 t + a 2 t 2 + a 3 t 3 Regression can be used to get the estimates for a 1, a 2 and a 3. We will use orthogonal polynomial model f(t) =β 0 + β 1 P 1 (t)+β 2 P 2 (t)+β 3 P 3 (t) where P 1 (t), P 2 (t) and P 3 (t) are pre-specified polynomials of order 1, 2 and 3, respectively. P 1 (t) is linear, P 2 (t) is quadratic and P 3 (t) is cubic. Let t 1, t 2,..., t a are the treatments (equally spaced), then the polynomials corresponds to the following contrasts: Page 14
15 t t 1 t 2 t a Contrasts D P 1 (t) P 1 (t 1 ) P 1 (t 2 ) P 1 (t a ) Γ 1 D 1 P 2 (t) P 2 (t 1 ) P 2 (t 2 ) P 2 (t a ) Γ 2 D 2 P 3 (t) P 3 (t 1 ) P 3 (t 2 ) P 3 (t a ) Γ 3 D 3 where D i = P i (t 1 ) 2 + P i (t 2 ) P i (t a ) 2 If Γ 1, Γ 2 and Γ 3 are orthogonal to each other, then we say P 1 (t), P 2 (t) and P 3 (t) are orthogonal polynomials. Coefficients β i can be estimated and tested by the contrasts Γ 1, Γ 2 and Γ 3. Predict f(t) when t is not a treatment used in the experiment. Page 15
16 tensile strength example: orthogonal polynomial effects Treatment levels t k : 15, 20, 25, 30, 35; Median: 25; Pace: 5 Orthogonal polynomials: let x =(t 25)/5. P 1 (t) =x P 2 (t) =x 2 2 P 3 (t) =5/6[x 3 17x/5] P 4 (t) =35/12[x 4 31x/7+72/35] Polynomial Contrasts and Effects t Contrast D Effect (Trend) P 1 (t) Γ 1 D 1 =10 linear P 2 (t) Γ 2 D 2 =14 quadratic P 3 (t) Γ 3 D 3 =10 cubic p 4 (t) Γ 4 D 4 =70 4th order The contrasts can be directly derived from Table IX or Table X. Page 16
17 Want to fit the model f(t) =β 0 + β 1 P 1 (t)+β 2 P 2 (t)+β 3 P 3 (t)+β 4 P 4 (t) Estimation and Testing β 1 :useγ 1, ˆβ 1 = c 11ȳ c 15 ȳ 5. D 1 Test: H 0 : β 1 =0, F 10 = SS C i MSE F 1,N 5. β 2 :useγ 2, ˆβ 2 = c 21ȳ c 25 ȳ 5. D 2 Test: H 0 : β 2 =0, F 20 = SS C i MSE F 1,N 5 Similar for β 3 and β 4 Question: what is the estimate for β 0? Page 17
18 General formulas for orthogonal polynomial of degrees 1-4 One factor of a levels l 1, l 2,..., l a, equally spaced. Let m be the median, δ be the difference between two consecutive levels: P 1 (t) =λ 1 ( t m ) δ P 2 (t) =λ 2 [( t m ) 2 a2 1 δ 12 ] P 4 (t) =λ 4 [( t m δ P 3 (t) =λ 3 [( t m δ ) 4 ( t m δ ) 3 ( t m δ )( 3a2 7 )] 20 ) 2 ( 3a2 13 )+ 3(a2 1)(a 2 9) (λ i ) are constants to make the polynomials have integer values at the treatment levels, they are available from Table IX or Table X. Tensile Strength Example: m=25, δ =5,(λ i )=(1, 1, 5/6, 35/12) Page 18
19 SAS tensile strength example data one; infile c:\saswork\data\tensile.dat ; input percent strength time; proc glm data=one; class percent; model strength=percent; estimate C1 percent ; estimate C2 percent ; estimate C3 percent ; estimate C4 percent ; contrast C1 percent ; contrast C2 percent ; contrast C3 percent ; contrast C4 percent ; run; Page 19
20 Output Dependent Variable: strength Sum of Source DF Squares Mean Square F Value Pr > F Model <.0001 Error Corrected Total Parameter Estimate Error t Value Pr > t C C <.0001 C C Contrast DF Contrast SS Mean Square F Value Pr > F C C <.0001 C C Page 20
21 Estimates Hence, ˆβ 1 =8.20/10 =.82; ˆβ2 = 31/14 = ˆβ 3 = 11.4/10 = 1.14; ˆβ4 = 21.8/70 = So the fitted functional relationship between tensile strength y and cotton percent (t) is y = ˆβ P 1 (t) 2.214P 2 (t) 1.14P 3 (t) 0.311P 4 (t), where P 1 (t),...,p 4 (t) are defined on Slide 16. Page 21
22 Testing Multiple Contrasts (Multiple Comparisons) Using Confidence Intervals One contrast: H 0 :Γ= c i µ i =Γ 0 vs H 1 :Γ Γ 0 at α 100(1-α) Confidence Interval (CI) for Γ: CI : c i ȳ i. ± t α/2,n a MS E c 2 i n i P (CI not contain L 0 H 0 )=α(= type I error) Decision Rule: Reject H 0 if CI does not contain Γ 0. Page 22
23 Multiple contrasts H 0 :Γ 1 =Γ 1 0,...Γ m =Γ m 0 vs H 1 : at least one does not hold If we construct CI 1,CI 2,..., CI m, each with 100(1-α) level, then for each CI i, P (CI i not containγ i 0 H 0 )=α, for i =1,...,m But the overall error rate (probability of type I error for H 0 vs H 1 ) is inflated and much larger than α, that is, P (at least one CI i not contain Γ i 0 H 0 ) >> α One way to achieve small overall error rate, we require much smaller error rate (α ) of each individual CI i. Page 23
24 Bonferroni Inequality Bonferroni Method for Testing Multiple Contrasts P ( at least one CI i not contain Γ i 0 H 0 ) = P (CI 1 not contain..or...or CI m not contain H 0 ) P ( CI 1 not H 0 )+ + P ( CI m not H 0 )=mα In order to control overall error rate (or, overall confidence level), let mα = α, we have,α = α/m Bonferroni CIs: CI i : c ij ȳ j. ± t α/2m (N a) MS E c 2 ij n j When m is large, Bonferroni CIs are too conservative ( overall type II error too large). Page 24
25 Scheffe s method for Testing All Contrasts Consider all possible contrasts: Γ= c i µ i Estimate: C = c i ȳ i., St. Error: S.E. C = MS E c 2 i n i Critical value: (a 1)Fα,a 1,N a Scheffe s simultaneous CI: C ± (a 1)F α,a 1,N a S.E. C Overall confidence level and error rate for m contrasts P (CIs contain true parameter for any contrast) 1 α P (at least one CI does not contain true parameter) α Remark: Scheffe s method is also conservative, too conservative when m is small Page 25
26 Methods for Pairwise Comparisons There are a(a 1)/2 possible pairs: µ i µ j (contrast for comparing µ i and µ j ). We may be interested in m pairs or all pairs. Standard Procedure: 1. Estimation: ȳ i. ȳ j. 2. Compute a Critical Difference (CD) (based on the method employed) 3. If ȳ i. ȳ j. > CD or equivalently if the interval (ȳ i. ȳ j. CD, ȳ i. ȳ j. +CD) does not contain zero, declare µ i µ j significant. Page 26
27 Methods for Calculating CD. Least significant difference (LSD): not control overall error rate Bonferroni method (for m pairs) CD = t α/2,n a MS E (1/n i +1/n j ) CD = t α/2m,n a MS E (1/n i +1/n j ) control overall error rate for the m comparisons. Tukey s method (for all possible pairs) CD = q α(a, N a) 2 MS E (1/n i +1/n j ) q α (a, N a) from studentized range distribution (Table VII or Table VIII). Control overall error rate (exact for balanced experiments). (Example 3.7). Page 27
28 Comparing treatments with control (Dunnett s method) 1. Assume µ 1 is a control, and µ 2,...,µ a are (new) treatments 2. Only interested in a 1 pairs: µ 2 µ 1,..., µ a µ 1 3. Compare ȳ i. ȳ 1. to CD = d α (a 1,N a) MS E (1/n i +1/n 1 ) where d α (p, f) from Table IX or Table VIII: critical values for Dunnett s test. 4. Remark: control overall error rate. Read Example 3-9 (or 3-10) For pairwise comparison, which method should be preferred? LSD, Bonferroni, Tukey, Dunnett or others? Page 28
29 SAS Code data one; infile c:\saswork\data\tensile.dat ; input percent strength time; proc glm data=one; class percent; model strength=precent; /* Construct CI for Treatment Means*/ means percent /alpha=.05 lsd clm; means percent / alpha=.05 bon clm; /* Pairwise Comparison*/ means percent /alpha=.05 lines lsd; means percent /alpha=.05 lines bon; means percent /alpha=.05 lines scheffe; means percent /alpha=.05 lines tukey; means percent /dunnett; run; Page 29
30 The GLM Procedure t Confidence Intervals for y Alpha 0.05 Error Degrees of Freedom 20 Error Mean Square 8.06 Critical Value of t Half Width of Confidence Interval % Confidence trt N Mean Limits Page 30
31 The GLM Procedure Bonferroni t Confidence Intervals for y Alpha 0.05 Error Degrees of Freedom 20 Error Mean Square 8.06 Critical Value of t Half Width of Confidence Interval Simultaneous 95% trt N Mean Confidence Limits Page 31
32 NOTE: This test controls the Type I comparisonwise error rate, not t t Tests (LSD) for y experimentwise error rate. Alpha 0.05 Error Degrees of Freedom 20 Error Mean Square 8.06 Critical Value of t Least Significant Difference Means with the same letter are not significantly different. t Grouping Mean N trt A B B C C Page 32
33 NOTE: This test controls the Type I experimentwise error rate, but it ge Bonferroni (Dunn) t Tests for y has a higher Type II error rate than REGWQ. Alpha 0.05 Error Degrees of Freedom 20 Error Mean Square 8.06 Critical Value of t Minimum Significant Difference Means with the same letter are not significantly different. Bon Grouping Mean N trt A B A B C C C Page 33
34 Scheffe s Test for y NOTE: This test controls the Type I experimentwise error rate. Alpha 0.05 Error Degrees of Freedom 20 Error Mean Square 8.06 Critical Value of F Minimum Significant Difference Means with the same letter are not significantly different. Scheffe Grouping Mean N trt A A B A B B C C C C C Page 34
35 NOTE: This test controls the Type I experimentwise error rate, but it ge Tukey s Studentized Range (HSD) Test for y has a higher Type II error rate than REGWQ. Alpha 0.05 Error Degrees of Freedom 20 Error Mean Square 8.06 Critical Value of Studentized Range Minimum Significant Difference Means with the same letter are not significantly different. Tukey Grouping Mean N trt A A B A B B C C D C D D Page 35
36 NOTE: This test controls the Type I experimentwise error for comparisons Dunnett s t Tests for y treatments against a control. Alpha 0.05 Error Degrees of Freedom 20 Error Mean Square 8.06 Critical Value of Dunnett s t Minimum Significant Difference Comparisons significant at the 0.05 level are indicated by ***. Difference trt Between Simultaneous 95% Comparison Means Confidence Limits *** *** *** Page 36
37 Determining Sample Size More replicates required to detect small treatment effects Operating Characteristic Curves for F tests Probability of type II error β = P ( accept H 0 H 0 is false) = P (F 0 <F α,a 1,N a H 1 is correct ) Under H 1, F 0 follows a noncentral F distribution with noncentrality λ and degrees of freedom, a 1 and N a. Let Φ 2 = n a i=1 τ 2 i aσ 2 OC curves of β vs n and Φ are included in Chart V for various α and a. Read Example 3.10 or Page 37
Linear Combinations. Comparison of treatment means. Bruce A Craig. Department of Statistics Purdue University. STAT 514 Topic 6 1
Linear Combinations Comparison of treatment means Bruce A Craig Department of Statistics Purdue University STAT 514 Topic 6 1 Linear Combinations of Means y ij = µ + τ i + ǫ ij = µ i + ǫ ij Often study
More informationLec 1: An Introduction to ANOVA
Ying Li Stockholm University October 31, 2011 Three end-aisle displays Which is the best? Design of the Experiment Identify the stores of the similar size and type. The displays are randomly assigned to
More informationOutline. Topic 19 - Inference. The Cell Means Model. Estimates. Inference for Means Differences in cell means Contrasts. STAT Fall 2013
Topic 19 - Inference - Fall 2013 Outline Inference for Means Differences in cell means Contrasts Multiplicity Topic 19 2 The Cell Means Model Expressed numerically Y ij = µ i + ε ij where µ i is the theoretical
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 informationLaboratory Topics 4 & 5
PLS205 Lab 3 January 23, 2014 Orthogonal contrasts Class comparisons in SAS Trend analysis in SAS Multiple mean comparisons Laboratory Topics 4 & 5 Orthogonal contrasts Planned, single degree-of-freedom
More informationThese are all actually contrasts (the coef sum to zero). What are these contrasts representing? What would make them large?
Lecture 12 Comparing treatment effects Orthogonal Contrasts What use are contrasts? Recall the Cotton data In this case, the treatment levels have an ordering to them this is not always the case) Consider
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 informationLinear Combinations of Group Means
Linear Combinations of Group Means Look at the handicap example on p. 150 of the text. proc means data=mth567.disability; class handicap; var score; proc sort data=mth567.disability; by handicap; proc
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 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 informationAnalysis of Variance II Bios 662
Analysis of Variance II Bios 662 Michael G. Hudgens, Ph.D. mhudgens@bios.unc.edu http://www.bios.unc.edu/ mhudgens 2008-10-24 17:21 BIOS 662 1 ANOVA II Outline Multiple Comparisons Scheffe Tukey Bonferroni
More information4.8 Alternate Analysis as a Oneway ANOVA
4.8 Alternate Analysis as a Oneway ANOVA Suppose we have data from a two-factor factorial design. The following method can be used to perform a multiple comparison test to compare treatment means as well
More informationT-test: means of Spock's judge versus all other judges 1 12:10 Wednesday, January 5, judge1 N Mean Std Dev Std Err Minimum Maximum
T-test: means of Spock's judge versus all other judges 1 The TTEST Procedure Variable: pcwomen judge1 N Mean Std Dev Std Err Minimum Maximum OTHER 37 29.4919 7.4308 1.2216 16.5000 48.9000 SPOCKS 9 14.6222
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 informationSTAT22200 Spring 2014 Chapter 5
STAT22200 Spring 2014 Chapter 5 Yibi Huang April 29, 2014 Chapter 5 Multiple Comparisons Chapter 5-1 Chapter 5 Multiple Comparisons Note the t-tests and C.I. s are constructed assuming we only do one test,
More informationOne-Way Analysis of Variance (ANOVA) There are two key differences regarding the explanatory variable X.
One-Way Analysis of Variance (ANOVA) Also called single factor ANOVA. The response variable Y is continuous (same as in regression). There are two key differences regarding the explanatory variable X.
More informationMore about Single Factor Experiments
More about Single Factor Experiments 1 2 3 0 / 23 1 2 3 1 / 23 Parameter estimation Effect Model (1): Y ij = µ + A i + ɛ ij, Ji A i = 0 Estimation: µ + A i = y i. ˆµ = y..  i = y i. y.. Effect Modell
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 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 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 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 informationBiological Applications of ANOVA - Examples and Readings
BIO 575 Biological Applications of ANOVA - Winter Quarter 2010 Page 1 ANOVA Pac Biological Applications of ANOVA - Examples and Readings One-factor Model I (Fixed Effects) This is the same example for
More information1 Tomato yield example.
ST706 - Linear Models II. Spring 2013 Two-way Analysis of Variance examples. Here we illustrate what happens analyzing two way data in proc glm in SAS. Similar issues come up with other software where
More informationLecture 3. Experiments with a Single Factor: ANOVA Montgomery Sections 3.1 through 3.5 and Section
Lecture 3. Experiments with a Single Factor: ANOVA Montgomery Sections 3.1 through 3.5 and Section 15.1.1 1 Lecture 3 Page 1 Tensile Strength Experiment Investigate the tensile strength of a new synthetic
More informationN J SS W /df W N - 1
One-Way ANOVA Source Table ANOVA MODEL: ij = µ* + α j + ε ij H 0 : µ = µ =... = µ j or H 0 : Σα j = 0 Source Sum of Squares df Mean Squares F J Between Groups nj( j * ) J - SS B /(J ) MS B /MS W = ( N
More informationThe Distribution of F
The Distribution of F It can be shown that F = SS Treat/(t 1) SS E /(N t) F t 1,N t,λ a noncentral F-distribution with t 1 and N t degrees of freedom and noncentrality parameter λ = t i=1 n i(µ i µ) 2
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 informationST505/S697R: Fall Homework 2 Solution.
ST505/S69R: Fall 2012. Homework 2 Solution. 1. 1a; problem 1.22 Below is the summary information (edited) from the regression (using R output); code at end of solution as is code and output for SAS. a)
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 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 informationSTAT 5200 Handout #7a Contrasts & Post hoc Means Comparisons (Ch. 4-5)
STAT 5200 Handout #7a Contrasts & Post hoc Means Comparisons Ch. 4-5) Recall CRD means and effects models: Y ij = µ i + ϵ ij = µ + α i + ϵ ij i = 1,..., g ; j = 1,..., n ; ϵ ij s iid N0, σ 2 ) If we reject
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 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 informationCh 3: Multiple Linear Regression
Ch 3: Multiple Linear Regression 1. Multiple Linear Regression Model Multiple regression model has more than one regressor. For example, we have one response variable and two regressor variables: 1. delivery
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 informationSTAT 705 Chapter 19: Two-way ANOVA
STAT 705 Chapter 19: Two-way ANOVA Adapted from Timothy Hanson Department of Statistics, University of South Carolina Stat 705: Data Analysis II 1 / 41 Two-way ANOVA This material is covered in Sections
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 informationEXST Regression Techniques Page 1. We can also test the hypothesis H :" œ 0 versus H :"
EXST704 - Regression Techniques Page 1 Using F tests instead of t-tests We can also test the hypothesis H :" œ 0 versus H :" Á 0 with an F test.! " " " F œ MSRegression MSError This test is mathematically
More informationIntroduction. Chapter 8
Chapter 8 Introduction In general, a researcher wants to compare one treatment against another. The analysis of variance (ANOVA) is a general test for comparing treatment means. When the null hypothesis
More informationLecture 9: Factorial Design Montgomery: chapter 5
Lecture 9: Factorial Design Montgomery: chapter 5 Page 1 Examples Example I. Two factors (A, B) each with two levels (, +) Page 2 Three Data for Example I Ex.I-Data 1 A B + + 27,33 51,51 18,22 39,41 EX.I-Data
More informationMultiple Testing. Tim Hanson. January, Modified from originals by Gary W. Oehlert. Department of Statistics University of South Carolina
Multiple Testing Tim Hanson Department of Statistics University of South Carolina January, 2017 Modified from originals by Gary W. Oehlert Type I error A Type I error is to wrongly reject the null hypothesis
More informationLecture 2: Basic Concepts and Simple Comparative Experiments Montgomery: Chapter 2
Lecture 2: Basic Concepts and Simple Comparative Experiments Montgomery: Chapter 2 Fall, 2013 Page 1 Random Variable and Probability Distribution Discrete random variable Y : Finite possible values {y
More informationA posteriori multiple comparison tests
A posteriori multiple comparison tests 11/15/16 1 Recall the Lakes experiment Source of variation SS DF MS F P Lakes 58.000 2 29.400 8.243 0.006 Error 42.800 12 3.567 Total 101.600 14 The ANOVA tells us
More informationMultiple Testing. Gary W. Oehlert. January 28, School of Statistics University of Minnesota
Multiple Testing Gary W. Oehlert School of Statistics University of Minnesota January 28, 2016 Background Suppose that you had a 20-sided die. Nineteen of the sides are labeled 0 and one of the sides is
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 informationWeek 14 Comparing k(> 2) Populations
Week 14 Comparing k(> 2) Populations Week 14 Objectives Methods associated with testing for the equality of k(> 2) means or proportions are presented. Post-testing concepts and analysis are introduced.
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 informationSTATISTICS 479 Exam II (100 points)
Name STATISTICS 79 Exam II (1 points) 1. A SAS data set was created using the following input statement: Answer parts(a) to (e) below. input State $ City $ Pop199 Income Housing Electric; (a) () Give the
More informationMultiple comparisons The problem with the one-pair-at-a-time approach is its error rate.
Multiple comparisons The problem with the one-pair-at-a-time approach is its error rate. Each confidence interval has a 95% probability of making a correct statement, and hence a 5% probability of making
More informationLinear models and their mathematical foundations: Simple linear regression
Linear models and their mathematical foundations: Simple linear regression Steffen Unkel Department of Medical Statistics University Medical Center Göttingen, Germany Winter term 2018/19 1/21 Introduction
More informationLecture 6 Multiple Linear Regression, cont.
Lecture 6 Multiple Linear Regression, cont. BIOST 515 January 22, 2004 BIOST 515, Lecture 6 Testing general linear hypotheses Suppose we are interested in testing linear combinations of the regression
More informationSTAT 705 Chapter 19: Two-way ANOVA
STAT 705 Chapter 19: Two-way ANOVA Timothy Hanson Department of Statistics, University of South Carolina Stat 705: Data Analysis II 1 / 38 Two-way ANOVA Material covered in Sections 19.2 19.4, but a bit
More informationMultiple Comparison Procedures Cohen Chapter 13. For EDUC/PSY 6600
Multiple Comparison Procedures Cohen Chapter 13 For EDUC/PSY 6600 1 We have to go to the deductions and the inferences, said Lestrade, winking at me. I find it hard enough to tackle facts, Holmes, without
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 informationOrthogonal contrasts and multiple comparisons
BIOL 933 Lab 4 Fall 2017 Orthogonal contrasts Class comparisons in R Trend analysis in R Multiple mean comparisons Orthogonal contrasts and multiple comparisons Orthogonal contrasts Planned, single degree-of-freedom
More informationThe One-Way Independent-Samples ANOVA. (For Between-Subjects Designs)
The One-Way Independent-Samples ANOVA (For Between-Subjects Designs) Computations for the ANOVA In computing the terms required for the F-statistic, we won t explicitly compute any sample variances or
More informationST430 Exam 1 with Answers
ST430 Exam 1 with Answers Date: October 5, 2015 Name: Guideline: You may use one-page (front and back of a standard A4 paper) of notes. No laptop or textook are permitted but you may use a calculator.
More informationSpecific Differences. Lukas Meier, Seminar für Statistik
Specific Differences Lukas Meier, Seminar für Statistik Problem with Global F-test Problem: Global F-test (aka omnibus F-test) is very unspecific. Typically: Want a more precise answer (or have a more
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 informationLecture 5: Determining Sample Size
Lecture 5: Determining Sample Size Montgomery: Section 3.7 and 13.4 1 Lecture 5 Page 1 Choice of Sample Size: Fixed Effects Can determine the sample size for OverallF test Contrasts of interest For simplicity,
More informationSTAT 350: Summer Semester Midterm 1: Solutions
Name: Student Number: STAT 350: Summer Semester 2008 Midterm 1: Solutions 9 June 2008 Instructor: Richard Lockhart Instructions: This is an open book test. You may use notes, text, other books and a calculator.
More informationChapter 6 Multiple Regression
STAT 525 FALL 2018 Chapter 6 Multiple Regression Professor Min Zhang The Data and Model Still have single response variable Y Now have multiple explanatory variables Examples: Blood Pressure vs Age, Weight,
More informationIntroduction to Analysis of Variance (ANOVA) Part 2
Introduction to Analysis of Variance (ANOVA) Part 2 Single factor Serpulid recruitment and biofilms Effect of biofilm type on number of recruiting serpulid worms in Port Phillip Bay Response variable:
More informationH0: Tested by k-grp ANOVA
Analyses of K-Group Designs : Omnibus F, Pairwise Comparisons & Trend Analyses ANOVA for multiple condition designs Pairwise comparisons and RH Testing Alpha inflation & Correction LSD & HSD procedures
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 informationPairwise multiple comparisons are easy to compute using SAS Proc GLM. The basic statement is:
Pairwise Multiple Comparisons in SAS Pairwise multiple comparisons are easy to compute using SAS Proc GLM. The basic statement is: means effects / options Here, means is the statement initiat, effects
More informationComparisons among means (or, the analysis of factor effects)
Comparisons among means (or, the analysis of factor effects) In carrying out our usual test that μ 1 = = μ r, we might be content to just reject this omnibus hypothesis but typically more is required:
More informationChapter Seven: Multi-Sample Methods 1/52
Chapter Seven: Multi-Sample Methods 1/52 7.1 Introduction 2/52 Introduction The independent samples t test and the independent samples Z test for a difference between proportions are designed to analyze
More informationLAMPIRAN. Ket. tn = tidak berbeda nyata
LAMPIRAN Lampiran 1.Tabel. Analisis Sidik Ragam Konsumsi ahan Kering analisis sidik ragam Sk Db JK KT Fhit Ftabel 0.05 0.01 Periode 3 4483.91 1494.63 61.23 tn 3.86 6.99 Sapi 3 15.66 5.22 0.21 tn 3.86 6.99
More informationDESIGNING EXPERIMENTS AND ANALYZING DATA A Model Comparison Perspective
DESIGNING EXPERIMENTS AND ANALYZING DATA A Model Comparison Perspective Second Edition Scott E. Maxwell Uniuersity of Notre Dame Harold D. Delaney Uniuersity of New Mexico J,t{,.?; LAWRENCE ERLBAUM ASSOCIATES,
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 informationSTA2601. Tutorial letter 203/2/2017. Applied Statistics II. Semester 2. Department of Statistics STA2601/203/2/2017. Solutions to Assignment 03
STA60/03//07 Tutorial letter 03//07 Applied Statistics II STA60 Semester Department of Statistics Solutions to Assignment 03 Define tomorrow. university of south africa QUESTION (a) (i) The normal quantile
More informationProblems. Suppose both models are fitted to the same data. Show that SS Res, A SS Res, B
Simple Linear Regression 35 Problems 1 Consider a set of data (x i, y i ), i =1, 2,,n, and the following two regression models: y i = β 0 + β 1 x i + ε, (i =1, 2,,n), Model A y i = γ 0 + γ 1 x i + γ 2
More informationUNIVERSITY EXAMINATIONS NJORO CAMPUS SECOND SEMESTER 2011/2012
UNIVERSITY EXAMINATIONS NJORO CAMPUS SECOND SEMESTER 2011/2012 THIRD YEAR EXAMINATION FOR THE AWARD BACHELOR OF SCIENCE IN AGRICULTURE AND BACHELOR OF SCIENCE IN FOOD TECHNOLOGY AGRO 391 AGRICULTURAL EXPERIMENTATION
More informationInference for Regression
Inference for Regression Section 9.4 Cathy Poliak, Ph.D. cathy@math.uh.edu Office in Fleming 11c Department of Mathematics University of Houston Lecture 13b - 3339 Cathy Poliak, Ph.D. cathy@math.uh.edu
More informationDETAILED CONTENTS PART I INTRODUCTION AND DESCRIPTIVE STATISTICS. 1. Introduction to Statistics
DETAILED CONTENTS About the Author Preface to the Instructor To the Student How to Use SPSS With This Book PART I INTRODUCTION AND DESCRIPTIVE STATISTICS 1. Introduction to Statistics 1.1 Descriptive and
More informationLecture 7: Latin Squares and Related Designs
Lecture 7: Latin Squares and Related Designs Montgomery: Section 4.2 and 4.3 1 Lecture 7 Page 1 Automobile Emission Experiment Four cars and four drivers are employed in a study of four gasoline additives(a,b,c,
More informationH0: Tested by k-grp ANOVA
Pairwise Comparisons ANOVA for multiple condition designs Pairwise comparisons and RH Testing Alpha inflation & Correction LSD & HSD procedures Alpha estimation reconsidered H0: Tested by k-grp ANOVA Regardless
More informationSTAT 525 Fall Final exam. Tuesday December 14, 2010
STAT 525 Fall 2010 Final exam Tuesday December 14, 2010 Time: 2 hours Name (please print): Show all your work and calculations. Partial credit will be given for work that is partially correct. Points will
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 information3 Variables: Cyberloafing Conscientiousness Age
title 'Cyberloafing, Mike Sage'; run; PROC CORR data=sage; var Cyberloafing Conscientiousness Age; run; quit; The CORR Procedure 3 Variables: Cyberloafing Conscientiousness Age Simple Statistics Variable
More information9 One-Way Analysis of Variance
9 One-Way Analysis of Variance SW Chapter 11 - all sections except 6. The one-way analysis of variance (ANOVA) is a generalization of the two sample t test to k 2 groups. Assume that the populations of
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 informationdata proc sort proc corr run proc reg run proc glm run proc glm run proc glm run proc reg CONMAIN CONINT run proc reg DUMMAIN DUMINT run proc reg
data one; input id Y group X; I1=0;I2=0;I3=0;if group=1 then I1=1;if group=2 then I2=1;if group=3 then I3=1; IINT1=I1*X;IINT2=I2*X;IINT3=I3*X; *************************************************************************;
More information4:3 LEC - PLANNED COMPARISONS AND REGRESSION ANALYSES
4:3 LEC - PLANNED COMPARISONS AND REGRESSION ANALYSES FOR SINGLE FACTOR BETWEEN-S DESIGNS Planned or A Priori Comparisons We previously showed various ways to test all possible pairwise comparisons for
More informationGeneral Linear Models (GLM) for Fixed Factors
Chapter 224 General Linear Models (GLM) for Fixed Factors Introduction This procedure performs analysis of variance (ANOVA) and analysis of covariance (ANCOVA) for factorial models that include fixed factors
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 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 informationChapter 2 Inferences in Simple Linear Regression
STAT 525 SPRING 2018 Chapter 2 Inferences in Simple Linear Regression Professor Min Zhang Testing for Linear Relationship Term β 1 X i defines linear relationship Will then test H 0 : β 1 = 0 Test requires
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 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 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 informationCOMPARING SEVERAL MEANS: ANOVA
LAST UPDATED: November 15, 2012 COMPARING SEVERAL MEANS: ANOVA Objectives 2 Basic principles of ANOVA Equations underlying one-way ANOVA Doing a one-way ANOVA in R Following up an ANOVA: Planned contrasts/comparisons
More informationOverview Scatter Plot Example
Overview Topic 22 - Linear Regression and Correlation STAT 5 Professor Bruce Craig Consider one population but two variables For each sampling unit observe X and Y Assume linear relationship between variables
More informationDescriptions of post-hoc tests
Experimental Statistics II Page 81 Descriptions of post-hoc tests Post-hoc or Post-ANOVA tests! Once you have found out some treatment(s) are different, how do you determine which one(s) are different?
More informationMAT3378 (Winter 2016)
MAT3378 (Winter 2016) Assignment 2 - SOLUTIONS Total number of points for Assignment 2: 12 The following questions will be marked: Q1, Q2, Q4 Q1. (4 points) Assume that Z 1,..., Z n are i.i.d. normal random
More informationMultiple Linear Regression
Multiple Linear Regression ST 430/514 Recall: a regression model describes how a dependent variable (or response) Y is affected, on average, by one or more independent variables (or factors, or covariates).
More informationTopic 25 - One-Way Random Effects Models. Outline. Random Effects vs Fixed Effects. Data for One-way Random Effects Model. One-way Random effects
Topic 5 - One-Way Random Effects Models One-way Random effects Outline Model Variance component estimation - Fall 013 Confidence intervals Topic 5 Random Effects vs Fixed Effects Consider factor with numerous
More informationLecture 10: 2 k Factorial Design Montgomery: Chapter 6
Lecture 10: 2 k Factorial Design Montgomery: Chapter 6 Page 1 2 k Factorial Design Involving k factors Each factor has two levels (often labeled + and ) Factor screening experiment (preliminary study)
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 information