R-如何汇集多重数据插补后SEM回归的结果

我正在R中对SEM进行多重数据插补。我正在测试一个特定的数据插补算法，这就是为什么我“手动”而不是使用MItools进行汇总

for (i in 1:m){ # m is number of imputation
 # Imputation
 df_imputed <-  myImputationAlgorithm(df)

 # Fit model
 fitted_model<- sem(model, data=df_imputed, se="bootstrap",bootstrap=100)     

 # Save parameters for pooling
 # https://rdrr.io/cran/lavaan/man/lavInspect.html
 betas <- append(betas, lavInspect(fitted_model, "coef"))
 vars <- append(vars,   lavInspect(fitted_model, "vcov"))

}
# Pooling
summary(MIcombine(betas, vars))

for（i in 1:m）{#m是插补数
#插补
df_插补我想我可以直接从
lavan
结果矩阵中提取东西

m<- 10 # Number of imputation

estimates <- as.data.frame(matrix(NA, nrow=29, ncol = m))   # Estimates
standErr <- as.data.frame(matrix(NA, nrow=29, ncol = m))    # Standard deviations
zvalue <- as.data.frame(matrix(NA, nrow=29, ncol = m))      # Z-value
pvalue <- as.data.frame(matrix(NA, nrow=29, ncol = m))      # P-value
tli <- as.data.frame(matrix(NA, nrow=1, ncol = m))          # TL 
cfi <- as.data.frame(matrix(NA, nrow=1, ncol = m))          # CFI
rmsea <- as.data.frame(matrix(NA, nrow=1, ncol = m))        # RMSEA 

for (i in 1:m){
  print(cat("Imputation #",i,"\n", sep= ""))
  df_imputed <-  myImputationAlgorithm() 

  # Estimation
  fitted_model<- sem(model, data=df_imputed, se="bootstrap",bootstrap=100)         

  # Extrcat results
  estimates[[i]] <- parameterEstimates(fitted_model)$est    # Estimate
  standErr[[i]]  <- parameterEstimates(fitted_model)$se     # Standard Error
  zvalue[[i]]    <- parameterEstimates(fitted_model)$z      # z-value  
  pvalue[[i]]    <- parameterEstimates(fitted_model)$pvalue # p-value
  tli[[i]]       <- inspect(fitted_model, "fit")["tli"]     # TLI
  cfi[[i]]       <- inspect(fitted_model, "fit")["cfi"]     # CFI
  rmsea[[i]]     <- inspect(fitted_model, "fit")["rmsea"]   # RMSEA


 }


  # Pooling
  mean_estimates  <- rowMeans(estimates)
  mean_standErr   <- rowMeans(standErr)
  mean_zvalue     <- rowMeans(zvalue)
  mean_pvalue     <- rowMeans(pvalue)
  mean_tli        <- rowMeans(tli)
  mean_cfi        <- rowMeans(cfi)
  mean_rmsea      <- rowMeans(rmsea)

m