Momento R do Dia – Reproduzindo um exemplo de Lettau & Ludvigson (2004)

Untitled 171 Uma das grandes qualidades dos bons pesquisadores é quando eles disponibilizam os resultados dos artigos e as bases de dados para que qualquer um possa reproduzir seu trabalho. É o caso do prof. Lettau que tem um artigo importante em Macroeconomia e Finanças que foi atualizado em 2004 por conta de diversas mudanças nos dados utilizados na primeira versão do artigo (a primeira versão é de 2001).

Pois bem, o prof. Lettau criou uma página para disponibilizar o artigo e os dados. Ela está aqui. O artigo é este: Understanding Trend and Cycle in Asset Values: Reevaluating the Wealth Effect on Consumption. O começo do artigo é um exemplo de boa didática e retórica, o que me faz reproduzi-lo aqui, só para a degustação do leitor.

Both textbook economics and common sense teach us that the value of household wealth should be related to consumer spending. Early academic work by Franco Modigliani (1971) suggested that a dollar increase in wealth (holding Ž fixed labor income) leads to an increase in consumer spending of about Žfive cents. Since then, the so-called “wealth effect” on consumption has increasingly crept into both mainstream and policy discussions of the macroeconomy.

Today, it is commonly presumed that signiŽficant movements in wealth will be associated with movements in consumer spending, either contemporaneously or subsequently. Quantitative estimates of roughly the magnitude reported by Modigliani are routinely cited in leading macroeconomic textbooks, and are important features of many contemporary macroeconomic models, including those still widely studied by both academic economists and practitioners.

In this paper, we reevaluate the empirical foundation for such estimates of the consumption-wealth link. Contrary to conventional wisdom, we Žfind that a surprisingly small fraction of the variation in household net worth is related to variation in aggregate consumer spending.

Legal, não? Tema situado na literatura com a citação do resultado clássico e essencial de Modigliani e mais uma fisgada no leitor com um resultado contra-intuitivo. Quando eu crescer, quero começar a escrever meus artigos assim (mas tenho que encontrar resultados legais para tal…).

Pois é. Duas coisas legais. No apêndice online dos dados, Lettau disponibiliza o vetor de cointegração estimado:

 cay=cc-0.621y -0.269a (demeaned)

Isso nos permite reproduzir seus resultado para, por exemplo, checar a qualidade do algoritmo de um programa específico. São três variáveis aí: consumo (cc), renda (y) e ativos (a). “cay” é o nome do termo de correção de erros (o famoso mecanismo de correção de erros).

A estimação, seguindo a leitura do artigo, retorna o seguinte vetor:

cay_R = cc -0.624y -0.263a .

A precisão é razoável, não? Até a segunda casa decimal e com os sinais corretos.

Eis o script.

# reprodução dos resultados de Lettau & Ludvigson (2004)

lettau<-read.table("C:/Users/cdshi_000/Documents/Meus Documentos/Meus Documentos/cursos/Econometria II/lettau_data.csv", 
                   header=TRUE, sep=",", na.strings="NA", dec=".", strip.white=TRUE)
head(lettau)
tail(lettau)

# http://faculty.haas.berkeley.edu/lettau/data/cay_pce_14Q3.txt
# Leia as observações sobre os dados

# transformando para séries de tempo.

lettau.ts<-ts(lettau, start=c(1951,4), freq=4)

head(lettau.ts)

# não precisamos da primeira coluna, nem da última

lettau.ts<-lettau.ts[,c(2,3,4)]

# conferindo...

head(lettau.ts)

# uma rápida olhadela nos dados...

plot(lettau.ts)

# para o exercício, vamos seguir a especificação do autor. Usamos o 
# pacote vars para tentar reproduzir seus resultados.

library(vars)

t<-ca.jo(lettau.ts, type=c("eigen"), K=2, ecdet=c("none"),spec=c("transitory"))

# feito o teste do maior autovalor, aproveito para reproduzir o VECM do autor
# com uma relação de cointegração. Ela é o "teste$beta" abaixo.

teste<-cajorls(t, r=1)
teste$beta

Legal, não?

Anúncios

Deixe um comentário

Preencha os seus dados abaixo ou clique em um ícone para log in:

Logotipo do WordPress.com

Você está comentando utilizando sua conta WordPress.com. Sair / Alterar )

Imagem do Twitter

Você está comentando utilizando sua conta Twitter. Sair / Alterar )

Foto do Facebook

Você está comentando utilizando sua conta Facebook. Sair / Alterar )

Foto do Google+

Você está comentando utilizando sua conta Google+. Sair / Alterar )

Conectando a %s