最终结果
image.png
直接再密度图上是无法实现分段设置颜色的，要变成geom_line才能行
因此核心是得到geom_density作图使用的值
方法是先做出图p
再利用ggplot_build算出画图时候实际使用的值
image.png
p <- ggplot(data,aes(x=log2FoldChange))+
geom_density()
pb <- ggplot_build(p)
data = pb[['data']][[1]]
ggplot(data,aes(x=x,y=y))+geom_line()
image.png
image.png
添加分组信息
image.png
image.png
能看到分组变量是group
颜色 分组信息都加上
p <- ggplot(data,aes(x = log2FoldChange, group=type,colour=type))+
geom_density()
pb <- ggplot_build(p)
Data = pb[['data']][[1]]
head(Data)
image.png
Data$colour=ifelse((Data$x<1 & Data$x> -1),'gray',Data$colour)
image.png
p <- ggplot(data,aes(x = log2FoldChange, group=type,colour=type))+
geom_density()
pb <- ggplot_build(p)
Data = pb[['data']][[1]]
head(Data)
Data$colour=ifelse((Data$x<1 & Data$x> -1),'#C0C0C0',Data$colour)
ggplot(Data,aes(x=x,y=y,group=group,colour=colour))+geom_line()+xlim(-5,5)+
geom_vline(xintercept=-1,linetype = "dashed",colour='#C0C0C0')+
geom_vline(xintercept=1,linetype = "dashed",colour='#C0C0C0')+
scale_colour_manual(values = c("#00BFC4","grey","#F8766D"))+theme_bw()+
theme(panel.grid.major = element_blank(),panel.grid.minor = element_blank())
image.png
