ANALISI BLOCK MAXIMA PER RX1 STAGIONALI

Serie stagionali ottenute dalle serie di valori giornalieri SCIA.

rm(list = objects())

library("evd")

# elimina warnings
options(warn = -1)


# funzioni per il calcolo del massimo annuale (RX1) block-maxima approach

analisiGEV <- function(col) {

    # data.frame per il livelli di ritorno
    idColumn <- rep(NA, length = length(tempi.ritorno))
    return.period <- rep(NA, length = length(tempi.ritorno))
    return.level <- rep(NA, length = length(tempi.ritorno))
    confidence.delta.lower <- rep(NA, length = length(tempi.ritorno))
    confidence.delta.upper <- rep(NA, length = length(tempi.ritorno))
    confidence.profile.lower <- rep(NA, length = length(tempi.ritorno))
    confidence.profile.upper <- rep(NA, length = length(tempi.ritorno))

    precReturn <- data.frame(idColumn, return.period, return.level, confidence.delta.lower, 
        confidence.delta.upper, confidence.profile.lower, confidence.profile.upper)

    ID <- names(dati)[col]

    cat(sprintf("Stazione: %s", ID))

    # associamo i dati
    x <- dati[, c(1, col)]
    names(x) <- c("yy", "prcp")

    # subset dei dati sulla base di quanto riportato nel file annoI
    file.annoI <- list.files(pattern = "^annoI.+txt$")
    if (length(file.annoI) != 1) 
        stop("Impossibile: ci deve essere un solo file")

    # cerchiamo anno inizio
    info.anni <- read.csv(file = file.annoI[1], sep = ";", head = T)
    index.id <- which(info.anni$idstaz == ID)

    if (!length(index.id)) 
        stop(paste0(ID, "stazione non trovata in file annoI"))
    YYYY <- info.anni[index.id, ]$annoI

    # YYYY: prima di questo anno dobbiamo eliminare la serie
    if (YYYY != 1951) 
        x <- subset(x, x >= YYYY)


    plot(x$yy, x$prcp, type = "l", xlab = "Years", ylab = "mm", main = paste0("Station -- ", 
        ID))
    points(x$yy, x$prcp, pch = 21, bg = "red")

    index <- which(is.na(x$prcp))

    # con pacchetto evd non sarebbe necessario
    if (length(index)) 
        x <- x[-index, ]


    annual.fit <- tryCatch({
        fgev(x$prcp)
    }, error = function(err) {
        print(paste(ID, "gev", sep = " "))
        print(err)
        precReturn <- data.frame(ID = ID, conf.level = NA, return.level = NA, 
            return.period = NA, confidence.delta.lower = NA, confidence.delta.upper = NA, 
            rl.confidence.profile.lower = NA, rl.confidence.profile.upperNA)
        my.shape <- c(NA, NA, NA)
        return(list(precReturn, my.shape))
    })


    plot(annual.fit, which = 1)
    plot(annual.fit, which = 2)
    plot(annual.fit, which = 3)
    plot(annual.fit, which = 4)
    parametri <- round(annual.fit$param, 2)

    cat(sprintf("LOcation Scale and Shape; PACKAGE evd: %s\n", parametri))

    lapply(1:length(tempi.ritorno), function(i) {

        print(sprintf("Calcolo livello di ritorno per tempo: %s anni", tempi.ritorno[i]))
        stima.evd <- fgev(x$prcp, prob = 1/(tempi.ritorno[i]))
        livello.evd <- round(stima.evd$param["quantile"], 2)
        cat(sprintf("RLevel; PACKAGE evd: %s\n", livello.evd))
        profile.object <- evd:::profile.evd(stima.evd, which = "quantile", conf = 0.95)
        plot(profile.object)
        profile.int <- confint(profile.object, level = 0.95, parm = "quantile")
        # intervallo mediante profile likelihood
        profLow <- profile.int[, "lower"]
        profUp <- profile.int[, "upper"]

        # intervalli di confidenza assumendo ipotesi normalità
        profLowW <- confint(stima.evd)[, 1][1]
        profUpW <- confint(stima.evd)[, 2][1]


        precReturn$idColumn[i] <<- ID
        precReturn$return.period[i] <<- tempi.ritorno[i]
        precReturn$return.level[i] <<- livello.evd
        precReturn$confidence.delta.lower[i] <<- profLowW
        precReturn$confidence.delta.upper[i] <<- profUpW

        precReturn$confidence.profile.lower[i] <<- profLow
        precReturn$confidence.profile.upper[i] <<- profUp

    })
    ## fine lapply plot dei tempi di ritorno e intervalli di confidenza (in
    ## rosso) con profile likelihood
    lineaL <- approx(seq(10, 100, 10), precReturn$confidence.profile.lower)
    lineaU <- approx(seq(10, 100, 10), precReturn$confidence.profile.upper)

    plot(annual.fit, which = 4)
    lines(lineaL$x, lineaL$y, col = "red")
    lines(lineaU$x, lineaU$y, col = "red")
    ########################### fine parte su ismev

    # delta intervals for shape
    shapeFi.delta <- confint(annual.fit, parm = "shape")

    # parametro shape e intervalli di confidenza mediante profile likelihood
    profile.object <- profile(annual.fit, which = "shape")
    shapeFi.prof <- confint(profile.object)

    my.shape <- c(ID = as.numeric(ID), xi = parametri[3], xi.delta.dwn = shapeFi.delta[1], 
        xi.delta.up = shapeFi.delta[2], xi.prof.dwn = shapeFi.prof[1], xi.prof.up = shapeFi.prof[2])

    list(precReturn, my.shape)


}  #fine analisiGEV



