venice: Sea Level Data

In marg: Approximate Marginal Inference for Regression-Scale Models

Description

The venice data frame has 51 rows and 2 columns.

Pirazzoli (1982) collected the ten largest values of sea levels in Venice (with a few exceptions) for the years 1887–1981. The venice data frame contains the maxima for the years 1931–1981.

Usage

data(venice)

Format

This data frame contains the following columns:

year

the years;

sea

the sea levels (in cm).

Source

The data were obtained from

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

References

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

Examples

data(venice)
attach(venice)
#
plot(sea ~ year, ylab = "sea level")
##
Year <- 1:51/51
venice.l <- rsm(sea ~ Year + I(Year^2), family = extreme)
lines(year, fitted(venice.l))
##
c11 <- cos(2*pi*1:51/11) ; s11 <- sin(2*pi*1:51/11)
c19 <- cos(2*pi*1:51/18.62) ; s19 <- sin(2*pi*1:51/18.62)
venice.p <- rsm(sea ~ Year + I(Year^2) + c11 + s11 + c19 + s19, 
                family = extreme)
lines(year, fitted(venice.p), col = "red")
##
detach()

Example output

Loading required package: statmod
Loading required package: survival

   Package "marg" 1.2-2 (2014-03-31) 
    Copyright (C) 2000-2014 A. R. Brazzale

 This is free software, and you are welcome to redistribute
 it and/or modify it under the terms of the GNU General
 Public License published by the Free Software Foundation.
 Package "marg" comes with ABSOLUTELY NO WARRANTY.

 type `help(package="marg")' for summary information

marg documentation built on May 2, 2019, 7:55 a.m.