evgam (version 1.0.2)

plot.geoevgam: Plot results of a geoevgam fit

Description

Plot results of a geoevgam fit

Usage

# S3 method for geoevgam
plot(
  x,
  type = c("threshold", "excess"),
  view = c("polar", "cartesian"),
  norm,
  margins = "original",
  ...
)

Arguments

x

An object of class geoevgam.

type

Character vector specifying plot types: "threshold", "excess", "gauge", or "qqplot".

view

Character specifying the coordinate system for gauge plots: "polar" or "cartesian".

norm

Numeric; the norm to use for the gauge plot. Defaults to the norm used in object.

margins

Character; margin type for QQ-plots. Defaults to "original".

...

Additional arguments passed to the underlying plot functions.

Examples

n <- 1e3
gauss_sim_data <- rmvnorm(n, numeric(2), matrix(c(1, .8, .8, 1), 2))
laplace_sim_data <- qlaplace(pnorm(gauss_sim_data))
polar1 <- polarise(t(laplace_sim_data), norm = 2)
data1 <- data.frame(radius = polar1$r, angle = polar1$phi[1, ])
fit1 <- geoevgam(data = data1)
plot(fit1)
plot(fit1, type = 'gauge', col = 2, lwd = 2, view = 'cartesian')
par(mfrow = c(2, 2))
plot(fit1, type = 'qqplot')
plot(fit1, type = 'qqplot2')
plot(fit1, type = 'qqplot', margins = 'uniform')
plot(fit1, type = 'qqplot2', margins = 'uniform')

# \donttest{
# example using L_1 norm
polar2 <- polarise(t(laplace_sim_data), norm = 1)
data2 <- data.frame(radius = polar2$r, angle = polar2$phi[1, ])
fit2 <- geoevgam(data = data2, norm = 1) 
# passing norm ensures plotted gauge is correct
plot(fit2, type = 'gauge', col = 2, lwd = 2, view = 'cartesian')
# }