不好意思,主要是有点着急了。
我刚开始照着用stata做地区发展能力的空间相关性, Moran I 值检验说明有数据空间相关性。文献指出此时要用极大似然法,我就按照找到的程序试着做了下,具体如下:
.program define jk
1.args lnf xb
2.quietly replace `lnf' = ln (normalden ($ML_y1,`xb')) %其实对这块不是很理解
3.end
.ml model lf jk ( r: lnjinjuan= )
. ml search
initial: log likelihood = -<inf> (could not be evaluated > )
feasible: log likelihood = -130.19333
improve: log likelihood = -130.19333
rescale: log likelihood = 21142.494
. ml maximize
initial: log likelihood = 21142.494
alternative: log likelihood = -130.19333
rescale: log likelihood = 21142.494
could not calculate numerical derivatives
discontinuous region with missing values encountered
结果如上。
谢谢斑竹!