stata空间杜宾,基于内生时空权重矩阵

use D:\w.dta,clear
spmat dta w var*,replace
spmat summarize w, links
spatwmat using D:\w.dta,name(w) standardize
use D:\7.5回归.dta
destring *, replace
xtset id year
xsmle d mmx_dra mmx_pergdp mmx_ava mmx_ale mmx_aad mmx_covid mmx_pattent mmx_farnumdra mmx_coopfarnum agrlevel ass main, model(sdm) wmat(w) type(both) nolog effects fe

这是我进行空间杜宾的命令,3年时间,450个地区,共1350个观测值,用的是反距离矩阵w.dat,现在把距离矩阵改为行列数都是1350的时空权重矩阵,请问如何修改命令。

在使用 spatwmat 命令时,将 using 选项改为你新准备的时空权重矩阵文件,例如:

spatwmat using D:\w_stmat.stmat, name(w) standardize

确保你的时空权重矩阵文件中包含了正确的时空权重信息。

xsmle 命令中,将 wmat 选项改为你新准备的时空权重矩阵文件,同时使用 type(both) 来指定同时考虑空间和时间权重,例如:

xsmle d mmx_ dra mmx_ _pergdp mmx_ _ava mmx_ ale mmx_ _aad mmx_ covid mmx_ pattent mmx_ farnumdra mmx_ coopfarnum agrlevel assmain, model(sdm) wmat(w) type(both) nolog effects fe

要将距离矩阵改为行列数都是1350的时空权重矩阵,您可以按照以下步骤修改命令:

1.使用spmat命令将权重矩阵加载为Stata的稀疏矩阵格式:
scss
Copy code
spmat dta w, from(D:\w.dta) replace
2.使用spmat命令获取权重矩阵的汇总统计信息:
Copy code
spmat summarize w, links
3.使用spatwmat命令加载时空权重矩阵:
arduino
Copy code
spatwmat using D:\w.dta, name(w) standardize
4.使用use命令加载数据文件"D:\7.5回归.dta":
mathematica
Copy code
use D:\7.5回归.dta
5.使用destring命令将变量转换为数值格式(如果需要):
Copy code
destring *, replace
6.使用xtset命令设置面板数据的索引变量:
bash
Copy code
xtset id year
7.使用xsmle命令进行时空面板数据分析,其中将模型类型设置为sdm(时空延迟模型),并将权重矩阵指定为时空权重矩阵wmat(w):
scss
Copy code
xsmle d mmx_dra mmx_pergdp mmx_ava mmx_ale mmx_aad mmx_covid mmx_pattent mmx_farnumdra mmx_coopfarnum agrlevel ass main, model(sdm) wmat(w) type(both) nolog effects fe
请注意,这只是对您给出的命令进行修改以适应时空权重矩阵的说明。确保您已经准备好正确的权重矩阵文件,并根据您的分析需求适当调整其他参数和选项。

以下答案参考newbing,回答由博主波罗歌编写:
要将距离矩阵改为行列数都是1350的时空权重矩阵,你需要进行以下步骤:

  1. 将这个时空权重矩阵保存为一个Stata格式的文件(.dta文件),确保它包含以下变量:

    • 权重变量(例如,w):表示每个观测值之间的空间权重。
    • ID变量:表示每个观测值的唯一标识符。
    • Year变量:表示每个观测值的时间。

    假设你已经将这个时空权重矩阵保存为了"wmat.dta"。

  2. 修改命令中的spmatuse部分,将距离矩阵改为时空权重矩阵。命令如下:

spmat dta w var*, replace
summarize w, link
use D:\wmat.dta, clear
  1. 确保你的数据集中包含以下变量(在你的命令中有列出):
    • ID变量(id):用于标识每个观测值。
    • Year变量(year):用于表示每个观测值的时间。
    • 要建模的自变量(d、mmx_dra、mmx_pergdp等)。
  2. 运行修改后的命令来进行空间杜宾回归分析。命令如下:
string *, replace
xtset id year
xsmle d mmx_dra mmx_pergdp mmx_ava mmx_ale mmx_aad mmx_covid mmx_pattent mmx_farnumdra mmx_coopfarnum agrlevel ass main, model(sdm) wmat(type(both)) nolog effects fe

以上代码应该可以满足将距离矩阵改为行列数都是1350的时空权重矩阵并进行杜宾回归分析的需求。请注意,代码中的变量名和文件路径可能需要根据你的实际情况进行修改。同时,请确保你的数据集和时空权重矩阵正确加载。
如果我的回答解决了您的问题,请采纳!

参考 https://praschky.github.io/assets/Lecture3_SpatialStata.pdf
错误参考 https://bbs.pinggu.org/thread-6396921-1-1.html

创建一个名为"w_temporal"的新变量,该变量将作为时空权重矩阵。使用sqrt操作将矩阵w中的元素转换为时空权重矩阵w_temporal中的元素。在进行杜宾模型的命令中,将wmat(w)替换为wmat(w_temporal)。

https://download.csdn.net/download/xiubo1288/10229189

stata实现经济生态的空间杜宾模型
非常详细,可以借鉴下