Я получаю сообщение об ошибке ниже:
Ошибка в solve.default (crossprod (attr (x, "отклонение)) [,! 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)
Я знаю, что ошибка может быть связана с сингулярной матрицей Гессе, но я хочу проверить, если онаприведет к решению, если я установлю аргумент tol=1e-33
в функции solve
.К сожалению, mlogit
, хотя я считаю, что solve
существует в пакете maxLik
, который интегрирован в пакет mlogit, он не распознает tol
.Как мне поступить, чтобы функция mlogit
распознала tol
?
Спасибо