i'm working set of results of inla package in r. these results stored in objects meaningful names can have, instance, model_a , model_b ... in current environment. each of these models i'd several processing tasks including extracting of data separate data frame, can used merge spatial data create map, etc. turning simpler, reproducible example let's assume 2 results ctl <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14) trt <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69) group <- gl(2, 10, 20, labels = c("ctl","trt")) weight <- c(ctl, trt) model_a <- lm(weight ~ group) model_b <- lm(weight ~ group - 1) i can handle steps individual model, instance: model_a_sum <- data.frame(var = character(1), model_a_value = numeric(1)) model_a_sum$var <- "intercept" model_a_sum$model_a_value <- model_a$coefficients[1] png("model_a_plot.png") plot(model_a, las = 1) dev.off() now, i'd reuse...