Ldiag: L-moment ratio diagram

In martindurocher/floodStat: FloodNet Regional Frequency Analysis

Description

Plot the L-skew and L-kurtosis on the l-moment ratio diagram.

Usage

Ldiag(x, pos = "bottomright", xlim = NULL, ylim = NULL, ...)

Arguments

x	Matrix (n x 2) containing the L-skew and L-kurtosis respectively.
pos	Position of the legend. See legend.
...	Other ploting parameters. See par.

Examples

library(lmomco)

## Replicate 20 times the simulation of a GEV distribution
Fsim <- function() rAmax(100,c(100,20,-.1),'gev')
x <- replicate(20,Fsim())

## Compute the L-ratios and display the diagram
FLcoef <- function(z) lmoms(z)$ratio[3:4]
ll <- apply(x,2, FLcoef )

Ldiag(t(ll), pch = 16, cex = .7)
points(t(rowMeans(ll)), col = 2, pch = 17, cex = 2)

martindurocher/floodStat documentation built on May 31, 2019, 12:42 a.m.