Я оцениваю модель того, как люди путешествуют. Есть 4 варианта: самолет, машина, автобус и поезд. Я хочу создать вложенную логит-модель, в которой путешественники сначала выбирают, лететь или нет, а затем, если они не летают, выбирают между оставшимися вариантами. Таким образом, мои гнезда это (воздух) и (автомобиль, автобус, поезд).
Эта статья предполагает, что эта процедура должна быть возможной.
Однако, когда я приступаю к оценке этой модели, я обнаруживаю, что она сингулярна, если только каждое гнездо не содержит по крайней мере два элемента. Кажется, что в гнезде не может быть только один предмет, в отличие от статьи выше. Почему это?
Пример:
#install.packages("AER")
data("TravelMode", package = "AER")
# Works fine
nl1 <- mlogit(choice ~ gcost + wait, TravelMode, shape = 'long', alt.var = 'mode',
nests = list(public = c('train', 'air'), other = c('bus', 'car')))
# Gives singularity, as does any permutation with only one object in a nest
nl2 <- mlogit(choice ~ gcost + wait, TravelMode, shape = 'long', alt.var = 'mode',
nests = list(public = c('air'), other = c('bus', 'car','train')))
