如何用误差条和 p 值绘制每天的治疗方法？

Question

所以我有在第 0 天、第 15 天和实验结束时 (WEAN) 收集的数据。这些是两次处理之间的动物重量。它是这样的：

PIG TRT DAY0 DAY15 WEAN 
1 SOY 1.48  5.22  8.35    
2 SOY 1.44  5.07  7.76   
3 SOY 1.66  5.30  9.46    
4 LIN 1.94  5.91 10.10   
5 LIN 1.54  4.96  9.01 
6 LIN 1.68  5.56  9.52

一共有61头猪。使用 ggplot2 或 ggpubr，我如何使用 x = 体重收集日期和 y = 权重在单线图中执行统计并绘制处理方法，并带有误差条和 p 值？

感谢您的帮助！

Answer 1

对于您的样本数据，我认为您的目的是绘制每个处理和日期的每头猪的平均值。你可以试试

library(tidyverse)

dummy %>%
  select(-PIG) %>%
  melt() %>%
  group_by(TRT, variable) %>%
  summarise(mean = mean(value),
            sd = sd(value),
            upper = mean + 1.96 * sd,
            lower = mean - 1.96 * sd) %>%
  ungroup %>%
  ggplot(aes(variable, mean, color = TRT, group = TRT)) +
  geom_line() + geom_point() + geom_errorbar(aes(ymin = lower, ymax = upper, width = 0.2), position = ) +
  xlab("date") + ylab("mean of pig's weights")

如何用误差条和 p 值绘制每天的治疗方法？

How to plot treatment means per day with error bars and p values?

r

mean

ggplot2

errorbar

ggpubr