###################### inizio programma

tempi.ritorno <- seq(10, 100, by = 10)

# dati di input: serie che cominciano dal 1951
dati <- read.csv(file = "dati1951.csv", sep = ";", head = T, check.names = F)

colonne <- ncol(dati)

risultati.gev <- lapply(2:colonne, analisiGEV)

## Stazione: 6220

## LOcation Scale and Shape; PACKAGE evd: 20.69
##  LOcation Scale and Shape; PACKAGE evd: 12.1
##  LOcation Scale and Shape; PACKAGE evd: 0.24
## [1] "Calcolo livello di ritorno per tempo: 10 anni"
## RLevel; PACKAGE evd: 57.41
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 20 anni"
## RLevel; PACKAGE evd: 73.05
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 30 anni"
## RLevel; PACKAGE evd: 83.69
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 40 anni"
## RLevel; PACKAGE evd: 91.49
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 50 anni"
## RLevel; PACKAGE evd: 97.83
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 60 anni"
## RLevel; PACKAGE evd: 103.62
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 70 anni"
## RLevel; PACKAGE evd: 105.06
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 80 anni"
## RLevel; PACKAGE evd: 110.96
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 90 anni"
## RLevel; PACKAGE evd: 115.59
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 100 anni"
## RLevel; PACKAGE evd: 116.38
## [1] "profiling quantile"

## [1] "profiling shape"
## Stazione: 6225

## LOcation Scale and Shape; PACKAGE evd: 11.17
##  LOcation Scale and Shape; PACKAGE evd: 9.8
##  LOcation Scale and Shape; PACKAGE evd: 0.36
## [1] "Calcolo livello di ritorno per tempo: 10 anni"
## RLevel; PACKAGE evd: 45.08
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 20 anni"
## RLevel; PACKAGE evd: 63.07
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 30 anni"
## RLevel; PACKAGE evd: 76.28
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 40 anni"
## RLevel; PACKAGE evd: 67.52
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 50 anni"
## RLevel; PACKAGE evd: 87.54
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 60 anni"
## RLevel; PACKAGE evd: 97.09
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 70 anni"
## RLevel; PACKAGE evd: 100.4
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 80 anni"
## RLevel; PACKAGE evd: 106.78
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 90 anni"
## RLevel; PACKAGE evd: 112.25
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 100 anni"
## RLevel; PACKAGE evd: 110.18
## [1] "profiling quantile"

## [1] "profiling shape"
## Stazione: 6259

## LOcation Scale and Shape; PACKAGE evd: 8.24
##  LOcation Scale and Shape; PACKAGE evd: 7.57
##  LOcation Scale and Shape; PACKAGE evd: 0.13
## [1] "Calcolo livello di ritorno per tempo: 10 anni"
## RLevel; PACKAGE evd: 28.15
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 20 anni"
## RLevel; PACKAGE evd: 35.89
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 30 anni"
## RLevel; PACKAGE evd: 40.7
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 40 anni"
## RLevel; PACKAGE evd: 44.26
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 50 anni"
## RLevel; PACKAGE evd: 47.11
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 60 anni"
## RLevel; PACKAGE evd: 49.46
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 70 anni"
## RLevel; PACKAGE evd: 51.57
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 80 anni"
## RLevel; PACKAGE evd: 53.13
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 90 anni"
## RLevel; PACKAGE evd: 55.05
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 100 anni"
## RLevel; PACKAGE evd: 56.02
## [1] "profiling quantile"

