Last data update: 2014.03.03
R: Return level plot of GEV
retlevel.gev.graph R 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
>