Last data update: 2014.03.03

R: Return level plot of GEV
retlevel.gev.graphR Documentation

Return level plot of GEV

Description

Return level plot of GEV distribution.

Usage

retlevel.gev.graph(vector)

Arguments

vector

a list object returned by posterior.gev

Value

The program returns return level plot from 1 until 100 periods in time, with 95% credibility intervals

Examples

# Return level plot for River nidd data
data(nidd.annual)
out=posterior.gev(nidd.annual,1,1000)
retlevel.gev.graph(out)
# Return level plot for ibovespa 15 day maxima
data(ibovespa)
postibv=posterior.gev(ibovespa[,4],15,300)
retlevel.gev.graph(postibv)

Results


R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(MCMC4Extremes)
Loading required package: evir
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/MCMC4Extremes/retlevel.gev.graph.Rd_%03d_medium.png", width=480, height=480)
> ### Name: retlevel.gev.graph
> ### Title: Return level plot of GEV
> ### Aliases: retlevel.gev.graph
> 
> ### ** Examples
> 
> # Return level plot for River nidd data
> data(nidd.annual)
> out=posterior.gev(nidd.annual,1,1000)
[1] 0.006666667
[1] 0.01333333
[1] 0.02
[1] 0.02666667
[1] 0.03333333
[1] 0.04
[1] 0.04666667
[1] 0.05333333
[1] 0.06
[1] 0.06666667
[1] 0.07333333
[1] 0.08
[1] 0.08666667
[1] 0.09333333
[1] 0.1
[1] 0.1066667
[1] 0.1133333
[1] 0.12
[1] 0.1266667
[1] 0.1333333
[1] 0.14
[1] 0.1466667
[1] 0.1533333
[1] 0.16
[1] 0.1666667
[1] 0.1733333
[1] 0.18
[1] 0.1866667
[1] 0.1933333
[1] 0.2
[1] 0.2066667
[1] 0.2133333
[1] 0.22
[1] 0.2266667
[1] 0.2333333
[1] 0.24
[1] 0.2466667
[1] 0.2533333
[1] 0.26
[1] 0.2666667
[1] 0.2733333
[1] 0.28
[1] 0.2866667
[1] 0.2933333
[1] 0.3
[1] 0.3066667
[1] 0.3133333
[1] 0.32
[1] 0.3266667
[1] 0.3333333
[1] 0.34
[1] 0.3466667
[1] 0.3533333
[1] 0.36
[1] 0.3666667
[1] 0.3733333
[1] 0.38
[1] 0.3866667
[1] 0.3933333
[1] 0.4
[1] 0.4066667
[1] 0.4133333
[1] 0.42
[1] 0.4266667
[1] 0.4333333
[1] 0.44
[1] 0.4466667
[1] 0.4533333
[1] 0.46
[1] 0.4666667
[1] 0.4733333
[1] 0.48
[1] 0.4866667
[1] 0.4933333
[1] 0.5
[1] 0.5066667
[1] 0.5133333
[1] 0.52
[1] 0.5266667
[1] 0.5333333
[1] 0.54
[1] 0.5466667
[1] 0.5533333
[1] 0.56
[1] 0.5666667
[1] 0.5733333
[1] 0.58
[1] 0.5866667
[1] 0.5933333
[1] 0.6
[1] 0.6066667
[1] 0.6133333
[1] 0.62
[1] 0.6266667
[1] 0.6333333
[1] 0.64
[1] 0.6466667
[1] 0.6533333
[1] 0.66
[1] 0.6666667
[1] 0.6733333
[1] 0.68
[1] 0.6866667
[1] 0.6933333
[1] 0.7
[1] 0.7066667
[1] 0.7133333
[1] 0.72
[1] 0.7266667
[1] 0.7333333
[1] 0.74
[1] 0.7466667
[1] 0.7533333
[1] 0.76
[1] 0.7666667
[1] 0.7733333
[1] 0.78
[1] 0.7866667
[1] 0.7933333
[1] 0.8
[1] 0.8066667
[1] 0.8133333
[1] 0.82
[1] 0.8266667
[1] 0.8333333
[1] 0.84
[1] 0.8466667
[1] 0.8533333
[1] 0.86
[1] 0.8666667
[1] 0.8733333
[1] 0.88
[1] 0.8866667
[1] 0.8933333
[1] 0.9
[1] 0.9066667
[1] 0.9133333
[1] 0.92
[1] 0.9266667
[1] 0.9333333
[1] 0.94
[1] 0.9466667
[1] 0.9533333
[1] 0.96
[1] 0.9666667
[1] 0.9733333
[1] 0.98
[1] 0.9866667
[1] 0.9933333
[1] 1
> retlevel.gev.graph(out)
  [1] -3382.59496    96.46255   115.90324   130.16895   140.89926   150.55411
  [7]   158.50652   165.06856   171.45652   176.80714   181.64919   186.15851
 [13]   190.10977   194.26781   198.12825   201.24069   204.57328   207.31969
 [19]   209.82301   212.22079   214.63637   216.77949   219.47470   222.06446
 [25]   223.47617   225.52520   227.59318   229.60800   231.56609   232.99886
 [31]   234.32743   235.73089   237.71571   239.04145   240.17853   241.71437
 [37]   243.44126   245.20054   246.92591   248.20884   249.12484   250.01779
 [43]   250.88879   251.93917   253.27586   254.60092   255.90642   257.18964
 [49]   258.08685   259.34028   260.64923   261.93854   262.85215   263.72203
 [55]   264.57678   265.41694   266.24302   267.05956   267.86599   268.65962
 [61]   269.44088   270.21016   270.77546   271.24599   271.70804   272.16188
 [67]   272.60779   273.04601   273.47680   273.90039   274.42770   275.14337
 [73]   275.84998   276.54777   277.23697   277.91780   278.59312   279.26547
 [79]   279.93000   280.58912   281.24307   281.88972   282.52925   283.16182
 [85]   283.78759   284.40535   284.84568   285.28090   285.91482   286.64630
 [91]   287.38936   287.96353   288.53201   289.09494   289.65243   290.20458
 [97]   290.75149   291.29328   291.83005   292.36227
