我正在尝试比较两组患者。一组服用药物A,另一组服用药物B。我在四个不同的时间点(BB1,BB2,BB3和BB4)中了解了他们的血压信息以及其他变量。我想使用中位数,第25和第75个百分位数在每个时间点比较组,并显示组之间差异的P值。
理想情况下,我想用不同的颜色和它们之间的偏移来区分两组,但在同一张图上显示的时间点可能是中位数由一条线相连。也就是说,将首先显示两个组的BB1,然后显示BB2等。
为了我的一生,我无法弄清楚。互联网智者,请帮助我!
这有点到达那里了,但不是真的。 https://community.rstudio.com/t/2-groups-multiple-variables/13151/3
这是一个不错的解决方案,但是我确实解决了。
library(tidyverse)
library(hrbrthemes)
library(viridis)
mmHg <- KetData$MAP1
MAP1Data <- data.frame(mmHg)
MAP1Data$Drug <- Data$Drug
MAP1Data$Time <- "1. timepoint"
mmHg <- KetData$MAP2
MAP2Data <- data.frame(mmHg)
MAP2Data$Drug <- Data$Drug
MAP2Data$Time <- "2. timepoint"
...
MAPData <- do.call("rbind", list(MAP1Data, MAP2Data, ...))
ggplot(MAPData, aes(x=Time, y=mmHg, fill=Drug)) +
geom_boxplot() +
scale_fill_viridis(discrete = TRUE, alpha=0.6) +
geom_jitter(color="black", size=0.4, alpha=0.9) +
theme_ipsum() +
theme(
legend.position="right",
plot.title = element_text(size=11)
) +
ggtitle("Mean arterial pressure at n timepoints") +
xlab("")