Как сохранить графики в списке как jpeg, используя lapply в R?

Question

У меня есть список из 10 графиков из model_list, для которых я использовал следующий код ниже. Я сохранил эти графики в списке var_list.

library(mixOmics)

var_list<-lapply(model_list, function(x) plotVar(x))

var_list содержит, таким образом, 10 графиков, например, под первым элементом списка:

> var_list[[1]]
                x          y Block  names pch cex     col font                  Overlap
TPI200 -0.6975577 -0.5582925     X TPI200   1   5 #388ECC    1 Correlation Circle Plots
TPI350 -0.8561514 -0.4101970     X TPI350   1   5 #388ECC    1 Correlation Circle Plots
TPI500 -0.9403552 -0.1074518     X TPI500   1   5 #388ECC    1 Correlation Circle Plots
TPI700 -0.9256605  0.3070954     X TPI700   1   5 #388ECC    1 Correlation Circle Plots
TPI900 -0.8697037  0.4699423     X TPI900   1   5 #388ECC    1 Correlation Circle Plots

Я хочу сохраните эти графики из этого списка как jpeg (в результате получится 10 различных jpeg). Я использовал следующий код, и R создает 10 изображений, но все изображения одинаковы (поэтому для остальных создается и дублируется только первый график).

lapply(1:length(model_list), function (x) {
  jpeg(paste0(names(model_list)[x], ".jpg"))
  lapply(model_list, function(x) plotVar(x))
  dev.off()
})

Я видел похожие вопросы, но я не могу найти правильное решение, чтобы иметь jpg для каждого графика для каждого фрейма данных в списке! Как я могу это решить? Заранее большое спасибо!

По этой ссылке вы можете найти dput(model_list[[1]]).

Answer 1

С данными, предоставленными вами в аналогичном сообщении, здесь можно найти возможное решение вашей проблемы. Лучше, если вы работаете с model_list, потому что при преобразовании в var_list все данные становятся графическими элементами. Следующий код содержит копию model_list с использованием datalist, но в вашей реальной проблеме он должен быть у вас, а также должен включать имена для каждого из компонентов списка:

library(mixOmics)
#Data
datalist <- list(df1 = structure(list(OID = c(-1, -1, -1, -1, -1, -1), POINTID = c(1, 
2, 3, 4, 5, 6), WETLAND = c("no wetl", "no wetl", "no wetl", 
"wetl", "wetl", "wetl"), TPI200 = c(70, 37, 45, 46, 58, 56), 
    TPI350 = c(67, 42, 55, 58, 55, 53), TPI500 = c(55, 35, 45, 
    51, 53, 51), TPI700 = c(50, 29, 39, 43, 49, 49), TPI900 = c(48, 
    32, 41, 46, 47, 46), TPI1000 = c(46, 16, 41, 36, 46, 46), 
    TPI2000 = c(53, 17, 53, 54, 54, 54), TPI3000 = c(47, 35, 
    47, 47, 47, 47), TPI4000 = c(49, 49, 49, 49, 49, 49), TPI5000 = c(63, 
    63, 63, 62, 62, 61), TPI2500 = c(48, 26, 48, 49, 49, 49)), row.names = c(NA, 
6L), class = "data.frame"), df2 = structure(list(OID = c(-1, 
-1, -1, -1, -1, -1), POINTID = c(1, 2, 3, 4, 5, 6), WETLAND = c("no wetl", 
"no wetl", "no wetl", "wetl", "wetl", "wetl"), TPI200 = c(70, 
37, 45, 46, 58, 56), TPI350 = c(67, 42, 55, 58, 55, 53), TPI500 = c(55, 
35, 45, 51, 53, 51), TPI700 = c(50, 29, 39, 43, 49, 49), TPI900 = c(48, 
32, 41, 46, 47, 46), TPI1000 = c(46, 16, 41, 36, 46, 46), TPI2000 = c(53, 
17, 53, 54, 54, 54), TPI3000 = c(47, 35, 47, 47, 47, 47), TPI4000 = c(49, 
49, 49, 49, 49, 49), TPI5000 = c(63, 63, 63, 62, 62, 61), TPI2500 = c(48, 
26, 48, 49, 49, 49)), row.names = c(NA, 6L), class = "data.frame")) 
#Function
custom_splsda <- function(datalist, ncomp, keepX, ..., Xcols, Ycol){
  Y <- datalist[[Ycol]]
  X <- datalist[Xcols]
  res <- splsda(X, Y, ncomp = ncomp, keepX = keepX, ...)
  res
}
#Create model_list, you must have the object
model_list <- lapply(datalist, custom_splsda,
                     ncomp = 2, keepX = c(5, 5),
                     Xcols = 4:8, Ycol = "WETLAND")

Далее l oop для графиков:

#Loop
for(i in 1:length(model_list))
{
  jpeg(paste0(names(model_list)[i], ".jpg"))
  plotVar(model_list[[i]],title = names(model_list)[i])
  dev.off()
}

Это создаст графики в вашей папке, как вы можете видеть здесь:

And also the plots that change (see titles):

введите описание изображения здесь