Is there an equivalent to R qqnorm and qqplot functions in some julia package?
There’s a gist here for how to do it with Gadfly https://gist.github.com/kleinschmidt/7ce8cec988a84ce73ba2
It would be fairly easy to make a recipe for this for StatPlots as well - we should definitely have that.
This should work with Plots:
using Distributions, RecipesBase @recipe function f(h::QQPair) seriestype --> :scatter h.qx, h.qy end @userplot QQPlot @recipe f(h::QQPlot) = qqbuild(h.args, h.args) @userplot QQNorm @recipe f(h::QQNorm) = qqbuild(Normal(), h.args)
Use like this:
x = rand(Normal(), 100) y = rand(Cauchy(), 100) using Plots qqplot(x,y) qqnorm(y)
And now there’s a PR on StatPlots:
@harven there’s a request for comment on the pull request, if you like.
That’s great, thank you! I will have a look at the RFC soon.
What’s the current status? How can we use it?
The package StatsPlots has the qq plot function available.
# Estimate distribution parameters: mw = mean(fix_Acid) StdAbw = std(fix_Acid) # Create QQ plot: qqplot(Normal(mw, StdAbw), fix_Acid, title = "QQ-Plot gebundene Säuren", ylabel = "Gehalt geb. Säuren")