如何将 LME4 输出提取到乳胶表？

Question

我通常使用 modelsummary() 或 stargazer() 将 table 提取到 R 中的乳胶中。但是，它们似乎不适用于 lme4。有人知道将这个可重现的示例提取到漂亮的乳胶中的方法吗table？

这是我的模型：

df <- tibble(
      y = rnorm(100000),
      x1 = rnorm(100000),
      x2 = rnorm(100000),
      school =sample.int(300,size=100000,replace=TRUE)-1,
      classes =sample.int(100,size=100000,replace=TRUE)-1
      )
    
    df$school = as.factor(df$school)
    df$classes = as.factor(df$classes)
    
    library(lme4)
    model1 <- lmer(y~ x1 + x2 +  
                  (x1 + x2 |classes) +
                  (x1 + x2 |school), data=df)

Answer 1

包 texreg 存在。

library(texreg)
texreg(model1)

或者，您可以使用 tab_model() 并将 html 输出转换为乳胶（未经测试）。

library(sjPlot)
model_html <- tab_model(model1)

library(rrtable)
HTMLcode2latex(model_html$page.content)

Answer 2

从版本 0.6.4（现在在 CRAN 上）开始，modelsummary 支持开箱即用的 lme4::lmer 模型。您只需要更新您的软件包并重试。

update.packages("modelsummary")

library(lme4)
library(modelsummary)

N <- 1000
df <- data.frame(
      y = rnorm(N),
      x1 = rnorm(N),
      x2 = rnorm(N),
      school =sample.int(300,size=N,replace=TRUE)-1,
      classes =sample.int(100,size=,replace=TRUE)-1)
    
df$school = as.factor(df$school)
df$classes = as.factor(df$classes)
    
model1 <- lmer(y~ x1 + x2 +  
              (x1 + x2 |classes) +
              (x1 + x2 |school), data=df)

modelsummary(model1, "markdown")

|            |  Model 1  |
|:-----------|:---------:|
|(Intercept) |   0.009   |
|            |  (0.036)  |
|x1          |  -0.037   |
|            |  (0.032)  |
|x2          |   0.042   |
|            |  (0.033)  |
|Num.Obs.    |   1000    |
|R2 Marg.    |   0.003   |
|R2 Cond.    |           |
|AIC         |  2858.0   |
|BIC         |  2936.6   |
|Log.Lik.    | -1413.013 |

如何将 LME4 输出提取到乳胶表？

How to extract LME4 output to latex tables?

r

lme4

stargazer

broom

modelsummary