使用texreg从强大的回归模型中获得重要性

时间:2018-01-17 23:29:08

标签: r stargazer texreg

默认情况下,

texreg至少不会使用来自rlm包的MASS函数拟合的鲁棒线性模型将有效星置于系数上。有什么方法可以做到吗?实际上,stargazer确实如以下MWE所示:

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

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

stargazer(mod)
texreg(mod)

stargazer(modR)
texreg(modR)

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

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

PPS:我认为我尝试使用lmRob软件包中的robust代替,但问题更简单:texreg有一个方法对于签名lmrob但不是lmRob的对象,这是我输出的内容。我假设robust已更新并更改了此对象名称,但texreg尚未找到。

1 个答案:

答案 0 :(得分:1)

如果您查看texreg的代码,请按照stars&#34;的跟踪进行操作。您最终会查看对未记录的ciforce - 强制函数(texreg:::ciforce)的调用,其ci.force - 参数默认为FALSE。将其设置为TRUE:

 texreg(modR, ci.force=TRUE)

\begin{table}
\begin{center}
\begin{tabular}{l c }
\hline
 & Model 1 \\
\hline
(Intercept) & $8.21^{*}$      \\
            & $[6.46;\ 9.97]$ \\
dist        & $0.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值。