next up previous

STAT 330 Lecture 29

Reading for Today's Lecture: 12.1, 12.2,12.3

Goals of Today's Lecture:

Today's notes

In the model

displaymath93

if the levels of tex2html_wrap_inline95 are the only levels of interest of Factor 1 we call Factor 1 (and the tex2html_wrap_inline97 ) a fixed effect. If on the other hand they are a sample of size I from a population of possible levels we refer to Factor 1 as a random effect. Often Randomized Blocks designs have blocks which are regarded as random. For instance in an experiment where 5 runs of some production process can be run on a single day we often treat DAY as a blocking factor and then pretend the days we tried are a sample of possible days.

We call the Model a Fixed Effects model if both factors are fixed, a Random Effects model if both are random and a Mixed model if we have one fixed and one random factor. For mixed models with replicates we get different F tests for main effects. Moreover, the injunction that we test main effects only when there are no interactions is no longer relevant.

END OF CHAPTER 11

Simple Linear Regression and Correlation

Here are two experimental designs used to investigate the relation between two continuous variables.

1: Controlled Experiment: A variable, X is set at values tex2html_wrap_inline105 and corresponding values tex2html_wrap_inline107 of a response variable are measured.

Example: Chapter 12, question 9. x is the "Burner area liberation rate" and Y is the tex2html_wrap_inline113 (nitrous oxides) emission rate.

2: A sample of n pairs: We sample a population of n pairs of numbers and get tex2html_wrap_inline119 .

Example: we sample 1074 families and measure the Father's height (X) and Son's height (Y) for each family.

In this section our goal is to predict Y from the value of X and not the other way around. We do not treat the variables symmetrically.

Regression Models:

We assume for each observation a model equation of the form

displaymath129

where

Assumptions:

  1. tex2html_wrap_inline143 are independent random variables with mean 0 and variance tex2html_wrap_inline145 .
  2. We often assume the tex2html_wrap_inline141 s are normally distributed.

Definition: The regression is called linear if tex2html_wrap_inline149 is a linear function of tex2html_wrap_inline151 . (This jargon is used also when each of X and tex2html_wrap_inline139 is a vector.)

Our example is Simple Linear Regression.

displaymath157

where the tex2html_wrap_inline159 are independent mean 0 homoscedastic errors. (Notice that the map tex2html_wrap_inline161 is a linear function of tex2html_wrap_inline163 . At the same time this model describes a straight line function of x.

Estimation

Estimation is based on least squares. We choose tex2html_wrap_inline167 to minimize

displaymath169

To minimize this we take the derivatives tex2html_wrap_inline171 and tex2html_wrap_inline173 and set them both equal to 0. We get

displaymath175

and

displaymath177

These two equations are called the normal equations usually written in the form

displaymath179

and

displaymath181

The solution is

displaymath183

and

displaymath185

There is an ANOVA table for this least squares analysis based on the identity

displaymath187

where tex2html_wrap_inline189 is the so-called fitted value namely tex2html_wrap_inline191 . The quantity tex2html_wrap_inline193 is called the Error Sum of Squares and the quantity tex2html_wrap_inline195 is called the Regression Sum of Squares. We get the following ANOVA table.

In this table the P value is used to test tex2html_wrap_inline199 . However, for simple linear regression it is usually better to use a technique which easily provides confidence intervals for tex2html_wrap_inline139 and can be used to test other values of tex2html_wrap_inline139 .

Let tex2html_wrap_inline205 .and note that because tex2html_wrap_inline207

displaymath209

If the errors tex2html_wrap_inline159 are normal so that the tex2html_wrap_inline213 s are normal then

displaymath215

is normal and we can compute the mean and variance of tex2html_wrap_inline217 as follows:

eqnarray61

So tex2html_wrap_inline217 is unbiased. Next:

eqnarray78


next up previous



Richard Lockhart
Mon Mar 9 16:22:10 PST 1998