运行分析相同的数据集并将系数提取到 table

Question

我有 198 个独特的数据集，每个数据集中的所有变量都是相同的。我想对每个数据集进行相同的分析（见下文）。从每次进行的分析中，我想从模型结果中提取系数并将它们存储在一个 table 中。为此，我编写了以下代码：

 one_table <- data.frame() #dataframe I will use to store results in later
for(i in 1:198){
  load(file = paste("Dataset_",i,".RData",sep = "")) #loads a dataset
  Dataset = get(paste0("Dataset_", i)) #renames the loaded Dataset_i to Dataset, important for model below
  model_result = try(summary(lme(A ~ B + C + B*C + D + E, random = ~1|company,na.action=na.omit,
                                 control=lmeControl(opt='optim'),method='REML',
                                 data=Dataset,correlation = corAR1())),silent=TRUE) #model
  coefficients = coef(model_result) #extracts coefficients from model results
  coefficients_df <- as.data.frame(coefficients) #puts coefficients into a table
  one_table <- rbind(one_table, coefficients_df) #updates table with each iteration of i

现在，代码运行到数据集编号 96，我收到以下错误：$ operator is invalid for atomic vectors.我不太确定如何继续，这个错误指的是什么？

编辑：我想我找到了问题所在。该模型在数据集编号 97 上没有运行，我收到以下与 B*C

相关的错误

[1] "Error in logLik.reStruct(object, conLin) : \n  NA/NaN/Inf in foreign function call (arg 3)\n"
attr(,"class")
[1] "try-error"
attr(,"condition")
<simpleError in logLik.reStruct(object, conLin): NA/NaN/Inf in foreign function call (arg 3)>

Answer 1

发生 B*C 错误是因为在数据集编号 97 中计算错误了一个变量 - 愚蠢的错误！

注意：您还可以将模型结果存储在列表中

list.res.H = list()
for(i in 1:198){
  load(file = paste("Dataset_",i,".RData",sep = "")) #loads a dataset
  Dataset = get(paste0("Dataset_", i)) #renames the loaded Dataset_i to Dataset, important for model below
  model_result = try(summary(lme(A ~ B + C + B*C + D + E, random = ~1|company,na.action=na.omit,
                                 control=lmeControl(opt='optim'),method='REML',
                                 data=Dataset,correlation = corAR1())),silent=TRUE) #model
  list.res.H[[i]] = coef(model_result) #extracts coefficients from model results
}

运行分析相同的数据集并将系数提取到 table

Running analysis on identical datasets and extracting coefficients into a table

r

for-loop

lme4

coefficients

运行 分析相同的数据集并将系数提取到 table

Running analysis on identical datasets and extracting coefficients into a table

r

for-loop

lme4

coefficients

运行分析相同的数据集并将系数提取到 table