R语言做空间杜宾模型代码

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

      用R自带数据集跑了一篇，关键是生成自变量的空间滞后项WX，此处用到了slag()。之后用自己的数据也跑了一次，成功了。如遇问题欢迎交流。

>data(Produc)
#Produc为R中自带数据集
>usa.lw<-mat2listw(usaww)
>Produc.pd<-pdata.frame(Produc, index=c("state", "year"))  
#对原数据格式进行一定的转换
>View(Produc.pd)
>Produc.pd$Wpcap<-slag(Produc.pd$pcap, usa.lw)
#生成自变量的滞后量即WX项，会发现数据集中多了一列Wpcap,下面的自变量作同样的处理
>Produc.pd$Wpc<-slag(Produc.pd$pc, usa.lw)
>Produc.pd$Wunemp<-slag(Produc.pd$unemp, usa.lw)
>View(Produc.pd)
>fm <- log(gsp)~log(pcap)+log(pc)+log(emp)+unemp+Wpcap+Wpc+Wunemp
>SPDM=spml(fm,data=Produc.pd,listw=usa.lw,model="within",effect = "individual",lag=TRUE,spatial.error = "none")
>summary(SPDM)
>effects.splm(SPDM)

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏14 回帖

关键词：空间杜宾模型 杜宾模型 R语言 Individual individua R语言做杜宾模型

当然，在这之前还要加载相关的空间分析包，如splm等。

注意：这是空间面板杜宾模型的代码。

>library(splm)
>library(plm)
#调用面板数据包对面板数据进行格式转换
>data(Produc)
#Produc为R中自带数据集
>usa.lw<-mat2listw(usaww)
>Produc.pd<-pdata.frame(Produc, index=c("state", "year"))  
#对原数据格式进行一定的转换
>View(Produc.pd)
>Produc.pdWpcap<−slag(Produc.pdpcap, usa.lw)
#生成自变量的滞后量即WX项，会发现数据集中多了一列Wpcap,下面的自变量作同样的处理
>Produc.pdWpc<−slag(Produc.pdpc, usa.lw)
>Produc.pdWunemp<−slag(Produc.pdunemp, usa.lw)
>View(Produc.pd)
>fm <- log(gsp)~log(pcap)+log(pc)+log(emp)+unemp+Wpcap+Wpc+Wunemp
>SPDM=spml(fm,data=Produc.pd,listw=usa.lw,model="within",effect = "individual",lag=TRUE,spatial.error = "none")
>summary(SPDM)
>effects.splm(SPDM)

wjj2111987 发表于 2019-9-2 20:47
>library(splm)
>library(plm)
#调用面板数据包对面板数据进行格式转换

以后面这个为准，一定要记得调一次plm包

谢谢大虾分享！

wjj2111987 发表于 2019-9-2 17:35
用R自带数据集跑了一篇，关键是生成自变量的空间滞后项WX，此处用到了slag()。之后用自己的数据也跑了 ...

学习了，感谢分享！

楼主，你好，固定效应模型有偏误的修正，不知你用R软件如何实现的。谢谢回答！！

您好，感谢您提出的方法，请问要怎么将总效应分解为直接效应和简介效应呢？

LeeHyeYoung 发表于 2020-1-28 22:33
您好，感谢您提出的方法，请问要怎么将总效应分解为直接效应和简介效应呢？

这个当时用R也没做出来，又转用STATA做了效应分解。发现R中有些功能还是不全，也可能是我水平有限吧。