R风格轴与Gadfly.jl对齐

时间:2015-03-09 14:01:24

标签: julia gadfly

我不喜欢Gadfly在绘图时选择轴限制的方式,例如我制作的一个图只有画布中心四分之一的数据。 MWE可以是:

plot(x=[2.9,8.01],y=[-0.01,0.81])

screenshot of above plot然后Gadfly为y轴选择[0,10]和[-0.5,1]的x轴范围,这对我来说似乎太宽了。这里的值显然是弥补的,但基本上是我的真实数据的边界框。

我更不希望拥有所有空白区域,例如R的默认4%模式(即par(xaxs='r',yaxs='r'))。我可以通过以下方式在Gadfly中获得类似的东西:

plot(x=[2.9,8.01],y=[-0.01,0.81]
  Guide.xticks(ticks=[3:8]),
  Guide.yticks(ticks=[0:0.2:0.8]))

即。 screenshot of second example

Gadfly中是否已存在类似的内容?鉴于我努力寻找Guide.[xy]ticks,我期待我需要为此编写一些代码......

指针赞赏!

1 个答案:

答案 0 :(得分:1)

作为一种解决方法,我有一个Heckbert's 1990's Graphics Gems code的更改版本,用于在给定的最小值/最大值内生成节拍。在(我的天真)朱莉娅看起来像这样:

# find a "nice" number approximately equal to x.
# round the number if round=true, take the ceiling if round=false
function nicenum{T<:FloatingPoint}(x::T, round::Bool)
    ex::T = 10^floor(log10(x))
    f::T = x/ex
    convert(T, if round
        if     f < 1.5; 1.
        elseif f < 3.;  2.
        elseif f < 7.;  5.
        else;           10.
        end
    else
        if     f <= 1.; 1.
        elseif f <= 2.; 2.
        elseif f <= 5.; 5.
        else;           10.
        end
    end)*ex
end

function pretty_inner{T<:FloatingPoint}(min::T, max::T, ntick::Int)
    if max < min
        error("min must be less than max")
    end
    if min == max
        return min:one(T):min
    end
    delta = nicenum((max-min)/(ntick-1),false)
    gmin =  ceil(min/delta)*delta
    gmax = floor(max/delta)*delta
    # shift max a bit in case of rounding errors
    gmin:delta:(gmax+delta*eps(T))
end

可以用作:

plot(x=[2.9,8.01],y=[-0.01,0.81],
  Guide.xticks(ticks=[pretty_inner(2.9,8.01,7)]),
  Guide.yticks(ticks=[pretty_inner(-0.01,0.81,7)]))

并且会得到与我从R获得的结果相同的结果。

如果范围可以自动拉出会很棒,但我无法在现有的Gadfly代码中看到如何做到这一点。