R语言空间面板数据统计分析

install.packages(“splm”)
library(splm)
library(spdep)
library(sp)
library(Matrix)
library(sf)
library(spData)
install.packages(“Ecdat”)
data(Produc,package = “Ecdat”)
library(Ecdat)
library(Ecfun)
View(Produc)
data(usaww)
View(usaww)
GM=spgm(log(gsp)~log(pcap)+log(pc)+log(emp)+unemp,data=Produc,listw = usaww,moments = “fullweights”,spatial.error = TRUE) #空间面板数据的广义矩估计
summary(GM)
fm=log(gsp)~log(pcap)+log(pc)+log(emp)+unemp#空间面板极大似然估计
#fixed effects panel with spatial errors
Fespaterr=spml(fm,data=Produc,listw=mat2listw(usaww),model=“within”,spatial.error=“b”,hess=FALSE)
summary(Fespaterr)
#伴有随机效应和序列误差相关的空间面板模型的极大似然估计，命令speml
#random effects panel with spatial lag and serial error correlation
#optimization method set to “BFGS”
Sarsrmod<-spreml(fm,data=Produc,w=usaww,error=“sr”,lag=TRUE,method=“BFGS”)
summary(Sarsrmod)
#模型检验
bsjktest，Baltigi,Song,Jung,and Koh LM test for spatial panels
#伴有随机效应和序列误差相关的空间面板模型的极大似然估计，命令speml
#random effects panel with spatial lag and serial error correlation
#optimization method set to “BFGS”
Sarsrmod<-spreml(fm,data=Froduc,w=usaww,error=“sr”,lag=TRUE,method=“BFGS”)
summary(Sarsrmod)
#模型检验
#bsjktest，Baltigi,Song,Jung,and Koh LM test for spatial panels
bsjktest(fm,data=Produc,listw=usaww,test=“C.1”)
#bsktest，Baltigi,Song and Koh LM test for spatial panels
bsktest(fm,data=Produc,listw=mat2listw(usaww),test=“LM1”,standardize = TRUE)
#Covariance extractor method for splm objects
sarremod=spml(fm,data=Produc,listw=mat2listw(usaww),model = “random”,lag = TRUE,spatial.error = “none”)
library(lmtest)
library(zoo)
coeftest(sarremod)
install.packages(“car”)
library(carData)
library(car)
lht(sarremod,“log(pcap)=log(pc)”)

R语言空间面板数据统计分析相关推荐

R语言做面板VAR例子
面板VAR步骤: (1)对各变量做平稳性检验(IPS.PP.ADF.LLC等方法检验) 是逐个变量检验??还是一起检验?? (2)面板数据的最优滞后阶数确定(AIC和SIC方法) (3)在PVAR系统 ...
R语言空间数据处理(part1)--基础数据操作与处理
学习笔记,仅供参考学习书目:<R语言空间数据处理与分析实践教程>–卢宾宾; 基础数据操作与处理设置工作路径,并导入包 workL = "F:/MyStudio/Rstudio ...
R语言空间数据处理(part2)--空间数据读写
学习笔记,仅供参考学习书目:<R语言空间数据处理与分析实践教程>–卢宾宾; 准备工作设置工作路径,并导包 workL = "F:/MyStudio/Rstudio/RSpac ...
R语言将向量数据按照行方式转化为矩阵数据（设置参数byrow为TRUE）、计算矩阵数据的特征值（eigenvalue）
R语言将向量数据按照行方式转化为矩阵数据(设置参数byrow为TRUE).计算矩阵数据的特征值(eigenvalue) 目录 R语言将向量数据按照行方式转化为矩阵数据(设置参数byrow为TRUE). ...
R语言处理缺失数据的5个常用包
R语言处理缺失数据的5个常用包 1.常用缺失数据处理包 2. MICE 包 2.1基本介绍 2.2 实例展示 3.Amelia包 3.1基本介绍 3.2实例展示 4.missForest包 4.1基本 ...
r语言中进行数据可视化_R中的数据可视化
r语言中进行数据可视化 R programming was developed in 1993 for making graphs and producing statistical results. ...
R语言将向量数据按照行方式转化为矩阵数据（设置参数byrow为TRUE）
R语言将向量数据按照行方式转化为矩阵数据(设置参数byrow为TRUE) 目录 R语言将向量数据按照行方式转化为矩阵数据(设置参数byrow为TRUE) R语言是解决什么问题的? R语言将向量数据按照 ...
R语言将向量数据按照行方式转化为矩阵数据（设置参数byrow为TRUE）、对矩阵进行转置操作
R语言将向量数据按照行方式转化为矩阵数据(设置参数byrow为TRUE).对矩阵进行转置操作目录 R语言将向量数据按照行方式转化为矩阵数据(设置参数byrow为TRUE).对矩阵进行转置操作 R语言 ...
R语言可以直接从数据库读取数据么？R语言读取数据库数据示例
R语言可以直接从数据库读取数据么?R语言读取数据库数据示例目录 R语言可以直接从数据库读取数据么?R语言读取数据库数据示例 R语言是解决什么问题的? R语言可以直接从数据库读取数据么?R语言读取数据 ...

R语言空间面板数据统计分析

R语言空间面板数据统计分析相关推荐

最新文章

热门文章