## [1] "profiling shape"
## Stazione: 6278

## LOcation Scale and Shape; PACKAGE evd: 14.82
##  LOcation Scale and Shape; PACKAGE evd: 10.63
##  LOcation Scale and Shape; PACKAGE evd: 0.26
## [1] "Calcolo livello di ritorno per tempo: 10 anni"
## RLevel; PACKAGE evd: 47.09
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 20 anni"
## RLevel; PACKAGE evd: 62.84
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 30 anni"
## RLevel; PACKAGE evd: 72.63
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 40 anni"
## RLevel; PACKAGE evd: 80.2
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 50 anni"
## RLevel; PACKAGE evd: 86.81
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 60 anni"
## RLevel; PACKAGE evd: 91.82
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 70 anni"
## RLevel; PACKAGE evd: 95.9
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 80 anni"
## RLevel; PACKAGE evd: 99.6
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 90 anni"
## RLevel; PACKAGE evd: 101.97
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 100 anni"
## RLevel; PACKAGE evd: 104.51
## [1] "profiling quantile"

## [1] "profiling shape"
## Stazione: 6429

## LOcation Scale and Shape; PACKAGE evd: 8.87
##  LOcation Scale and Shape; PACKAGE evd: 8.16
##  LOcation Scale and Shape; PACKAGE evd: 0.32
## [1] "Calcolo livello di ritorno per tempo: 10 anni"
## RLevel; PACKAGE evd: 35.6
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 20 anni"
## RLevel; PACKAGE evd: 48.99
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 30 anni"
## RLevel; PACKAGE evd: 58.21
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 40 anni"
## RLevel; PACKAGE evd: 65.33
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 50 anni"
## RLevel; PACKAGE evd: 70.93
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 60 anni"
## RLevel; PACKAGE evd: 63.18
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 70 anni"
## RLevel; PACKAGE evd: 81.22
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 80 anni"
## RLevel; PACKAGE evd: 66.92
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 90 anni"
## RLevel; PACKAGE evd: 84.51
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 100 anni"
## RLevel; PACKAGE evd: 93.24
## [1] "profiling quantile"

## [1] "profiling shape"
## Stazione: 6709

## LOcation Scale and Shape; PACKAGE evd: 34.04
##  LOcation Scale and Shape; PACKAGE evd: 15.11
##  LOcation Scale and Shape; PACKAGE evd: -0.02
## [1] "Calcolo livello di ritorno per tempo: 10 anni"
## RLevel; PACKAGE evd: 67.3
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 20 anni"
## RLevel; PACKAGE evd: 77.71
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 30 anni"
## RLevel; PACKAGE evd: 83.72
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 40 anni"
## RLevel; PACKAGE evd: 87.85
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 50 anni"
## RLevel; PACKAGE evd: 90.99
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 60 anni"
## RLevel; PACKAGE evd: 93.62
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 70 anni"
## RLevel; PACKAGE evd: 95.97
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 80 anni"
## RLevel; PACKAGE evd: 97.82
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 90 anni"
## RLevel; PACKAGE evd: 99.47
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 100 anni"
## RLevel; PACKAGE evd: 101
## [1] "profiling quantile"

## [1] "profiling shape"
## Stazione: 6720

## LOcation Scale and Shape; PACKAGE evd: 29.3
##  LOcation Scale and Shape; PACKAGE evd: 15.53
##  LOcation Scale and Shape; PACKAGE evd: 0.22
## [1] "Calcolo livello di ritorno per tempo: 10 anni"
## RLevel; PACKAGE evd: 74.37
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 20 anni"
## RLevel; PACKAGE evd: 94.06
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 30 anni"
## RLevel; PACKAGE evd: 106.86
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 40 anni"
## RLevel; PACKAGE evd: 112.44
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 50 anni"
## RLevel; PACKAGE evd: 117.42
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 60 anni"
## RLevel; PACKAGE evd: 121.46
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 70 anni"
## RLevel; PACKAGE evd: 135.73
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 80 anni"
## RLevel; PACKAGE evd: 139.36
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 90 anni"
## RLevel; PACKAGE evd: 130.57
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 100 anni"
## RLevel; PACKAGE evd: 144.01
## [1] "profiling quantile"

## [1] "profiling shape"
## Stazione: 6721

