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
On Dec 4, 2006, at 12:05 PM, Eric Thompson wrote:
> 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 <[EMAIL PROTECTED]> 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:
>> [EMAIL PROTECTED]
>> [EMAIL PROTECTED]
>> [EMAIL PROTECTED]
>>
>> [[alternative HTML version deleted]]
>>
>> __
>> R-help@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@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@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.
