mixed为混合模型提供tidiers。您也可以使用purrr::map_dfr()代替lapply()来获取一个 Dataframe 中的所有系数。

library(lme4)
library(lmerTest)
library(tidyverse)
library(broom.mixed)
set.seed(13)

allModels <- map_dfr(
  set_names(depVarList), 
  \(x) {
    tidy(lmer(
       formula = paste0(x, " ~ visit + (1| id)"),
       data = df_rand0, 
       na.action = na.omit
   ))
  },
  .id = "dv"
)

allModels

# A tibble: 12 × 9
   dv    effect   group    term          estim…¹ std.e…² stati…³    df   p.value
   <chr> <chr>    <chr>    <chr>           <dbl>   <dbl>   <dbl> <dbl>     <dbl>
 1 x1    fixed    <NA>     (Intercept)   6372.     286.   22.3    32.9  2.00e-21
 2 x1    fixed    <NA>     visit          229.     278.    0.821  21.0  4.21e- 1
 3 x1    ran_pars id       sd__(Interce…  973.      NA    NA      NA   NA       
 4 x1    ran_pars Residual sd__Observat…  923.      NA    NA      NA   NA       
 5 x2    fixed    <NA>     (Intercept)   2278.     123.   18.5    42.0  8.84e-22
 6 x2    fixed    <NA>     visit          -19.2    174.   -0.110  42.0  9.13e- 1
 7 x2    ran_pars id       sd__(Interce…    0       NA    NA      NA   NA       
 8 x2    ran_pars Residual sd__Observat…  578.      NA    NA      NA   NA       
 9 x3    fixed    <NA>     (Intercept)    314.      12.1  26.0    42.0  1.63e-27
10 x3    fixed    <NA>     visit           -8.82    17.1  -0.516  42.0  6.09e- 1
11 x3    ran_pars id       sd__(Interce…    0       NA    NA      NA   NA       
12 x3    ran_pars Residual sd__Observat…   56.7     NA    NA      NA   NA       
# … with abbreviated variable names ¹​estimate, ²​std.error, ³​statistic