//打开自己要用的矩阵数据(这个矩阵是截面矩阵,就是指:n个id,权重矩阵是n*n)
use "D:\stata16\实证\jingjijuzhen01.dta"
//扩展权重矩阵{是指把n*n的矩阵变成【(t*n)*(t*n)】矩阵,t是指一共是t年的数据}。V1-V31是指这个权重矩阵的变量名是从v1-v31,标红的根据自己的矩阵进行对应修改。
spcs2xt V1-V31,matrix(W) time(t)
//导入新的权重矩阵(是扩展了的)
spatwmat using Wxt,name(W) standardize
//导入待用面板数据
use 路径\dta.
//面板数据设置
xtset id year
//对所有变量回归
reg y x1 x2 x3
//LM检验
spatdiag,weights(W)
• 贴个结果图。欢迎大家多多评论~