#
#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