gev.d.diag: Diagnostic Plots for d-gev Models

In IDF: Estimation and Plotting of IDF Curves

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)

IDF documentation built on July 20, 2022, 5:06 p.m.