21.05.2014 Views

Modelos para Dados de Contagem com Estrutura Temporal

Modelos para Dados de Contagem com Estrutura Temporal

Modelos para Dados de Contagem com Estrutura Temporal

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.

em que f t = F ′ ta t = E(log(λ t ) | D t−1 ) e q t = F ′ tR t F t = V ar(log(λ t ) | D t−1 ).<br />

Da análise conjugada, temos que<br />

E(λ t | D t−1 ) = α t<br />

β t<br />

(B.7)<br />

e<br />

V ar(λ t | D t−1 ) = α t<br />

. (B.8)<br />

βt<br />

2<br />

Po<strong>de</strong> ser mostrado que<br />

E(log(λ t ) | D t−1 ) = γ(α t ) − log(β t )<br />

(B.9)<br />

e<br />

V ar(log(λ t ) | D t−1 ) = ˙γ(α t ),<br />

(B.10)<br />

em que γ(·) e ˙γ(·) são, respectivamente, a função digama e função trigama. Entretanto,<br />

na prática, utilizar os valores reais da esperança em (B.9) e da variância (B.10) po<strong>de</strong><br />

gerar um gran<strong>de</strong> custo <strong>com</strong>putacional. Portanto, uma alternativa simples é aproximar<br />

E(log(λ t ) | D t−1 ) e V ar(log(λ t ) | D t−1 ) através <strong>de</strong> uma expansão em Taylor até primeira<br />

or<strong>de</strong>m <strong>com</strong> centro em E(λ t | D t−1 ). Assim, temos que<br />

e<br />

E(log(λ t ) | D t−1 ) ∼ = log(E(λ t | D t−1 )) = log<br />

(<br />

αt<br />

)<br />

= log(α t ) − log(β t ) (B.11)<br />

β t<br />

V ar(log(λ t ) | D t−1 ) ∼ = V ar(λ t | D t−1 )<br />

E(λ t | D t−1 ) 2 = 1 α t<br />

. (B.12)<br />

Resolvendo agora um sistema simples, encontramos os parâmetros da distribuição a priori<br />

conjugada <strong>para</strong> λ t da forma<br />

α t = 1 q t<br />

(B.13)<br />

e<br />

β t =<br />

1<br />

q t exp(f t ) .<br />

(B.14)<br />

148

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

Saved successfully!

Ooh no, something went wrong!