multilevLCA 2.1.0运行报错pinv(): svd failed 求解决方案

问题背景 #

升级multilevLCA包至2.1.0版本后，运行原有代码或官方示例代码均抛出Error: pinv(): svd failed错误，降级回2.0.1版本、切换至Posit Cloud/Google Colab环境均无法解决该问题。

复现代码 #

library(multilevLCA)
data("dataIEA")
data = dataIEA
Y = colnames(dataIEA)[5:16]
iT = 3
out = multiLCA(data, Y, iT)

错误输出 #

Discarding rows with missing values for any indicator...

Fitting measurement model...

Error: pinv(): svd failed

临时解决方法 #

• 检查数据矩阵奇异性
  SVD失败多因矩阵秩不足（不可逆），可先验证指标矩阵的秩：

  indicator_data <- na.omit(dataIEA[, Y])
  cat("指标矩阵秩：", qr(indicator_data)$rank, "\n")
  cat("指标矩阵列数：", ncol(indicator_data), "\n")
  

  若秩小于列数，说明存在多重共线性，需移除高度相关的指标后重新拟合模型。

• 调整分类数量
  尝试修改iT参数（分类数），比如改为2或4，可能绕过SVD计算的异常：

  # 尝试2分类
  out <- multiLCA(data, Y, iT = 2)
  # 尝试4分类
  out <- multiLCA(data, Y, iT = 4)
  

• 降级R版本
  当前使用的R 4.5.1为较新版本，可能与multilevLCA存在兼容性问题，可尝试降级至R 4.4.x版本，再搭配multilevLCA 2.0.1使用。

• 反馈给包开发者
  官方示例代码也出现该错误，属于包版本bug，建议直接向包开发者提交问题，附上会话信息、错误日志及复现代码，以便快速修复。

本地会话信息 #

R version 4.5.1 (2025-06-13)
Platform: aarch64-apple-darwin20
Running under: macOS Sequoia 15.6.1

内容的提问来源于stack exchange，提问作者Robin Schädler