xx<-rnorm(500)
parms<-JohnsonFit(xx)
sJohnson(parms)
plot(function(xx)dJohnson(xx,parms),-2,2)
pJohnson(1,parms)
parms2<-JohnsonFit(rexp(50))
qJohnson(p=0.5,list(parms,parms2))
## JohnsonFit with moments=TRUE is not always possible, and even
## when possible, may produce odd results.
## parms<-JohnsonFit(x,moment=TRUE)
sJohnson(parms)
# Fit illustration
data(cars)
xx<-cars$speed
parms<-JohnsonFit(xx)
hist(xx,freq=FALSE)
plot(function(x)dJohnson(x,parms),0,25,add=TRUE)
Run the code above in your browser using DataCamp Workspace