2.用老师分享的love 2015程序包分析 1)直接用原变量去做pvar
程序:xtset region year
pvarsoc lnpgdp urb_1 urb_2, maxlag(3) pvaropts(instl(1/4))*选最优滞后为滞后1期
pvar lnpgdp urb_1 urb_2, instl(1/4) gmmstyle
pvarstable, graph
pvarfevd, mc(200) save("fevd.dta")
pvarirf, mc(200) oirf byopt(yrescale)
由于模型不稳定,脉冲发散 遂放弃 2)加入 Ryan Decker的helm包对数据先进行helmert变换
脉冲同样发散 3)一阶差分处理后再helmert变化
gen Dlnpgdp = lnpgdp- L.lnpgdp
gen Durb_1 =urb_1-L.urb_1
gen Durb_2 =urb_2-L.urb_2
rename region id
helm id year
之后程序同上
可能我的时间序列太短 又经过多次差分,所以出现老师说的shocks die very quickly,
并且进行格兰杰因果分析也几乎没有因果关系
我的疑问:
1.我认为数据平稳不存在单位根的认识是错误的吗?毕竟后来用原数据做脉冲都是发散的
2.love博士给的例子没有进行helm变化,fd指令是与之相关的,然后解释说是默认会进行差分(fodspecifies that the panel fixed-effects be removed using forward orthogonaldeviation or Helmert transformation. By default, the first # lags of transformed depvarlist in the model are instrumented by the same lags in level),这个我不是很理解
3.shocks die very quickly 该怎么处理呢?增加数据量吗?
4.变量之间大多不存在格兰杰因果,这与现实不太相符,是我的模型不能做格兰杰因果吗?应该协整检验,建立ECM吗?