venice {VGAM} | R Documentation |
The maximum heights of sea levels recorded at Venice, Italy, between 1931 and 1981.
data(venice)
A data frame with 51 observations on the following 11 variables.
The units are cm. In 1935 only the top six values were recorded.
Pirazzoli, P. (1982) Maree estreme a Venezia (periodo 1872–1981). Acqua Aria, 10, 1023–1039.
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.
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)