执行误差项空间自相关LM检验时遇“非一致参数”错误求助

问题背景 #

• 面板数据为pdata.frame格式，涵盖RN地区166个市镇4年观测数据，共664条记录
• 基于RN市镇SHP文件构建空间权重矩阵：

  nblist <- poly2nb(rn, queen = TRUE)
  w.rn.queen <- spdep::nb2listw(nblist, style="W")
  

• 执行空间滞后模型的LM检验时触发错误：

  splm::slmtest(within_model, data = rn_df, listw = w.rn.queen, test='lme')
  

  错误提示：Error in X %*% coef(mod) : non-conformable arguments

排查与解决建议 #

1. 核对维度一致性 #

• 确认权重矩阵w.rn.queen对应的空间单元数（166个市镇），与within_model拟合时使用的个体数量完全匹配。若rn_df存在缺失市镇、重复个体，会导致模型矩阵与权重矩阵维度不兼容。
• 用以下代码快速验证：

  # 权重矩阵对应的空间单元数
  length(w.rn.queen$neighbours)
  # 模型中个体的唯一数量
  length(unique(index(within_model)[[1]]))
  

2. 确保模型与数据匹配 #

• within_model必须是基于rn_df拟合的固定效应模型，且个体标识（如市镇ID）的顺序与权重矩阵的空间要素顺序完全一致。若shp文件的市镇顺序和rn_df的个体顺序不同，会导致空间滞后计算时维度不匹配。
• 可将shp文件的ID与rn_df的个体ID做匹配，重新排序后再构建权重矩阵。

3. 检查模型参数提取 #

• 固定效应模型（within模型）会吸收个体固定效应，coef(within_model)仅返回自变量的系数。需确保模型矩阵X的列数与系数向量的长度一致，避免因模型设定错误（如包含未识别的变量）导致维度不匹配。
• 用model.matrix(within_model)查看模型矩阵的结构，确认其列数等于length(coef(within_model))。

4. 重新构建权重矩阵 #

• 验证poly2nb生成的邻接列表是否正确：检查nblist中每个市镇的邻接数量是否合理，是否存在孤立单元（可能导致标准化后的权重矩阵异常）。
• 重新生成权重矩阵时，可添加检查步骤：

  # 检查邻接列表是否有缺失
  table(card(nblist))
  # 转换为矩阵确认维度
  w_mat <- nb2mat(nblist, style="W")
  dim(w_mat) # 应返回166 166
  

内容的提问来源于stack exchange，提问作者Alex Gois