decluster: Decluster Point Process

In evir: Extreme Values in R

Description

Declusters clustered point process data so that Poisson assumption is more tenable over a high threshold.

Usage

decluster(series, run = NA, picture = TRUE)

Arguments

series	a numeric vector of threshold exceedances with a times attribute which should be a numeric vector containing either the indices or the times/dates of each exceedance (if times/dates, the attribute should be an object of class "POSIXct" or an object that can be converted to that class; see as.POSIXct)
run	parameter to be used in the runs method; any two consecutive threshold exceedances separated by more than this number of observations/days are considered to belong to different clusters
picture	whether or not a picture of declustering should be drawn

Value

The declustered object.

References

Embrechts, P., Klueppelberg, C., Mikosch, T. (1997). Modelling Extremal Events. Springer. Chapter 8, 413–429.

See Also

pot, exindex, as.POSIXct

Examples

# decluster the 200 exceedances of a particular threshold in 
# the negative BMW data
data(bmw)
out <- pot(-bmw, ne = 200) 
decluster(out$data, 30) 

Example output

Declustering picture...
Data reduced from 200 to 70 
 [1] 0.05525992 0.05084288 0.06689143 0.06887800 0.04688708 0.10617520
 [7] 0.05865478 0.04071008 0.02896766 0.02913723 0.04122546 0.03094040
[13] 0.02659780 0.03465775 0.02606819 0.02710781 0.06955422 0.03859640
[19] 0.03845126 0.06807869 0.02732490 0.04194781 0.02563799 0.04957226
[25] 0.03744841 0.02652559 0.03021099 0.04541566 0.02695412 0.03051167
[31] 0.03984488 0.02648610 0.04525858 0.03938106 0.03236528 0.05561317
[37] 0.03191673 0.04609111 0.06519183 0.03668970 0.03540193 0.04491698
[43] 0.04880244 0.04565154 0.04148352 0.02619353 0.10852160 0.02954286
[49] 0.03225168 0.14061565 0.02895132 0.02834847 0.03360120 0.07529138
[55] 0.04328683 0.03405556 0.02677080 0.03221330 0.10577463 0.02764203
[61] 0.05806183 0.03004121 0.03973532 0.02892858 0.02688352 0.04001846
[67] 0.03412405 0.03232889 0.03123302 0.03490287
attr(,"times")
 [1] "1973-02-05 UTC" "1973-07-26 UTC" "1973-11-14 UTC" "1973-12-31 UTC"
 [5] "1974-05-09 UTC" "1974-07-05 UTC" "1974-09-30 UTC" "1975-05-27 UTC"
 [9] "1976-04-06 UTC" "1976-07-28 UTC" "1976-10-12 UTC" "1976-11-15 UTC"
[13] "1977-02-18 UTC" "1977-07-08 UTC" "1977-12-06 UTC" "1978-06-21 UTC"
[17] "1979-07-11 UTC" "1980-01-21 UTC" "1980-03-18 UTC" "1980-06-20 UTC"
[21] "1980-09-08 UTC" "1981-01-27 UTC" "1981-04-29 UTC" "1981-06-26 UTC"
[25] "1981-08-25 UTC" "1982-04-27 UTC" "1982-08-17 UTC" "1982-09-27 UTC"
[29] "1983-01-24 UTC" "1983-07-12 UTC" "1984-02-09 UTC" "1984-03-12 UTC"
[33] "1984-07-05 UTC" "1984-10-19 UTC" "1985-03-25 UTC" "1985-07-10 UTC"
[37] "1985-09-13 UTC" "1985-11-05 UTC" "1986-02-27 UTC" "1986-04-28 UTC"
[41] "1986-06-03 UTC" "1986-07-04 UTC" "1986-09-05 UTC" "1987-02-04 UTC"
[45] "1987-04-27 UTC" "1987-06-02 UTC" "1987-11-09 UTC" "1988-03-25 UTC"
[49] "1988-05-11 UTC" "1989-10-16 UTC" "1990-01-04 UTC" "1990-02-26 UTC"
[53] "1990-04-23 UTC" "1990-09-25 UTC" "1991-01-14 UTC" "1991-02-26 UTC"
[57] "1991-05-17 UTC" "1991-06-28 UTC" "1991-08-19 UTC" "1992-05-13 UTC"
[61] "1992-09-24 UTC" "1993-01-25 UTC" "1993-05-14 UTC" "1993-07-27 UTC"
[65] "1993-11-22 UTC" "1994-06-20 UTC" "1994-11-23 UTC" "1995-03-09 UTC"
[69] "1995-09-22 UTC" "1995-10-23 UTC"

evir documentation built on May 2, 2019, 12:21 p.m.