15.01.2013 Views

an introduction to generalized linear models - GDM@FUDAN ...

an introduction to generalized linear models - GDM@FUDAN ...

an introduction to generalized linear models - GDM@FUDAN ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

We also need expressions for the expected value <strong>an</strong>d vari<strong>an</strong>ce of the derivatives<br />

of the log-likelihood function. From (3.3), the log-likelihood function for<br />

a distribution in the exponential family is<br />

l(θ; y) =a(y)b(θ)+c(θ)+d(y).<br />

The derivative of l(θ; y) with respect <strong>to</strong> θ is<br />

dl(θ; y)<br />

U(θ; y) =<br />

dθ = a(y)b′ (θ)+c ′ (θ).<br />

The function U is called the score statistic <strong>an</strong>d, as it depends on y, it c<strong>an</strong><br />

be regarded as a r<strong>an</strong>dom variable, that is<br />

Its expected value is<br />

U = a(Y )b ′ (θ)+c ′ (θ). (3.13)<br />

E(U) =b ′ (θ)E[a(Y )] + c ′ (θ).<br />

From (3.9)<br />

E(U) =b ′ �<br />

(θ) − c′ (θ)<br />

b ′ �<br />

+ c<br />

(θ)<br />

′ (θ) =0. (3.14)<br />

The vari<strong>an</strong>ce of U is called the information <strong>an</strong>d will be denoted by I. Using<br />

the formula for the vari<strong>an</strong>ce of a <strong>linear</strong> tr<strong>an</strong>sformation of r<strong>an</strong>dom variables<br />

(see (1.3) <strong>an</strong>d (3.13))<br />

Substituting (3.12) gives<br />

I = var(U) = � b ′ (θ) 2� var[a(Y )].<br />

var(U) = b′′ (θ)c ′ (θ)<br />

b ′ (θ) − c′′ (θ). (3.15)<br />

The score statistic U is used for inference about parameter values in <strong>generalized</strong><br />

<strong>linear</strong> <strong>models</strong> (see Chapter 5).<br />

Another property of U which will be used later is<br />

var(U) =E(U 2 )=−E(U ′ ). (3.16)<br />

The first equality follows from the general result<br />

var(X) =E(X 2 ) − [E(X)] 2<br />

for <strong>an</strong>y r<strong>an</strong>dom variable, <strong>an</strong>d the fact that E(U) = 0 from (3.14). To obtain<br />

the second equality, we differentiate U with respect <strong>to</strong> θ; from (3.13)<br />

U ′ = dU<br />

dθ = a(Y )b′′ (θ)+c ′′ (θ).<br />

Therefore the expected value of U ′ is<br />

© 2002 by Chapm<strong>an</strong> & Hall/CRC<br />

E(U ′ ) = b ′′ (θ)E[a(Y )] + c ′′ (θ)<br />

= b ′′ �<br />

(θ) − c′ (θ)<br />

b ′ �<br />

+ c<br />

(θ)<br />

′′ (θ) (3.17)<br />

= −var(U) =−I

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!