在R mlogit中处理计算奇异系统

时间:2018-10-09 21:11:53

标签: r mlogit singular

我收到以下错误:

  

solve.default(crossprod(attr(x,“ gradi”)[,!fixed]))中的错误:     系统在计算上是奇异的:倒数条件数= 5.95531e-32

通过运行mlogit函数:

MIXL_no.socdemTR <- mlogit(formula = route~ cost + total_time.Egnatia + total_time.PEO |1 , data = dataTR_MIXL, rpar=c(cost="ln", total_time.Egnatia="ln", total_time.PEO="ln"), panel=T, R=200, halton=list("primes"= c(3,5,7), drop= rep(10,3)), method="bfgs", start = c(-1,0.1,0.5,0.2,0.1,0.3,0.2), print.init = T)

我知道该错误可能是由于奇异的Hessian矩阵引起的,但是如果要在tol=1e-33函数中设置参数solve,我想检查是否会导致解决方案。不幸的是,mlogit尽管我相信solve存在于集成在mlogit软件包中的maxLik软件包中,但是它无法识别tol。我应该怎么做才能让mlogit函数识别tol

谢谢

0 个答案:

没有答案