我有一个看起来像这样的数据:
ENSG00000211521 MIR665 x 89
....
ENSG00000207793 MIR432 y 50
....
我想做的是制作一个条形图并在每个条形上添加“人口百分比”。例如,值为 100 的值类别 'y' 的百分比为 45.6 (132/289),这是因为 100 中有 132 个,总人口为 289 (对于每个 "x" 和 "y")。
最后,我想要一个大致如下所示的情节:
但我坚持使用以下代码。正确的做法是什么?
library(ggplot2)
dat.m <- read.delim("http://dpaste.com/1269939/plain/",sep="")
colnames(dat.m) <- c("ensg","mirna_hgc","variable","value")
qplot(value,data=dat.m, geom="bar", binwidth=1, origin=-0.05, xlim=c(50,100),ylim=c(0,75), facets=variable~.,main="")+
xlab("Value")+
ylab("Frequency")+
theme(legend.position="none")
更新:计算百分比
可以使用此代码获得上图中的百分比。但不知何故,我找不到将它们包含到 qplot 中的方法:
dat.m <- read.delim("http://dpaste.com/1269939/plain/",sep="")
colnames(dat.m) <- c("ensg","mirna_hgc","variable","value")
# the following steps can be applied for "x"
y <- subset(dat.m,dat.m$variable=="y")
y.df <- data.frame(table(y$value))
y.df$percentage <- ((y.df$Freq)/sum(y.df$Freq) * 100)
y.df