【问题标题】:How can I display significant stars using Zelig package如何使用 Zelig 包显示重要的星星
【发布时间】:2019-12-17 08:58:21
【问题描述】:

我想在 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

系数右侧没有星号。

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

【问题讨论】:

    标签: r regression non-linear-regression r-zelig


    【解决方案1】:

    即使有点晚了:

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

    【讨论】:

      猜你喜欢
      • 1970-01-01
      • 2023-03-31
      • 2023-03-17
      • 1970-01-01
      • 2023-03-02
      • 2011-04-20
      • 2010-11-10
      • 2020-10-23
      • 2014-05-14
      相关资源
      最近更新 更多