采用plm包,固定、随机效应模型,Hausman过度识别检验(原假设是两个模型一致);
加入截面变系数,即是按个体分别进行回归分析。
> #载入面板数据分析包
library("plm", lib.loc="D:/R/R/R-3.1.2/library")
p_csy = plm.data(csy_zero, indexes = c("industry", "year"))
#fix_e = plm(data = p_csy, export ~ human, model = "within")
#summary(fix_e)
#取对数
p_csy$export = log(p_csy$export)
p_csy$human = log(p_csy$human)