如何用R循环实现多数据框的Lavaan CFA重复分析？

我太懂你这种重复复制粘贴改编号的痛苦了——8个波手动操作不仅容易手抖写错，后期要调整模型参数或者分析设置时，还要挨个改8遍，简直是浪费时间。之前尝试for循环没成功？大概率是动态调用数据框或者结果存储的方式没搞对，下面给你两个亲测好用的方案，彻底告别复制粘贴：

方案1：用列表+lapply批量处理（最简洁） #

这种方式代码更紧凑，R处理列表的效率也很高，步骤如下：

第一步：把所有时间波的数据框放进一个列表 #

先把分散的bds_1到bds_8打包成一个列表，方便统一遍历：

# 创建数据列表，按顺序放入8个波的数据框
bds_waves <- list(bds_1, bds_2, bds_3, bds_4, bds_5, bds_6, bds_7, bds_8)
# 给列表命名，后续看结果时能清楚对应哪个波
names(bds_waves) <- paste0("wave_", 1:8)

第二步：定义你的CFA模型 #

和你之前的模型定义完全一样，不用改：

model_ADI_aff <- "aff =~ bds89 + bds39 + bds50 + bds29 + bds84 + bds49 + bds70 + bds88 + bds11 + bds28 "

第三步：批量运行CFA并存储结果 #

用lapply遍历列表里的每个数据框，自动运行CFA，结果会存在一个新的列表里：

# 批量执行CFA，返回包含所有波拟合结果的列表
fit_list <- lapply(bds_waves, function(current_data) {
  cfa(model = model_ADI_aff,
      data = current_data,
      missing = "fiml",
      estimator = "MLR",
      se = "robust.huber.white",
      test = "yuan.bentler")
})

第四步：批量输出结果 #

写个小循环，一次性输出所有波的summary和修正指数：

# 遍历结果列表，逐个输出
for (wave_name in names(fit_list)) {
  cat("\n====================================\n")
  cat("===== 分析结果：", wave_name, " =====\n")
  cat("====================================\n")
  # 输出标准化结果和拟合指标
  summary(fit_list[[wave_name]], standardized = TRUE, fit.measures = TRUE)
  
  cat("\n------------------------------------\n")
  cat("----- 修正指数（MI≥3）：", wave_name, " -----\n")
  cat("------------------------------------\n")
  # 输出排序后的修正指数
  print(modindices(fit_list[[wave_name]], sort.=TRUE, minimum.value=3))
}

方案2：用for循环（更直观，适合新手） #

如果你觉得lapply有点抽象，用for循环也能实现，关键是用get()函数动态获取数据框：

# 初始化空列表存储所有拟合结果
fit_list <- list()

# 遍历8个时间波
for (i in 1:8) {
  # 动态获取当前波的数据框（比如i=1时，get("bds_1")就拿到bds_1）
  current_data <- get(paste0("bds_", i))
  # 运行CFA并把结果存入列表
  fit_list[[paste0("wave_", i)]] <- cfa(model_ADI_aff,
                                        data = current_data,
                                        missing = "fiml",
                                        estimator = "MLR",
                                        se = "robust.huber.white",
                                        test = "yuan.bentler")
  
  # 实时输出当前波的结果
  cat("\n====================================\n")
  cat("===== 分析结果：Wave", i, " =====\n")
  cat("====================================\n")
  summary(fit_list[[paste0("wave_", i)]], standardized = TRUE, fit.measures = TRUE)
  
  cat("\n------------------------------------\n")
  cat("----- 修正指数（MI≥3）：Wave", i, " -----\n")
  cat("------------------------------------\n")
  print(modindices(fit_list[[paste0("wave_", i)]], sort.=TRUE, minimum.value=3))
}

为什么之前的for循环可能失败？ #

大概率是这两个问题：

• 直接写bds_i这种代码，R会把它当成一个叫bds_i的变量，而不是bds_1、bds_2，必须用get(paste0("bds_", i))来动态调用数据框；
• 没有把结果存入列表，而是每次循环都覆盖同一个变量（比如fit_ADI_aff），最后只保留了最后一个波的结果。

小提示 #

把所有结果存在fit_list里之后，后续提取单个波的结果也很方便，比如想单独看第3波的拟合情况，直接写fit_list[["wave_3"]]就行，比单独存8个变量整洁太多了！

内容的提问来源于stack exchange，提问作者Eva Billen