Я пытаюсь запустить hmftest в пакете mlogit.В Интернете есть пример, но я действительно его не понимаю.
Скажем, у меня есть 3 категории: 1,2,3 и две непрерывные переменные: A, B. Мне просто нужно запустить тест Хаусмана-Макфаддена, чтобы проверить, верно ли предположение о независимости в моей полиномиальной логистической регрессии,Я не понимаю пример главным образом потому, что я не знаю природу переменных в этом примере.
## from Greene's Econometric Analysis p. 731
data("TravelMode",package="AER")
TravelMode <- mlogit.data(TravelMode,choice="choice",shape="long",
alt.var="mode",chid.var="individual")
## Create a variable of income only for the air mode
TravelMode$avinc <- with(TravelMode,(mode=='air')*income)
## Estimate the model on all alternatives, with car as the base level
## like in Greene's book.
#x <- mlogit(choice~wait+gcost+avinc,TravelMode,reflevel="car")
x <- mlogit(choice~wait+gcost+avinc,TravelMode)
## Estimate the same model for ground modes only (the variable avinc
## must be dropped because it is 0 for every observation
g <- mlogit(choice~wait+gcost,TravelMode,reflevel="car",
alt.subset=c("car","bus","train"))
## Compute the test
hmftest(x,g)
Может ли кто-нибудь, кто понимает приведенный выше пример, оказать некоторую помощь?Большое спасибо,