繁体   English   中英

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

[英]Getting significance stars from a robust regression model using texreg

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

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 : rlm 和意义-在-RLM-质量包-R

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

PPS:我以为我会尝试从robust软件包中使用lmRob ,但是还有一个更愚蠢的问题: texreg有一个用于签名lmrob的对象的方法,但不是lmRob ,这就是我要输出的方法。 我假设robust更新并更改了该对象名称,但texreg尚未赶上。

如果您查看texreg的代码并遵循“ stars踪迹”,您最终将最终看到对未记录的ciforce -force函数( texreg:::ciforce )的调用,该函数的ci.force -parameter默认为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值。

暂无
暂无

声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.

 
粤ICP备18138465号  © 2020-2024 STACKOOM.COM