R中的累积死亡率图答案

【问题标题】：Cumulative Mortality Graph in RR中的累积死亡率图
【发布时间】：2017-12-25 17:56:00
【问题描述】：

我需要一张图表，其中包含 4 种不同盐度下两个牡蛎种群的平均累积死亡率。

Date = c("5/29/2017","5/31/2017","6/2/2017", "5/29/2017","5/31/2017","6/2/2017","5/29/2017","5/31/2017","6/2/2017","5/29/2017","5/31/2017","6/2/2017","5/29/2017","5/31/2017 ","6/2/2017","5/29/2017 ","5/31/2017","6/2/2017", "5/29/2017","5/31/2017","6/2/2017","5/29/2017", "5/31/2017","6/2/2017")
Mortality =  c(0, 0.4, 0.5, 0.1, 0.3, 0.4, 0.2, 0.3, 0.5, 0.1, 0.2, 0.6, 0, 0.7, 0.8, 0.1, 0.6, 0.9, 0.1, 0.5, 0.9, 0.1, 0.6, 0.8)
Tank = c("1A", "1A", "1A", "1B", "1B", "1B", "5A", "5A", "5A", "5B", "5B", "5B", "1A", "1A", "1A", "1B", "1B", "1B", "5A", "5A", "5A", "5B", "5B", "5B")
Population = c("SL", "SL", "SL", "SL", "SL", "SL", "SL", "SL", "SL", "SL", "SL", "SL", "LC", "LC", "LC", "LC", "LC", "LC", "LC", "LC", "LC", "LC", "LC", "LC") = c (SL SL SL SL SL SL SL SL SL SL SL SL LC LC LC LC LC LC LC LC LC LC LC LC)

我需要平均这两个值（例如 1A 和 1B），然后为日期绘制该值。我只想将计算累积死亡率的日期显示在 x 轴上。

我在 Excel 中制作了以下图表，但不是我需要的。

【问题讨论】：

那么，你到底需要什么？
另外，您发布的代码不可重现...您应该将字符串括在引号中（单引号或双引号）并用逗号分隔，即：c('1A,'1A'...) 而不是 c(1A 1A...)

标签： r graph cumulative-line-chart

【解决方案1】：

这是我使用ggplot2 和dplyr 的解决方案

您可以在绘图之前使用其他软件包来汇总您的数据。这里我做了两组（1 = 1A 和 1B，2 = 其他坦克）。

只是出于好奇，您为什么不对这些数据使用 Kapplan-Meier 图？

Date = c("5/29/2017","5/31/2017","6/2/2017", "5/29/2017","5/31/2017","6/2/2017","5/29/2017","5/31/2017","6/2/2017","5/29/2017","5/31/2017","6/2/2017","5/29/2017","5/31/2017 ","6/2/2017","5/29/2017 ","5/31/2017","6/2/2017", "5/29/2017","5/31/2017","6/2/2017","5/29/2017", "5/31/2017","6/2/2017")
Mortality =  c(0, 0.4, 0.5, 0.1, 0.3, 0.4, 0.2, 0.3, 0.5, 0.1, 0.2, 0.6, 0, 0.7, 0.8, 0.1, 0.6, 0.9, 0.1, 0.5, 0.9, 0.1, 0.6, 0.8)
Tank = c("1A", "1A", "1A", "1B", "1B", "1B", "5A", "5A", "5A", "5B", "5B", "5B", "1A", "1A", "1A", "1B", "1B", "1B", "5A", "5A", "5A", "5B", "5B", "5B")
Population = c("SL", "SL", "SL", "SL", "SL", "SL", "SL", "SL", "SL", "SL", "SL", "SL", "LC", "LC", "LC", "LC", "LC", "LC", "LC", "LC", "LC", "LC", "LC", "LC")


df <- data.frame(Date, Mortality, Tank, Population)

library(dplyr)
library(ggplot2)

dfs <- df %>% mutate(TankGroup =  if_else(Tank %in% c('1A','1B'),1,2)) %>% 
  group_by(Date, TankGroup) %>% 
  summarise(Mortality = mean(Mortality))

ggplot(dfs, aes(x = Date, y = Mortality, group = factor(TankGroup), color = factor(TankGroup))) +
  geom_point() + geom_line() + theme_bw()

【讨论】：