plot.mevd {mevd}R Documentation

Plot return levels of MEVD fit

Description

Plot return levels of an mevd object.

Usage

     ## S3 method for class 'mevd'
plot(x, q = c(2,5,10,20,50,100,200), 
     ci = FALSE, type=c("all","rl","qq"), ...)

Arguments

x

mevd object.

q

vector of return periods, q > 1.

ci

if ci = TRUE, confidence intervals will be computed.

type

if omitted a return level plot type = 'rl' and a qq-plot type = 'qq' are shown.

...

Further parameters may also be supplied as arguments. See e.g. plot.

Examples

w <- 0.7
c <- 48
data <- matrix(NA,365,50)
for(y in 1:50){
     data[,y] <- rweibull(365,shape=w,scale=c)
}

# mean number of wet days per year
threshold <- 20
n <- apply(data,2,function(x){
  length(which(x>threshold))
})

fit <- fmev(data,n)
plot(fit)


[Package mevd version 0.1-0 Index]