r - 在 for 循环中使用 mlogit 给我一个错误

Question

for(i in 1:ncol(mldata))
{
  mlogit.mydata=mlogit(severity ~ 1|mldata[,i],data=mldata)
  mlogit[i]=summary(mlogit.mydata)$lratio$p.value
}

它给出以下错误：solve.default(H, g[!fixed]) 中的错误：Lapack 例程 dgesv：系统完全是奇异的：U[57,57] = 0

标签： rmlogit

Answer 1

请添加 mldata 详细信息。的第一个参数mlogit应该是 type formula。试试下面的代码

for(i in 1:ncol(mldata))
   {
      formula <- as.formula(paste0("y ~ 1 | ", colnames(mldata)[i])
      mlogit.mydata=mlogit(formula,data=mldata)
      mlogit[i]=summary(mlogit.mydata)$lratio$p.value
   }