我收到以下错误:
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
?
谢谢