我已经制作了一个图表,facet_grid用于显示每天每个治疗组中每个组中的litium百分比.
library(ggplot2)
library(Rmisc)
library(plyr)
mus2 <- summarySE(mus, measurevar="litium",
groupvars=c("treatment", "group", "day"), na.rm = TRUE)
mus2
mus3 <- mus2
mus3$group <- factor(mus3$group)
ms.chl<- ggplot(mus3, aes(x=group, y=litium, fill=treatment)) +
geom_bar(stat="identity", colour="black") + facet_grid(~day) + theme_bw()
ms.chl
Run Code Online (Sandbox Code Playgroud)
结果如下:
为此我有两个问题:
我无法为litium内容PER GROUP制作正确的误差线.我试过这个,但每次治疗我只得到误差棒.
ms.chl + geom_errorbar(aes(ymin=litium-se, ymax=litium+se), size=0.5,
width=.25,
position=position_dodge(.9)) +
facet_grid(~day)
Run Code Online (Sandbox Code Playgroud)
我想从每组的总数中得到误差线
之后,我的第二个问题是:是否有可能代表每组的绝对值和每种治疗的百分比?
数据集(mus):
litium group treatment day
0.009439528 1 Control day1
0.005115057 1 Control day1
0.009742297 1 Control day1
0.016515625 2 Control day1
0.01074537 2 Control day1
0.016300836 2 Control day1 …Run Code Online (Sandbox Code Playgroud)