venice {VGAM}R Documentation

Venice Maximum Sea Levels

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
a numeric vector; the second highest recorded value.
r3
a numeric vector; the third highest recorded value.
r4
a numeric vector; the fourth highest recorded value.
r5
a numeric vector; the fifth highest recorded value.
r6
a numeric vector; the sixth highest recorded value.
r7
a numeric vector; the seventh highest recorded value.
r8
a numeric vector; the eighth highest recorded value.
r9
a numeric vector; the ninth highest recorded value.
r10
a numeric vector; 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

data(venice)
## 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)
fitted(fit1)[1:4,]

## 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-1 Index]