venice {VGAM} | R Documentation

## Venice Maximum Sea Levels Data

### Description

The maximum heights of sea levels recorded at Venice, Italy, between 1931 and 1981.

### Usage

    data(venice)

### Format

A data frame with 51 observations on the following 11 variables.

year
  a numeric vector.
r1
  a numeric vector; the highest recorded value.
r2,r3,r4,r5,r6,r7,r8,r9,r10
  numeric vectors; the second highest recorded value down to the tenth highest recorded value.

### Details

The units are cm. In 1935 only the top six values were recorded.

### Source

Pirazzoli, P. (1982) Maree estreme a Venezia (periodo 1872–1981). Acqua Aria, 10, 1023–1039.

### References

Smith, R. L. (1986) Extreme value theory based on the r largest annual events. Journal of Hydrology, 86, 27–43.

Rosen, O. and Cohen, A. (1996) Extreme Percentile Regression. In: Haerdle, W. and Schimek, M. G. (eds.) Statistical Theory and Computational Aspects of Smoothing: Proceedings of the COMPSTAT '94 Satellite Meeting held in Semmering, Austria, 27–28 August 1994, Physica-Verlag: Heidelberg, pp.200–214.

### Examples

    ## Not run:
    matplot(venice[["year"]], venice[,-1], xlab="Year", ylab="Sea level (cm)")
    ## End(Not run)

    y = as.matrix(venice[,paste("r",1:10,sep="")])
    fit1 = vgam(y ~ s(year, df=3), gumbel(R=365, mpv=TRUE),
                data=venice, trace=TRUE, na.action=na.pass)
    head(fitted(fit1))

    ## Not run:
    par(mfrow=c(2,1), xpd=TRUE)
    plot(fit1, se=TRUE, lcol="blue", llwd=2, slty="dashed")

    par(mfrow=c(1,1), bty="l", xpd=TRUE, las=1)
    qtplot(fit1, mpv=TRUE, lcol=c(1,2,5), tcol=c(1,2,5),
           llwd=2, pcol="blue", tadj=0.1)
    ## End(Not run)

[Package VGAM version 0.7-9 Index]