## LOcation Scale and Shape; PACKAGE evd: 37.52
##  LOcation Scale and Shape; PACKAGE evd: 15.55
##  LOcation Scale and Shape; PACKAGE evd: -0.02
## [1] "Calcolo livello di ritorno per tempo: 10 anni"
## RLevel; PACKAGE evd: 71.89
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 20 anni"
## RLevel; PACKAGE evd: 82.68
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 30 anni"
## RLevel; PACKAGE evd: 88.82
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 40 anni"
## RLevel; PACKAGE evd: 93.1
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 50 anni"
## RLevel; PACKAGE evd: 96.39
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 60 anni"
## RLevel; PACKAGE evd: 99.1
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 70 anni"
## RLevel; PACKAGE evd: 101.41
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 80 anni"
## RLevel; PACKAGE evd: 103.51
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 90 anni"
## RLevel; PACKAGE evd: 105.25
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 100 anni"
## RLevel; PACKAGE evd: 106.83
## [1] "profiling quantile"

## [1] "profiling shape"
## Stazione: 6723

## LOcation Scale and Shape; PACKAGE evd: 35.6
##  LOcation Scale and Shape; PACKAGE evd: 13.38
##  LOcation Scale and Shape; PACKAGE evd: 0.13
## [1] "Calcolo livello di ritorno per tempo: 10 anni"
## RLevel; PACKAGE evd: 70.36
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 20 anni"
## RLevel; PACKAGE evd: 83.75
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 30 anni"
## RLevel; PACKAGE evd: 92.12
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 40 anni"
## RLevel; PACKAGE evd: 98.01
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 50 anni"
## RLevel; PACKAGE evd: 103.12
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 60 anni"
## RLevel; PACKAGE evd: 107.24
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 70 anni"
## RLevel; PACKAGE evd: 110.79
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 80 anni"
## RLevel; PACKAGE evd: 113.81
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 90 anni"
## RLevel; PACKAGE evd: 116.61
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 100 anni"
## RLevel; PACKAGE evd: 119.35
## [1] "profiling quantile"

## [1] "profiling shape"
## Stazione: 6729

## LOcation Scale and Shape; PACKAGE evd: 34.45
##  LOcation Scale and Shape; PACKAGE evd: 12.47
##  LOcation Scale and Shape; PACKAGE evd: 0.01
## [1] "Calcolo livello di ritorno per tempo: 10 anni"
## RLevel; PACKAGE evd: 62.97
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 20 anni"
## RLevel; PACKAGE evd: 72.25
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 30 anni"
## RLevel; PACKAGE evd: 77.67
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 40 anni"
## RLevel; PACKAGE evd: 81.5
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 50 anni"
## RLevel; PACKAGE evd: 84.47
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 60 anni"
## RLevel; PACKAGE evd: 86.84
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 70 anni"
## RLevel; PACKAGE evd: 88.74
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 80 anni"
## RLevel; PACKAGE evd: 90.63
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 90 anni"
## RLevel; PACKAGE evd: 92.12
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 100 anni"
## RLevel; PACKAGE evd: 93.76
## [1] "profiling quantile"

## [1] "profiling shape"
## Stazione: 6736

## LOcation Scale and Shape; PACKAGE evd: 37.08
##  LOcation Scale and Shape; PACKAGE evd: 15.39
##  LOcation Scale and Shape; PACKAGE evd: 0.14
## [1] "Calcolo livello di ritorno per tempo: 10 anni"
## RLevel; PACKAGE evd: 77.74
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 20 anni"
## RLevel; PACKAGE evd: 93.87
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 30 anni"
## RLevel; PACKAGE evd: 103.58
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 40 anni"
## RLevel; PACKAGE evd: 111.67
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 50 anni"
## RLevel; PACKAGE evd: 116.78
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 60 anni"
## RLevel; PACKAGE evd: 127.98
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 70 anni"
## RLevel; PACKAGE evd: 126.2
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 80 anni"
## RLevel; PACKAGE evd: 130.02
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 90 anni"
## RLevel; PACKAGE evd: 137.18
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 100 anni"
## RLevel; PACKAGE evd: 139.53
## [1] "profiling quantile"

## [1] "profiling shape"
## Stazione: 6769

