R调用rankFD运行非参数方差分析时报错svd(X)存在无穷或缺失值

你遇到的svd(X)报错并非由数据中的显性NA/Inf值导致，核心原因是五因素全交互模型的设计矩阵秩亏，无法完成奇异值分解，具体排查和解决方法如下：

• 首先排查因子组合空单元格：全交互项Jobs*Stages*ShopCondition*PD*CV会生成五阶交叉的所有水平组合，只要存在任意一组水平组合没有对应观测值，就会导致设计矩阵奇异，触发该报错。你可以运行代码ftable(Jobs, Stages, ShopCondition, PD, CV, data = Datos)统计所有交叉组合的观测频数，确认是否存在频数为0的空单元格。
• 调整模型复杂度：非参数秩方差分析对高维全交互模型的支持度较差，无需一开始就拟合五阶全交互。可以先尝试拟合主效应模型验证函数运行逻辑：

  # 先拟合仅包含主效应的模型
  rankFD(Tardiness~Jobs + Stages + ShopCondition + PD + CV, data = Datos)
  

  若主效应模型可正常运行，再逐步加入二阶、三阶交互项，定位到导致矩阵奇异的交互层级后再做针对性调整。
• 统一自变量数据类型：若你把分类自变量作为数值型变量传入模型，多变量交互时也容易出现共线问题。运行下方代码将所有分类自变量转换为因子类型后再重试：

  Datos[,c("Jobs","Stages","ShopCondition","PD","CV")] <- lapply(Datos[,c("Jobs","Stages","ShopCondition","PD","CV")], as.factor)
  

• 若确实需要保留全交互结构，可以通过合并因子水平、使用基于 Wald 统计量的非参数检验替代方案，或者补充空单元格对应观测的方式解决秩亏问题。

内容的提问来源于stack exchange，提问作者Diego Rodas