PPmatG: Plotting positions of a matrix of data

In ilapros/ilaprosUtils: A Miscellaneous Collection of Functions I Tend to Use

Plotting positions of a matrix of data

Description

This can be useful when building DDF models or when comparing seasonal maxima of the same record. It takes each columns of a matrix and plots them as Gringorten plotting positions. The x-axis is shown using the Gumbel variate -log(-log(f)), with f the non exceedance probability If a distribution is given to fit.lines, it adds an estimated return level curve based on L-moments estimate

Usage

PPmatG(
  mat,
  cols = NULL,
  fit.lines = "gev",
  ff = seq(0.005, 0.995, by = 0.0025),
  yll = NULL,
  timeGrid = TRUE,
  vgrid = c(1.2, 2, 5, 10, 50, 100, 200),
  ...
)

Arguments

mat	data matrix, each column will be analysed separately.
cols	colour in which the data points of each column (and estimated line) should be displayed
fit.lines	a character string indicating the distribution to be used to draw the fitted line. Current options are "gev" (default), "glo", "gamma", "gpa" and "none" which results in no lines
ff	the non-exceedance probabilities for which the lines are estimated. It also affects the width of the x-axis. Default is seq(0.005,0.995,by=0.0025).
yll	ylimits (optional)
timeGrid	logical. If TRUE a time grid with return periods expressed in years corresponding to values in vgrid will be displayed.
vgrid	return periods expressed in years shown as a grid
...	other parameters passed into plot

Value

A plot with the plotting positions of each column and a matrix of the estimated parameters

Examples

x <- matrix(c(rgev(20, 10, 3,   -0.2), 
              rgev(20, 15, 4.5, -0.2), 
              rgev(20, 40, 6.5, -0.2)),byrow = FALSE, ncol=3)
PPmatG(x, cols = c(4,5,2), fit.lines = "gev")

ilapros/ilaprosUtils documentation built on April 6, 2023, 4:44 a.m.