全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
17273 21
2019-09-02
      用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)

二维码

扫码加我 拉你入群

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

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

全部回复
2019-9-2 17:49:10
当然,在这之前还要加载相关的空间分析包,如splm等。
二维码

扫码加我 拉你入群

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

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

2019-9-2 20:16:01
注意:这是空间面板杜宾模型的代码。
二维码

扫码加我 拉你入群

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

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

2019-9-2 20:47:06
>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)
二维码

扫码加我 拉你入群

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

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

2019-9-2 20:47:50
wjj2111987 发表于 2019-9-2 20:47
>library(splm)
>library(plm)
#调用面板数据包对面板数据进行格式转换
以后面这个为准,一定要记得调一次plm包
二维码

扫码加我 拉你入群

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

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

2019-9-3 17:38:32
谢谢大虾分享!
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群