powered by
Creates a plot showing how the estimate of shape varies with threshold or number of extremes.
shape(data, models = 30, start = 15, end = 500, reverse = TRUE, ci = 0.95, auto.scale = TRUE, labels = TRUE, …)
numeric vector of data
number of consecutive gpd models to be fitted
lowest number of exceedances to be considered
maximum number of exceedances to be considered
should plot be by increasing threshold (TRUE) or number of extremes (FALSE)
TRUE
FALSE
probability for asymptotic confidence band; for no confidence band set to zero
whether or not plot should be automatically scaled; if not, xlim and ylim graphical parameters may be entered
xlim
ylim
whether or not axes should be labelled
other graphics parameters
A table of results is returned invisibly.
For every model gpd is called. Evaluation may be slow.
gpd
gpd, plot.gpd, hill
plot.gpd
hill
# NOT RUN { data(danish) # } # NOT RUN { shape(danish) # } # NOT RUN { # Shape plot of heavy-tailed Danish fire insurance data # }
Run the code above in your browser using DataLab