求助关于空间面板计量

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

求问各位牛人，①如何将geoda生成的权重矩阵导入stata里面，或者如何将geoda生成的权重矩阵转换成【0，1】形式的？②如何在stata里面生成空间面板模型的权重矩阵？非常谢谢

扫码加我 拉你入群

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

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

分享0 收藏5 回帖

关键词：空间面板 空间面板模型 Stata 权重矩阵 tata 空间 模型 如何

看看能不能转化成excel格式，然后用xlsread（）命令读入？

zjying2000 发表于 2012-8-22 12:35
看看能不能转化成excel格式，然后用xlsread（）命令读入？

非常谢谢，可是用记事本读出来的矩阵是这个样式，不知道如何转变？
0 76 zhejiang BOUNT_
1625 2
1665 1649
1641 1
1647
1647 4
1734 1721 1680 1641
1649 4
1701 1680 1665 1625
1665 5
1756 1649 1625 1701 1737
1666 0

1668 0

1680 5
1737 1701 1647 1734 1649
1701 4
1737 1680 1649 1665
1721 2
1734 1647
1729 0

1734 5
1768 1737 1721 1647 1680
1737 8
1821 1770 1768 1756 1734 1680 1701 1665
1755 0

1756 5
1856 1841 1821 1737 1665
1768 7
1866 1821 1770 1826 1828 1734 1737
1770 3
1821 1768 1737
1786 1
1801
1801 4
1901 1826 1877 1786
1813 0

1821 6
1841 1768 1866 1756 1770 1737
1826 4
1901 1828 1801 1768
1828 5
1901 1866 1850 1826 1768
1841 6
1934 1918 1856 1866 1821 1756
1849 0

1850 1
1828
1856 6
2030 2007 1979 1918 1841 1756
1866 8
1964 1934 1901 1982 1828 1768 1841 1821
1877 1
1801
1901 6
1982 1962 1801 1826 1828 1866
1916 1
1972
1918 6
2020 2007 1934 1988 1841 1856
1934 5
1988 1866 1964 1841 1918
1943 1
1972
1962 4
2012 2006 1982 1901
1964 6
2053 2013 1988 1982 1866 1934
1972 3
2041 1943 1916
1979 2
2030 1856
1982 6
2053 2012 1901 1962 1964 1866
1988 6
2027 2020 2013 1964 1934 1918
2006 3
2071 2012 1962
2007 7
2083 2063 2030 2020 2098 1918 1856
2012 6
2076 2071 2053 2006 1982 1962
2013 6
2098 2061 2027 2020 1964 1988
2020 5
2013 2007 2098 1988 1918
2027 2
1988 2013
2030 4
2083 2007 1979 1856
2041 2
2058 1972
2053 5
2076 2061 2012 1982 1964
2058 2
2104 2041
2061 6
2135 2114 2098 2053 2076 2013
2063 1
2007
2071 4
2134 2076 2006 2012
2076 7
2153 2114 2071 2134 2053 2012 2061
2083 3
2098 2007 2030
2098 7
2175 2135 2061 2013 2083 2020 2007
2104 3
2141 2140 2058
2114 5
2199 2135 2153 2076 2061
2131 0

2134 4
2206 2153 2071 2076
2135 4
2175 2114 2098 2061
2140 2
2141 2104
2141 3
2178 2104 2140
2153 5
2232 2199 2134 2076 2114
2175 4
2251 2199 2135 2098
2178 1
2141
2199 6
2254 2251 2232 2153 2114 2175
2200 1
2230
2206 1
2134
2226 0

2230 2
2266 2200
2232 3
2254 2153 2199
2251 2
2199 2175
2254 2
2232 2199
2266 2
2309 2230
2309 1
2266

用R吧，splm包可以做，你这个权重文件是geoda生成的gal格式文件，splm是认可的，数据和权重文件准备好后，一个命令就行了。网上搜一下《splm:spatial panel data models in R》,仔细看看说明。

淘宝网橙迷橙橙 发表于 2012-8-24 20:42
用R吧，splm包可以做，你这个权重文件是geoda生成的gal格式文件，splm是认可的，数据和权重文件准备好后，一 ...

非常谢谢

淘宝网橙迷橙橙 发表于 2012-8-24 20:42
用R吧，splm包可以做，你这个权重文件是geoda生成的gal格式文件，splm是认可的，数据和权重文件准备好后，一 ...

你好，我现在用geoda生成的gal文件有了，然后stata的相关数据也有，求问gal文件如何导入到R，相关的ID如何设定（我尝试w<- read.gal("E:/splm/zj.gal") ，但是显示read.gal("E:/splm/zj.gal") : GAL file IDs and region.id differ）？有没有相关的程序？非常谢谢。
PS.这是我草拟的程序，但是进行到引入权重时就进行不下去了。
library(foreign)
mydata<-read.dta("E:/splm/merge.dta")
library(sp)
library(boot)
library(lattice)
library(MASS)
library(nlme)
library(Matrix)
library(deldir)
library(coda)
library(maptools)
library(spdep)
library(bdsmatrix)
library(Formula)
library(zoo)
library(sandwich)
library(plm)
library(spam)
library(ibdreg)
library(nnet)
library(car)
library(lmtest)
library(Ecdat)
library(splm)
###引入权重###
library(spdep)
w<- read.gal("E:/splm/zj.gal")
*****计算***
fm <-rincome ~ ishare  
reg1<- spml(fm, data =mydata, listw = mat2listw(w), model="within", spatial.error="b")
summary(reg1)
reg2<- spml(fm, data =mydata , listw = mat2listw(w), model="random", spatial.error="none", lag=TRUE)
summary(reg2)

wujunqian 发表于 2012-8-28 16:33
你好，我现在用geoda生成的gal文件有了，然后stata的相关数据也有，求问gal文件如何导入到R，相关的ID如何 ...

去看Anselin的《Spatial Regression Analysis in R: A  Workbook》以及《splm-Spatial Panel Data Models in R》

gal文件读入用read.gal函数。如果不懂看一下相关的帮助文件。
数据读入用read.csv函数。要保证数据是用逗号分隔的文本文件。
以上两种文件最好用geoda统一生成，以保证空间权重文件与数据中观测的一致性。

淘宝网橙迷橙橙 发表于 2012-8-28 20:19
gal文件读入用read.gal函数。如果不懂看一下相关的帮助文件。
数据读入用read.csv函数。要保证数据是用逗号 ...

非常谢谢你，我已经尝试做出来了。
library(spdep)
w<- read.gal("E:/splm/zj.GAL",override.id=TRUE)
ww<-nb2listw(w,zero.policy=TRUE)

wujunqian 发表于 2012-8-30 16:17
非常谢谢你，我已经尝试做出来了。
library(spdep)
w

你好，请问一下你最后一个语句是什么意思，我按照你的方法尝试怎么没成功呢？