我有一个数据框,其中包含不同的组(“label”列)。对于每个标签,我想绘制一个从bootstrapping获得的null分布(值在“null”列中),并在顶部绘制真实的统计量(值在“sc”列中)。理想情况下,我希望统计量后面的区域具有不同的颜色,以标记这是我的p值。这是否可以用stat_density_ridges
来完成?
下面是一个R代码示例:
library(ggplot2)
library(tidyverse)
library(ggridges)
df <- data.frame()
for (label in LETTERS) {
mean=rnorm(1,0.5,0.2)
null = rnorm(1000,mean,0.1);
sc = rnorm(1,0.5,0.2)
df <- rbind(df, data.frame(label=label, null=null, sc=sc))
}
df <- df %>%
mutate(label=as.factor(label))
ggplot(df, aes(x = null, y = label)) +
stat_density_ridges(scale=1.2,alpha = 1, size=1)+
scale_x_continuous(limits=c(0,1),breaks=seq(0,1,0.2)) +
geom_segment(aes(x=sc, xend=sc, y=as.numeric(label)-0.1, yend=as.numeric(label)+0.5), size=1) +
coord_flip()
结果图如下:
但理想情况下,我希望每个脊更像这样:
随着sc值后颜色的变化。这是可能的吗?谢谢:)
1条答案
按热度按时间6l7fqoea1#
您可以使用
fill
和..x..
在图的固定x值处创建不同的颜色。因此阴影区域在所有图中都是相同的。您可以通过使用ggplot_build
和一个单独的 Dataframe 来修改这一点,该 Dataframe 具有作为阈值的p_values。因此,使用这些阈值,您可以有条件地更改图层中的color
。下面是一些可复制的代码:创建于2023-03-28带有reprex v2.0.2
正如你在最新的图中看到的,阴影区域现在是根据每组 Dataframe 的sc值。