use w1.dta,clear
spatwmat using w1.dta,name(w) standardize
spmat dta w c1-c30, norm(row) replace //c1-c30是变量名
drop c1-c30
set matsize 330 //30*11
mat TMAT=I(11) //11年数据
mat Wxt=TMAT#w //30*30的截面矩阵与TMAT合并
svmat Wxt
save Wxt.dta,replace //保存LM检验所需空间面板权重矩阵
spatwmat using Wxt.dta,name(w2) standardize