循环遍历国家波并保存地块

Question

我有一个国家/地区随时间变化的数据集（波），结构如下：

country wave income
1       1     1000
1       2     1400
1       3     900
..     ..     .. 
1       8     40000
2       1     3000
2       2     5432

等等 15 个国家和 8 个波浪。我想在家庭总净收入上运行 foreach country 和 foreach wave qqplot 以首先检测右尾的 Pareto 分布，然后保存该图。我尝试使用以下内容列出国家和波浪

countries<-unique(final_echp$country) 
waves<-unique(final_echp$wave)

然后使用我运行的 foreach 包：

foreach(c=countries)
  foreach(w=waves) %do% { 
    qqnorm(final_echp$hi100_eur)
    qqline(final_echp$hi100_eur)
    } 

但是它崩溃了，因为它运行了几分钟而没有任何结果，我必须打破它。如何进行循环并使用相应的国家波标题保存地块？

谢谢

Answer 1

据我了解，您想为一堆国家/地区绘制一个参数并保存图以便以后检查它们。为什么不将它们保存为pdf？让我告诉你怎么做。我打算用ggplot2:

library(ggplot2)
library(dplyr)

df <- data.frame(country = c(rep('USA',20), rep('Canada',20), rep('Mexico',20)),
                 wave = c(1:20, 1:20, 1:20),
                 par = c(1:20 + 5*runif(20), 21:40 + 10*runif(20), 1:20 + 15*runif(20)))

countries <- unique(df$country)

plot_list <- list()
i <- 1

for (c in countries){
  pl <- ggplot(data = df %>% filter(country == c)) +
    geom_point(aes(wave, par), size = 3, color = 'red') +
    labs(title = as.character(c), x = 'wave', y = 'value') +
    theme_bw(base_size = 16)

  plot_list[[i]] <- pl
  i <- i + 1
}

pdf('path/to/pdf')
pdf.options(width = 9, height = 7)
for (i in 1:length(plot_list)){
  print(plot_list[[i]])
}
dev.off()

生成的 pdf 可以在 path/to/pdf 找到。