Another particle physicist who likes step histograms here.

I’ve googled and came across this:

https://github.com/JuliaStats/StatsBase.jl/pull/243

I very much like the plots shown here, but this doesn’t seem to work (and I think this is unrelated to the fact that this pull request isn’t in a released version). Plots.jl does weird things when trying to plot a StatsBase.Histogram.

It also does weird things when doing things like `histogram(randn(1000), bins=100, seriestype=:steps)`

So my question is: What’s the recommended way to plot a StatsBase.Histogram?

There is a companion PR (https://github.com/JuliaPlots/Plots.jl/pull/713) on Plots that depends on this PR (EDIT in fact that PR was already merged onto the `new_histogram`

branch). If you checkout the linked PR and the `new_histogram`

branch in Plots it should work. It’d be cool if you could test it.

Until this gets merged, there is functionality in StatPlots to plot StatsBase.Histogram, but you won’t have the `steps`

seriestype.

2 Likes

Checked out the JuliaPlots branch, and applied the pull request to StatsBase.

Non-uniform binning is no longer broken. (try, for example to histogram a normal distribution with the bins [-2.5, -2, -0.5, 0.1, 1.9, 2.5]. Using `histogram`

, the bars overlap.

With the new version I can just do

```
plot(h)```
or
`plot(h, seriestype=:steps)` and everything looks as expected.
That's as easy as it should be. Admittedly this is not a particularly thorough test, but I hope this gets merged soon.
```

2 Likes