McFadden R^2 및 Likelihood Ratio Test가 mLogit에 표시되지 않았습니다.

Question

제 문제는 이해하기가 쉽지만 솔루션을 찾을 수 없습니다. 나는 mLogit 라이브러리를 사용하여 R에 다음 코드를 사용하고

:

Call: 
mlogit(formula = Choice ~ price + distance + inveh + onoff + 
    prob | 0 , data = ExpData, method = "nr", print.level = 0) 

Frequencies of alternatives: 
    alt1 alt2 
0.51431 0.48569 

nr method 
4 iterations, 0h:0m:0s 
g''(-H)^-1g = 1.55E-07 
gradient close to zero 

Coefficients : 
      Estimate Std. Error t-value Pr(>|t|)  
price -7.3472e-01 3.1842e-02 -23.0735 < 2.2e-16 *** 
distance -5.8012e-04 6.6842e-05 -8.6790 < 2.2e-16 *** 
inveh -1.0994e-02 4.5466e-03 -2.4180 0.0156048 * 
onoff  1.1858e-01 3.4718e-02 3.4157 0.0006363 *** 
prob  5.6877e-01 8.2690e-02 6.8784 6.053e-12 *** 

--- 
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 

Log-Likelihood: -2912.3 

내가 맥파든과 가능성의 비율을 좀하고 싶습니다 :

library("mlogit") 
dat = read.csv("ExpeData.csv", header = TRUE) 
ExpData<- mlogit.data(dat,shape="wide", varying = 3:14, choice = "Choice",sep=".") 
wrf<- mlogit(Choice ~ price+distance+inveh+onoff+prob|0, ExpData) 
summary(wrf) 

내가 얻을 출력은 다음입니다 테스트

무엇이 잘못 되었습니까 ??

Answer 1

내 자신의 질문에 대답하기 위해, 내가 MC Fadden의 사용 R^2 테스트를 이해하지 않았다 아니었다. 내 문제는 R^2가 모델 요약에 표시되지 않았다는 것입니다.

내 R 버전은 2. *입니다. 최근에 내 컴퓨터를 업그레이드하고 내 문제를 해결 한 3.1.3 버전을 얻었습니다. 이제 모델 요약의 결과는 다음과 같습니다

Log-Likelihood: -7205.8 McFadden R^2: 0.067533 Likelihood ratio test : chisq = 1043.7 (p.value = < 2.22e-16)

나는 손으로 R^2을 추정 할 필요가 없습니다.

Answer 2

나는 McFadden에 대한 검색을 맡기고 있습니다. 게시하기 전에 그렇게 했어야합니다. LRT를하려면 관심이있는 공변량 모델과 모델이없는 작은 모델을 비교해야합니다. ? lrtest의 예를 수정 :

library("mlogit") 
data("TravelMode", package = "AER") 
ml <- mlogit(choice ~ wait + travel + vcost, TravelMode, 
      shape = "long", chid.var = "individual", alt.var = "mode") 
ml0 <- mlogit(choice ~ 1, TravelMode, 
       shape = "long", chid.var = "individual", alt.var = "mode") 
lrtest(ml,ml0) 
#--------------------- 
Likelihood ratio test 

Model 1: choice ~ wait + travel + vcost 
Model 2: choice ~ 1 
    #Df LogLik Df Chisq Pr(>Chisq)  
1 6 -192.89       
2 3 -283.76 -3 181.74 < 2.2e-16 *** 
--- 
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1