next up previous

STAT 330 Lecture 14

Reading for Today's Lecture: 6.2, 8.1, 8.2

Goals of Today's Lecture:

Today's notes

Two general methods to develop estimators:

  1. Method of moments: set sample mean equal to population mean, sample second moment equal to population second moment and so on. Solve these equations for tex2html_wrap_inline121 . Need 1 equation for each component of tex2html_wrap_inline121 .
  2. Maximum likelihood: Choose tex2html_wrap_inline125 to be the value of tex2html_wrap_inline121 which maximizes

    displaymath129

    This is easier to understand for discrete data.

General framework for MLEs:

If tex2html_wrap_inline131 are independent and tex2html_wrap_inline133 is the probability mass function for an individual tex2html_wrap_inline135 and if we actually observe tex2html_wrap_inline137 then the likelihood function is

eqnarray16

Example: tex2html_wrap_inline131 independent Poisson tex2html_wrap_inline141 rvs. So

displaymath143

The likelihood is

displaymath145

It is easier to maximize the logarithm of this:

displaymath147

To do the maximization set the derivative with respect to tex2html_wrap_inline149 equal to 0. Get

displaymath151

from which we get

displaymath153

so that tex2html_wrap_inline155 is the mle (maximum likelihood estimate) of tex2html_wrap_inline149 .

Extension to Continuous Distributions: If tex2html_wrap_inline131 are independent and tex2html_wrap_inline161 is the probability density of an individual tex2html_wrap_inline135 (so that the X's are continuous rvs) then we can interpret the density as follows. Let tex2html_wrap_inline167 denote a small positive number. Then

eqnarray39

(The equation is the definition of density and the approximation is the meaning of integration over very small intervals -- it expresses the idea that integration is the opposite of differentiation.)

This prompts us to define

displaymath169

and

displaymath171

The MLE of tex2html_wrap_inline173 maximizes the Likelihood or equivalently the log likelihood.

Example: tex2html_wrap_inline131 a sample from tex2html_wrap_inline177 population.

displaymath179

eqnarray53

We need to set tex2html_wrap_inline181 and tex2html_wrap_inline183 to find estimates of tex2html_wrap_inline185 and of tex2html_wrap_inline187 .

We find:

displaymath189

and

displaymath191

Set tex2html_wrap_inline181 to get

displaymath195

so that tex2html_wrap_inline197 . Then put this solution in the second equation to get

displaymath199

which produces the solution

displaymath201

Notice that in the denominator there is an n and not n-1. The MLE is not quite the usual estimate of tex2html_wrap_inline207 .

Property of MLE

Suppose a different statistician used tex2html_wrap_inline209 where tex2html_wrap_inline211 as the parameters? What would tex2html_wrap_inline213 and tex2html_wrap_inline215 be? The log likelihood is now

displaymath217

When you take the derivatives and set them equal to 0 the tex2html_wrap_inline181 equation is unchanged and tex2html_wrap_inline197 . The derivative with respect to tex2html_wrap_inline223 gives the estimate (I leave you to do the algebra)

displaymath225

so that tex2html_wrap_inline227 . This is a general principle. If

displaymath229

is a transformation (for some function like say tex2html_wrap_inline231 or any other function) of the parameters then the mles satisfy

displaymath233

Power, tex2html_wrap_inline119 , Sample Size

Definition: the power function of a hypothesis test procedure is

displaymath237

Recall: Type I error is incorrect rejection. Type 2 error is incorrect acceptance.

Definition: The level of a test is the largest probability of rejection for tex2html_wrap_inline121 in the null hypothesis:

displaymath241

Typically if the null hypothesis is like tex2html_wrap_inline243 the value of tex2html_wrap_inline245 is actually tex2html_wrap_inline247 , that is, the edge of the null hypothesis gives the highest risk of incorrect rejection.

Definition: The probability of a type two error is a function tex2html_wrap_inline249 defined for tex2html_wrap_inline251 by

displaymath253

Thus for tex2html_wrap_inline121 not in the Null:

displaymath257


next up previous



Richard Lockhart
Thu Jan 29 16:11:58 PST 1998