如何在R中实现对空间面板数据的LM检验

Sophiagxf7955 发表于 2015-10-29 10:13
楼主解决这个问题了吗？我现在也在做类似的内容。。。求告知具体程序。。谢谢。。

没有呢  最近落下了

楼主问题解决了吗？我也遇到同样的问题

elangshanren 发表于 2015-7-28 20:43
呵呵，你好好看看这篇文章里的介绍吧
http://www.jstatsoft.org/v47/i01/paper
要是还有问题可以再问 ...

这里面的LM检验也是针对模型做的估计呢，怎么对单个变量进行空间相关性检验呢？这么久了不知道大家都解决了问题没有~

请问楼主首页截图里面的各种检验是用lm.LMtests做出来的吗？还是做模型的时候用summary显示出来的？谢谢！

你说的东西貌似  在这里面有，
在最后一部分，关于检验这块的第一个检验，          是不是 备择假设是存在 误差相关，P值没有通过检验，说明存在空间误差相关性啊？
求讨论

蓝色 发表于 2015-6-28 11:17
随便执行一个，命令自带的例子：你可以看到各种检验的结果。
但这些命令式非官方命令，不知道准确性如何。 ...

问一下 spregxt y x1 x2,nc() wmfile() model(sar) mfx() stand tests 中的wmfile和mfx括号里分别表示是什么文件呀？初次用stata接触空间计量 都不是特别懂 看了帮助 没找到  

我也是被这个问题困扰好久了！终于找到了，slmtest这个函数，空间面板模型有很多种，用LM-error，LM-lag，RLM-error，RLM-lag这几个统计量可以确定应该用哪种模型，slmtest就是算的这几个统计量。楼主应该是想算的这个吧，我在R的函数介绍里看到的，应该没有错，大家康康，讨论一下~~

以下是R里的例子-------

data(Produc, package="plm")
data(usaww)
fm <- log(gsp)~log(pcap)+log(pc)+log(emp)+unemp
## robust LM test for spatial error sub spatial lag
## model on original data, pooling hypothesis
slmtest(fm, data=Produc, listw = usaww, test="rlme")
## model on within-transformed (time-demeaned) data,
## eliminates individual effects
slmtest(fm, data=Produc, listw = usaww, test="rlme",
  model="within")

找了两天终于找到了，开熏~~
另外截面数据用lm.lmtests
（我一开始搞混了希望没有误导人）

空气污染排放与经济增长的关系研究 ——基于中国省际面板数据的空间计量模型 艾小青 这篇文章中有R代码

最近也在用这两个包。spdep适合截面数据分析；splm是通过最大似然和glm估计系数，但是不能直接得出R方，还要其他命令。不当之处还请多多指正