2017-03-29 55 views
0

我写了一个函数,可以生成一个直方图两个垂直条,用于指示一个值的范围。我想修改这个函数,使得指定范围内的条纹是不同的颜色。R - ggplot直方图,基于范围的颜色

继承人我的功能和快速演示:

require(ggplot2) 
niceHist <- function(data, cutpoint1, cutpoint2, title = "Supply a title, genius") { 
    temp_dat = data.frame(Data = data, Col = 0) 
    temp_dat = temp_dat[! is.na(temp_dat$Data),] 
    temp_dat[temp_dat$Data >= cutpoint1 & temp_dat$Data <= cutpoint2,]$Col = 1 
    my_hist = qplot(data) + 
    geom_histogram(fill = "forestgreen") + 
    geom_vline(xintercept = cutpoint1) + 
    geom_vline(xintercept = cutpoint2) + 
    ggtitle(paste(title)) + 
    theme_minimal() + 
    theme(text = element_text(size = 16), axis.line.y = element_line(color = "black", size = 0.5), axis.line.x = element_line(color = "black", size = 0.5)) 
    my_hist 
} 
u = rnorm(100) 
c1 = mean(u) - sd(u) 
c2 = mean(u) + sd(u) 
niceHist(u, c1, c2) 

我见过similar question,因为我想保持原来的直方图的形状,其接受的解决方案是不是很适合我的需要。我也不希望更改垃圾箱的数量,并且如果可能的话,应用颜色差异,使得直方图中的单个条纹可以是两种颜色,如果垂直的黑色线条恰好平分它的话。

*我的主要目标是清楚地显示在提供的范围内捕捉了多少分布,而不改变分布的形状。 *因此,一个替代但不太理想的解决方案是简单地将背景与所提供的范围相关联。另外,我需要我的函数来返回一个ggplot对象,因为它偶尔会使用ggplot语法进一步修改。

UPDATE: 在评论的建议,我一直在使用scale_fill_gradientn尝试,但这个不起作用:

niceHist <- function(data, cutpoint1, cutpoint2, title = "Supply a title, genius") { 
     temp_dat = data.frame(Data = data, Col = 0) 
     temp_dat = temp_dat[! is.na(temp_dat$Data),] 
     temp_dat[temp_dat$Data >= cutpoint1 & temp_dat$Data <= cutpoint2,]$Col = 1 
     my_hist = qplot(data) + 
     geom_histogram() + 
     scale_fill_gradientn(colours = c("blue", "red", "red", "blue"), values = c(min(data, na.rm = TRUE), cutpoint1, cutpoint2, max(data, na.rm = TRUE))) + 
     geom_vline(xintercept = cutpoint1) + 
     geom_vline(xintercept = cutpoint2) + 
     ggtitle(paste(title)) + 
     theme_minimal() + 
     theme(text = element_text(size = 16), axis.line.y = element_line(color = "black", size = 0.5), axis.line.x = element_line(color = "black", size = 0.5)) 
     my_hist 
    } 
+1

你看过scale_fill_gradient – user2510479

回答

0

我的解决办法是建立计data.frame并添加系数来表示bin的区域。例如:

df <- ggplot_build(niceHist(u,c1,c2))$data[[1]] #recreate the count df 
require(dplyr) 
df <- df %>% mutate(col=cut(x,c(min(x)-0.001,c1,c2,max(x)+0.001))) 
ggplot(df,aes(x,count))+ geom_col(aes(fill=col)) + geom_vline(xintercept = c(c1,c2))