Я планирую использовать пакет R, mlo git, для расчета модели смешанного эффекта lo git.
> library("mlogit")
> fut<-mlogit.data(after,choice = "choice",shape="long",alt.var = "mode",chid.var = "tripid")
> randomfuturemodel<-mlogit(choice ~ time0+price|as.factor(trippurpose)+as.factor(age)+gender+as.factor(educationdegree)
+as.factor(adultnumber)+as.factor(employedadults)+as.factor(vehicle)
+as.factor(homeincome)+as.factor(personalincome)+ Ometro+Dmetro+Dtransit+Dpark+Ddensity,halton = NA,
R = 100,print.level = 1,fut,reflevel = "transit")
summary(randomfuturemodel)
Результат может быть хорошо рассчитан. Но это представляет что-то неправильно, когда я пытаюсь вычислить случайные эффекты. Тип указанных выше переменных цифра c в R. Но в моей работе некоторые переменные являются категориальными переменными, такими как 0,1,2,3. Я нарисовал QQplot, чтобы представить распределение соответственно. Эти две переменные, время и цена, являются логарифмически нормальными распределениями.
> alt<-mlogit(choice ~ time0+price|as.factor(trippurpose)+as.factor(age)+gender+as.factor(educationdegree)
+as.factor(adultnumber)+as.factor(employedadults)+as.factor(vehicle)
+as.factor(homeincome)+as.factor(personalincome)+ Ometro+Dmetro+Dtransit+Dpark+Ddensity,halton = NA,
R = 100,print.level = 1,fut,reflevel = "transit",rpar = c(time0="ln",price="ln"))
Ошибка в solve.default (crossprod (attr (x, "отклонение)) [,! Исправлено ])): система вычислительно единственная: число взаимных условий = 3.39628e-22