我不喜欢Gadfly在绘图时选择轴限制的方式,例如我制作的一个图只有画布中心四分之一的数据。 MWE可以是:
plot(x=[2.9,8.01],y=[-0.01,0.81])
然后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]))
即。
Gadfly中是否已存在类似的内容?鉴于我努力寻找Guide.[xy]ticks
,我期待我需要为此编写一些代码......
答案 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代码中看到如何做到这一点。