Help for L-moment Ratio Diagram

Amina,
If you can follow this code, you will see that the black dots scatter  
around an L-skew of zero and L-kurtosis of 0.122602, which are the  
theoretical L-moments of the standard normal distribution.

library(lmomco)
lmrdia <- lmrdia()
plotlmrdia(lmrdia,ylim=c(-.2,.2),xlim=c(-.2,.2))
lms <- lmoms(rnorm(500))
points(lms$ratios[3],lms$ratios[4],pch=16)
lms <- lmoms(rnorm(500))
points(lms$ratios[3],lms$ratios[4],pch=16)
lms <- lmoms(rnorm(500))
points(lms$ratios[3],lms$ratios[4],pch=16)
lms <- lmoms(rnorm(500))
points(lms$ratios[3],lms$ratios[4],pch=16)

Further interpretion of the diagram is well documented in Hosking's  
1990  paper and the references shown on the help page.

William

I think there are a few packages for doing this. I have used "lmomco":
the function lmom.ub() will calculate the sample lmoments, and
lmrdia() gives theoretical lmoments for different distributions.

Hope this is helpful.

Eric


On 12/4/06, amna khan <amnakhan493 at gmail.com> wrote:

Respected Sir

I have to select a probability distribution using L-moment Ratio  
Diagram. I
am not understanding how to plot sample TAU3 and TUA4 on L-moment  
Ratio
Diagram.
Please Guide me
Best Regards

--
AMINA SHAHZADI
Department of Statistics
GC University Lahore, Pakistan.
Email:
amnakhan493 at gmail.com
amna_989 at hotmail.com
amna_989 at yahoo.com

        [[alternative HTML version deleted]]

______________________________________________
R-help at stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting- 
guide.html
and provide commented, minimal, self-contained, reproducible code.

______________________________________________
R-help at stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting- 
guide.html
and provide commented, minimal, self-contained, reproducible code.