我在一个包含列表列表的数据集上运行lmer函数时遇到了困难。

lapply(list, function(x) lapply(x, function(x) as.data.frame(x)))

lapply(list, function(x) lapply(x, function(x) 
      lmer(x[,4:7] ~ var + (1| id), x)))

model.frame.default（data = x，drop.unused.levels = TRUE，formula = x[，：变量“x”的类型（列表）无效
有没有人能建议通过lapply运行lmer的方法？

复制

library(tidyverse)
library(readxl)
require(writexl)
write_xlsx(mtcars, 'PQ_ele_com.xlsx')
        
dir = c('com', 'set', 'rit')
init = c('PQ', 'MG')
inpst = c('_ele_', '_mel_', '_col_')
           
for (i in dir) {
   for (j in init) {
      for (k in inpst){
        write_xlsx(mtcars, paste0(j, k, i, '.xlsx')) 
      }
   }
}
    
files = list.files(recursive= FALSE, full.names= FALSE)
     
init = c('PQ', 'MG')
dir = c('com', 'set', 'rit')
     
list3 = NULL
for (j in init){
   for (k in dir){
      df=c()
      for (i in files){
           if(str_detect(i, pattern = j) & str_detect(i, pattern = k)) {
             df=rbind(df,read_excel(i))}
         }
         list3[[j]][[k]] = df
      }
 }

让我们假设我想通过每个子列表拟合lmer模型。例如：

lapply(list3, function(x) lapply(x, 
    function(x) lmer(x[,3:7]  ~ vs + (1| cyl), x)))

它返回报告的错误。你认为有可能以某种方式强制它吗？