ggplot2 密度图 分段设置线段颜色

最终结果

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