evir (version 1.7-4)

plot.potd: Plot Fitted POT Model

Description

The plot method plot.potd provides seven different plots for assessing fitted POT model. The user selects the plot type from a menu. See the examples below.

Usage

# S3 method for potd
plot(x, …)

Arguments

x

an object returned by the function pot

…

other graphics parameters

Value

If plot 1 or 2 from the GPD plots is selected as the final plot (i.e. option 8 is selected, followed by option 1 or 2), a list object containing details of the plot is returned invisibly. This object should be used as the first argument of gpd.q or gpd.sfall to add quantile estimates or expected shortfall estimates to the plot.

Details

Plot 1 displays the exceedance process of the chosen threshold. Plots 2-4 assess the Poisson nature of the exceedance process by looking at the scaled gaps between exceedances, which should be iid unit exponentially distributed. Plots 5-6 assess the GPD nature of the excesses by looking at suitably defined residuals, which should again be iid unit exponentially distributed. Option 8 allows the user to call GPD plotting functions.

See Also

gpd, pot, plot.gpd

Examples

# NOT RUN {
data(danish)
out <- pot(danish,10) 
# }
# NOT RUN {
plot(out)
# }
# NOT RUN {
  
# }
# NOT RUN {
Make a plot selection (or 0 to exit):
# }
# NOT RUN {
1: plot: Point Process of Exceedances
# }
# NOT RUN {
2: plot: Scatterplot of Gaps
# }
# NOT RUN {
3: plot: Qplot of Gaps
# }
# NOT RUN {
4: plot: ACF of Gaps
# }
# NOT RUN {
5: plot: Scatterplot of Residuals
# }
# NOT RUN {
6: plot: Qplot of Residuals
# }
# NOT RUN {
7: plot: ACF of Residuals
# }
# NOT RUN {
8: plot: Go to GPD Plots
# }