if (FALSE) {
par1 <- vec2par(c(-100, 30), type="sla")
X <- rlmomco(500, par1)
lmr <- TLmoms(X, trim=1)
par2 <- parsla(lmr)
F <- seq(0.001,.999, by=0.001)
plot(qnorm(pp(X)), sort(X), pch=21, col=8,
xlab="STANDARD NORMAL VARIATE",
ylab="QUANTILE")
lines(qnorm(F), quasla(F, par1), lwd=3)
lines(qnorm(F), quasla(F, par2), col=2)
}
Run the code above in your browser using DataLab