我面临以下问题:一些极端值支配着我的
geom_raster
图的颜色比例尺。一个示例可能更清楚(请注意,此示例仅适用于最近的ggplot2版本,我使用0.9.2.1):library(ggplot2)
library(reshape)
theme_set(theme_bw())
m_small_sd = melt(matrix(rnorm(10000), 100, 100))
m_big_sd = melt(matrix(rnorm(100, sd = 10), 10, 10))
new_xy = m_small_sd[sample(nrow(m_small_sd), nrow(m_big_sd)), c("X1","X2")]
m_big_sd[c("X1","X2")] = new_xy
m = data.frame(rbind(m_small_sd, m_big_sd))
names(m) = c("x", "y", "fill")
ggplot(m, aes_auto(m)) + geom_raster() + scale_fill_gradient2()
现在我通过将超过某个分位数的值设置为该分位数来解决这个问题:
qn = quantile(m$fill, c(0.01, 0.99), na.rm = TRUE)
m = within(m, { fill = ifelse(fill < qn[1], qn[1], fill)
fill = ifelse(fill > qn[2], qn[2], fill)})
这不是最优解决方案。我想要做的是将颜色非线性映射到值范围,即在具有更多观测值的区域中存在更多颜色。在spplot
中,我可以使用classInt
包中的classIntervals
计算适当的类边界:
library(sp)
library(classInt)
gridded(m) = ~x+y
col = c("#EDF8B1", "#C7E9B4", "#7FCDBB", "#41B6C4",
"#1D91C0", "#225EA8", "#0C2C84", "#5A005A")
at = classIntervals(m$fill, n = length(col) + 1)$brks
spplot(m, at = at, col.regions = col)
spplot
中那样硬编码颜色到类间的映射。我可以转换fill
轴,但由于fill
变量中存在负值,这种方法行不通。因此,我的问题是:是否有使用ggplot2解决此问题的方法?