如何在texreg表中包含赔率比？

Question

所以我尝试在我的texreg表中显示赔率比：

library(texreg)

df <- mtcars

logit <- glm(
  formula = am ~ cyl,
  data = df,
  family = binomial(link = "logit")
)

odds_ratio <- exp(coef(logit))

screenreg(logit, odds_ratio)

不幸的是，这会产生以下错误：

Error in screenreg(logit, odds_ratio) : 
  The 'file' argument must be a character string.

如何包含赔率比？

常规模型(我在上面称为logit )显示赔率的对数赔率/对数，这不是正确的吗？

Answer 1

要将对数赔率替换为赔率，您可以执行以下操作：

library("texreg")
logit <- glm(formula = am ~ cyl,
             data = mtcars,
             family = binomial(link = "logit"))
screenreg(logit, override.coef = exp(coef(logit)))

override.coef参数只是将系数替换为您提供的向量。如果你需要标准错误，p值等，也有类似的论据。

结果是：

=========================
                Model 1  
-------------------------
(Intercept)      43.71 * 
                 (1.55)  
cyl               0.50 **
                 (0.25)  
-------------------------
AIC              37.95   
BIC              40.88   
Log Likelihood  -16.98   
Deviance         33.95   
Num. obs.        32      
=========================
*** p < 0.001; ** p < 0.01; * p < 0.05

请注意，在上面的示例中没有转换标准错误。首先转换它们是没有意义的，因为它们需要是不对称的。因此，您不妨将它们从表中删除。为此，您可以将结果保存到中间texreg对象中，操作此对象，然后将其移交给screenreg函数，如下所示：

tr <- extract(logit)
tr@coef <- exp(tr@coef)
tr@se <- numeric()
screenreg(tr)

结果：

=========================
                Model 1  
-------------------------
(Intercept)      43.71 * 
cyl               0.50 **
-------------------------
AIC              37.95   
BIC              40.88   
Log Likelihood  -16.98   
Deviance         33.95   
Num. obs.        32      
=========================
*** p < 0.001; ** p < 0.01; * p < 0.05