来自 vglm 的 Tobit 模型输出 - 哪个截距是哪个

编程问答 2022-05-24

问题描述

我有一个 tobit 模型

  / Y_ if Y_>0
Y=                     with Y_=X^TB+e
  \ 0 if Y_<= 0

和似然函数的对数 log l(beta,sigma|Y)=

enter image description here

我有一个结果 Y=lneconaid 和预测变量 X=(l.physint,l.poliity2,l.lneconaid,l.lnrgdp,l.lnpop,Colony)^T。我在 VGAM 包中使用 vglm 拟合了一个 tobit 模型，发现如下输出：

> m=vglm(lneconaid ~ l.physint + l.polity2 + l.lneconaid + l.lnrgdp + l.lnpop + colony,tobit(Lower=0),dat=dat)
> summary(m)

Call:
vglm(formula = lneconaid ~ l.physint + l.polity2 + l.lneconaid + 
    l.lnrgdp + l.lnpop + colony,family = tobit(Lower = 0),data = dat)

Coefficients: 
              Estimate Std. Error z value Pr(>|z|)    
(Intercept):1  0.62107    0.73535   0.845   0.3983    
(Intercept):2  1.08276    0.02071  52.288  < 2e-16 ***
l.physint      0.07140    0.03791   1.884   0.0596 .  
l.polity2      0.01395    0.01203   1.160   0.2461    
l.lneconaid    1.21495    0.02460  49.392  < 2e-16 ***
l.lnrgdp      -0.52425    0.07335  -7.147 8.84e-13 ***
l.lnpop        0.11769    0.05221   2.254   0.0242 *  
colony        -0.12038    0.36169  -0.333   0.7393    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Names of linear predictors: mu,loglink(sd)

Log-likelihood: -3879.845 on 6250 degrees of freedom

Number of Fisher scoring iterations: 8 

No Hauck-Donner effect found in any of the estimates

很明显，命名系数是 Beta 的 MLE 估计值。但我不清楚 sigma 和 beta0 的 MLE 估计是否分别对应于截距 1 和截距 2 中的哪一个。谢谢！

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

glmlog-likelihoodrr