我们正在尝试使用ggplot2绘制区域图,其中x轴上方的正区域是一种颜色,负区域是另一种颜色。
给定此数据集,我希望一个面积图在轴的每一侧上阴影不同的颜色。
我可以看到一种将数据集划分为两个子集的方法,一个为正数,其中所有负值均为零,一个为负数,所有正值均为零,然后将它们分别绘制在同一轴上,但似乎会有一个更多类似ggplot的方式来做到这一点。
张贴在this question上的解决方案未提供准确的结果(请参见下文)。
示例数据准确显示为条形图
通过此代码生成:
# create some fake data with zero-crossings
yvals=c(2,2,-1,2,2,2,0,-1,-2,2,-2)
test = data.frame(x=seq(1,length(yvals)),y=yvals)
# generate the bar plot
ggplot(data=test,aes(x=x,y=y))
+ geom_bar(data=test[test$y>0,],aes(y=y), fill="blue",stat="identity", width=.5)
+ geom_bar(data=test[test$y<0,],aes(y=y), fill="red",stat="identity", width=.5)
RLE方法不通用
当应用到我们的数据集时,在the other question上提出的RLE方法会产生与零交叉相关的伪像:
由以下代码生成(不使用):
# set up grouping function
rle.grp <- function(x) {
xx <- rle(x)
xx$values = seq_along(xx$values)
inverse.rle(xx) }
# generate ribbon plot
ggplot(test, aes(x=x,y=y,group = factor(rle.grp(sign(y))))) +
geom_ribbon(aes(ymax = pmax(0,y),ymin = pmin(0,y),
fill = factor(sign(y), levels = c(-1,0,1), labels = c('-','0','+'))))
+ scale_fill_brewer(name = 'sign', palette = 'RdBu')
参见@baptiste和Kohske建议的以下最终答案。
最佳答案
根据@baptiste的评论(已删除),我想这是最好的答案。它基于Kohske的this post。它将在零交叉处向数据集添加新的x-y对,并生成以下图:
# create some fake data with zero-crossings
yvals = c(2,2,-1,2,2,2,0,-1,-2,2,-2)
d = data.frame(x=seq(1,length(yvals)),y=yvals)
rx <- do.call("rbind",
sapply(1:(nrow(d)-1), function(i){
f <- lm(x~y, d[i:(i+1),])
if (f$qr$rank < 2) return(NULL)
r <- predict(f, newdata=data.frame(y=0))
if(d[i,]$x < r & r < d[i+1,]$x)
return(data.frame(x=r,y=0))
else return(NULL)
}))
d2 <- rbind(d,rx)
ggplot(d2,aes(x,y)) + geom_area(data=subset(d2, y<=0), fill="pink")
+ geom_area(data=subset(d2, y>=0), fill="lightblue") + geom_point()
生成以下输出:
关于r - 基于值的曲线下的填充面积,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/17959817/