gev.d.diag {IDF}

Diagnostic Plots for d-gev Models

Description

Produces diagnostic plots for d-gev models using the output of the function gev.d.fit. Values for different durations can be plotted in different colors of with different symbols.

Usage

gev.d.diag(
  fit,
  subset = NULL,
  cols = NULL,
  pch = NULL,
  which = "both",
  mfrow = c(1, 2),
  legend = TRUE,
  title = c("Residual Probability Plot", "Residual Quantile Plot"),
  emp.lab = "Empirical",
  mod.lab = "Model",
  ci = FALSE,
  ...
)

Arguments

fit

object returned by gev.d.fit

subset

an optional vector specifying a subset of observations to be used in the plot

cols

optional either one value or vector of same length as unique(fit$ds) to specify the colors of plotting points. The default uses the rainbow function.

pch

optional either one value or vector of same length as unique(fit$ds) containing integers or symbols to specify the plotting points.

which

string containing 'both', 'pp' or 'qq' to specify, which plots should be produced.

mfrow

vector specifying layout of plots. If both plots should be produced separately, set to c(1,1).

legend

logical indicating if legends should be plotted

title

character vector of length 2, giving the titles for the pp- and the qq-plot

emp.lab, mod.lab

character string containing names for empirical and model axis

ci

logical indicating whether 0.95 confidence intervals should be plotted

...

additional parameters passed on to the plotting function

Examples

data('example',package ='IDF')

fit <- gev.d.fit(xdat=example$dat,ds = example$d,ydat=as.matrix(example[,c('cov1','cov2')])
                 ,mutl=c(1,2),sigma0l=1)
# diagnostic plots for complete data
gev.d.diag(fit,pch=1,ci = TRUE)
# diagnostic plots for subset of data (e.g. one station)
gev.d.diag(fit,subset = example$cov1==1,pch=1,ci = TRUE)