空间面板模型-空间分析.docx

文档介绍：空间面板模型

我们生活在时间和空间中，每个事件都在一定的时间和地点发生，因而可以标度出时间和空间坐标，这样的数据可以称为空间面板数据，它是指一定空间单元的时间序列观测。在研究实际问题时，空间面板数据本身具有更大自由度32




110101
201108








110101
201109








110101
201110








表8-2各变量的含义
变量名称
变量含义
单位
code
地理编码
无
date
年月信息
rate
发病率
popuDen
人口密度
千人/平方公里
PerGdp
当月人均国内生产总值
千元
aveTem
平均温度
摄氏度
precp
降水量
毫米
relHum
相对湿度
%
sunShn
日照强度
瓦/平方米
wndspd
风速
千米/时
此外，还需要的数据是包含该地区所有区县的地图文件(dbf文件和shp文件)。该文件为存放于光盘中的""",另外，该文件也可以通过ArcGis软件从全国各区县地图中选择生成。
(2)采用R语言建立空间面板模型第一步，加载如下程序包。代码如下：library(splm)
library(spdep)
library(sp)
library(maptools)
library(foreign)
第二步，导入所需的数据，代码如下：
hData<-("C:/")#导入发病率和影响因素的数据
dbf<-("C:/")#导入地图的数据(dbf格式)
第三步，将导入的两组数据合并，代码如下：
hData<-merge(hData,dbf,="code",="CNTY_CODE",=T)
第四步，生成新变量：人口密度(PopuDen)、人均GDP(PerGdp),代码如下：
hData<-cbind(hData,PopuDen=hData$population/hData$area*10A6,PerGdp=hData$gdp/hData$population)
第五步，导入shp格式的地图文件。由于dbf文件无法存储各地理单元的位置关系，导致无法生成空间权重矩阵，因此还需要导入shp格式的地图文件，代码如下：
<-readShapePoly("C:/")
第六步，生成空间权重矩阵。改步使用poly2nb函数实现，需要两个输入参数。第一个参数是导入的shp格式的地图数据名称。第二个参数是空间权重的确定方法(逻辑变量queen)