Я пытаюсь понять модель логита в R. Я использовал набор данных в R, вызовите TravelMode в пакете AER.
Цель состоит в том, чтобы найти оптимальный режим для каждого человека.
в наборе данных 4 режима: воздух, поезд, автобус, машина
и 6 переменных: ожидание, стоимость, путешествия, стоимость, доход, размер
Я запустил Mlogit и получил следующее лето:
Call:
mlogit(formula = choice ~ wait | vcost | travel | gcost | income |
size, data = TM, method = "nr")
Frequencies of alternatives:
air train bus car
0.27619 0.30000 0.14286 0.28095
nr method
5 iterations, 0h:0m:1s
g'(-H)^-1g = 0.000153
successive function values within tolerance limits
Coefficients :
Estimate Std. Error z-value Pr(>|z|)
train:(intercept) -0.3618062 0.5295741 -0.6832 0.49448
bus:(intercept) -1.3173150 0.8146869 -1.6170 0.10589
car:(intercept) -4.3062975 0.7187265 -5.9916 2.078e-09 ***
wait -0.0944647 0.0103010 -9.1705 < 2.2e-16 ***
train:vcost -0.0375975 0.0083331 -4.5118 6.428e-06 ***
bus:vcost -0.0412963 0.0212837 -1.9403 0.05235 .
car:vcost -0.0642957 0.0162907 -3.9468 7.921e-05 ***
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Log-Likelihood: -189.77
McFadden R^2: 0.33124
Likelihood ratio test : chisq = 187.98 (p.value = < 2.22e-16)
Теперь у меня есть следующие вопросы:
1- почему нет воздушного перехвата?
2- почему нет оценочной стоимости для командировки, гкоста, дохода и размера?
3 - почему есть 3 vcost для каждого режима? Я думал, что в утилите мы должны найти вектор B (U = Bx + e) и B одинаков для всех режимов (кроме перехвата)