71

Click here to load reader

GMM, VI , MQ2S e Equações Simultâneas

Embed Size (px)

Citation preview

Page 1: GMM, VI , MQ2S e Equações Simultâneas

Método dos Momentos Generalizados (GMM), Variáveis Instrumentais (VI), Mínimos quadrados 2

Estágios (MQ2S) e Equações Simultâneas

Disciplina: Econometria I

Professor: Marcelo Savino Portugal

Elaborado por: Enrique, Fernanda Victor, Julio Cesar Junior

1

Page 2: GMM, VI , MQ2S e Equações Simultâneas

O método generalizado dos momentos (GMM)

• O GMM foi “melhorado” por Lars Hansen (1982), e trata-se de igualar um momento a um determinado valor. Em sua origem o método se associava diretamente aos momentos de uma distribuição. Porém, o método foi estendido para os parâmetros de momentos decorrentes de condições econômicas.

• Dito isso, “igualar um momento a um determinado valor não é mais do que satisfazer uma média amostral”(Bueno, 2008).

2Hansen, Lars, P. Large sample proprieties of generalized method of moments estimators. Econométrica, v.50, n.4, p.1029-1054, 1982

Page 3: GMM, VI , MQ2S e Equações Simultâneas

O método generalizado dos momentos (GMM)

• Exemplo simples:

(1)• Onde:

• β = é um escalar a ser estimado por algum método;• = pode representar, por exemplo, alguma defasagem de ou

algum termo exógeno que represente .

• A matriz de regressores tem dimensão n x K, onde n é o número de observações.

• Dado como um ruído branco, IID (idêntica e independentemente distribuído).

ttt xy 1

tx tyty

t

3

tx

Page 4: GMM, VI , MQ2S e Equações Simultâneas

O método generalizado dos momentos (GMM)

• Dadas as suposições anteriores é natural supormos que e são independentes entre si.

Formalizando: E( )=0. (2)

A equação (2) é exatamente uma condição de momento. A forma de obtê-la usando uma amostra finita de observações é calculando a seguinte média:

(3)

txttx

t

01

T

xT

t tt4

Page 5: GMM, VI , MQ2S e Equações Simultâneas

O método generalizado dos momentos (GMM)

• A equação (3) é a condição de momento para estimar . Além de propor o método para estimar esse parâmetro, Hansen (1982) derivou a distribuição de de tal modo que se pudessem fazer inferências estatísticas.

• Utilizando (3) e (2) temos:

=> (4)

1

1

