为了方便使用,我把变量用大写绿色字母表示了。保姆级教程,每一步出现问题了怎么办也都有。
S1输入数据并清洗数据data cleaning
*将经费变量进行CPI换算,以2010年为基期,1978=100的CPI为调节因子
foreach v of X-Z {
gen cpi`v'=(`v'*536.1/cpi
}
变量X Y Z变为cpiX cpiY cpiZ
order id, before(year) //将地区变量移到时间之前这步骤是必须的
xtset id year //对面板数据中的个体变量和时间变量进行设定
xtdes //显示面板数据的结构
*S2描述性统计
xtsum cpiX-cpiZ //显示面板数据的统计特征
*S3单位根检验,三种检验方法
*单位根检验包括检验同质单位根的LLC 和Breitung,检验异质单位根的IPS、ADF-Fisher和PP-Fisher五种方法。
*LLC检验
//xtunitroot llc lX, trend demean lags(bic 12) #包含线性时间趋势项又包含个体固定效应项,是条件最苛刻的一种检验
//xtunitroot llc lX, demean lags(bic 12) #仅含个体固定效应项
//xtunitroot llc lX , noconstant demean lags(bic 12) #不包含线性时间趋势项,含个体固定效应项,是条件最宽松的一种检验
*IPS检验
//xtunitroot ips lX trend demean lags(bic 12) #包含线性时间趋势项又包含个体固定效应项,是条件最苛刻的一种检验
//xtunitroot ips lX demean lags(bic 12) #仅含个体固定效应项
//xtunitroot ips lX, noconstant demean lags(bic 12) #不包含线性时间趋势项,也不包含个体固定效应项,是条件最宽松的一种检验
*ADF-Fisher检验
//xtunitroot fisher lX ,dfuller lags(12) trend demean
//xtunitroot fisher lX,dfuller lags(12) demean
//xtunitroot fisher lX ,dfuller lags(12) noconstant demean
*结果全部平稳啦,进行下一步,我用的LLC,结果怎么看论坛里有
*原始数据如果不平稳做ln处理,重复S2-S4
foreach v of cpiX-cpiZ {
gen l`v'=ln(`v' )
}
变量cpiX cpiY cpiZ变为lcpiX lcpiY lcpiZ
*如果还不平稳,把不平稳的变量一阶差分,意义是变化量
gen dX=d.X
*我的就已经平稳了,如果你的还不平稳理论上可以做二阶差分,但是我不介意做。直接去做协整分析。
*S4 最优滞后阶数
pvar2 lcpiX lcpiY lcpiZ,lag(5) soc //看结果哪个星星多选哪个,我的结果是滞后4阶
*S5:格兰杰因果检验
pvar2 lcpiX lcpiY lcpiZ, lag(4) granger //滞后4阶
*S6:脉冲响应和方差分解(平稳的情况,不平稳协整)
pvar2 lcpiX lcpiY lcpiZ,lag(4) irf(10) //滞后一阶
pvar2lcpiX lcpiY lcpiZ,lag(4) decomp(30)
完结撒花,需要pvar2安装包的我也有哦,直接放进来了
安装方法,不管是windows还是mac都是直接把解压后的文件存放到第一个子目录下(我去,论坛已经有了,不能上传)。


雷达卡



[em23]

京公网安备 11010802022788号







