使用 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 值。
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 值。