根据计数改变线图的厚度

Question

数据框“id”具有“year”、“id”和“matriline”列，其中每一行都是一个事件。我想按母系计算每年的事件数量，所以我这样做了：

events.bymatr = 
id %>%
group_by(year, matr, .drop = FALSE) %>%
dplyr::summarise(n = n()) %>%
ungroup()
events.bymatr

我绘制了按母系随时间变化的事件数量的折线图。

ggplot(events.bymatr, aes(x=year, y=n, group=matr)) + geom_line(aes(color=matr))

我的问题是双重的：

1. 有没有一种方法可以重新创建这个折线图，其中线条的粗细由每个母系有多少个 ID 决定？我想这将涉及重塑我上面的数据，但是当我尝试 group_by(year,matr,id,.drop=FALSE) 时，我的数据变得很不稳定。

2. 我想更改调色板，使每种颜色都非常独特 - 如何附加新调色板？我尝试使用 这个 c25 调色板 与此代码，但它使我所有线条消失。 ggplot(events.bymatr, aes(x=year, y=n, group=matr)) + geom_line(aes(color=c25))

提前非常感谢！

“id”的输出（缩短为每列前五行）：

> dput(id)
    structure(list(date = structure(c(8243, 8243, 8243, 8248, 8947,
    class = "Date"), year = c(1992L, 1992L, 1992L, 1992L, 1994L), 
    event.id = c(8L, 8L, 8L, 10L, 11L), id = structure(c(51L, 55L, 59L, 
    46L, 51L), .Label = c("J11", "J16", "J17", "J2", "J22"),
    class = "factor"), sex = structure(c(1L, 2L, 2L, 1L, 1L),
    .Label = c("0", "1"), class = "factor"), age = c(28L, 12L, 6L, 42L, 
    30L), matr = structure(c(20L, 20L, 20L, 11L, 20L), .Label = c("J2", 
    "J4", "J7", "J9", "K11"), class = "factor"),
    matralive = structure(c(2L, 2L, 2L, 2L, 2L),
    .Label = c("0", "1"), class = "factor"), pod = structure(c(3L, 3L, 
    3L, 3L, 3L), .Label = c("J", "K", "L"), class = "factor")),
    row.names = c(NA, -134L), class = c("tbl_df", "tbl", "data.frame"))

events.bymatr 的输出：

> dput(events.bymatr)
      structure(list(year = c(1992L, 1992L, 1992L, 1992L, 1992L),
      matr = structure(c(1L, 2L, 3L, 4L, 5L), .Label = c("J2", "J4", 
      "J7", "J9", "K11"), class = "factor"), n = c(0L, 0L, 0L, 0L, 0L)), 
      row.names = c(NA, -380L), class = c("tbl_df", "tbl", 
      "data.frame"))

Answer 1

正如 @r2evans 所指出的，要清楚地区分多种颜色是非常困难的。我在这里使用了一个 20 色标尺的示例，效果相当不错，但即便如此，仍有一些色标很难区分。这是使用 dplyr 附带的 storms 数据集的尝试。

library(dplyr)
storms %>%
  group_by(name, year) %>%
  summarize(n = n(), .groups = "drop") %>%   # = number of name per year View
  tidyr::complete(name, year = 1975:2015, fill = list(n = 0)) %>% 
  group_by(name) %>%
  mutate(total = sum(n)) %>%                 # = number of name overall
  ungroup()  %>% 
  filter(total %% 12 == 0) %>% # Arbitrary, to reduce scope of data for example
  ggplot(aes(year, n, color = name, size = total, group = name)) +
  geom_line() +
  guides(color = guide_legend(override.aes = list(size = 3))) +
  ggthemes::scale_color_tableau(palette = "Tableau 20")