tests/test-zvalue.R

library(RTDE)


#####
# (1) small example

omega <- 1/2
m <- 10
n <- 100
obs <- cbind(rupareto(n), rupareto(n)) + rupareto(n)

#unit Pareto transform
zvalueRTDE(obs, omega, output="orig")

relexcess(zvalueRTDE(obs, omega, output="orig"), m)
zvalueRTDE(obs, omega, nbpoint=m, output="relexcess")

#unit Frechet transform
zvalueRTDE(obs, omega, output="orig", marg="ufrechet")

relexcess(zvalueRTDE(obs, omega, output="orig", marg="ufrechet"), m)
zvalueRTDE(obs, omega, nbpoint=m, output="relexcess", marg="ufrechet")

Try the RTDE package in your browser

Any scripts or data that you put into this service are public.

RTDE documentation built on Jan. 8, 2020, 5:09 p.m.