1 item foi adicionado ao seu carrinho.

>

ESTATA DOORDENADOR

Modelos de comutação Markov

>

>

>

Luzes altas

    >

  • Modelos de comutação Markov-modelo de transição
    • Modelo auto-regressivo
    • Modelo de regressão dinâmica
  • Parâmetros de regressão dependentes do estado
  • Parâmetros de regressão dependentes do estadoparâmetros de variância dependente
  • Tabelas de
    • Probabilidades de transição
    • Durações esperadas do estado
  • Previsões
    • Valores esperados da variável dependente
    • Probabilidades de estar num estado
    • Estático (umpasso)
    • Dinâmico (multistep)
    • RMSEs de previsões

De que se trata?

Por vezes, os processos evoluem ao longo do tempo com alterações discretas nos resultados.

Pense em recessões e expansões económicas. No início de uma recessão, a produção e o emprego caem e permanecem baixos e depois, mais tarde, a produção e o emprego aumentam. Pense nas desordens bipolares em que há períodos maníacos seguidos de períodos depressivos, e o processo se repete. Estatisticamente, as médias, variâncias e outros parâmetros estão mudando ao longo dos episódios (regimes). O nosso problema é estimar quando os regimes mudam e os valores dos parâmetros associados a cada regime. Perguntar quando os regimes mudam é equivalente a perguntar por quanto tempo os regimes persistem.

Nos modelos de Markov-transição, além de estimar as médias, variâncias, etc. de cada regime, estimamos também a probabilidade de mudança de regime. As probabilidades estimadas de transição para algum problema podem ser, as seguintes:


de/para
>estado 1 2
1 0.82 0,18
2 0,75 0,25

Iniciar no estado 1. A probabilidade de transitar do estado 1 para o estado 1 é de 0,82. Dito de forma diferente, uma vez no estado 1, o processo tende a permanecer lá. Com probabilidade 0,18, entretanto, o processo transita para o estado 2. O estado 2 não é tão persistente. Com probabilidade 0,75, o processo reverte do estado 2 para o estado 1 no próximo período de tempo.

Os modelos de comutação de Markov não estão limitados a dois regimes, embora modelos de dois regimes sejam comuns.

No exemplo acima, descrevemos a comutação como sendo abrupta; a probabilidade mudou instantaneamente. Tais modelos de Markov são chamados modelos dinâmicos. Os modelos de Markov também podem acomodar mudanças mais suaves modelando as probabilidades de transição como um processo autoregressivo.

Assim, a comutação pode ser suave ou abrupta.

Vejamos como funciona

Vejamos as mudanças médias entre regimes. Em particular, vamos analisar a Taxa Federal de Fundos. A Taxa Federal de Fundos é a taxa de juros que o banco central dos EUA cobra dos bancos comerciais por empréstimos overnight. Vamos analisar as mudanças na taxa de fundos federais de 1954 até o final de 2010. Aqui estão os dados:

graph

Temos dados trimestrais. As altas taxas de juros parecem caracterizar os anos setenta e oitenta. Vamos assumir que existe outro regime para taxas de juro mais baixas que parecem caracterizar as outras décadas.

Para encaixar um modelo de comutação dinâmica (mudança abrupta) com dois regimes, digitamos

. mswitch dr fedfundsPerforming EM optimizaton:Performing gradient-based optimization:
Iteração 0: probabilidade log = -508.66031
Iteração 1: verosimilhança log = -508.6382
Iteração 2: verosimilhança log = -508.63592
Iteração 3: verosimilhança log = -508.63592

Regressão dinâmica de comutação Markov Amostra: 1954q3 – 2010q4 Não. de obs = 226Número de estados = 2 AIC = 4.5455Probabilidades incondicionais: HQIC de transição = 4.5760 SBIC = 4.6211Probabilidade de log = -508,63592

fedfunds > Coef. Std. Err. z P>||z|
State1
_cons 3.70877 .1767083 20.99 0.000 3.362428 4.055112
State2
_cons 9,556793 .2999889 31,86 0,000 8,968826 10,14476
2,107562 .1008692 1,918851 2.314831
p11 .9820939 .0104002 .9450805 .9943119
p21 .0503587 .0268434 .0173432 .1374344

Reportados na saída acima são

  • as médias dos dois estados (_cons);
  • um único desvio padrão para todo o processo (sigma); e
  • as probabilidades de transição para o estado 1 para 1 e estado 2 para 1 (p11 e p21).

State1 é o estado de taxa moderada (média de 3,71%).

State2 é o estado de taxa alta (média de 9,56%).


de/para
estado 1 2
1 0.98 1 – 0,98
2 0,05 1 – 0,05

Ambos os estados são incrivelmente persistentes (1->1 e 2->2 probabilidades de 0,98 e 0,95).

Entre as coisas que você pode prever após a estimativa está a probabilidade de estar nos vários estados. Temos apenas dois estados, e assim a probabilidade de estar no (digamos) estado 2 diz-nos a probabilidade de estar nos dois estados. Podemos obter a probabilidade prevista e fazer um gráfico junto com os dados originais:

. predict prfed, pr

graph

O modelo tem pouca incerteza quanto ao regime em cada momento. Vemos três períodos de estados de taxa elevada e quatro períodos de estados de taxa moderada.

Vejamos o seu funcionamento

Vejamos um exemplo de surto de doença, nomeadamente papeira por 10.000 residentes na cidade de Nova Iorque entre 1929 e 1972. Você pode pensar que surtos correspondem a mudanças significativas, mas o que vemos nos dados é uma mudança ainda maior na variância:

graph

Pegamos a variável S12.mumpspc, significando casos de papeira per capita sazonalmente diferentes durante um período de 12 meses, e vamos analisar a S12.mumpspc.

Vamos assumir dois regimes em que a média e a variância da S12.mumpspc mudam. Para ajustar um modelo dinâmico (mudança brusca), vamos digitar

. mswitch dr S12.mumpspc, varswitch switch(LS12.mumpspc, noconstant)Performing EM optimizaton:Performing gradient-based optimization:
Iteração 0: verosimilhança log = 110,9372 (não côncava)
Iteração 1: verosimilhança log = 120.68028
Iteração 2: verosimilhança log = 123,23244
Iteração 3: verosimilhança log = 131,47084
Iteração 3: verosimilhança log = 131.72182
Iteração 3: verosimilhança log = 131,7225
Iteração 3: verosimilhança log = 131,7225

Regressão dinâmica de comutação Markov Amostra: 1929m2 – 1972m6 Não. de obs = 521Número de estados = 2 AIC = -0,4826Probabilidades incondicionais: HQIC de transição = -0,4634 SBIC = -0,4336Probabilidade de log = 131,7225

S12.mumspc Coef. Std. Err. z P>||z|
State1
mumpspc
LS12. .4202751 .0167461 25.10 0.000 .3874533 .4530968
State2
mumpspc
LS12. .9847369 .0258383 38.11 0.000 .9340947 1.035379
sigma1 .0562405 .0050954 .0470901 .067169
sigma2 .2611362 .0111191 .2402278 .2838644
p11 .762733 .0362619 .6846007 .8264175
p12 .1473767 .0257599 .1036675 .205294

>

Reportados são

  • os meios dos dois estados de S12.mumpspc (0,42 e 0,98);
  • os desvios padrão dos dois estados (0,06 e 0,26); e
  • as probabilidades de transição para o estado 1 para 1 e estado 2 para 1 (0,76 e 0,15).

o estado 1 é o estado de baixa variância.

O conjunto completo de probabilidades de transição é o seguinte:


Deixe uma resposta

O seu endereço de email não será publicado.