> # Return level plot for ibovespa 15 day maxima
> data(ibovespa)
> postibv=posterior.gev(ibovespa[,4],15,300)
[1] 0.02222222
[1] 0.04444444
[1] 0.06666667
[1] 0.08888889
[1] 0.1111111
[1] 0.1333333
[1] 0.1555556
[1] 0.1777778
[1] 0.2
[1] 0.2222222
[1] 0.2444444
[1] 0.2666667
[1] 0.2888889
[1] 0.3111111
[1] 0.3333333
[1] 0.3555556
[1] 0.3777778
[1] 0.4
[1] 0.4222222
[1] 0.4444444
[1] 0.4666667
[1] 0.4888889
[1] 0.5111111
[1] 0.5333333
[1] 0.5555556
[1] 0.5777778
[1] 0.6
[1] 0.6222222
[1] 0.6444444
[1] 0.6666667
[1] 0.6888889
[1] 0.7111111
[1] 0.7333333
[1] 0.7555556
[1] 0.7777778
[1] 0.8
[1] 0.8222222
[1] 0.8444444
[1] 0.8666667
[1] 0.8888889
[1] 0.9111111
[1] 0.9333333
[1] 0.9555556
[1] 0.9777778
[1] 1
> retlevel.gev.graph(postibv)
  [1] -0.06593151  0.02814844  0.03344719  0.03716782  0.03992601  0.04210985
  [7]  0.04404546  0.04581900  0.04741653  0.04874576  0.04996950  0.05115331
 [13]  0.05225276  0.05328273  0.05425233  0.05516894  0.05603297  0.05683067
 [19]  0.05758184  0.05828184  0.05894941  0.05958484  0.06019513  0.06080114
 [25]  0.06140789  0.06194925  0.06246679  0.06298472  0.06348630  0.06397149
 [31]  0.06444276  0.06490095  0.06539550  0.06588549  0.06639462  0.06689413
 [37]  0.06738297  0.06786167  0.06833070  0.06879048  0.06924145  0.06968396
 [43]  0.07009301  0.07044856  0.07079629  0.07113744  0.07147226  0.07179716
 [49]  0.07210497  0.07240725  0.07270420  0.07299604  0.07328293  0.07356506
 [55]  0.07384259  0.07411568  0.07438448  0.07464913  0.07490976  0.07516651
 [61]  0.07541949  0.07566837  0.07591364  0.07615549  0.07639401  0.07662930
 [67]  0.07686144  0.07709054  0.07731796  0.07755501  0.07778916  0.07802050
 [73]  0.07824909  0.07847481  0.07869732  0.07891730  0.07913480  0.07934989
 [79]  0.07956262  0.07977305  0.07998123  0.08018721  0.08039104  0.08059278
 [85]  0.08079246  0.08099013  0.08118584  0.08137962  0.08157152  0.08176158
 [91]  0.08194984  0.08213632  0.08232107  0.08250412  0.08268551  0.08286527
 [97]  0.08304343  0.08322001  0.08339506  0.08356859
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>