paneldata<-read.csv('panel002.csv',header=T) #读入数据csv格式
> form<-inps~alr+car+fat+enf+ens+eny+emp+enn+lai+tei+cai+cur+cii+are+ssr+fae+gov+tar+ope+pro
> #1.判断随机效应or固定效应or混合效应
> m1<-plm(form,data=paneldata,index=c("code","year"),effect="individual",model="random")
summary(m1)
#1.判断随机效应or固定效应or混合效应
m2<-plm(form,data=paneldata,index=c("code","year"),effect="individual",model="within")
summary(m1)
m1<-plm(form,data=paneldata,index=c("code","year"),model="random")
m2<-plm(form,data=paneldata,index=c("code","year"),model="within")
m3<-plm(form,data=paneldata,index=c("code","year"),model="pooling")
summary(m1)
summary(m2)
summary(m3)
phtest(m1,m2) #hausman检验
|