McFadden R^2和似然比測試未在mLogit中顯示R

我的問題相當直接瞭解，但我無法找到解決方案。McFadden R^2和似然比測試未在mLogit中顯示R

我使用mLogit庫使用以下代碼中R：

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)

輸出I得到的是以下內容：

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

我想獲得的麥克法登和似然比測試

什麼是錯？

來源

2014-11-23 Manos C

輕軌需要兩種車型。你可以建立一個空模型並運行'lrtest'。要找到McFadden測試，你只需要用'??「McFadden」'在R控制檯上搜索 – 2014-11-23 16:09:19

回答我自己的問題，這並不是說我不明白MC法登的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。

來源

2015-03-12 13:26:03

我要離開搜尋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

來源

2014-11-23 16:16:42

McFadden R^2和似然比測試未在mLogit中顯示R

回答

相關問題