Analysis of Variance and Design of Experiments-II

Similar documents
Chapter 7 Partially Balanced Incomplete Block Design (PBIBD)

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur

where I = (n x n) diagonal identity matrix with diagonal elements = 1 and off-diagonal elements = 0; and σ 2 e = variance of (Y X).

Bose (1942) showed b t r 1 is a necessary condition. PROOF (Murty 1961): Assume t is a multiple of k, i.e. t nk, where n is an integer.

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

LINEAR REGRESSION ANALYSIS. MODULE VIII Lecture Indicator Variables

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)

[The following data appear in Wooldridge Q2.3.] The table below contains the ACT score and college GPA for eight college students.

Chapter 12 Analysis of Covariance

Chapter 13: Multiple Regression

Polynomial Regression Models

PARTIALLY BALANCED INCOMPLETE BLOCK DESIGNS

F statistic = s2 1 s 2 ( F for Fisher )

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur

ECONOMICS 351*-A Mid-Term Exam -- Fall Term 2000 Page 1 of 13 pages. QUEEN'S UNIVERSITY AT KINGSTON Department of Economics

Chapter 8 Indicator Variables

Chapter 6. Supplemental Text Material

Two-factor model. Statistical Models. Least Squares estimation in LM two-factor model. Rats

Lecture 2: Prelude to the big shrink

Chapter 4 Experimental Design and Their Analysis

x i1 =1 for all i (the constant ).

17 Nested and Higher Order Designs

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur

Chapter 11: Simple Linear Regression and Correlation

# c i. INFERENCE FOR CONTRASTS (Chapter 4) It's unbiased: Recall: A contrast is a linear combination of effects with coefficients summing to zero:

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Lecture 9: Linear regression: centering, hypothesis testing, multiple covariates, and confounding

Lecture 9: Linear regression: centering, hypothesis testing, multiple covariates, and confounding

International Journal of Mathematical Archive-9(3), 2018, Available online through ISSN

Economics 130. Lecture 4 Simple Linear Regression Continued

Tests of Single Linear Coefficient Restrictions: t-tests and F-tests. 1. Basic Rules. 2. Testing Single Linear Coefficient Restrictions

Regression. The Simple Linear Regression Model

Interval Estimation in the Classical Normal Linear Regression Model. 1. Introduction

BLOCK DESIGNS WITH FACTORIAL STRUCTURE

UNIVERSITY OF TORONTO Faculty of Arts and Science. December 2005 Examinations STA437H1F/STA1005HF. Duration - 3 hours

BOOTSTRAP METHOD FOR TESTING OF EQUALITY OF SEVERAL MEANS. M. Krishna Reddy, B. Naveen Kumar and Y. Ramu

STAT 3008 Applied Regression Analysis

Topic 23 - Randomized Complete Block Designs (RCBD)

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis

Lecture 6: Introduction to Linear Regression

Lecture 6 More on Complete Randomized Block Design (RBD)

ANOVA. The Observations y ij

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

a. (All your answers should be in the letter!

Correlation and Regression. Correlation 9.1. Correlation. Chapter 9

Homework 10 Stat 547. Problem ) Z D!

Statistics II Final Exam 26/6/18

Department of Statistics University of Toronto STA305H1S / 1004 HS Design and Analysis of Experiments Term Test - Winter Solution

x yi In chapter 14, we want to perform inference (i.e. calculate confidence intervals and perform tests of significance) in this setting.

6.3.7 Example with Runga Kutta 4 th order method

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests

NANYANG TECHNOLOGICAL UNIVERSITY SEMESTER I EXAMINATION MTH352/MH3510 Regression Analysis

Comparison of Regression Lines

Selected Student Solutions for Chapter 2

1. Inference on Regression Parameters a. Finding Mean, s.d and covariance amongst estimates. 2. Confidence Intervals and Working Hotelling Bands

x = , so that calculated

A Comparative Study for Estimation Parameters in Panel Data Model

Lecture 4 Hypothesis Testing

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Y = β 0 + β 1 X 1 + β 2 X β k X k + ε

Joint Statistical Meetings - Biopharmaceutical Section

Durban Watson for Testing the Lack-of-Fit of Polynomial Regression Models without Replications

Systems of Equations (SUR, GMM, and 3SLS)

18. SIMPLE LINEAR REGRESSION III