## LOcation Scale and Shape; PACKAGE evd: 31.01
##  LOcation Scale and Shape; PACKAGE evd: 15.04
##  LOcation Scale and Shape; PACKAGE evd: 0.07
## [1] "Calcolo livello di ritorno per tempo: 10 anni"
## RLevel; PACKAGE evd: 67.87
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 20 anni"
## RLevel; PACKAGE evd: 80.98
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 30 anni"
## RLevel; PACKAGE evd: 88.7
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 40 anni"
## RLevel; PACKAGE evd: 94.42
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 50 anni"
## RLevel; PACKAGE evd: 98.83
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 60 anni"
## RLevel; PACKAGE evd: 102.64
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 70 anni"
## RLevel; PACKAGE evd: 103.04
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 80 anni"
## RLevel; PACKAGE evd: 105.26
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 90 anni"
## RLevel; PACKAGE evd: 110.85
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 100 anni"
## RLevel; PACKAGE evd: 113
## [1] "profiling quantile"

## [1] "profiling shape"
## Stazione: 6773

## LOcation Scale and Shape; PACKAGE evd: 24.26
##  LOcation Scale and Shape; PACKAGE evd: 12.74
##  LOcation Scale and Shape; PACKAGE evd: 0.19
## [1] "Calcolo livello di ritorno per tempo: 10 anni"
## RLevel; PACKAGE evd: 60.53
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 20 anni"
## RLevel; PACKAGE evd: 75.14
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 30 anni"
## RLevel; PACKAGE evd: 84.6
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 40 anni"
## RLevel; PACKAGE evd: 91.77
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 50 anni"
## RLevel; PACKAGE evd: 97.58
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 60 anni"
## RLevel; PACKAGE evd: 102.6
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 70 anni"
## RLevel; PACKAGE evd: 107.04
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 80 anni"
## RLevel; PACKAGE evd: 110.67
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 90 anni"
## RLevel; PACKAGE evd: 113.71
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 100 anni"
## RLevel; PACKAGE evd: 115.66
## [1] "profiling quantile"

## [1] "profiling shape"
## Stazione: 6779

## LOcation Scale and Shape; PACKAGE evd: 25.06
##  LOcation Scale and Shape; PACKAGE evd: 13.06
##  LOcation Scale and Shape; PACKAGE evd: 0
## [1] "Calcolo livello di ritorno per tempo: 10 anni"
## RLevel; PACKAGE evd: 54.3
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 20 anni"
## RLevel; PACKAGE evd: 63.57
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 30 anni"
## RLevel; PACKAGE evd: 68.89
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 40 anni"
## RLevel; PACKAGE evd: 72.55
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 50 anni"
## RLevel; PACKAGE evd: 75.46
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 60 anni"
## RLevel; PACKAGE evd: 77.83
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 70 anni"
## RLevel; PACKAGE evd: 79.83
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 80 anni"
## RLevel; PACKAGE evd: 81.57
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 90 anni"
## RLevel; PACKAGE evd: 83.09
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 100 anni"
## RLevel; PACKAGE evd: 84.46
## [1] "profiling quantile"

## [1] "profiling shape"
## Stazione: 6782

## LOcation Scale and Shape; PACKAGE evd: 34.24
##  LOcation Scale and Shape; PACKAGE evd: 15.36
##  LOcation Scale and Shape; PACKAGE evd: 0.09
## [1] "Calcolo livello di ritorno per tempo: 10 anni"
## RLevel; PACKAGE evd: 72.66
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 20 anni"
## RLevel; PACKAGE evd: 86.59
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 30 anni"
## RLevel; PACKAGE evd: 94.94
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 40 anni"
## RLevel; PACKAGE evd: 100.92
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 50 anni"
## RLevel; PACKAGE evd: 105.98
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 60 anni"
## RLevel; PACKAGE evd: 109.5
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 70 anni"
## RLevel; PACKAGE evd: 112.99
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 80 anni"
## RLevel; PACKAGE evd: 115.16
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 90 anni"
## RLevel; PACKAGE evd: 117.69
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 100 anni"
## RLevel; PACKAGE evd: 120.83
## [1] "profiling quantile"

## [1] "profiling shape"
## Stazione: 6794

## LOcation Scale and Shape; PACKAGE evd: 22.21
##  LOcation Scale and Shape; PACKAGE evd: 11.81
##  LOcation Scale and Shape; PACKAGE evd: 0.09
## [1] "Calcolo livello di ritorno per tempo: 10 anni"
## RLevel; PACKAGE evd: 51.76
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 20 anni"
## RLevel; PACKAGE evd: 62.02
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 30 anni"
## RLevel; PACKAGE evd: 69.1
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 40 anni"
## RLevel; PACKAGE evd: 73.94
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 50 anni"
## RLevel; PACKAGE evd: 77.69
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 60 anni"
## RLevel; PACKAGE evd: 80.67
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 70 anni"
## RLevel; PACKAGE evd: 83.21
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 80 anni"
## RLevel; PACKAGE evd: 85.73
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 90 anni"
## RLevel; PACKAGE evd: 87.74
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 100 anni"
## RLevel; PACKAGE evd: 89.59
## [1] "profiling quantile"

