Add points or lines to an L-moment ratio diagram
lmrdpoints adds points,
and lmrdlines adds connected line segments,
to an L-moment ratio diagram.
lmrdpoints(x, y=NULL, type="p", ...) lmrdlines(x, y=NULL, type="l", ...)
x |
Numeric vector of L-skewness values. |
y |
Numeric vector of L-kurtosis values. May be omitted: see “Details” below. |
type |
Character indicating the type of plotting.
Can be any valid value for the |
... |
Further arguments (graphics parameters),
passed to |
The functions lmrdpoints and lmrdlines are equivalent to
points and lines respectively,
except that if argument y is omitted, x is assumed to be
an object that contains both L-skewness and L-kurtosis values.
As in lmrd, it can be a vector with elements named
"t_3" and "t_4" (or "tau_3" and "tau_4"),
a matrix or data frame with columns named
"t_3" and "t_4" (or "tau_3" and "tau_4"),
or an object of class "regdata" (as defined in package lmomRFA).
J. R. M. Hosking jrmhosking@gmail.com
# Plot L-moment ratio diagram of Wind from the airquality data set data(airquality) lmrd(samlmu(airquality$Wind), xlim=c(-0.2, 0.2)) # Sample L-moments of each month's data ( lmom.monthly <- with(airquality, t(sapply(5:9, function(mo) samlmu(Wind[Month==mo])))) ) # Add the monthly values to the plot lmrdpoints(lmom.monthly, pch=19, col="blue") # Draw an L-moment ratio diagram and add a line for the # Weibull distribution lmrd(xaxs="i", yaxs="i", las=1) weimom <- sapply( seq(0, 0.9, by=0.01), function(tau3) lmrwei(pelwei(c(0,1,tau3)), nmom=4) ) lmrdlines(t(weimom), col='darkgreen', lwd=2)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.