STAT 3340 Assignment 1 solutions. 1. Find the equation of the line which passes through the points (1,1) and (4,5).

First Year Examination Department of Statistics, University of Florida

28. SIMPLE LINEAR REGRESSION III

Analysis of Covariance of Reinforced Balanced Incomplete Block Designs With a Single Explanatory Variable

Chapter 14 Simple Linear Regression

ECON 351* -- Note 23: Tests for Coefficient Differences: Examples Introduction. Sample data: A random sample of 534 paid employees.

Statistics for Economics & Business

EE215 FUNDAMENTALS OF ELECTRICAL ENGINEERING

ANSWERS CHAPTER 9. TIO 9.2: If the values are the same, the difference is 0, therefore the null hypothesis cannot be rejected.

Linearity. If kx is applied to the element, the output must be ky. kx ky. 2. additivity property. x 1 y 1, x 2 y 2

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential

Statistics Chapter 4

4 Analysis of Variance (ANOVA) 5 ANOVA. 5.1 Introduction. 5.2 Fixed Effects ANOVA

STAT 511 FINAL EXAM NAME Spring 2001

A New Refinement of Jacobi Method for Solution of Linear System Equations AX=b

2016 Wiley. Study Session 2: Ethical and Professional Standards Application

Estimation: Part 2. Chapter GREG estimation

MMA and GCMMA two methods for nonlinear optimization

Answers Problem Set 2 Chem 314A Williamsen Spring 2000

