全部版块 我的主页
论坛 经济学论坛 三区 区域经济学
5172 9
2012-08-22
求问各位牛人,①如何将geoda生成的权重矩阵导入stata里面,或者如何将geoda生成的权重矩阵转换成【0,1】形式的?②如何在stata里面生成空间面板模型的权重矩阵?非常谢谢
二维码

扫码加我 拉你入群

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

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

全部回复
2012-8-22 12:35:43
看看能不能转化成excel格式,然后用xlsread()命令读入?
二维码

扫码加我 拉你入群

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

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

2012-8-22 15:32:21
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
二维码

扫码加我 拉你入群

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

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

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

扫码加我 拉你入群

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

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

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

扫码加我 拉你入群

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

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

2012-8-28 16:33:03
淘宝网橙迷橙橙 发表于 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)
二维码

扫码加我 拉你入群

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

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

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

说点什么

分享

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