我想在 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
系数右侧没有星号。
有谁知道如何显示重要的星星?