next up previous

STAT 330 Lecture 13

Reading for Today's Lecture: 4.4, 6.1, 6.2

Goals of Today's Lecture:

Today's notes

The Gamma distribution and gamma function:

Definition:

displaymath135

Properties:

displaymath137

and, if n is an integer

displaymath141

Note:

eqnarray15

Stirling's formula

displaymath143

The function

displaymath145

is called the Gamma density with scale parameter tex2html_wrap_inline147 and shape parameter tex2html_wrap_inline149 . (Note that tex2html_wrap_inline151 and tex2html_wrap_inline153 are required.)

Some sketches of Gamma densities ( tex2html_wrap_inline155 )

Exponential density tex2html_wrap_inline157 .

Gamma density with tex2html_wrap_inline159 , tex2html_wrap_inline161 .

Gamma density with tex2html_wrap_inline163 .

Gamma density with large tex2html_wrap_inline149 (approximately normal).

Method of moments example: example 6.12 in text. (Note: most examples lead either to trivial equations or to very difficult equations to solve or to very hard integrals to do; that's why I am duplicating the example in the text.)

Suppose tex2html_wrap_inline167 are iid with density

displaymath145

Population moments

First moment

displaymath171

Make the substitution tex2html_wrap_inline173 and tex2html_wrap_inline175 to get

eqnarray59

where we have used the recursive property of the Gamma function tex2html_wrap_inline177 .

Second moment

displaymath179

Make the substitution tex2html_wrap_inline173 and tex2html_wrap_inline175 to get

eqnarray76

To get method of moments estimates of tex2html_wrap_inline149 and tex2html_wrap_inline147 we solve the equations:

displaymath189

and

displaymath191

Divide the second equation by the square of the first to get

displaymath193

so that

displaymath195

Then

displaymath197

Notice the use of the bar over the top of tex2html_wrap_inline199 to indicate the average of the squares; this is quite different than averaging and then squaring.

Maximum Likelihood Estimation

New treatment for disease tried until fourth success. Resulting data: SFFSSS.

Problem: estimate tex2html_wrap_inline201 .

Likelihood:

displaymath203

which depends on the unknown parameter p.

displaymath207

Idea: to estimate p make L(p) as large as possible, that is choose as your estimate of p the one which assigns the largest possible probability to the data you observed.

Solution: set

displaymath215

and get the equation

displaymath217

or

displaymath219

or

displaymath221

so that

displaymath223

The Maximum Likelihood Estimate (MLE) of p is

displaymath227


next up previous



Richard Lockhart
Fri Jan 30 14:07:19 PST 1998