模型为：log(Y)=a+blog(WY)+u，u 是正态分布，Y是已知的样本向量，W为空间邻接矩阵。

其中待估参数为a,b及u的方差西格玛方

寻求大神帮忙写代码估计模型参数，最好是R代码，谢谢啊！