# Sample L-moments of Ozone from the airquality data
data(airquality)
lmom <- samlmu(airquality$Ozone,nmom=6)
distrib <- c("exp","gam","gev","glo","gpa","gno","gum","kap","ln3","nor",
"pe3","wak","wei")
para_list <- pel(distrib=distrib,lmom=lmom)
f <- (1:10)/10
qua_list <- qua(para=para_list,f=f)
qua_gam <- qua(para=para_list$gam,f=f)
qua_gam2 <- qua(para=as.vector(para_list$gam),f=f,distrib="gam")
if (any(qua_gam!=qua_gam2,na.rm=TRUE)) stop("Any possible errors after 0.6.3 package updates!")
Run the code above in your browser using DataLab