However, since P is a symmetric idempotent matrix, of P are either 0 or 1 [Eigen-values

STATISTICS QUESTIONS. Step by Step Solutions.

β0 + β1xi. You are interested in estimating the unknown parameters β

Econometrics of Panel Data

Number of cases Number of factors Number of covariates Number of levels of factor i. Value of the dependent variable for case k

Here is the rationale: If X and y have a strong positive relationship to one another, then ( x x) will tend to be positive when ( y y)

Now we relax this assumption and allow that the error variance depends on the independent variables, i.e., heteroskedasticity

Analysis of Variance and Design of Experiments-II

PHYS 215C: Quantum Mechanics (Spring 2017) Problem Set 3 Solutions

ISQS 6348 Final Open notes, no books. Points out of 100 in parentheses. Y 1 ε 2

Question 1 carries a weight of 25%; question 2 carries 20%; question 3 carries 25%; and question 4 carries 30%.

Department of Quantitative Methods & Information Systems. Time Series and Their Components QMIS 320. Chapter 6

Regression Analysis. Regression Analysis

Transcription:

Analyss of Varance and Desgn of Experments-II MODULE - III LECTURE - 8 PARTIALLY BALANCED INCOMPLETE BLOCK DESIGN (PBIBD) Dr Shalah Department of Mathematcs & Statstcs Indan Insttute of Technology Kanpur

Intralock analyss of PBIBD wth two assocates Consder a PBIBD under two assocates scheme The parameters of ths scheme are,, r, k, λ, λ, n, n, p, p, p, p, p p and y µ + β + τ + ε ;,,,,,,,, The lnear model nolng lock and treatment effects s µ s the general mean effect; β τ ε m s the fxed addte th lock effect satsfyng s the fxed addte th treatment effect satsfyng s the d random error wth ε N σ m ~ (0, ) β 0; r τ 0 and The PBIBD s a nary, proper and equreplcate desgn So n ths case, the alues of the parameters ecome n 0 or, k k for all,,, and r r for all,,, There can e two types of null hypothess whch can e consdered one for the equalty of treatment effects and another for the equalty of lock effects We are consderng here the ntralock analyss, so we consder the null hypothess related to the treatment effects only As done earler n the case of BIBD, the lock effects can e consdered under the nterlock analyss of PBIBD and the recoery of nterlock nformaton The null hypothess of nterest s H0: τ τ τ aganst alternate hypothess H : at least one par of τ s dfferent The null hypothess related to the lock effects s of not much practcal releance and can e treated smlarly

3 Ths was llustrated earler n the case of BIBD In order to otan the least squares estmates of the sum of squares due to resduals ( y µ β τ ) wth respect to µβ, and τ Ths results n the three normal equatons, whch can e unfed usng the matrx notatons The set of reduced normal equatons n matrx notaton after elmnatng the lock effects are expressed as Q Cτ C R N K N, ', Q V N K B R ri K ki ', Then the dagonal elements of C are otaned as n rk ( ) c r, (,,, ), k k the off-dagonal elements of C are otaned as µβ, and τ, we mnmze c λ f treatment and ' are the frst assocates k nn f treatment and ' are the second assocates ( ',,, ) k ' ' k λ

4 and the th alue n Q s [ th Q V Sum of lock totals n whch treatment occurs ] k rk ( ) τ nn ' τ k '( ') Next, we attempt to smplfy the expresson of Q S S Let e the sum of all treatments whch are the frst assocates of th treatment and e the sum of all treatments whch are the second assocates of th treatment Then τ + S + S τ Thus the equaton n Q usng ths relatonshp for,,, ecomes [ rk ] kq r( k ) τ ( λs + λs ) Imposng the sde condton rk ( ) τ λs λ τ τ S ( ) + λ τ + ( λ λ ) S λ τ a rk ( ) + λ and a λ λ τ 0, we hae kq r( k ) τ + λ τ + ( λ λ ) S a τ + S,,,,

These equatons are used to otan the adusted treatment sum of squares Let Q denote the adusted sum of Q ' soer the set of those treatment,s whch are the frst assocate of th treatment We note that when we add the terms for all, then occurs n tmes n the sum, eery frst assocate of occurs tmes n the sum and eery second assocate of occurs tmes n the sum wth p + p n Then usng K ki and we hae [ ] S S r( k ) + λ S + ( λ λ ) nτ + p S + p S rk ( ) + λ + ( λ λ )( p p S + ( λ λ ) pτ p kq k Q Sum of Q 's whch are the frst assocates of treatment a S + a τ ( λ λ ) p rk ( ) + λ + ( λ λ )( p p ) 5 p τ 0, Now we hae followng two equatons n kq and kq : kq a τ + S kq a τ + S Now solng these two equatons, the estmate of τˆ s otaned as k [ Q Q ] ˆ τ, (,,, ) a a

6 We see that Q Q 0, so ˆ τ 0 Thus τˆ s a soluton of reduced normal equaton The analyss of arance can e carred out y otanng the unadusted lock sum of squares as Block ( unad) B G k k, the adusted sum of squares due to treatment as ˆ τ Q Treat( ad) G y The sum of squares due to error as Error () t Total Block ( unad) Treat( ad) Total G y k

7 A test for H : τ τ τ 0 s then ased on the statstc F Tr Treat( ad) /( ) /( k + ) Error( unad) F Tr > F α ;, k + 0 If then H s reected The ntralock analyss of arance for testng the sgnfcance of treatment effects s gen n the followng tale: Source Sum of Degrees of Mean squares F squares freedom Between treatments (adusted) Treat( ad) ˆ τ Q dftreat MS Treat Treat( ad) df Treat MS Treat MSE Between locks (unadusted) Block ( unad) B k G k dfblock Intralock error Total Error () t ( By sustracton) Total G y k df ET k + dft k MSE df Error ET

8 Note that n the ntralock analyss of PBIBD analyss, at the step we otaned and elmnated to otan the followng equaton [ + ] + kq r( k ) λ τ ( λ λ ) S λ τ, S S another posslty s to elmnate S nstead of S If we elmnate S nstead of (as we approached), then the S soluton has less work noled n the summng of f < Q n n >, S Q Q If n n then elmnatng wll nole less work n otanng n otanng denotes the adusted sum of Q ' s s oer the set of those treatments whch are the second assocate of th treatment When we do so, the followng estmate of the treatment effect s otaned : ˆ τ k Q a Q a a rk ( ) + λ a λ λ ( λ λ ) p rk ( ) + λ + ( λ λ )( p p ) The analyss of arance s then ased on τˆ and can e carred out smlarly

9 The arance of the elementary contrasts of estmates of treatments (n case of n < n ) ( kq kq ) ( kq kq ) ˆ τ ˆ τ ' ' ' a a s k ( + ) a Var( ˆ τ ˆ τ ' ) k a a a f treatment and ' are the frst assocates f treatment and ' are the second assocates We osere that the arance of ˆ τ ˆ τ depends on the nature of and ' n the sense that whether they are the frst or ' second assocates So desgn s not (arance) alanced But arance of any elementary contrast are equal under a gen order of assocaton, z, frst or second That s why the desgn s sad to e partally alanced n ths sense