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的时空权重矩阵,你需要进行以下步骤:
将这个时空权重矩阵保存为一个Stata格式的文件(.dta文件),确保它包含以下变量:
假设你已经将这个时空权重矩阵保存为了"wmat.dta"。
修改命令中的spmat
和use
部分,将距离矩阵改为时空权重矩阵。命令如下:
spmat dta w var*, replace
summarize w, link
use D:\wmat.dta, clear
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)。