首页 > 解决方案 > 如何使用 Zelig 包显示重要的星星

问题描述

我想在 Zelig 回归的结果中显示重要的星星,并为 tobit 模型提供稳健的标准误差。代码是这样的。但是没有显着的星星,比如使用 AER 包的正常 tobit 回归(tobit 命令)

library(Zelig)
# input
dat <- read.csv("http://www.omori.e.u-tokyo.ac.jp/STATA/Sample/select.csv")
head(dat)

# robust tobit
rb.tobit <- zelig(HOUR~CHILD+AGE+EDU+WAGE+HINC, below=0, above=Inf, model="tobit", data=dat, robust=T)
summary(rb.tobit)
Model: 

Call:
z5$zelig(formula = HOUR ~ CHILD + AGE + EDU + WAGE + HINC, below = 0, 
    above = Inf, robust = T, data = dat)

Observations:
         Total  Left-censored     Uncensored Right-censored 
           753            325            428              0 

Coefficients:
              Estimate Std. Error z value Pr(>|z|)
(Intercept) 1130.70774  443.05762   2.552  0.01071
CHILD       -806.15058  245.32541  -3.286  0.00102
AGE          -24.10593    7.68759  -3.136  0.00171
EDU           22.86477   31.14134   0.734  0.46281
WAGE         201.39910   34.23250   5.883 4.02e-09
HINC         -53.92795   23.04316  -2.340  0.01927
Log(scale)     7.04989    0.04425 159.316  < 2e-16

Scale: 1153 

Gaussian distribution
Number of Newton-Raphson Iterations: 4 
Log-likelihood: -3812 on 7 Df
Wald-statistic: 75.14 on 5 Df, p-value: 8.6858e-15 

Next step: Use 'setx' method

系数右侧没有星号。

有谁知道如何显示重要的星星?

标签: rregressionnon-linear-regressionr-zelig

解决方案


即使有点晚了:

-> 输入:summary(rb.tobit, signif.stars=TRUE ) 那么它应该可以工作。


推荐阅读