## [1] "profiling shape"
## Stazione: 6812

## LOcation Scale and Shape; PACKAGE evd: 12.26
##  LOcation Scale and Shape; PACKAGE evd: 8.63
##  LOcation Scale and Shape; PACKAGE evd: 0.25
## [1] "Calcolo livello di ritorno per tempo: 10 anni"
## RLevel; PACKAGE evd: 38.36
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 20 anni"
## RLevel; PACKAGE evd: 50.4
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 30 anni"
## RLevel; PACKAGE evd: 58.33
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 40 anni"
## RLevel; PACKAGE evd: 64.35
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 50 anni"
## RLevel; PACKAGE evd: 69.56
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 60 anni"
## RLevel; PACKAGE evd: 73.05
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 70 anni"
## RLevel; PACKAGE evd: 77.63
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 80 anni"
## RLevel; PACKAGE evd: 80.23
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 90 anni"
## RLevel; PACKAGE evd: 82.17
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 100 anni"
## RLevel; PACKAGE evd: 83.51
## [1] "profiling quantile"

## [1] "profiling shape"
## Stazione: 6817

## LOcation Scale and Shape; PACKAGE evd: 20.14
##  LOcation Scale and Shape; PACKAGE evd: 13.02
##  LOcation Scale and Shape; PACKAGE evd: 0.04
## [1] "Calcolo livello di ritorno per tempo: 10 anni"
## RLevel; PACKAGE evd: 50.64
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 20 anni"
## RLevel; PACKAGE evd: 60.88
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 30 anni"
## RLevel; PACKAGE evd: 66.91
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 40 anni"
## RLevel; PACKAGE evd: 71.2
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 50 anni"
## RLevel; PACKAGE evd: 74.48
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 60 anni"
## RLevel; PACKAGE evd: 77.16
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 70 anni"
## RLevel; PACKAGE evd: 79.42
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 80 anni"
## RLevel; PACKAGE evd: 81.36
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 90 anni"
## RLevel; PACKAGE evd: 83.21
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 100 anni"
## RLevel; PACKAGE evd: 84.88
## [1] "profiling quantile"

## [1] "profiling shape"
## Stazione: 6820

## LOcation Scale and Shape; PACKAGE evd: 18.68
##  LOcation Scale and Shape; PACKAGE evd: 12.42
##  LOcation Scale and Shape; PACKAGE evd: 0.19
## [1] "Calcolo livello di ritorno per tempo: 10 anni"
## RLevel; PACKAGE evd: 53.55
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 20 anni"
## RLevel; PACKAGE evd: 68.15
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 30 anni"
## RLevel; PACKAGE evd: 77.81
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 40 anni"
## RLevel; PACKAGE evd: 84.5
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 50 anni"
## RLevel; PACKAGE evd: 89.98
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 60 anni"
## RLevel; PACKAGE evd: 95.1
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 70 anni"
## RLevel; PACKAGE evd: 98.97
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 80 anni"
## RLevel; PACKAGE evd: 102.67
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 90 anni"
## RLevel; PACKAGE evd: 106.52
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 100 anni"
## RLevel; PACKAGE evd: 109.07
## [1] "profiling quantile"

## [1] "profiling shape"
## Stazione: 6822

## LOcation Scale and Shape; PACKAGE evd: 13.16
##  LOcation Scale and Shape; PACKAGE evd: 10.95
##  LOcation Scale and Shape; PACKAGE evd: 0.27
## [1] "Calcolo livello di ritorno per tempo: 10 anni"
## RLevel; PACKAGE evd: 47.14
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 20 anni"
## RLevel; PACKAGE evd: 63.23
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 30 anni"
## RLevel; PACKAGE evd: 73.45
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 40 anni"
## RLevel; PACKAGE evd: 81.49
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 50 anni"
## RLevel; PACKAGE evd: 88.51
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 60 anni"
## RLevel; PACKAGE evd: 93.26
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 70 anni"
## RLevel; PACKAGE evd: 96.8
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 80 anni"
## RLevel; PACKAGE evd: 80.36
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 90 anni"
## RLevel; PACKAGE evd: 105.95
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 100 anni"
## RLevel; PACKAGE evd: 105.64
## [1] "profiling quantile"

## [1] "profiling shape"
## Stazione: 6826

