Question

我有一个名为my_list的列表，其中包含17个数据框，每个数据框都包含一个“HEIGHT”列。我想比较

样式（大致）中每个数据帧的HEIGHT密度图

https://github.com/halhen/viz-pub/blob/master/sports-time-of-day/out.png

目前，我已经想出了如何制作一个看起来很糟糕的所有17个密度图的ggplot：

library(tidyverse)

PR2016_demo_data <- read.csv("DRAFT_PR2016_demo_data.csv")
PR2016_demo_data <- arrange(PR2016_demo_data, CORAL_CD)

species_vector <- c("ACR CERV", "ACR PALM", "ORB ANNU","ORB FRAN", "ORB FAVE", "MYC FERO", "DEN CYLI", "STE INTE", "SID SIDE", "POR ASTR", "MON CAVE", "MEA MEAN", "COL NATA", "PSE STRI", "DIP LABY", "AGA AGAR", "DIC STOK")

my_list <- list()

for (i in 1:17)
{
  my_list[[i]] <- subset(PR2016_demo_data, CORAL_CD == species_vector[i])
}

ggplot(data = my_list[[1]], mapping = aes(x = HEIGHT)) +
  geom_density() + 
  geom_density(data = my_list[[2]], mapping = aes(x = HEIGHT)) +
  geom_density(data = my_list[[3]], mapping = aes(x = HEIGHT)) +
  geom_density(data = my_list[[4]], mapping = aes(x = HEIGHT)) +
  geom_density(data = my_list[[5]], mapping = aes(x = HEIGHT)) +
  geom_density(data = my_list[[6]], mapping = aes(x = HEIGHT)) +
  geom_density(data = my_list[[7]], mapping = aes(x = HEIGHT)) +
  geom_density(data = my_list[[8]], mapping = aes(x = HEIGHT)) +
  geom_density(data = my_list[[9]], mapping = aes(x = HEIGHT)) +
  geom_density(data = my_list[[10]], mapping = aes(x = HEIGHT)) +
  geom_density(data = my_list[[11]], mapping = aes(x = HEIGHT)) +
  geom_density(data = my_list[[12]], mapping = aes(x = HEIGHT)) +
  geom_density(data = my_list[[13]], mapping = aes(x = HEIGHT)) +
  geom_density(data = my_list[[14]], mapping = aes(x = HEIGHT)) +
  geom_density(data = my_list[[15]], mapping = aes(x = HEIGHT)) +
  geom_density(data = my_list[[16]], mapping = aes(x = HEIGHT)) +
  geom_density(data = my_list[[17]], mapping = aes(x = HEIGHT))

但我不明白halhen所做的所有花哨的东西，所以我很感激帮助。

比较R中的17个密度图

0 个答案: