R Under development (unstable) (2023-11-25 r85635 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2023 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64

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(HyperbolicDist)
> testVG <- FALSE
> 
> ### Test ghyp
> m1 <- momIntegrated("ghyp", order = 1, param = c(1/2,3,1,1,0), about = 0)
> m1
[1] 0.4785534
> ghypMean(c(1/2,3,1,1,0))
[1] 0.4785534
> momIntegrated("ghyp", order = 1, param = c(1/2,3,1,1,0), about = m1)
[1] -1.387779e-17
> m2 <- momIntegrated("ghyp", order = 2, param = c(1/2,3,1,1,0), about = 0)
> m2
[1] 0.7830109
> m2 - m1^2
[1] 0.5539976
> momIntegrated("ghyp", order = 2, param = c(1/2,3,1,1,0), about = m1)
[1] 0.5539976
> momIntegrated("generalized hyperbolic", order = 2, param = c(1/2,3,1,1,0),
+               about = m1)
[1] 0.5539976
> ghypVar(c(1/2,3,1,1,0))
[1] 0.5539976
> 
> 
> ### Test hyperb
> m1 <- momIntegrated("hyperb", order = 1, param = c(2,1,1,0), about = 0)
> m1
[1] 5.398968
> hyperbMean(c(2,1,1,0))
[1] 5.398968
> momIntegrated("hyperb", order = 1, param = c(2,1,1,0), about = m1)
[1] -3.932415e-15
> m2 <- momIntegrated("hyperb", order = 2, param = c(2,1,1,0), about = 0)
> m2
[1] 49.89123
> m2 - m1^2
[1] 20.74237
> momIntegrated("hyperb", order = 2, param = c(2,1,1,0), about = m1)
[1] 20.74237
> momIntegrated("hyperbolic", order = 2, param = c(2,1,1,0),
+               about = m1)
[1] 20.74237
> hyperbVar(c(2,1,1,0))
[1] 20.74237
> 
> ### Test gig
> m1 <- momIntegrated("gig", order = 1, param = c(1,2,3), about = 0)
> m1
[1] 1.353536
> gigMean(c(1,2,3))
[1] 1.353536
> gigRawMom(1, c(1,2,3))
[1] 1.353536
> momIntegrated("gig", order = 1, param = c(1,2,3), about = m1)
[1] -6.471603e-16
> 
> m2 <- momIntegrated("gig", order = 2, param = c(1,2,3), about = 0)
> m2
[1] 2.471381
> m2 - m1^2
[1] 0.6393221
> momIntegrated("gig", order = 2, param = c(1,2,3), about = m1)
[1] 0.6393221
> momIntegrated("generalized inverse Gaussian", order = 2, param = c(1,2,3),
+               about = m1)
[1] 0.6393221
> gigVar(c(1,2,3))
[1] 0.6393221
> gigMom(2, c(1,2,3), about = m1)
[1] 0.6393221
> 
> ### Test gamma
> m1 <- momIntegrated("gamma", order = 1, param = c(2,3), about = 0)
> m1
[1] 0.6666667
> shape <- 2
> rate <- 3
> gigMom(1, c(shape, 0, 2*rate))
[1] 0.6666667
> momIntegrated("gamma", order = 1, param = c(2,3), about = m1)
[1] -4.510281e-17
> m2 <- momIntegrated("gamma", order = 2, param = c(2,3), about = 0)
> m2
[1] 0.6666667
> m2 - m1^2
[1] 0.2222222
> momIntegrated("gamma", order = 2, param = c(2,3), about = m1)
[1] 0.2222222
> gigMom(2, c(shape, 0, 2*rate), about = m1)
[1] 0.2222222
> 
> ### Test inverse gamma
> m1 <- momIntegrated("invgamma", order = 1, param = c(3,4), about = 0)
> m1
[1] 0.125
> shape <- 3
> rate <- 1/4
> gigMom(order = -1, c(shape, 0, 2*rate))
[1] 0.125
> momIntegrated("invgamma", order = 1, param = c(3,4), about = m1)
[1] -5.809155e-16
> m2 <- momIntegrated("invgamma", order = 2, param = c(3,4), about = 0)
> m2
[1] 0.03125
> gigMom(order = -2, c(shape, 0, 2*rate))
[1] 0.03125
> m2 - m1^2
[1] 0.015625
> momIntegrated("invgamma", order = 2, param = c(3,4), about = m1)
[1] 0.015625
> 
> 
> proc.time()
   user  system elapsed 
   0.18    0.00    0.18