#
#Example: find the percentage points for an SL distribution
#    with mean 1, standard deviation 1, skewness 3
#    corresponding to observed values 1, 2, 3, 4, 5
ans <- FitJohnsonDistribution(1, 1, 3, -1)
GAMMA <- ans["GAMMA"]
DELTA <- ans["DELTA"]
XLAM <- ans["XLAM"]
XI <- ans["XI"]
ITYPE <- 1
y <- 1:5
Z <- zJohnsonDistribution(y, ITYPE, GAMMA, DELTA, XLAM, XI)
pnorm(Z)
Run the code above in your browser using DataLab