我有一系列看起来像这样的大型数据框:
Year Variable Class
1999 0.4504 a
1999 0.4244 a
1999 0.8980 b
1999 0.8901 b
2000 0.2552 a
2000 0.6504 a
2000 0.7830 b
2000 0.2030 b
... ... ...
... ... ...
2018 0.6520 a
2018 0.1403 a
2018 0.4502 b
2018 0.1240 b
有 15 个不同的班级,每年和班级至少有 500 次观察。我想做的是为每年和每个班级创建一个子集(所以 20 年 * 15 个班级 = 300 个子集),并制作一个 QQ 图来查看它们的分布,看看它们与正常情况有多大不同。
问题是,即使我可以使用循环创建子集,我也不知道如何存储这些图以便以后查看它们。
这个想法是这样的:
library(data.table)
library(ggpubr)
dataframe_list <- list.files(path="B:/Data/", pattern=".csv", full.names=TRUE)
years <- c(1999:2018)
for (i in 1:13){
db <- fread(dataframe_list[i])
for (j in 1:length(years)){
db_year <- subset(db, db$Year==years[j])
class <- levels(as.factor(db_year$Class))
for (k in 1:length(class)){
db_year_class <- subset(db_year, db_year$Class==class[k])
}
}
}
使用以下方法在第三个循环中绘制图:
ggqqplot(db_year_class$Variable)
并以某种方式存储它,直到我把它们都放在循环的末尾。由于列表中的每个数据框至少有 300 个,因此将每个数据框的图存储在一个文件中会很棒。
编辑好的,这就是我改编它的方式:
library(data.table)
library(ggpubr)
dataframe_list <- list.files(path="B:/Data/", pattern=".csv", full.names=TRUE)
df <- c(A, B, C, D, E, F, G, H, I, J, K, L, M)
years <- c(1999:2018)
plot_list <- list()
for (i in 1:13){
db <- fread(dataframe_list[i])
for (j in 1:length(years)){
db_year <- subset(db, db$Year==years[j])
class <- levels(as.factor(db_year$Class))
for (k in 1:length(class)){
db_year_class <- subset(db_year, db_year$Class==class[k])
plot_list[[k]] <- ggqqplot(db_year_class$Variable)+ggtitle(paste0("Code_0_", df[i], "_", years[j], "_", class[k]))
outfiles <- paste0("B:/Results/Plot/Code_0_", df[i], "_", years[j], ".pdf")
pdf(file=outfiles, onefile = TRUE)
for(z in seq_along(list)) {
print(plot_list[[z]])
}
dev.off()
}
}
}
结果是每年的 pdf,这不是那么方便,但否则使用列表对象将不起作用:列表将在循环结束时被覆盖。