0)(1

1

t

T

ttt xyx

T

t t

T

t tt

x

yx

1

2

11

^

5Hansen, Lars, P. Large sample proprieties of generalized method of moments estimators. Econométrica, v.50, n.4, p.1029-1054, 1982

Page 6: GMM, VI , MQ2S e Equações Simultâneas

O método generalizado dos momentos (GMM)

• O resultado de (4) é o estimador de MQO também, com isso mostramos que MQO é um caso especial de GMM.

• Contudo, empiricamente alguns dos regressores são endógenos, de forma que E(xtεt) ≠0, isso implicaria na impossibilidade de estimação via MQO, pois geraria viés. Porém, o método dos momentos permitem tal estimação, desde que exista outra variável correlacionada com xt e não correlacionada com εt

Segundo Bueno (2009), “O método dos momentos apresenta a vantagem de permitir trabalhar com total inexistência de hipóteses sobre a distribuição dos erros”. Em contrapartida, tem a desvantagem de possuir mal desempenho para pequenas amostras.

6

Page 7: GMM, VI , MQ2S e Equações Simultâneas

O método generalizado dos momentos (GMM)

• Se existe , tal que E( )≠0 e E( )=0, pode-se estimar β usando a seguinte média amostral:

• => =>

• (5)

7

tz tx t

0)(1

1

t

T

ttt xyz

T

t tt

T

t tt

xz

yz

1

11

^

01

T

zT

t tt

tz tz

Page 8: GMM, VI , MQ2S e Equações Simultâneas

O método generalizado dos momentos (GMM)

• Com o resultado anterior obtemos o estimador de variáveis instrumentais, IV, também conhecido como mínimos quadrados a dois ou três estágios, que também é uma forma especial do GMM.

• Nada impede de se usar mais de um instrumento para estimar β. Para fins de exposição, vamos supor o que chamaremos de “Caso 1”, com duas variáveis e , suas condições de momento são:

(6)8

tztm

0

0

tt

tt

z

mE

Page 9: GMM, VI , MQ2S e Equações Simultâneas

O método generalizado dos momentos (GMM)

• Utilizando a idéia de obter médias amostrais usando variáveis instrumentais,

• (7)

• O problema nesse caso é ter duas equações (condições de momento) e apenas uma incógnita (parâmetro).Em casos com mais de um instrumento, quantos devemos utilizar? 1, 2,3, ..., T?

• Não é razoável deixar de utilizar uma boa variável instrumental pois com isso perdemos informações. Então como utilizar mais de uma?

• Para solução desse tipo de problema Hansen (1982) propôs ponderar cada um dos momentos de modo a obter os parâmetros que minimizam uma função quadrática.

9

0)(1

1

t

T

ttt xyz 0)(

11

t

T

ttt xym

Page 10: GMM, VI , MQ2S e Equações Simultâneas

O método generalizado dos momentos (GMM)

• Chamamos o Caso 1 de sistema sobreidentificado, pois existem mais parâmetros que momentos.

• A idéia para solucionar o Caso 1 é mais ou menos a idéia de estimação de MQO, porém no caso dos momentos é a matriz de covariância dos momentos que serve como ponderador, ao invés da matriz de variância dos parâmetros.

• Essa ponderação tem uma intuição simples: os momentos de maior variância recebem o menor peso (ocorre naturalmente quando tomamos o inverso da matriz de covariância dos momentos).

10

Page 11: GMM, VI , MQ2S e Equações Simultâneas

O método generalizado dos momentos (GMM)

• A idéia então é achar uma matriz W, simétrica e positiva semi definida de forma a minimizar:

• (8)

• Em termos de médias amostrais significa:

• (9)11

][

tt

tttttt z

mWzmE

)(

)())()((min

11

11

11

11

1

t

T

ttt

t

T

ttt

t

T

tttt

T

ttt

xyz

xymWxyzxym

Page 12: GMM, VI , MQ2S e Equações Simultâneas

O método generalizado dos momentos (GMM)

• Supondo W seja uma matriz identidade (significando que os momentos teriam os mesmos pesos) o problema seria:

=>

(10)

Resolvendo a minimização (10), resulta em:

(11) 12

)(

)(

10

01))()((min

11

11

11

11

1

t

T

ttt

t

T

ttt

t

T

tttt

T

ttt

xyz

xymxyzxym

2

11

2

11 )()(min

1

t

T

tttt

T

ttt xyzxym

)(

)(

1 1

1 11

^

jtjt

T

t j

T

t t

T

t jtjtj

T

t t

zzmmxx

zzmmxy

Page 13: GMM, VI , MQ2S e Equações Simultâneas

(GMM) - EspecificaçãoPara formalização da especificação de GMM em termos de seus momentos

populacionais e contrapartida amostral, consideremos a seguinte definição:

Definição (Condição de momento populacional) Sejam: um vetor de variáveis aleatórias um vetor de parâmetros k x 1g (.) um vetor m x 1 de funções reais, a condição de momentos populacionais é

definida como:

(12)

A contrapartida amostral dessa condição de momento é:

(13)

twt

0[ ( , )] 0tE g w

01

0

( , )( , )

T

tt

T

g wg w

T

Page 14: GMM, VI , MQ2S e Equações Simultâneas

(GMM) - Especificação

• O estimador GMM será definido como o valor de que resolve o seguinte problema :

• (14)

Em que W é uma matriz m X m positiva semidefinida. Naturalmente o número de momentos deve ser superior ao número de parâmetros (sobreidentificados). Os sistemas podem ser classificados como:indentificados -> m = k, (m = momentos, k = parâmetros), Subidentificados -> k > mSobreidentificados -> m > k

0 0min ( ) ( , ) ' ( , )T T T TJ g w W g w

arg min ( )GMM

TJ

Page 15: GMM, VI , MQ2S e Equações Simultâneas

(GMM) - Estimação

• Diferenciar a função objetivo (14) em relação ao parâmetro • com isso temos a condição de primeira ordem (15)

• (15)

• A matriz , de dimensão m X k, é o jacobiano das condições de momento tomadas em relação a ϴ

• (16)

0 02 ( , ) ' ( , ) 0T T TG w W g w

1 10

( , )( , )

'

T

tt

T

g wG w T

t

),( wGt

Page 16: GMM, VI , MQ2S e Equações Simultâneas

(GMM) - Estimação

• As condições de primeira ordem podem ser não-lineares. Logo, o problema pode não ter solução analítica, precisando de métodos numéricos para solução. Obter um mínimo global pode ser computacionalmente intensivo.

• Uma pergunta persiste: qual a matriz de pesos ideal?

• Defina a matriz de covariância de longo prazo dos momentos:

• (17)1

20lim var( ( , ))TT

S T g w

Page 17: GMM, VI , MQ2S e Equações Simultâneas

(GMM) - Estimação• Assim, a matriz ideal de pesos é:

• A intuição é: Momentos com variância maior devem ter peso menor. Por isso a inversa.

• Dado que a obtenção das covariâncias depende do parâmetro que é, na maioria das vezes, desconhecido. Com isso, será preciso usar um sistema de estimação com múltiplos estágios.

• Primeiro estágio: impõe-se uma matriz W;

• Seleciona-se o parâmetro resultante do estágio 1. Depois calcula-se a matriz de covariância . É comum selecionar uma matriz identidade para iniciar o processo.

1TW S

t

t1

^

S

Page 18: GMM, VI , MQ2S e Equações Simultâneas

(GMM) - Estimação• Admitindo a existência de um método consistente para S, há três

métodos de para estimação GMM.

• A. método a dois estágios: impõe-se uma matriz identidade no primeiro passo e, no passo seguinte, encontram-se os parâmetros do modelo.

• B. método seqüencial: é a repetição do método de dois estágios, em que cada novo vetor de parâmetros é usado para estimar a nova matriz de covariância.

• C. método de atualização contínua dos parâmetros: a matriz de covariância é simultaneamente estimada com o vetor de parâmetros. 18

Page 19: GMM, VI , MQ2S e Equações Simultâneas

(GMM) - Estimação• No caso do método A (dois estágios) teremos

o qual será utilizado no segundo estágio.• O caso B (procedimento seqüencial ), tende a ser mais robusto em

pequenas amostras. (caso o teste a dois estágios seja robusto o procedimento seqüencial também será).

• O caso C (atualização contínua), o modelo encontra um ϴ que satisfaz

0 0( , ) ' ( , ) 0T TG w g w

0 0

10

( , ) ( , ) '( , )

T

t tt

g w g wS w

T

10 0min ( ) ( , ) ' ( , )T T TJ g w S g w

Page 20: GMM, VI , MQ2S e Equações Simultâneas

(GMM) - Estimação

• O caso C é de complexa solução. Se o estimador converge, trata-se de um método cujas propriedades de pequenas amostras são melhores. Com isso, Shepard (2006) sugere utilizar o estimador seqüencial para encontrar os valores iniciais da função a ser minimizada e depois utilizar o estimador de atualização contínua.

20

Page 21: GMM, VI , MQ2S e Equações Simultâneas

(GMM) – Propriedades do estimador

• Consistente e assintoticamente Normal se

• São ergódicos e possuem poucos momentos. • Outro ponto: para que o problema de otimização possua

solução única é necessário que a matriz W seja positiva definida. Caso W seja eficiente:

21

0( , )g w

1 10( ) (0, ( ' ) )dT N G S G

Page 22: GMM, VI , MQ2S e Equações Simultâneas

Testes: GMM• O GMM minimiza uma função representando as condições de

momentos devidamente ponderadas. Se as condições estiverem corretas, vão ter média zero. Assim, o teste de superidentificação será:

• Observe que o grau de liberdade é igual ao numero de excesso de momentos em relação aos parâmetros estimados.

1 20 0( , ) ' ( , ) ~T T m kJ Tg w S g w

Page 23: GMM, VI , MQ2S e Equações Simultâneas

Testes: GMM

• O que significa a rejeição de superidentificação? Que existem momentos que não são estatisticamente iguais a zero. Ou seja, rejeita-se o modelo, pois a condição de momento não é válida.

• Ou ainda, podemos interpretar da seguinte forma: se a inclusão de mais um momento não implica em rejeição da hipótese nula, esse momento é valido e contribui para estimar os parâmetros.

23

Page 24: GMM, VI , MQ2S e Equações Simultâneas

Variáveis Instrumentais e o Estimador IV

Tal como o MQO, o estimador IV (variáveis instrumentais) pertence à classe de estimadores GMM. O IV é um GMM linear e

a sua utilização justifica-se principalmente em modeloscom endogeneidade:

0)(),( iijiij uxEuxCov (1)

Para algum j = 1,....., k e i = 1, .... , n

onde se mantém a hipótese de E(u) = 0. Na estimação ao IV é necessário recorrer ao uso de

variáveis instrumentais a(s) qual(is) estão associadas a um particular regressor endógeno.

Page 25: GMM, VI , MQ2S e Equações Simultâneas

Definição 1: z é uma variável instrumental para x se:

(i) z é exógeno em relação ao erro:

(ii) z está (fortemente) correlacionado com o regressor x:

Variáveis Instrumentais e o Estimador IV

0)(),( uzEuzCov

0),( xzCov (3)

(2)

Em relação à primeira, e como não observamos u, partimos apenas de pressupostos da teoria econômica ou senso comum. A análise da correlação amostral entre z e os resíduos é um tanto

quanto simplista e não rigorosa.

Page 26: GMM, VI , MQ2S e Equações Simultâneas

26

Para verificar a segunda condição e qual a intensidade da relação entre z e x, uma das possibilidades é usar um teste

t na regressão auxiliar:

)(

),(; 221 zV

zxCovuzx

Após a definição da variável instrumental é muito simples a expressão do estimador IV. No modelo de Regressão Linear

Simples:

Variáveis Instrumentais e o Estimador IV

(4)

0)(),(;21 uxExuCovuxy (5)

Page 27: GMM, VI , MQ2S e Equações Simultâneas

Para uma amostra de dimensão n e um único instrumento para x,

Variáveis Instrumentais e o Estimador IV

))((

))((

1

1,2

xxzz

yyzz

ii

n

i

ii

n

iIV (6)

Naturalmente:

xy IVIV ,2,1

Ao contrário do MQO, o IV não resulta na minimização dos quadrados dos erros. A sua expressão pode no entanto ser

justificada pelo método dos momentos.

(7)

Page 28: GMM, VI , MQ2S e Equações Simultâneas

28

Variáveis Instrumentais e o Estimador IV

Da expressão (5):

),(

),(),(

),(),(

),(),(

22

2

21

zxCov

zyCovzxCov

zuCovzxCov

zuxCovzyCov

Alternativamente, resulta dos momentos:

0)( uE 0)( uzEe

(8)

Page 29: GMM, VI , MQ2S e Equações Simultâneas

29

Quando a variável instrumental é válida, o estimador IV é consistente1lim

IVn

p

Ou seja, num modelo com endogeneidade, 0),( xuCov

22,

22,

),(

),(lim

),(lim

x

uIVs

n

x

uMQOs

n

zxCorr

zuCorrp

xuCorrp

Variáveis Instrumentais e o Estimador IV

Em um instrumento válido, essa correlação é zeroEm um instrumento válido, essa correlação é diferente de zero

(9)

Page 30: GMM, VI , MQ2S e Equações Simultâneas

30

O IV tem um menor viés assimptótico do que o MQO se:

Variáveis Instrumentais e o Estimador IV

),(),(

),(uxCorr

xzCorr

uzCorr (10)

Esta condição não é fácil de testar a não ser de uma forma intuitiva em que se calculam as correlações amostrais onde u é

representado pelos resíduos de .

u

Quando essa condição não é satisfeita, dizemos que Z é uma IV pobre.

Page 31: GMM, VI , MQ2S e Equações Simultâneas

31

Variáveis Instrumentais e o Estimador IV

Sob a hipótese de que : )()/( 22 uVzuE

22

2

,2 ),(),/(

zxCorrnzxV

xVI

2,

2

2

,2 ),/(

zxx

VI

Rn

zxV

Que pode ser consistentemente estimada por:

(11)

(12)

2

1

2

2

1

un

niOnde: e são os resíduos de IV.

u

Page 32: GMM, VI , MQ2S e Equações Simultâneas

32

Variáveis Instrumentais e o Estimador IVIsso significa que se X e Z são linearmente dependentes:

E o estimador de VI será igual ao estimador inconsistente por MQO.

Em geral, isso não acontece, e porque:

O IV é menos eficiente (maior que que MQO). Mas, o importante é que ao contrário do IV, o MQO é inconsistente sob a hipótese de endogeneidade. Quanto maior é a correlação entre X e Z (melhor a qualidade do instrumento), maior é a precisão na estimação do IVdo modelo.

2,

2

. zxx RSQT

10 2, zxR

Igual a 1 (13)

Page 33: GMM, VI , MQ2S e Equações Simultâneas

EXEMPLO: ESTIMAÇÃO DO RETORNO DO SALÁRIO PARA MULHERES CASADAS

Primeiro obtemos as estimativas por MQO:

educsalário 109,0185,0)log(

(0,185) (0,014)118,0,428 2 Rn

Em seguida, usamos e educação do pai (educp) como variável instrumental de educ. Temos que sustentar que educp é não correlacionada com o erro, e testamos

se ela é correlacionada com educ usando uma regressão simples:

(20)

educpeduc 269,024,10

(0,28) (0,029)

173,0,428 2 Rn

(21)

eaptideducsalário 210)log(

eeducsalário 10)log(

Sendo a estatística t de educp 9,28, indica que educp tem ma correlação positiva estatisticamente significativa com educ.

Page 34: GMM, VI , MQ2S e Equações Simultâneas

EXEMPLO: ESTIMAÇÃO DO RETORNO DO SALÁRIO PARA MULHERES CASADAS

Utilizar educp como uma VI de educ produz:

educsalário 059,0441,0)log(

(0,446) (0,035)093,0,428 2 Rn

A estimativa de VI do retorno da educação é 5,9%, que é cerca de metade da estimativa MQO. Isso sugere que a estimativa de MQO é alta demais e é

consistente com o viés de aptidão omitida. Entretanto, devemos lembrar que essas estimativas são de apenas uma amostra: nunca podemos saber se 0,109 está acima do verdadeiro retorno da educação, ou se 0,059 está mas próximo

do verdadeiro retorno da educação. Além disso, o erro-padrão da estimativa VI é duas vezes e meia maior que o erro-padrão de MQO.

O intervalo de confiança de 95% de , utilizando MQO, é muito mais justo do que utilizando VI; de fato, o intervalo de confiança de VI, na realidade, contém a estimativa de MQO. Mais adiante será demonstrado como saber se a diferença

entre os resultados é estatisticamente significante.

(22)

1

Page 35: GMM, VI , MQ2S e Equações Simultâneas

Mínimos Quadrados de dois Estágios MQ2E

Considere o seguinte modelo estrutural:

0;0; 43241321 vzzzx

0),(;0),(;321 zuCovxuCovuzxy

Em que x é endógeno, z é exógeno, e z1, z2 são instrumentos válidos para x. Num primeiro passo, estima-se consistentemente por MQO:

(14)

(15)

Em que x tem como regressores todas as variáveis exógenas (instrumentos e do modelo estrutural). Como a combinação linear de

todas as variáveis exógenas:

24132121 ),,/( zzzxzzzxE

(16)

Page 36: GMM, VI , MQ2S e Equações Simultâneas

36

Mínimos Quadrados de dois Estágios MQ2E

No segundo passo, estima-se consistentemente por MQO o modelo estrutural (14)

uzxy

321 (17)

onde é determinado no primeiro passo e o estimador resultante para β coincidente com o IV.

x

Podemos interpretar essa regressão da seguinte forma: o valor estimado é a versão estimada de X , que por sua vez é não

correlacionado com u. Portanto, primeiro o MQ2E expurga X de sua correlação com u antes de fazer a regressão por MQO.

x

Page 37: GMM, VI , MQ2S e Equações Simultâneas

Teste de endogeneidadeTorna-se necessário nesta discussão apresentar testes quanto à

exogeneidade dos regressores. No caso de existir mais de um instrumento por regressor endógeno, também podemos testar a validade de algum dos instrumentos no sentido de não estarem

correlacionados com os erros.

O teste de Hausman é o procedimento de referência para inferir sobre a

exogeneidade dos regressores. Em resumo, o teste de Hausman procura

comparar estatisticamente dois estimadores , para o mesmo modelo e

vetor de parâmetros β. O estimador é consistente sob ambas as

hipóteses , enquanto que o estimador apenas é

consistente sob a hipótese nula. Por outro lado, sob a hipótese nula,

é (assimptoticamente) mais eficiente do que o .

VI

MQO

10;HH

MQO

0H

VI

Page 38: GMM, VI , MQ2S e Equações Simultâneas

Se cada z é não-correlacionado com o erro, x será não-correlacionado com o erro se, e somente se, v não o for. Logo, u e v são não-

correlacionados se, e somente se,

uvzxy

321(18)

(19)0

Teste de endogeneidadeVejamos o teste de Hausman com o recurso a regressão auxiliar.

Dado o fato de que: , os erros u e v dos modelos (14) e (15) estão auto-correlacionados, ou seja:

O teste à exogeneidade é simplesmente o teste t ao coeficiente δ na regressão auxiliar (aumentada em relação ao modelo estrutural

em ):

0),( xuCov

0),( vuCov

v

Page 39: GMM, VI , MQ2S e Equações Simultâneas

Teste de endogeneidadeSe existirem mais do que uma variável endógena (por exemplo, 5) então usa-se o teste F da significância conjunta dos 5 coeficientes

associados aos 5 resíduos do primeiro passo. Para testar se um particular regressor x, que tem associado o

parâmetro β2 no modelo, não é endógeno:

01

lim0)(),(:1

0

xiui

n

in npuxExuCovH

21

,2,2

2,2,2

)()(

)(

d

MQOIV

MQOIV

VV

pode-se utilizar a estatística:

(20)

(21)

Page 40: GMM, VI , MQ2S e Equações Simultâneas

Teste de restrições sobreidentificadasPara finalizar, apresentamos um teste quanto à auto-correlação de

erros e instrumentos onde o número de instrumentos é pelo menos dois. Num contexto IV ou GMM linear, procura-se testar a

hipótese nula de:

em que os m > 1 instrumentos para x não estão auto-correlacionados com u. Se a hipótese nula não é aceita, concluímos

que pelo menos um dos instrumentos não é válido e está correlacionado com os erros u. O modelo diz-se sobre-identificado

pois m>1. No caso de exata identificação, m = 1, não é possível desenvolver este teste.

10 ,0)(: mxZuZEH (22)

Page 41: GMM, VI , MQ2S e Equações Simultâneas

Teste de restrições sobreidentificadas

ezzzzu mm

......221121

O teste é similar ao anterior de Hausman em que se recorre a

resíduos. Neste caso, no entanto, a estatística de teste é:

em que o resulta do modelo auxiliar.

onde são os resíduos IV do modelo estrutural, e z1, z2,

..., zm são os (possíveis) instrumentos para x. A estatística é

assimptoticamente distribuída como: 21m

u

2nR 2R

(23)

Page 42: GMM, VI , MQ2S e Equações Simultâneas

Modelos de Equações Simultâneas - MESCaracterísticas:

-Modelo onde existem duas ou mais variáveis dependentes em vez de apenas uma.

-Diferem da maioria dos modelos econométricos porque eles consistem de um conjunto de equações.

-Procedimento de estimação por MQO não é mais apropriado nesses modelos para obter estimativas confiáveis dos parâmetros econômicos.

Page 43: GMM, VI , MQ2S e Equações Simultâneas

A oferta e a demanda conjuntamente determinam o preço de mercado de um bem e a quantidade que é vendida desse produto.

Um modelo econométrico que explique o preço de mercado e a quantidade vendida deve ter duas equações, uma para a S e outra para D.

1 2Demanda: dq p y e

1Oferta: sq p e

p e q são chamadas de variáveis endógenas A variável renda y é uma variável exógena.

MES – Exemplo da Oferta e Demanda

Page 44: GMM, VI , MQ2S e Equações Simultâneas

Acrescentamos Erros aleatórios

2

2

( ) 0, var( )

( ) 0, var( )

cov( , ) 0

d d d

s s s

d s

E e e

E e e

e e

O fato de que p é aleatório significa que do lado direito das equações de oferta e demanda nós temos uma variável explanatória

que é aleatória.

p, , ed e es, são correlacionados, fazendo com que o estimador de MQO seja viesado e inconsistente.

MES – Exemplo da Oferta e Demanda

Page 45: GMM, VI , MQ2S e Equações Simultâneas

s

d

epq

eypq

1

21

Resolvendo as equações estruturais

1 1 2s dp e p y e

Igualando as equações de demanda e oferta, encontramos a forma reduzida, simultaneamente para p e q.

2

1 1 1 1

1 1

d se ep y

y v

MES – Equações de forma reduzida

Page 46: GMM, VI , MQ2S e Equações Simultâneas

Substituindo p na equação de oferta

21 1

1 1 1 1

1 2 1 1

1 1 1 1

2 2

d ss s

d s

e eq p e y e

e ey

y v

1 e 2 são os parâmetros de forma reduzida.

v1 e v2 são os erros de forma reduzida

MES – Equações de forma reduzida

Page 47: GMM, VI , MQ2S e Equações Simultâneas

As equações de forma reduzida (1 e 2. ) podem ser estimadas consistentemente por MQO.

As equações de forma reduzida são importantes para análise econômica.Essas equações relacionam os valores de equilíbrio das variáveis endógenas com as variáveis exógenas. Assim, se existir um aumento na renda y, é o aumento esperado no preço, depois dos ajustamentos de mercado levarem ao novo

equilíbrio para p e q.

Em segundo lugar, as equações de forma reduzida podem ser utilizadas para prever valores de equilíbrio do preço e quantidade

para diferentes níveis de renda.

1

MES – Equações de forma reduzida

Page 48: GMM, VI , MQ2S e Equações Simultâneas

MES – Falha de estimação de MQO

Porque não se deve usar os estimadores de MQO para estimar a equação de Oferta no MES?

Porque o termo p e correlacionado com o erro e

Page 49: GMM, VI , MQ2S e Equações Simultâneas

EXPLICAÇÃO: Suponha que exista uma pequena variação no termo de erro es, digamos es. A variação es no termo de erro é diretamente transmitida para o valor de equilíbrio de p.

Observando a equação reduzida vemos que dada um es ocorre um efeito direto linear sobre p.

Como e , se es > 0 , então .

Assim, toda vez que existe uma variação em es, existe uma mudança associada em p na direção oposta. Conseqüentemente, p

e es são negativamente correlacionados.

1 0 1 0 0p

)()( 1111

2

sd eeyp

MES – Falha de estimação de MQO

Page 50: GMM, VI , MQ2S e Equações Simultâneas

EXPLICAÇÃO: A estimação por mínimos quadrados da relação entre q e p atribui ao preço o crédito pelo efeito de variações nas

perturbações. Em amostras grandes, o estimador de mínimos quadrados tendem

a ser negativamente viesados. O viés persiste mesmo quando o tamanho da amostra é grande e,

assim, o estimador de mínimos quadrados é inconsistente.

O estimador de mínimos quadrados dos parâmetros na equação simultânea estrutural é viesado e inconsistente em

virtude da correlação entre o erro aleatório e a variável endógena no lado direito da equação.

MES – Falha de estimação de MQO

Page 51: GMM, VI , MQ2S e Equações Simultâneas

Uma Condição Necessária para Identificação: Em um sistema de M equações simultâneas, que conjuntamente determinam os valores de M variáveis endógenas, pelo menos M1 variáveis devem estar ausentes de uma equação para que a estimação de seus parâmetros seja possível. Quando a estimação dos parâmetros de uma equação é possível de ser feita, então se diz que a equação é identificável e seus parâmetros podem ser consistentemente estimados. Se menos do que M1 variáveis estiverem fora da equação, então se diz que ela não é identificável e seus parâmetros não podem ser estimados consistentemente.

MES - Identificação

Page 52: GMM, VI , MQ2S e Equações Simultâneas

M=2 equações e 3 (três) variáveis: p, q e y.

Nenhuma variável é omitida Não é identificável e seus parâmetros não podem ser estimados consistentemente.

M1=1 variável, renda é omitida; a curva de oferta é identificável e seus parâmetros podem ser estimados.

A condição de identificação deve ser verificada antes de se tentar estimar uma equação.

1 2Demanda: dq p y e

1Oferta: sq p e

MES - Identificação

Page 53: GMM, VI , MQ2S e Equações Simultâneas

Procedimento de Estimação por MQ2E

1 sq p e A variável p é composta de uma parte sistemática, que é o seu valor esperado E(p), e de uma parte aleatória, que é o erro aleatório da forma reduzida v1.

1 1 1( )p E p v y v

Método mais utilizado para estimar os parâmetros de uma equação estrutural identificada.

1 *ˆ ˆq p e

se conhecermos o valor de 1 então podemos substituir p por seu valor

1 1

1 1 1

1 *

[ ( ) ]

( ) ( )

( )

s

s

q E p v e

E p v e

E p e

Page 54: GMM, VI , MQ2S e Equações Simultâneas

Regras de Identificação 1

Modelo com 2 variáveis endógenasNenhuma variável predeterminada M-1 = 1Nenhuma equação identificada

ttt

ttt

uPQOferta

uPQDemanda

210

110

Page 55: GMM, VI , MQ2S e Equações Simultâneas

Regras de Identificação 2

Modelo com 2 variáveis endógenas (P,Q)I é exógenaFunção oferta esta identificada

ttt

tttt

uPQOferta

uIPQDemanda

210

1210

Page 56: GMM, VI , MQ2S e Equações Simultâneas

Regras de Identificação 3

Modelo com 2 variáveis endógenas (P,Q)I e Pt-1 são exógenosFunção oferta e demanda estão identificadas portanto o modelo está identificado

tttt

tttt

uPPQOferta

uIPQDemanda

21210

1210

Page 57: GMM, VI , MQ2S e Equações Simultâneas

Regras de Identificação 4

Modelo com 2 variáveis endógenas (P,Q)I, R e Pt-1 são exógenos

Função oferta e demanda estão identificadas portanto o modelo está sobre-identificado em excesso

tttt

ttttt

uPPQOferta

uRIPQDemanda

21210

13210

Page 58: GMM, VI , MQ2S e Equações Simultâneas

•Como E(p) não esta correlacionada com e podemos estimar por mínimos quadrados 1.

•Nós podemos estimar 1 utilizando da equação de forma reduzida para p. •Um estimador consistente para E(p) é

1ˆ ˆ p y• Utilizando como um substituto para E(p)

1 *ˆ ˆq p e p̂

Procedimento de Estimação por MQ2E

Page 59: GMM, VI , MQ2S e Equações Simultâneas

Mínimos Quadrados de Dois Estágios

Resumindo os dois estágios do processo de estimação

1. Estimação por MQO da equação reduzida de p, e o calculo de seu valor predito

2. a estimação por MQO da equação estrutural em que a variável endógena do lado direito p é substituída por seu valor predito

Page 60: GMM, VI , MQ2S e Equações Simultâneas

Equações Simultâneas, MQ2E e Variáveis Instrumentais - EXEMPLO

• Um modelo de oferta e demanda para trufas:

- q é a quantidade de trufas comercializada - p é o preço de mercado das trufas, - ps é o preço de mercado de um substituto- di é renda per capita disponível. - Assumimos que p e q são variáveis endógenas. - As variáveis exógenas são ps, di, pf e a variável de

intercepto.

stttt

dttttt

epfpqOferta

edipspqDemanda

321

4321

Page 61: GMM, VI , MQ2S e Equações Simultâneas

A regra para identificar uma equação é que, em um sistema de M equações, pelo menos M1 variáveis devem estar

omitidas de cada equação para ela ser identificável.

Na equação de demanda, a variável pf não está incluída Na equação de oferta, ps e di estão ausentes

Nós concluímos que cada equação nesse sistema é identificável e pode, assim, ser estimada por mínimos

quadrados em dois estágios.

MES, MQ2E e VI – Exemplo - Identificação

Page 62: GMM, VI , MQ2S e Equações Simultâneas

As equações de forma reduzida expressam cada variável endógena, p e q, em termos das variáveis exógenas ps, di, pf e

a variável de intercepto, mais um termo de erro.

11 21 31 41 1

12 22 32 42 2

t t t t t

t t t t t

q ps di pf v

p ps di pf v

MES, MQ2E e VI – Exemplo – Equações da Forma Reduzida

Page 63: GMM, VI , MQ2S e Equações Simultâneas

pfdipsq 507,0217,0656,0895,7

MES, MQ2E e VI – Exemplo – Equação na forma reduzida para a quantidade (q)

Page 64: GMM, VI , MQ2S e Equações Simultâneas

MES, MQ2E e VI – Exemplo – Equação na forma reduzida para o preço (p)

pfdipsp 0.4513 0.2534 0.5693 -10.8374

Page 65: GMM, VI , MQ2S e Equações Simultâneas

As equações de forma reduzida são utilizadas para obter q e p predito, que serão usados no lugar de q e p no lado direito das equações de oferta e demanda no segundo estágio dos mínimos quadrados em dois estágios.

MES, MQ2E e VI – Exemplo – Equações de forma Reduzida

Page 66: GMM, VI , MQ2S e Equações Simultâneas

MES, MQ2E e VI – Exemplo – Estimando a Oferta por MQ2E

Page 67: GMM, VI , MQ2S e Equações Simultâneas

MES, MQ2E e VI – Exemplo – Estimando a Demanda por MQ2E

Page 68: GMM, VI , MQ2S e Equações Simultâneas

INTERPRETAÇÃO DA EQUAÇÃO:

Aumentos no preço elevam a quantidade ofertada,

Aumentos na taxa de aluguel dos porcos farejadores, que é um aumento no custo de produção, reduzem a oferta.

Ambas as variáveis têm estimativas de coeficiente estatisticamente significantes.

MES, MQ2E e VI – Exemplo – Interpretação da Oferta

pfpq 0009,10139,10328,20 Oferta

Page 69: GMM, VI , MQ2S e Equações Simultâneas

INTERPRETAÇÃO DA EQUAÇÃO:- O coeficiente do preço é negativo, indicando que à medida que o preço aumenta, a quantidade demandada de trufas cai, como predito pela lei de demanda. - O valor-p indica que a inclinação estimada da curva de demanda é significativamente diferente de zero. - Aumentos nos preços dos substitutos para as trufas elevam a

demanda por trufas, que é uma característica dos bens substitutos.

- Finalmente, o efeito da renda é positivo, indicando que as trufas são um bem normal.

MES, MQ2E e VI – Exemplo – Interpretação da oferta e demanda

dipspq 5016,02965,11234,1279,4 Demanda

Page 70: GMM, VI , MQ2S e Equações Simultâneas

ReferênciasBUENO, R. L (2008). Econometria de Séries Temporais. São Paulo: Cengage Learning. Cap 5.

CARTER HILL, R.; GRIFITHS, Willian E.; JUDGE, George G. (2003). Econometria. 2. ed. São Paulo: Saraiva.

DAVIDSON, R. e MACKINNON, J. G. (2004). Econometric Theory qand Methods. Oxford University Press, Oxford. Capítulo 9.

FLORES JR., Renato G. Notas de aula da Fundação Getulio Vargas – 2003.

FIGUEIREDO, Erik A. Notas de aula da Universidade Federal da Paraíba – 2010.

GREENE, William. H. Econometric Analysis. (ANO) (5th ed). Cap 18.

HAYASHI, F. (2000). Econometrics. New Jersey: Princeton University Press. Cap 3 e 4.

NEDER, Henrique D. Notas de aula da Universidade Federal de Uberlândia – 2009.

WOOLDRIDGE, J. M. (2003). Introductory Econometrics: A Modern Approach. 2nd ed. New York: Thomson Learning. 70

Page 71: GMM, VI , MQ2S e Equações Simultâneas

Evidências empíricas

• Issler & Piqueira. Estimating Relative Risk Aversion, the Discount Rate, and the Intertemporal Elasticity of Substitution in Consumption for Brazil using Three Types of Utility Function. BRE, 2000.

• Diversas funções de utilidade: função utilidade potência (CRRA), utilidade com hábito externo, e aversão ao desapontamento (Kreps-Porteus).

• Outras matrizes de ponderação: Newey e West (1987), com janela fixa, de Newey e West (1994) com janela variável, e o método de Andrews (1991), com janela variável;