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)
> 
> ### Test gigRawMom
> Theta <- c(-0.5,5,2.5)
> gigMean(Theta)
[1] 1.414214
> gigRawMom(1, Theta)
[1] 1.414214
> momIntegrated("gig", order = 1, param = Theta, about = 0)
[1] 1.414214
> gigRawMom(2, Theta)
[1] 2.565685
> momIntegrated("gig", order = 2, param = Theta, about = 0)
[1] 2.565685
> gigRawMom(10, Theta)
[1] 306061
> momIntegrated("gig", order = 10, param = Theta, about = 0)
[1] 306061
> gigRawMom(2.5, Theta)
[1] 3.782698
> momIntegrated("gig", order = 2.5, param = Theta, about = 0, absolute = TRUE)
[1] 3.782698
> 
> ### Test gigMom
> Theta <- c(-0.5,5,2.5)
> m1 <- gigRawMom(1, Theta)
> m1
[1] 1.414214
> gigMom(1, Theta)
[1] 1.414214
> gigMom(2, Theta, m1)
[1] 0.5656854
> m2 <- momIntegrated("gig", order = 2, param = Theta, about = m1)
> m2
[1] 0.5656854
> gigMom(1, Theta, m1)
[1] 0
> gigVar(Theta)
[1] 0.5656854
> gigMom(3, Theta, m1)
[1] 0.6788225
> momIntegrated("gig", order = 3, param = Theta, about = m1)
[1] 0.6788225
> #gigMom(2.5, Theta, m1)
> 
> ### Test gammaRawMom
> shape <- 2
> rate <- 3
> Theta <- c(shape, rate)
> gammaRawMom(1, shape, rate)
[1] 0.6666667
> momIntegrated("gamma", order = 1, param = Theta, about = 0)
[1] 0.6666667
> gammaRawMom(2, shape, rate)
[1] 0.6666667
> momIntegrated("gamma", order = 2, param = Theta, about = 0)
[1] 0.6666667
> gammaRawMom(10, shape, rate)
[1] 675.9945
> momIntegrated("gamma", order = 10, param = Theta, about = 0)
[1] 675.9945
> 
> ### Test gigMom for gamma case
> shape <- 2
> rate <- 3
> Theta <- c(shape, rate)
> m1 <- gammaRawMom(1, shape, rate)
> m1
[1] 0.6666667
> gigMom(1, c(shape, 0, 2*rate))
[1] 0.6666667
> gigMom(2, c(shape, 0, 2*rate), about = m1)
[1] 0.2222222
> m2 <- momIntegrated("gamma", order = 2, param = Theta, about = m1)
> m2
[1] 0.2222222
> gigMom(1, c(shape, 0, 2*rate), about = m1)
[1] 0
> gigMom(3, c(shape, 0, 2*rate), about = m1)
[1] 0.1481481
> momIntegrated("gamma", order = 3, param = Theta, about = m1)
[1] 0.1481481
> 
> ### Test gigRawMom for gamma case
> library(actuar)

Attaching package: 'actuar'

The following objects are masked from 'package:stats':

    sd, var

The following object is masked from 'package:grDevices':

    cm

> Theta <- c(0.5,0,2.5)
> gigRawMom(2, Theta)
[1] 0.48
> mgamma(2, Theta[1], Theta[3]/2)
[1] 0.48
> momIntegrated("gamma", order = 2,
+               param = c(Theta[1],Theta[3]/2), about = 0)
[1] 0.48
> ### Infinite moment
> gigRawMom(-2, Theta)
[1] Inf
> mgamma(-2, Theta[1], Theta[3]/2)
[1] Inf
> #momIntegrated("gamma", order = -2,
> #              param = c(Theta[1],Theta[3]/2), about = 0)
> ### Test gigRawMom for inverse gamma case
> Theta <- c(-0.5,5,0)
> ### Infinite moment
> gigRawMom(2, Theta)
[1] Inf
> minvgamma(2, -Theta[1], Theta[2]/2)
[1] Inf
> #momIntegrated("invgamma", order = 2,
> #              param = c(-Theta[1],Theta[2]/2), about = 0)
> gigRawMom(-2, Theta)
[1] 4.6875
> minvgamma(-2, -Theta[1], Theta[2]/2)
[1] 4.6875
> momIntegrated("invgamma", order = -2,
+               param = c(-Theta[1],Theta[2]/2), about = 0)
[1] 4.6875
> 
> ### Test gigMom for gamma case
> Theta <- c(0.5,0,2.5)
> m1 <- gigMom(1, Theta)
> m1
[1] 0.4
> mgamma(1, Theta[1], Theta[3]/2)
[1] 0.4
> momIntegrated("gamma", order = 1,
+               param = c(Theta[1],Theta[3]/2), about = 0)
[1] 0.4
> gigMom(2, Theta, m1)
[1] 0.32
> momIntegrated("gamma", order = 2,
+               param = c(Theta[1],Theta[3]/2), about = m1)
[1] 0.32
> ### Infinite moments
> gigMom(-2, Theta)
[1] Inf
> #momIntegrated("gamma", order = -2,
> #              param = c(Theta[1],Theta[3]/2), about = 0)
> #gigMom(-2, Theta, m1)
> #momIntegrated("gamma", order = -2,
> #              param = c(Theta[1],Theta[3]/2), about = m1)
> 
> ### Test gigRawMom for inverse gamma case
> Theta <- c(-0.5,5,0)
> ### Infinite moments
> m1 <- gigMom(1, Theta)
> m1
[1] Inf
> minvgamma(1, -Theta[1], Theta[2]/2)
[1] Inf
> #momIntegrated("invgamma", order = 1,
> #              param = c(-Theta[1],Theta[2]/2), about = 0)
> gigMom(2, Theta, m1)
[1] NaN
> #momIntegrated("invgamma", order = 2,
> #              param = c(-Theta[1],Theta[2]/2), about = m1)
> 
> gigMom(-2, Theta)
[1] 4.6875
> momIntegrated("invgamma", order = -2,
+               param = c(-Theta[1],Theta[2]/2), about = 0)
[1] 4.6875
> #gigMom(-2, Theta, m1)
> momIntegrated("invgamma", order = -2,
+               param = c(-Theta[1],Theta[2]/2), about = m1)
[1] 0
> 
> 
> 
> proc.time()
   user  system elapsed 
   0.14    0.06    0.18