Ho una domanda riguardante il pacchetto mnlogit in R che chiederò su StackOverflow poiché è correlato a una lingua e una libreria specifiche, tuttavia non sarò offeso se qualcuno decide di spostarlo su Cross Validated (era una scelta difficile di quale sito di StackExchange era più appropriato).Multinomial Logit Choice Modello in R con mnlogit()
Sto solo cercando di capire come funziona, dal momento che la documentazione non è troppo dettagliata, vedo alcune voci di blog confuse/in conflitto su questo pacchetto e perché sono trascorsi un paio di anni da quando ho affrontato i modelli di scelta.
L'esempio di utilizzo di mnlogit è fondamentalmente come segue:
> require(mnlogit)
> data(Fish)
> head(Fish)
mode income alt price catch chid
1.beach FALSE 7083.332 beach 157.930 0.0678 1
1.boat FALSE 7083.332 boat 157.930 0.2601 1
1.charter TRUE 7083.332 charter 182.930 0.5391 1
1.pier FALSE 7083.332 pier 157.930 0.0503 1
2.beach FALSE 1250.000 beach 15.114 0.1049 2
2.boat FALSE 1250.000 boat 10.534 0.1574 2
> fm <- formula(mode ~ price | income | catch)
> result <- mnlogit(fm, Fish, "alt", ncores = 2)
Quello che trovo confuso è mode e alt. Avrei pensato che la variabile dipendente sarebbe stata la scelta multinomiale, che sembra essere alt (spiaggia, barca, charter o molo).
Invece è questa variabile logica, mode. Cos'è mode?
Per provare a capire meglio ho letto uno R-Bloggers article sull'argomento. In realtà reso le cose più confuse ripetendo l'esempio Fish, ma con un preambolo che descrive il set di dati in questo modo:
Un frame di dati contenente:
mode - The choice set: beach, pier, boat, and charter price - price for a mode for an individual catch - fish catch rate for a mode for an individual income - monthly income of the individual decision-maker chid - decision maker ID
Descrive mode come se fosse alt quindi non menziona alt. Qualcuno può spiegarmelo?
Se mode davvero ci volle sui valori di scelta discreta, come descritto nell'articolo R-blogger, allora il modello formula(mode ~ price | income | catch) avrebbe senso per me ... ma mode è una variabile logica, quindi sono completamente confuso.