How to generate distributions given, mean, SD, skew and kurtosis in R?

by

There is a Johnson distribution in the SuppDists package. Johnson will give you a distribution that matches either moments or quantiles. Others comments are correct that 4 moments does not a distribution make. But Johnson will certainly try.

Here’s an example of fitting a Johnson to some sample data:

require(SuppDists)

## make a weird dist with Kurtosis and Skew
a <- rnorm( 5000, 0, 2 )
b <- rnorm( 1000, -2, 4 )
c <- rnorm( 3000,  4, 4 )
babyGotKurtosis <- c( a, b, c )
hist( babyGotKurtosis , freq=FALSE)

## Fit a Johnson distribution to the data
## TODO: Insert Johnson joke here
parms<-JohnsonFit(babyGotKurtosis, moment="find")

## Print out the parameters 
sJohnson(parms)

## add the Johnson function to the histogram
plot(function(x)dJohnson(x,parms), -20, 20, add=TRUE, col="red")

The final plot looks like this:

enter image description here

You can see a bit of the issue that others point out about how 4 moments do not fully capture a distribution.

Good luck!

EDIT
As Hadley pointed out in the comments, the Johnson fit looks off. I did a quick test and fit the Johnson distribution using moment="quant" which fits the Johnson distribution using 5 quantiles instead of the 4 moments. The results look much better:

parms<-JohnsonFit(babyGotKurtosis, moment="quant")
plot(function(x)dJohnson(x,parms), -20, 20, add=TRUE, col="red")

Which produces the following:

enter image description here

Anyone have any ideas why Johnson seems biased when fit using moments?

More Related Contents:

  1. In R, how to turn characters (grades) into a number and put in a separate column
  2. How to find the statistical mode?
  3. Workflow for statistical analysis and report writing
  4. How to make execution pause, sleep, wait for X seconds in R?
  5. R tick data : merging date and time into a single object
  6. Why am I getting “algorithm did not converge” and “fitted prob numerically 0 or 1” warnings with glm?
  7. Linear Regression loop for each independent variable individually against dependent
  8. Function to calculate R2 (R-squared) in R
  9. Changing whisker definition in geom_boxplot
  10. Geometric Mean: is there a built-in?
  11. Using cbind on an arbitrarily long list of objects
  12. Plotting a 3D surface plot with contour map overlay, using R
  13. Exact number of bins in Histogram in R
  14. Interpretation of ordered and non-ordered factors, vs. numerical predictors in model summary
  15. Standard Deviation in R Seems to be Returning the Wrong Answer – Am I Doing Something Wrong?
  16. Stepwise regression using p-values to drop variables with nonsignificant p-values
  17. predict.lm() in a loop. warning: prediction from a rank-deficient fit may be misleading
  18. Add error bars to show standard deviation on a plot in R
  19. Calculating percentile of dataset column
  20. How to shade a region under a curve using ggplot2
  21. Efficient calculation of matrix cumulative standard deviation in r
  22. How to use the ‘sweep’ function
  23. R optimization with equality and inequality constraints
  24. What does the period mean when used with ~ (in a formula)?
  25. Controlling number of decimal digits in print output in R
  26. ggplot2 change axis limits for each individual facet panel
  27. ignore NA in dplyr row sum
  28. How do I change the background color of a plot made with ggplot2
  29. How to split a data frame into multiple dataframes with each two columns as a new dataframe?
  30. R – How to make a click on webpage using rvest or rcurl

Leave a Comment