R语言交叉试验重复测量的描述性统计与箱线图绘制问题

问题修复方案 #

1. 分组汇总返回单行结果 #

你当前代码中的median_dbp=(diastolic_bp)没有调用中位数计算函数，直接返回了原始的舒张压字段向量，导致summarise没有完成聚合，仍然保留所有原始行。修正代码如下：

mydata %>% 
  group_by(timepoint, drug) %>% 
  summarise(mean_dbp = mean(diastolic_bp, na.rm = TRUE), 
            sd_dbp = sd(diastolic_bp, na.rm = TRUE), 
            median_dbp = median(diastolic_bp), 
            p25_dbp = quantile(diastolic_bp, probs = 0.25), 
            p75_dbp = quantile(diastolic_bp, probs = 0.75))

修正后输出仅包含4行结果，对应timepoint和drug的4种组合，符合分组汇总的预期。

2. 箱线图按drug分组 #

fill是ggplot2的美学映射参数，需要放在aes()内部才能关联数据集的drug字段，你当前写在了aes()外部，不会触发分组逻辑。修正代码如下：

ggplot(data = mydata, aes(x = timepoint, y = diastolic_bp, fill = drug)) + 
  geom_boxplot()

运行后每个时间点会生成两个分别对应Drug A、Drug B的箱线图，以不同填充色区分。

内容的提问来源于stack exchange，提问作者tcvdb1992