mantao/Plots.jl
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

fix-2dhistogram-bins

Browse Source
This commit is contained in:
Michael K. Borregaard 2017-09-05 14:01:59 +02:00
parent d5e5d81ee2
commit d6e5b57a08
1 changed files with 7 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
src/recipes.jl
View File

@ -506,7 +506,8 @@ function _auto_binning_nbins{N}(vs::NTuple{N,AbstractVector}, dim::Integer; mode
n_samples = length(linearindices(first(vs))) n_samples = length(linearindices(first(vs)))
# Estimator for number of samples in one row/column of bins along each axis: # Estimator for number of samples in one row/column of bins along each axis:
n = max(1, n_samples^(1/N)) nd = n_samples^(1/(2+N))
v = vs[dim] v = vs[dim]
@ -515,15 +516,15 @@ function _auto_binning_nbins{N}(vs::NTuple{N,AbstractVector}, dim::Integer; mode
end end
if mode == :sqrt # Square-root choice if mode == :sqrt # Square-root choice
_cl(sqrt(n)) _cl(sqrt(n_samples))
elseif mode == :sturges # Sturges' formula elseif mode == :sturges # Sturges' formula
_cl(log2(n)) + 1 _cl(log2(n_samples)) + 1
elseif mode == :rice # Rice Rule elseif mode == :rice # Rice Rule
_cl(2 * n^(1/3)) _cl(2 * nd)
elseif mode == :scott # Scott's normal reference rule elseif mode == :scott # Scott's normal reference rule
_cl(_span(v) / (3.5 * std(v) / n^(1/3))) _cl(_span(v) / (3.5 * std(v) / nd))
elseif mode == :fd # Freedman–Diaconis rule elseif mode == :fd # Freedman–Diaconis rule
_cl(_span(v) / (2 * _iqr(v) / n^(1/3))) _cl(_span(v) / (2 * _iqr(v) / nd))
elseif mode == :wand elseif mode == :wand
wand_edges(v) # this makes this function not type stable, but the type instability does not propagate wand_edges(v) # this makes this function not type stable, but the type instability does not propagate
else else