标签： r matrix spatial
我正在通过R中的Gibbs采样器估计二元空间自回归模型。当我使用邻域矩阵执行此操作时，一切都非常好，但是当我切换到距离矩阵时，程序几乎立即停止，并输出臭名昭著：
＆＃34; ginv（A）中的错误：＆＃39; X＆＃39;必须是数字或复杂矩阵＆＃34;
令我感到奇怪的是，当我尝试自己计算广义逆时，一切都顺利进行，并且所有矩阵似乎都是数字的。您可以在下面找到可复制示例的数据和代码。
Data Link
代码（执行最后一行后发生错误，这也是通过bsar函数估算的玩具模型）：
Code Link
有什么想法吗？