package mlogit R: À propos de Logit imbriqué

Question

Dans les modèles logit imbriqués, vous pouvez définir des régresseurs à chaque niveau de l'arborescence. Dans tous les exemples que j'ai lus dans le manuel et dans d'autres exemples, les régresseurs sont définis uniquement pour le dernier niveau. Je vais utiliser un exemple communément discuté, le mode de pêche.

Nids:

Shore = { Beach, Pier} , Boat = { Charter, Private} 

dire que j'ai régresseurs Price, CatchRate et Income. Comment puis-je utiliser Price et CatchRate pour expliquer le dernier niveau et Income pour expliquer le premier.

En R je peux faire:

mlogit(choice~price+catch,nests=list(shore=c("pier","beach"),boat=c("charter","private"))) 

mais je ne sais pas où coller la income variable.

Answer 1

L'exemple que vous utilisez et la réponse à votre question sont à la fois dans la documentation:

## model with charter as the reference level 

m <- mlogit(mode ~ price+ catch | income, data = Fish, reflevel = "charter") 

## same model with a subset of alternatives : charter, pier, beach 

m <- mlogit(mode ~ price+ catch | income, data = Fish, 
      alt.subset = c("charter", "pier", "beach")) 

## a pure "multinomial model" 

summary(mlogit(mode ~ 0 | income, data = Fish)) 

## which can also be estimated using multinom (package nnet) 

library("nnet") 
summary(multinom(mode ~ income, data = Fishing)) 

## a "mixed" model 

m <- mlogit(mode ~ price+ catch | income, data = Fish) 
summary(m)