plot.mevd {mevd} | R Documentation |
Plot return levels of an mevd object.
## S3 method for class 'mevd' plot(x, q = c(2,5,10,20,50,100,200), ci = FALSE, type=c("all","rl","qq"), ...)
x |
mevd object. |
q |
vector of return periods, q > 1. |
ci |
if |
type |
if omitted a return level plot |
... |
Further parameters may also be supplied as arguments. See e.g. |
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)