要让ggplot
在对数刻度上正确地绘制小的中断,我必须这样做:
faceplant1 <- function(x) {
return (c(x[1]*10^.25, x[2]/10^.25))
}
faceplant2 <- function(x) {
return (rep(seq(1,9),5)*rep(10^seq(-6,-2), each=9))
}
ggplot(mydata, aes(x=myseries)) +
geom_density() +
scale_x_log10(limits=c(1e-6, 1e-1),
breaks=10^seq(-6,-1),
minor_breaks=trans_breaks(faceplant1, faceplant2, n=45))
有没有更简单的方法来实现这个目标?
最终结果如下:
答案 0 :(得分:10)
以下是我对该问题的解决方案:
library(ggplot2)
log10_minor_break = function (...){
function(x) {
minx = floor(min(log10(x), na.rm=T))-1;
maxx = ceiling(max(log10(x), na.rm=T))+1;
n_major = maxx-minx+1;
major_breaks = seq(minx, maxx, by=1)
minor_breaks =
rep(log10(seq(1, 9, by=1)), times = n_major)+
rep(major_breaks, each = 9)
return(10^(minor_breaks))
}
}
mydata = data.frame(myseries = 10^(rnorm(1e4, mean=0, sd=0.5)))
myplot =
ggplot(mydata, aes(x=myseries))+
geom_density()+
scale_x_log10(minor_breaks=log10_minor_break())+
theme(panel.grid.major.x = element_line(size=1.0),
panel.grid.minor.x = element_line(size=2))
myplot
它与您已经完成的工作非常相似,但通常适用。 还有一个小改进:在你的例子中,它会将次要突破扩展到1e-6以下和1e-1以下。
我从查看函数trans_break
开始,并将其缩减为最基本的元素。
还值得考虑annotation_logticks()函数:
myplot+annotation_logticks(side="b")
答案 1 :(得分:4)
多年以后:我已经将@ gvrocha的答案改编为(甚至更多)通用。
log_breaks = function(maj, radix=10) {
function(x) {
minx = floor(min(logb(x,radix), na.rm=T)) - 1
maxx = ceiling(max(logb(x,radix), na.rm=T)) + 1
n_major = maxx - minx + 1
major_breaks = seq(minx, maxx, by=1)
if (maj) {
breaks = major_breaks
} else {
steps = logb(1:(radix-1),radix)
breaks = rep(steps, times=n_major) +
rep(major_breaks, each=radix-1)
}
radix^breaks
}
}
scale_x_log_eng = function(..., radix=10) {
scale_x_continuous(...,
trans=log_trans(radix),
breaks=log_breaks(TRUE, radix),
minor_breaks=log_breaks(FALSE, radix))
}
scale_y_log_eng = function(..., radix=10) {
scale_y_continuous(...,
trans=log_trans(radix),
breaks=log_breaks(TRUE, radix),
minor_breaks=log_breaks(FALSE, radix))
}
然后使用更优雅:
ggplot(...) +
geom_line() +
scale_x_log_eng() +
scale_y_log_eng()
答案 2 :(得分:0)
我想说的最简单的解决方案是使用注解_logticks()并将短,中和长参数布置为 mic 网格线,然后删除网格线并使用适当的比例来绘制数据。例如,
annotation_logticks(short = unit(1, "npc"), mid = unit(1, "npc"), long = unit(1,"npc"), ...) +
theme(panel.grid = element_blank()) +
scale_x_log10(breaks = 10^(seq(-6,-1,1)), limits = c(10^-6, 0.1),...)