Nothing
#############################################################################
# Copyright (c) 2014 Mathieu Ribatet
#
# This program is free software; you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation; either version 2 of the License, or
# (at your option) any later version.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with this program; if not, write to the
# Free Software Foundation, Inc.,
# 59 Temple Place, Suite 330, Boston, MA 02111-1307, USA
#
#############################################################################
## This file contains several functions to plot Peaks Over a Threshold.
plot.uvpot <- function(x, npy, main, which = 1:4,
ask = nb.fig < length(which) &&
dev.interactive(),ci = TRUE, ...){
if(!inherits(x, "uvpot"))
stop("Use only with 'uvpot' objects")
if (!is.numeric(which) || any(which < 1) || any(which > 4))
stop("`which' must be in 1:4")
show <- rep(FALSE, 4)
show[which] <- TRUE
nb.fig <- prod(par("mfcol"))
if (ask) {
op <- par(ask = TRUE)
on.exit(par(op))
}
if (show[1]) {
pp.uvpot(x, ci = ci, main, xlim = c(0, 1), ylim = c(0, 1),
...)
}
if (show[2]) {
qq.uvpot(x, ci = ci, main, ...)
}
if (show[3]) {
dens.uvpot(x, main, ...)
}
if (show[4]) {
retlev.uvpot(x, npy, main, ...)
}
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.