## LOcation Scale and Shape; PACKAGE evd: 21.25
##  LOcation Scale and Shape; PACKAGE evd: 9.94
##  LOcation Scale and Shape; PACKAGE evd: 0.11
## [1] "Calcolo livello di ritorno per tempo: 10 anni"
## RLevel; PACKAGE evd: 46.77
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 20 anni"
## RLevel; PACKAGE evd: 56.59
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 30 anni"
## RLevel; PACKAGE evd: 61.51
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 40 anni"
## RLevel; PACKAGE evd: 64.95
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 50 anni"
## RLevel; PACKAGE evd: 67.61
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 60 anni"
## RLevel; PACKAGE evd: 69.77
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 70 anni"
## RLevel; PACKAGE evd: 71.6
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 80 anni"
## RLevel; PACKAGE evd: 73.18
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 90 anni"
## RLevel; PACKAGE evd: 74.57
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 100 anni"
## RLevel; PACKAGE evd: 81.42
## [1] "profiling quantile"

## [1] "profiling shape"
## Stazione: 6827

## LOcation Scale and Shape; PACKAGE evd: 18.59
##  LOcation Scale and Shape; PACKAGE evd: 12.61
##  LOcation Scale and Shape; PACKAGE evd: 0.08
## [1] "Calcolo livello di ritorno per tempo: 10 anni"
## RLevel; PACKAGE evd: 49.85
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 20 anni"
## RLevel; PACKAGE evd: 60.38
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 30 anni"
## RLevel; PACKAGE evd: 67.85
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 40 anni"
## RLevel; PACKAGE evd: 72.79
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 50 anni"
## RLevel; PACKAGE evd: 76.73
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 60 anni"
## RLevel; PACKAGE evd: 79.83
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 70 anni"
## RLevel; PACKAGE evd: 82.57
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 80 anni"
## RLevel; PACKAGE evd: 84.94
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 90 anni"
## RLevel; PACKAGE evd: 87.19
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 100 anni"
## RLevel; PACKAGE evd: 88.92
## [1] "profiling quantile"

## [1] "profiling shape"
## Stazione: 6841

## LOcation Scale and Shape; PACKAGE evd: 12.04
##  LOcation Scale and Shape; PACKAGE evd: 10.42
##  LOcation Scale and Shape; PACKAGE evd: 0.16
## [1] "Calcolo livello di ritorno per tempo: 10 anni"
## RLevel; PACKAGE evd: 40.26
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 20 anni"
## RLevel; PACKAGE evd: 51.71
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 30 anni"
## RLevel; PACKAGE evd: 58.74
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 40 anni"
## RLevel; PACKAGE evd: 63.96
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 50 anni"
## RLevel; PACKAGE evd: 68.53
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 60 anni"
## RLevel; PACKAGE evd: 72.15
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 70 anni"
## RLevel; PACKAGE evd: 75.07
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 80 anni"
## RLevel; PACKAGE evd: 76.78
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 90 anni"
## RLevel; PACKAGE evd: 80.46
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 100 anni"
## RLevel; PACKAGE evd: 82.25
## [1] "profiling quantile"

## [1] "profiling shape"
## Stazione: 6850

## LOcation Scale and Shape; PACKAGE evd: 14.32
##  LOcation Scale and Shape; PACKAGE evd: 9.46
##  LOcation Scale and Shape; PACKAGE evd: 0.01
## [1] "Calcolo livello di ritorno per tempo: 10 anni"
## RLevel; PACKAGE evd: 35.93
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 20 anni"
## RLevel; PACKAGE evd: 42.98
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 30 anni"
## RLevel; PACKAGE evd: 47.07
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 40 anni"
## RLevel; PACKAGE evd: 49.96
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 50 anni"
## RLevel; PACKAGE evd: 52.16
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 60 anni"
## RLevel; PACKAGE evd: 53.99
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 70 anni"
## RLevel; PACKAGE evd: 55.55
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 80 anni"
## RLevel; PACKAGE evd: 56.95
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 90 anni"
## RLevel; PACKAGE evd: 58.02
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 100 anni"
## RLevel; PACKAGE evd: 59.14
## [1] "profiling quantile"

## [1] "profiling shape"
## Stazione: 6851

