使用 texreg 从稳健的回归模型中获取重要的星星

Getting significance stars from a robust regression model using texreg

texreg,至少在默认情况下,不会在使用 MASS 包中的 rlm 函数拟合的稳健线性模型的系数上放置显着性星号。有什么办法吗? stargazer 实际上,正如以下 MWE 所示:

library('MASS','stargazer','texreg')

mod <- lm(speed ~ dist, cars)
modR <- rlm(speed ~ dist, cars)

stargazer(mod)
texreg(mod)

stargazer(modR)
texreg(modR)

P.S.:这可能与 MASS 的作者显然只是不喜欢 P 值有关,因为它们在 rlm 输出中被明显省略: https://stats.stackexchange.com/questions/205614/p-values-and-significance-in-rlm-mass-package-r

然而 stargazer 可以让星星变得很好。大多数人似乎认为 texreg 更好,这是我问的主要原因。

P.P.S.: 我想我会尝试使用 robust 包中的 lmRob 来代替,但还有一个更愚蠢的问题:texreg有一个签名对象的方法 lmrob,但不是 lmRob,这是我正在输出的。我假设 robust 更新并更改了这个对象名称,但 texreg 没有跟上。

如果您查看 texreg 的代码并遵循 "trail of stars",您最终会看到对未记录的 ciforce-force 函数 (texreg:::ciforce) 的调用,其 ci.force-参数默认为 FALSE。将其设置为真:

 texreg(modR, ci.force=TRUE)

\begin{table}
\begin{center}
\begin{tabular}{l c }
\hline
 & Model 1 \
\hline
(Intercept) & .21^{*}$      \
            & $[6.46;\ 9.97]$ \
dist        & [=10=].17^{*}$      \
            & $[0.13;\ 0.20]$ \
\hline
Num. obs.   & 50              \
\hline
\multicolumn{2}{l}{\scriptsize{$^*$ 0 outside the confidence interval}}
\end{tabular}
\caption{Statistical models}
\label{table:coefficients}
\end{center}
\end{table}

...然后看星星出来。

请注意,summary(modR)(来自 MASS)中的标准误差被声明用于构造,而不是 p 值,而是 t 值。