1

我的问题很容易理解，但我找不到解决方案。

我使用 mLLogit 库在 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)

我得到的输出如下：

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

我想获得McFadden 和似然比检验

怎么了？？

4

3 回答 3

2

我将把寻找麦克法登的工作留给你。您应该在发布之前完成此操作。要进行 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-23T16:16:42.513 回答
2

回答我自己的问题，并不是我不了解 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。

于 2015-03-12T13:26:03.423 回答
0

对于可能遇到类似问题的其他人。如果您在 mlogit 中运行随机参数模型（例如 MIXL），则不会报告这些统计数据（McFadden 和 LRT）。在这种情况下不估计空对数似然，并且不存在空 LL 值。您可以如上所述手动执行此操作（即仅使用截距估计空模型，然后是完整模型）。

于 2019-06-13T20:55:28.357 回答