## LOcation Scale and Shape; PACKAGE evd: 19.41
##  LOcation Scale and Shape; PACKAGE evd: 11.74
##  LOcation Scale and Shape; PACKAGE evd: 0.23
## [1] "Calcolo livello di ritorno per tempo: 10 anni"
## RLevel; PACKAGE evd: 54.31
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 20 anni"
## RLevel; PACKAGE evd: 69.28
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 30 anni"
## RLevel; PACKAGE evd: 78.69
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 40 anni"
## RLevel; PACKAGE evd: 86.92
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 50 anni"
## RLevel; PACKAGE evd: 92.45
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 60 anni"
## RLevel; PACKAGE evd: 97.15
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 70 anni"
## RLevel; PACKAGE evd: 102.11
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 80 anni"
## RLevel; PACKAGE evd: 104.52
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 90 anni"
## RLevel; PACKAGE evd: 107.5
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 100 anni"
## RLevel; PACKAGE evd: 108.83
## [1] "profiling quantile"

## [1] "profiling shape"
## Stazione: 6856

## LOcation Scale and Shape; PACKAGE evd: 11.41
##  LOcation Scale and Shape; PACKAGE evd: 9.34
##  LOcation Scale and Shape; PACKAGE evd: 0.34
## [1] "Calcolo livello di ritorno per tempo: 10 anni"
## RLevel; PACKAGE evd: 42.97
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 20 anni"
## RLevel; PACKAGE evd: 59.44
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 30 anni"
## RLevel; PACKAGE evd: 70.2
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 40 anni"
## RLevel; PACKAGE evd: 79.2
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 50 anni"
## RLevel; PACKAGE evd: 85.12
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 60 anni"
## RLevel; PACKAGE evd: 89.88
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 70 anni"
## RLevel; PACKAGE evd: 94.11
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 80 anni"
## RLevel; PACKAGE evd: 99.24
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 90 anni"
## RLevel; PACKAGE evd: 101
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 100 anni"
## RLevel; PACKAGE evd: 103.82
## [1] "profiling quantile"

## [1] "profiling shape"
## Stazione: 6869

## LOcation Scale and Shape; PACKAGE evd: 11.78
##  LOcation Scale and Shape; PACKAGE evd: 10.41
##  LOcation Scale and Shape; PACKAGE evd: 0.23
## [1] "Calcolo livello di ritorno per tempo: 10 anni"
## RLevel; PACKAGE evd: 42.62
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 20 anni"
## RLevel; PACKAGE evd: 56.48
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 30 anni"
## RLevel; PACKAGE evd: 65.42
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 40 anni"
## RLevel; PACKAGE evd: 72.03
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 50 anni"
## RLevel; PACKAGE evd: 78.14
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 60 anni"
## RLevel; PACKAGE evd: 82.71
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 70 anni"
## RLevel; PACKAGE evd: 86.09
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 80 anni"
## RLevel; PACKAGE evd: 90.31
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 90 anni"
## RLevel; PACKAGE evd: 92.79
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 100 anni"
## RLevel; PACKAGE evd: 92.85
## [1] "profiling quantile"

## [1] "profiling shape"
## Stazione: 6874

## LOcation Scale and Shape; PACKAGE evd: 5.48
##  LOcation Scale and Shape; PACKAGE evd: 5.74
##  LOcation Scale and Shape; PACKAGE evd: 0.5
## [1] "Calcolo livello di ritorno per tempo: 10 anni"
## RLevel; PACKAGE evd: 29.24
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 20 anni"
## RLevel; PACKAGE evd: 44.45
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 30 anni"
## RLevel; PACKAGE evd: 55.75
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 40 anni"
## RLevel; PACKAGE evd: 62.39
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 50 anni"
## RLevel; PACKAGE evd: 66.48
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 60 anni"
## RLevel; PACKAGE evd: 65.91
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 70 anni"
## RLevel; PACKAGE evd: 68.99
## [1] "profiling quantile"

## [1] "Calcolo livello di ritorno per tempo: 80 anni"

## Error: observed information matrix is singular; use std.err = FALSE

rl.gevDataFrame <- do.call("rbind", lapply(1:(colonne - 1), function(i) {
    risultati.gev[[i]][[1]]
}))

## Error: object 'risultati.gev' not found

xi.gevDataFrame <- do.call("rbind", lapply(1:(colonne - 1), function(i) {
    risultati.gev[[i]][[2]]
}))

## Error: object 'risultati.gev' not found

write.table(rl.gevDataFrame, file = "rl.gevDataFrame1951.csv", sep = ";", col.names = T, 
    row.names = F, quote = F)

## Error: object 'rl.gevDataFrame' not found

write.table(xi.gevDataFrame, file = "xi.gevDataFrame1951.csv", sep = ";", col.names = T, 
    row.names = F, quote = F)

## Error: object 'xi.gevDataFrame' not found