下面的代码对我来说很有用：

# create a new list by extracting only models from the list sfa_out
newlist <- as.list(sfa_out$models)

# define a function which gets the required output from each model

outputs <- function(newlist) {
  
# create an empty data.frame
  eff_out <- data.frame()
  
# loop through each model in the newlist
  for(i in 1:length(newlist)) {

      eff <- as.data.frame(efficiencies(newlist[[i]])) %>%
      mutate(col = row.names(efficiencies(newlist[[i]]))) %>%
      pivot_longer(cols = 1:23,
                   names_to = "year",
                   values_to = "efficiency") %>%
      drop_na()

# current model output  
   model_output <- data.frame(eff)
  
# append output to empty dataset we created before
   eff_out <- rbind(eff_out, model_output)
  }
  return(eff_out)
}
# apply the function to the list containing models
sfa_eff_outputs <- outputs(newlist)