X1 <- rcauchy(20)
lmomcau( parcau( TLmoms(X1, trim=1) ) )
alpha <- 30
tlmr <- theoTLmoms(vec2par(c(100, alpha), type="cau"), nmom=6, trim=1)
print( c(tlmr$lambdas[2] / alpha, tlmr$ratios[c(4,6)]), 8 )
Run the code above in your browser using DataLab