VGAM (version 0.8-2)

venice: Venice Maximum Sea Levels Data

Description

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

Usage

data(venice)

source

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

Details

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

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

matplot(venice[["year"]], venice[,-1], xlab="Year", ylab="Sea level (cm)")

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))

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)