带有每次访问每个序列的置信区间的意大利面条图

时间:2019-07-30 10:14:02

标签: r ggplot2

我正在尝试绘制一些纵向数据(浓度与时间的关系),用于不同的治疗序列和随访。我还试图为每个系列添加95%的置信区间。

到目前为止,我只能为每个时间点绘制一条具有置信区间的曲线,并且每次绘制一次(使用构面网格每次访问),如图所示。my plot

我的数据框的格式如下:

  • LBTPTNUM(对应于浓度的时间点)

  • 平均值:是相对于时间绘制的浓度

  • lower / upper:这是95%CI的边界
  • USUBJID:每个主题的ID
  • VISITNUM:造访次数

我尝试了以下代码:

ggplot(data = CI_spag(), 
       aes(x = LBTPTNUM, 
           y = mean,
           colour = as.factor(TRTSEQP))) + 
  geom_line(aes(linetype = factor(TRTSEQP))) + 
  facet_grid(~ VISITNUM , labeller = label_both)  + 
  geom_segment(data = CI_spag(), 
               aes(x = LBTPTNUM ,y = lower, 
                   xend = LBTPTNUM, yend = upper, 
                   group = LBTPTNUM, color = "green", 
                   linetype = "dashed" ))

但是它并没有返回我所期望的东西,就像这样: what I expect

我用geom_segment表示CI,因为我还没有找到其他解决方案。

有什么帮助吗?

谢谢

0 个答案:

没有答案