SkewHyperbolic
The Skew Hyperbolic Student t-Distribution

Package index
Search the SkewHyperbolic package
Vignettes
Functions
51
Source code
19
Man pages
16
Home
/
R-Forge
/
SkewHyperbolic
/
tests/testqskewhyp.R

tests/testqskewhyp.R
In SkewHyperbolic: The Skew Hyperbolic Student t-Distribution 

require(SkewHyperbolic)
## source("../R/skewhypCalcRange.R")
## source("../R/dskewhyp.R")
## source("../data/skewhypParam.R")
## library(RUnit)

options(digits=20)
param <- c(0,1,0,10)
q <- c(-Inf,-1,0,1,Inf)
pskewhyp(q, param = param)
pskewhyp(q, param = param, lower.tail = FALSE)
pskewhyp(q, param = param, valueOnly = FALSE)
pskewhyp(q, param = param, valueOnly = FALSE, intTol = 10^(-12))

x <- rskewhyp(1, param = param)
x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-10)),
            param = param) - x
qskewhyp(pskewhyp(x, param = param),
            param = param, uniTol = 10^(-10)) - x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-8)),
            param = param, uniTol = 10^(-8)) - x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-10)),
            param = param, uniTol = 10^(-10)) - x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-12)),
            param = param, uniTol = 10^(-12)) - x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-10)),
         param = param, method = "integrate") - x
qskewhyp(pskewhyp(x, param = param),
            param = param, uniTol = 10^(-10)) - x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-8)),
         param = param, uniTol = 10^(-8), method = "integrate") - x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-10)),
         param = param, uniTol = 10^(-10), method = "integrate") - x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-12)),
         param = param, uniTol = 10^(-12), method = "integrate") - x
qskewhyp(pskewhyp(10, param = param, intTol = 10^(-12)),
         param = param, uniTol = 10^(-12), method = "integrate") - 10
qskewhyp(pskewhyp(-10, param = param, intTol = 10^(-12)),
         param = param, uniTol = 10^(-12), method = "integrate") + 10


param <- c(0,1,10,20)
q <- c(-Inf,-1,0,1,Inf)
pskewhyp(q, param = param)
pskewhyp(q, param = param, lower.tail = FALSE)
pskewhyp(q, param = param, valueOnly = FALSE)
pskewhyp(q, param = param, valueOnly = FALSE, intTol = 10^(-12))

x <- rskewhyp(1, param = param)
x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-10)),
            param = param) - x
qskewhyp(pskewhyp(x, param = param),
            param = param, uniTol = 10^(-10)) - x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-8)),
            param = param, uniTol = 10^(-8)) - x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-10)),
            param = param, uniTol = 10^(-10)) - x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-12)),
            param = param, uniTol = 10^(-12)) - x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-10)),
         param = param, method = "integrate") - x
qskewhyp(pskewhyp(x, param = param),
            param = param, uniTol = 10^(-10)) - x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-8)),
         param = param, uniTol = 10^(-8), method = "integrate") - x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-10)),
         param = param, uniTol = 10^(-10), method = "integrate") - x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-12)),
         param = param, uniTol = 10^(-12), method = "integrate") - x
qskewhyp(pskewhyp(10, param = param, intTol = 10^(-12)),
         param = param, uniTol = 10^(-12), method = "integrate") - 10
qskewhyp(pskewhyp(-10, param = param, intTol = 10^(-12)),
         param = param, uniTol = 10^(-12), method = "integrate") + 10

param <- c(0,1,1,1)
q <- c(-Inf,-1,0,1,Inf)
pskewhyp(q, param = param)
pskewhyp(q, param = param, lower.tail = FALSE)
pskewhyp(q, param = param, valueOnly = FALSE)
pskewhyp(q, param = param, valueOnly = FALSE, intTol = 10^(-12))


param <- c(0,1,-10,5)
q <- c(-Inf,-1,0,1,Inf)
pskewhyp(q, param = param)
pskewhyp(q, param = param, lower.tail = FALSE)
pskewhyp(q, param = param, valueOnly = FALSE)
pskewhyp(q, param = param, valueOnly = FALSE, intTol = 10^(-12))



param <- c(0,1,5,5)
q <- c(-Inf,-1,0,1,Inf)
pskewhyp(q, param = param)
pskewhyp(q, param = param, lower.tail = FALSE)
pskewhyp(q, param = param, valueOnly = FALSE)
pskewhyp(q, param = param, valueOnly = FALSE, intTol = 10^(-12))

x <- rskewhyp(1, param = param)
x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-10)),
            param = param) - x
qskewhyp(pskewhyp(x, param = param),
            param = param, uniTol = 10^(-10)) - x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-8)),
            param = param, uniTol = 10^(-8)) - x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-10)),
            param = param, uniTol = 10^(-10)) - x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-12)),
            param = param, uniTol = 10^(-12)) - x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-10)),
         param = param, method = "integrate") - x
qskewhyp(pskewhyp(x, param = param),
            param = param, uniTol = 10^(-10)) - x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-8)),
         param = param, uniTol = 10^(-8), method = "integrate") - x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-10)),
         param = param, uniTol = 10^(-10), method = "integrate") - x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-12)),
         param = param, uniTol = 10^(-12), method = "integrate") - x
qskewhyp(pskewhyp(10, param = param, intTol = 10^(-12)),
         param = param, uniTol = 10^(-12), method = "integrate") - 10
qskewhyp(pskewhyp(-10, param = param, intTol = 10^(-12)),
         param = param, uniTol = 10^(-12), method = "integrate") + 10

param <- c(1,2,20,10)
q <- c(-Inf,-1,0,1,Inf)
pskewhyp(q, param = param)
pskewhyp(q, param = param, lower.tail = FALSE)
pskewhyp(q, param = param, valueOnly = FALSE)
pskewhyp(q, param = param, valueOnly = FALSE, intTol = 10^(-12))

x <- rskewhyp(1, param = param)
x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-10)),
            param = param) - x
qskewhyp(pskewhyp(x, param = param),
            param = param, uniTol = 10^(-10)) - x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-8)),
            param = param, uniTol = 10^(-8)) - x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-10)),
            param = param, uniTol = 10^(-10)) - x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-12)),
            param = param, uniTol = 10^(-12)) - x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-10)),
         param = param, method = "integrate") - x
qskewhyp(pskewhyp(x, param = param),
            param = param, uniTol = 10^(-10)) - x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-8)),
         param = param, uniTol = 10^(-8), method = "integrate") - x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-10)),
         param = param, uniTol = 10^(-10), method = "integrate") - x
qskewhyp(pskewhyp(x, param = param, intTol = 10^(-12)),
         param = param, uniTol = 10^(-12), method = "integrate") - x
qskewhyp(pskewhyp(10, param = param, intTol = 10^(-12)),
         param = param, uniTol = 10^(-12), method = "integrate") - 10
qskewhyp(pskewhyp(-10, param = param, intTol = 10^(-12)),
         param = param, uniTol = 10^(-12), method = "integrate") + 10

Try the SkewHyperbolic package in your browser

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

SkewHyperbolic documentation built on Nov. 26, 2023, 3 p.m.

R Package Documentation

rdrr.io home R language documentation Run R code online

Browse R Packages

CRAN packages Bioconductor packages R-Forge packages GitHub packages

We want your feedback!

Note that we can't provide technical support on individual packages. You should contact the package authors for that.
GitHub issue tracker
ian@mutexlabs.com
Personal blog

 
Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.

Close