R 每行具有多个值的密度图

我遇到了一个似乎无法解决的问题。我有一个数据框架，可以总结如下：

peptide_sequence <- c("YKYTGFTG", "YNHRPDVRF", "MNALHHPPCS")

start_position <- c(5, 33, 79)

df <- data.frame(peptide_sequence, start_position)

df$peptide_length <- str_count(df$peptide_sequence)

df$end_position <- (df$start_position + df$peptide_length - 1)

上面的情节基本上就是我希望看到的。我正在使用的数据帧有许多行的肽，因此我无法像上面那样手动转换它，因此我正在寻找一种可以完成转换的方法，或者一种重新处理数据的替代方法，也可以生成相同类型的图。 任何帮助或建议都将不胜感激

---

谢谢

创建一个列表列，它是从开始位置到结束位置的顺序。取消此列的测试将产生所需的结果

library(tidyverse)

df %>%
  mutate(position = map2(start_position, end_position, seq)) %>%
  unnest(position) %>%
  select(peptide_sequence, position)
#> # A tibble: 27 x 2
#>   peptide_sequence position
#>   <fct>               <int>
#> 1 YKYTGFTG                5
#> 2 YKYTGFTG                6
#> 3 YKYTGFTG                7
#> 4 YKYTGFTG                8
#> 5 YKYTGFTG                9
#> 6 YKYTGFTG               10
#> 7 YKYTGFTG               11
#> 8 YKYTGFTG               12
#> 9 YNHRPDVRF              33
#> 10 YNHRPDVRF             34
#> # ... with 17 more rows

---

创建一个列表列，该列是从开始位置到结束位置的顺序。取消此列的测试将产生所需的结果

library(tidyverse)

df %>%
  mutate(position = map2(start_position, end_position, seq)) %>%
  unnest(position) %>%
  select(peptide_sequence, position)
#> # A tibble: 27 x 2
#>   peptide_sequence position
#>   <fct>               <int>
#> 1 YKYTGFTG                5
#> 2 YKYTGFTG                6
#> 3 YKYTGFTG                7
#> 4 YKYTGFTG                8
#> 5 YKYTGFTG                9
#> 6 YKYTGFTG               10
#> 7 YKYTGFTG               11
#> 8 YKYTGFTG               12
#> 9 YNHRPDVRF              33
#> 10 YNHRPDVRF             34
#> # ... with 17 more rows

---

您能否添加一个示例，说明您正试图精确获取的绘图类型，或者至少说明当前绘图代码能够获取的绘图距离？我不明白这将如何转化为密度图，我假设你指的是Hi@camille，你是对的，像链接的图基本上就是我希望得到的，x轴是总蛋白质的长度。我已经更新了这个问题，加入了一个示例，说明如何将示例数据重新转换为生成我所希望的密度图的格式。您是否可以添加一个示例，说明您正试图精确获取的图类型，或者至少您当前的绘图代码能够为您提供多大的帮助？我不明白这将如何转化为密度图，我假设你指的是Hi@camille，你是对的，像链接的图基本上就是我希望得到的，x轴是总蛋白质的长度。我已经更新了这个问题，加入了一个示例，说明如何将示例数据重新转换为生成我所希望的密度图的格式。

library(tidyverse)

df %>%
  mutate(position = map2(start_position, end_position, seq)) %>%
  unnest(position) %>%
  select(peptide_sequence, position)
#> # A tibble: 27 x 2
#>   peptide_sequence position
#>   <fct>               <int>
#> 1 YKYTGFTG                5
#> 2 YKYTGFTG                6
#> 3 YKYTGFTG                7
#> 4 YKYTGFTG                8
#> 5 YKYTGFTG                9
#> 6 YKYTGFTG               10
#> 7 YKYTGFTG               11
#> 8 YKYTGFTG               12
#> 9 YNHRPDVRF              33
#> 10 YNHRPDVRF             34
#> # ... with 17 more rows