11.04.2013 Views

Modelos Lineares Generalizados em Experimentação Agronômica

Modelos Lineares Generalizados em Experimentação Agronômica

Modelos Lineares Generalizados em Experimentação Agronômica

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.

44 Clarice G.B. D<strong>em</strong>étrio<br />

2<br />

sendo que d i mede a diferença dos logaritmos das funções de verossimilhanças observada e<br />

ajustada, para a observação correspondente e é chamado componente da deviance. A soma<br />

deles mede a discrepância total entre as duas funções de verossimilhanças. É portanto, uma<br />

medida da distância dos valores ajustados ˆµ 's<br />

<strong>em</strong> relação aos dados observados y’s, ou de<br />

forma equivalente, do modelo corrente <strong>em</strong> relação ao modelo saturado. Verifica-se que a<br />

deviance eqüivale a uma constante menos duas vezes o máximo da função de verossimilhança<br />

para o modelo corrente, isto é,<br />

S = 2lˆ − 2lˆ<br />

= constante − 2lˆ<br />

.<br />

p<br />

n<br />

p<br />

2<br />

Ex<strong>em</strong>plo 11: Seja Y 1,<br />

Y2<br />

, ..., Yn<br />

uma amostra aleatória de uma distribuição N(<br />

µ i , σ ) ,<br />

2<br />

sendo que x β e σ > 0 , conhecida. Considerando como função de ligação a<br />

identidade, isto é, η = µ , t<strong>em</strong>-se que<br />

T<br />

i i = µ<br />

Logo,<br />

S<br />

p<br />

1<br />

= 2<br />

σ<br />

1<br />

2<br />

σ<br />

n<br />

∑<br />

i=<br />

1<br />

n<br />

= ∑<br />

i=<br />

1<br />

i<br />

⎧<br />

2⎨y<br />

⎩<br />

i<br />

2<br />

φ= σ ; wi = 1;<br />

θ i = µ i e<br />

i<br />

y<br />

µ ˆ<br />

⎫<br />

1<br />

σ<br />

2 2<br />

θ i µ i<br />

b(<br />

θ i ) = = .<br />

2 2<br />

2 2<br />

n<br />

i i<br />

2<br />

2 2<br />

[ y − ] − + ⎬ = ∑(<br />

− − + )<br />

i µ ˆ i<br />

2y<br />

i 2µ<br />

ˆ i yi<br />

yi<br />

µ ˆ<br />

2<br />

i<br />

2<br />

SQRes<br />

=<br />

σ<br />

2<br />

( y − ˆ ) i µ i<br />

2<br />

2<br />

⎭<br />

2<br />

que coincide com a estatística clássica SQRes com (n-p) graus de liberdade dividida por σ .<br />

Ex<strong>em</strong>plo 12: Sejam Y variáveis aleatórias representando contagens de sucessos <strong>em</strong> amostras<br />

i<br />

independentes de tamanhos m i . Supondo que Yi ~ Bin ( mi<br />

, πi<br />

) , então,<br />

e<br />

Logo,<br />

S<br />

ou ainda,<br />

p<br />

φ = 1;<br />

w = 1;<br />

i<br />

θ<br />

i<br />

i=<br />

1<br />

⎛ π ⎞ ⎛ µ ⎞<br />

i<br />

i<br />

= l n ⎜<br />

⎟ = ln<br />

⎜<br />

⎟<br />

⎝1<br />

- πi<br />

⎠ ⎝ mi<br />

- µ i ⎠<br />

θ<br />

⎛ m − µ ⎞<br />

i<br />

i i<br />

b(<br />

θ = + = − − π = − ⎜<br />

⎟<br />

i ) mi l n ( 1 e ) mi<br />

ln<br />

( 1 i ) mi<br />

ln<br />

.<br />

⎝ mi<br />

⎠<br />

⎪⎧<br />

i<br />

i<br />

i i<br />

i i<br />

∑ 2⎨y ⎢ ⎜<br />

⎟ −<br />

⎜<br />

⎟<br />

⎟⎥<br />

+<br />

⎜<br />

⎟ −<br />

⎜<br />

⎟<br />

i l n ln<br />

mi<br />

ln<br />

mi<br />

ln<br />

⎬<br />

i= 1 mi<br />

− yi<br />

mi<br />

− µ ˆ i<br />

mi<br />

mi<br />

⎪⎭<br />

= n<br />

⎪⎩<br />

⎡<br />

⎣<br />

⎛<br />

⎝<br />

S<br />

p<br />

y<br />

⎞<br />

⎠<br />

⎛<br />

⎝<br />

µ ˆ<br />

⎞⎤<br />

⎠⎦<br />

⎝<br />

p<br />

⎛ m − y<br />

⎡ ⎛ y ⎞<br />

⎛ m − y ⎞⎤<br />

= 2 l l .<br />

n<br />

i<br />

i i<br />

∑ ⎢y<br />

⎜<br />

⎟ + −<br />

⎜<br />

⎟<br />

i n ( mi<br />

yi<br />

) n ⎥<br />

i= 1 ⎢⎣<br />

⎝ µ ˆ i ⎠<br />

⎝ mi<br />

− µ ˆ i ⎠⎥⎦<br />

⎞<br />

⎠<br />

⎛ m − µ ˆ ⎞⎪⎫<br />

⎝<br />

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

Saved successfully!

Ooh no, something went wrong!