在整理转录组下游的时候,看到中科新生命的报告中的基因表达水平分布部分有这么一个图
从图中可以非常直观的看出来不同样本在不同表达区间的分布情况。由于报告没有给出源代码,我们模仿的画一画。...[2]
# 载入R包
library(reshape)
# 宽变长
longdata <- melt(data = express_cpm)
# 将数据划分成6个区间
cut_data <- cut...(as.numeric(longdata$value), breaks = c(0,1,5,10,30,50,Inf),right = F,include.lowest = TRUE)
# 将分组结果添加到...cut_data) <- c("0<=CPM<1", "1<=CPM<5","5<=CPM<10","10<=CPM<30","30<=CPM<50","50<=CPM<+∞")
画图
画图需要用到R包ggplot2...[3]
# 载入R包
library(ggplot2)
# 画图
ggplot(longdata, aes(x = X2, fill = group, y = value)) +
geom_bar