lmr <- lmoms(rnorm(20))
para <- lmom2par(lmr,type="nor")
# The lmom2par() calls will error if trim != 1.
X <- rcauchy(20)
cauchy <- lmom2par(TLmoms(X, trim=1), type="cau")
slash <- lmom2par(TLmoms(X, trim=1), type="sla")
if (FALSE) {
plot(pp(X), sort(X), xlab="PROBABILITY", ylab="CAUCHY")
lines(nonexceeds(), par2qua(nonexceeds(), cauchy))
lines(nonexceeds(), par2qua(nonexceeds(), slash), col=2)
}
Run the code above in your browser using DataLab