Last updated on 2024-03-28 17:49:29 CET.
Flavor | Version | Tinstall | Tcheck | Ttotal | Status | Flags |
---|---|---|---|---|---|---|
r-devel-linux-x86_64-debian-clang | 0.8 | 3.04 | 23.96 | 27.00 | OK | |
r-devel-linux-x86_64-debian-gcc | 0.8 | 2.17 | 18.89 | 21.06 | OK | |
r-devel-linux-x86_64-fedora-clang | 0.8 | 34.15 | NOTE | |||
r-devel-linux-x86_64-fedora-gcc | 0.8 | 32.93 | OK | |||
r-devel-windows-x86_64 | 0.8 | 6.00 | 215.00 | 221.00 | ERROR | |
r-patched-linux-x86_64 | 0.8 | 2.76 | 23.51 | 26.27 | OK | |
r-release-linux-x86_64 | 0.8 | 2.73 | 23.30 | 26.03 | OK | |
r-release-macos-arm64 | 0.8 | 17.00 | OK | |||
r-release-macos-x86_64 | 0.8 | 23.00 | OK | |||
r-release-windows-x86_64 | 0.8 | 8.00 | 44.00 | 52.00 | OK | |
r-oldrel-macos-arm64 | 0.8 | 20.00 | OK | |||
r-oldrel-windows-x86_64 | 0.8 | 16.00 | 49.00 | 65.00 | OK |
Version: 0.8
Check: Rd cross-references
Result: NOTE
Undeclared package ‘Runuran’ in Rd xrefs
Flavor: r-devel-linux-x86_64-fedora-clang
Version: 0.8
Check: tests
Result: ERROR
Running 'examples.R' [172s]
Running the tests in 'tests/examples.R' failed.
Complete output:
>
> library(GIGrvg)
>
> ## --- auxiliary functions ---------------------------------------------------
>
> ## check whether error has been thrown
> iserror <- function (expr) { is(try(expr), "try-error") }
>
> ## mode of GIG distribution
> gig.mode <- function (lambda, chi, psi) {
+ omega = sqrt(psi*chi)
+
+ ## lambda can be < 0 !!
+ if (lambda >= 1.)
+ mode <- (sqrt((lambda-1)*(lambda-1) + psi*chi)+(lambda-1)) / psi
+ else
+ mode <- chi / (sqrt((1-lambda)*(1-lambda) + psi*chi)+(1-lambda))
+
+ mode
+ }
>
> ## check GIG parameters
> check.params <- function (lambda, chi, psi) {
+ valid <- FALSE
+ if ((chi >= 0 && psi >= 0) &&
+ (chi != 0 || lambda > 0) &&
+ (psi != 0 || lambda < 0)) {
+ valid <- TRUE
+ }
+ valid
+ }
>
> ## show warnings immediately
> options(warn=1)
>
> ##############################################################################
> ##
> ## dgig()
> ##
> ##############################################################################
>
> ## --- check handling of invalid parameters ----------------------------------
>
> xvals <- c(1,NA,NaN)
> lambda <- c(1,NA,NaN,Inf)
> chi <- c(1, NA, NaN, Inf)
> psi <- c(1, NA, NaN, Inf)
>
> for (x in xvals) { for (l in lambda) { for (c in chi) { for (p in psi) {
+ f <- dgig(x,l,c,p)
+ lf <- dgig(x,l,c,p,TRUE)
+ a <- c(x,l,c,p)
+ if (length(a[!is.finite(a)])>0 && !(is.na(f) && is.na(lf))) {
+ cat ("f=",f,"; logf=",lf,"\t",
+ "; x=",x,"; lambda=",l,"; chi=",c,"; psi=",p,"\n", sep="")
+ stop ("NA or NaN expected")
+ }
+ }}}}
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=1, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=1, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=1, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=1, psi=nan
Warning in dgig(x, l, c, p) : NaNs produced
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=1, psi=inf
Warning in dgig(x, l, c, p, TRUE) : NaNs produced
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=1, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=inf
Warning in dgig(x, l, c, p) : NaNs produced
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=inf, psi=1
Warning in dgig(x, l, c, p, TRUE) : NaNs produced
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=inf, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=inf, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=inf, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=inf, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=inf, psi=nan
Warning in dgig(x, l, c, p) : NaNs produced
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=inf, psi=inf
Warning in dgig(x, l, c, p, TRUE) : NaNs produced
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=inf, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=inf
Warning in dgig(x, l, c, p) : NaNs produced
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=1, psi=1
Warning in dgig(x, l, c, p, TRUE) : NaNs produced
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=1, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=1, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=1, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=1, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=1, psi=nan
Warning in dgig(x, l, c, p) : NaNs produced
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=1, psi=inf
Warning in dgig(x, l, c, p, TRUE) : NaNs produced
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=1, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=inf
Warning in dgig(x, l, c, p) : NaNs produced
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=inf, psi=1
Warning in dgig(x, l, c, p, TRUE) : NaNs produced
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=inf, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=inf, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=inf, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=inf, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=inf, psi=nan
Warning in dgig(x, l, c, p) : NaNs produced
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=inf, psi=inf
Warning in dgig(x, l, c, p, TRUE) : NaNs produced
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=inf, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=1, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=1, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=1, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=1, psi=nan
Warning in dgig(x, l, c, p) : NaNs produced
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=1, psi=inf
Warning in dgig(x, l, c, p, TRUE) : NaNs produced
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=1, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=inf
Warning in dgig(x, l, c, p) : NaNs produced
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=inf, psi=1
Warning in dgig(x, l, c, p, TRUE) : NaNs produced
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=inf, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=inf, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=inf, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=inf, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=inf, psi=nan
Warning in dgig(x, l, c, p) : NaNs produced
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=inf, psi=inf
Warning in dgig(x, l, c, p, TRUE) : NaNs produced
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=inf, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=inf
Warning in dgig(x, l, c, p) : NaNs produced
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=1, psi=1
Warning in dgig(x, l, c, p, TRUE) : NaNs produced
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=1, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=1, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=1, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=1, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=1, psi=nan
Warning in dgig(x, l, c, p) : NaNs produced
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=1, psi=inf
Warning in dgig(x, l, c, p, TRUE) : NaNs produced
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=1, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=inf
Warning in dgig(x, l, c, p) : NaNs produced
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=inf, psi=1
Warning in dgig(x, l, c, p, TRUE) : NaNs produced
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=inf, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=inf, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=inf, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=inf, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=inf, psi=nan
Warning in dgig(x, l, c, p) : NaNs produced
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=inf, psi=inf
Warning in dgig(x, l, c, p, TRUE) : NaNs produced
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=inf, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=1, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=1, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=1, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=1, psi=nan
Warning in dgig(x, l, c, p) : NaNs produced
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=1, psi=inf
Warning in dgig(x, l, c, p, TRUE) : NaNs produced
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=1, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=nan, psi=inf
Warning in dgig(x, l, c, p) : NaNs produced
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=inf, psi=1
Warning in dgig(x, l, c, p, TRUE) : NaNs produced
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=inf, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=inf, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=inf, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=inf, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=inf, psi=nan
Warning in dgig(x, l, c, p) : NaNs produced
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=1, chi=inf, psi=inf
Warning in dgig(x, l, c, p, TRUE) : NaNs produced
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=1, chi=inf, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=1, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=nan, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=nan, chi=inf, psi=inf
Warning in dgig(x, l, c, p) : NaNs produced
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=1, psi=1
Warning in dgig(x, l, c, p, TRUE) : NaNs produced
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=1, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=1, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=1, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=1, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=1, psi=nan
Warning in dgig(x, l, c, p) : NaNs produced
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=1, psi=inf
Warning in dgig(x, l, c, p, TRUE) : NaNs produced
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=1, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=inf
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=1
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=inf
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=nan, psi=inf
Warning in dgig(x, l, c, p) : NaNs produced
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=inf, psi=1
Warning in dgig(x, l, c, p, TRUE) : NaNs produced
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=inf, psi=1
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=inf, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=inf, psi=nan
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=inf, psi=nan
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=inf, psi=nan
Warning in dgig(x, l, c, p) : NaNs produced
Warning in dgig(x, l, c, p) :
invalid parameters for GIG distribution: lambda=inf, chi=inf, psi=inf
Warning in dgig(x, l, c, p, TRUE) : NaNs produced
Warning in dgig(x, l, c, p, TRUE) :
invalid parameters for GIG distribution: lambda=inf, chi=inf, psi=inf
>
> if (! (is.nan(dgig(1, -1,-1, 1)) &&
+ is.nan(dgig(1, -1, 0, 1)) &&
+ is.nan(dgig(1, -1, 1,-1)) &&
+ is.nan(dgig(1, 0,-1, 1)) &&
+ is.nan(dgig(1, 0, 1, 0)) &&
+ is.nan(dgig(1, 0, 0, 1)) &&
+ is.nan(dgig(1, 0, 1,-1)) &&
+ is.nan(dgig(1, 1,-1, 1)) &&
+ is.nan(dgig(1, 1, 1, 0)) &&
+ is.nan(dgig(1, 1, 1,-1)) ) ) {
+ stop ("NA or NaN expected")
+ }
Warning in dgig(1, -1, -1, 1) : NaNs produced
Warning in dgig(1, -1, -1, 1) :
invalid parameters for GIG distribution: lambda=-1, chi=-1, psi=1
Warning in dgig(1, -1, 0, 1) : NaNs produced
Warning in dgig(1, -1, 0, 1) :
invalid parameters for GIG distribution: lambda=-1, chi=0, psi=1
Warning in dgig(1, -1, 1, -1) : NaNs produced
Warning in dgig(1, -1, 1, -1) :
invalid parameters for GIG distribution: lambda=-1, chi=1, psi=-1
Warning in dgig(1, 0, -1, 1) : NaNs produced
Warning in dgig(1, 0, -1, 1) :
invalid parameters for GIG distribution: lambda=0, chi=-1, psi=1
Warning in dgig(1, 0, 1, 0) : NaNs produced
Warning in dgig(1, 0, 1, 0) :
invalid parameters for GIG distribution: lambda=0, chi=1, psi=0
Warning in dgig(1, 0, 0, 1) : NaNs produced
Warning in dgig(1, 0, 0, 1) :
invalid parameters for GIG distribution: lambda=0, chi=0, psi=1
Warning in dgig(1, 0, 1, -1) : NaNs produced
Warning in dgig(1, 0, 1, -1) :
invalid parameters for GIG distribution: lambda=0, chi=1, psi=-1
Warning in dgig(1, 1, -1, 1) : NaNs produced
Warning in dgig(1, 1, -1, 1) :
invalid parameters for GIG distribution: lambda=1, chi=-1, psi=1
Warning in dgig(1, 1, 1, 0) : NaNs produced
Warning in dgig(1, 1, 1, 0) :
invalid parameters for GIG distribution: lambda=1, chi=1, psi=0
Warning in dgig(1, 1, 1, -1) : NaNs produced
Warning in dgig(1, 1, 1, -1) :
invalid parameters for GIG distribution: lambda=1, chi=1, psi=-1
>
> ## --- verify finite (non-NaN, non-Inf) results: evaluate at mode ------------
>
> lambda <- c(-100, -10, -1, -0.1, 0, 0.1, 1, 10, 100)
> chi <- c(0, 0.01, 0.1, 1, 10, 100)
> psi <- c(0, 0.01, 0.1, 1, 10, 100)
>
> for (l in lambda) { for (c in chi) { for (p in psi) {
+ if (!check.params(l,c,p)) next
+ m <- gig.mode(l,c,p)
+ if (m==0.0) m <- 1
+ f <- signif(dgig(m,l,c,p),5)
+ lf <- signif(dgig(m,l,c,p,TRUE),5)
+ cat ("f=",f,"; logf=",lf,"\t",
+ "; lambda=",l,"; chi=",c,"; psi=",p,";mode=",m,"\n", sep="")
+ if (! isTRUE(is.finite(f)) || ! isTRUE(is.finite(lf)))
+ stop ("result not finite")
+ }}}
f=80120; logf=11.291 ; lambda=-100; chi=0.01; psi=0;mode=4.950495e-05
f=80120; logf=11.291 ; lambda=-100; chi=0.01; psi=0.01;mode=4.950495e-05
f=80120; logf=11.291 ; lambda=-100; chi=0.01; psi=0.1;mode=4.950495e-05
f=80120; logf=11.291 ; lambda=-100; chi=0.01; psi=1;mode=4.950494e-05
f=80121; logf=11.291 ; lambda=-100; chi=0.01; psi=10;mode=4.950483e-05
f=80124; logf=11.291 ; lambda=-100; chi=0.01; psi=100;mode=4.950374e-05
f=8012; logf=8.9887 ; lambda=-100; chi=0.1; psi=0;mode=0.0004950495
f=8012; logf=8.9887 ; lambda=-100; chi=0.1; psi=0.01;mode=0.0004950495
f=8012; logf=8.9887 ; lambda=-100; chi=0.1; psi=0.1;mode=0.0004950494
f=8012.1; logf=8.9887 ; lambda=-100; chi=0.1; psi=1;mode=0.0004950483
f=8012.4; logf=8.9887 ; lambda=-100; chi=0.1; psi=10;mode=0.0004950374
f=8016; logf=8.9892 ; lambda=-100; chi=0.1; psi=100;mode=0.0004949282
f=801.2; logf=6.6861 ; lambda=-100; chi=1; psi=0;mode=0.004950495
f=801.2; logf=6.6861 ; lambda=-100; chi=1; psi=0.01;mode=0.004950494
f=801.21; logf=6.6861 ; lambda=-100; chi=1; psi=0.1;mode=0.004950483
f=801.24; logf=6.6862 ; lambda=-100; chi=1; psi=1;mode=0.004950374
f=801.6; logf=6.6866 ; lambda=-100; chi=1; psi=10;mode=0.004949282
f=805.2; logf=6.6911 ; lambda=-100; chi=1; psi=100;mode=0.004938422
f=80.12; logf=4.3835 ; lambda=-100; chi=10; psi=0;mode=0.04950495
f=80.121; logf=4.3835 ; lambda=-100; chi=10; psi=0.01;mode=0.04950483
f=80.124; logf=4.3836 ; lambda=-100; chi=10; psi=0.1;mode=0.04950374
f=80.16; logf=4.384 ; lambda=-100; chi=10; psi=1;mode=0.04949282
f=80.52; logf=4.3885 ; lambda=-100; chi=10; psi=10;mode=0.04938422
f=84.054; logf=4.4315 ; lambda=-100; chi=10; psi=100;mode=0.04834777
f=8.012; logf=2.0809 ; lambda=-100; chi=100; psi=0;mode=0.4950495
f=8.0124; logf=2.081 ; lambda=-100; chi=100; psi=0.01;mode=0.4950374
f=8.016; logf=2.0814 ; lambda=-100; chi=100; psi=0.1;mode=0.4949282
f=8.052; logf=2.0859 ; lambda=-100; chi=100; psi=1;mode=0.4938422
f=8.4054; logf=2.1289 ; lambda=-100; chi=100; psi=10;mode=0.4834777
f=11.496; logf=2.442 ; lambda=-100; chi=100; psi=100;mode=0.4113022
f=2626.3; logf=7.8733 ; lambda=-10; chi=0.01; psi=0;mode=0.0004545455
f=2626.3; logf=7.8733 ; lambda=-10; chi=0.01; psi=0.01;mode=0.0004545454
f=2626.3; logf=7.8733 ; lambda=-10; chi=0.01; psi=0.1;mode=0.0004545445
f=2626.4; logf=7.8734 ; lambda=-10; chi=0.01; psi=1;mode=0.0004545361
f=2627.6; logf=7.8738 ; lambda=-10; chi=0.01; psi=10;mode=0.0004544516
f=2639.6; logf=7.8784 ; lambda=-10; chi=0.01; psi=100;mode=0.0004536102
f=262.63; logf=5.5708 ; lambda=-10; chi=0.1; psi=0;mode=0.004545455
f=262.63; logf=5.5708 ; lambda=-10; chi=0.1; psi=0.01;mode=0.004545445
f=262.64; logf=5.5708 ; lambda=-10; chi=0.1; psi=0.1;mode=0.004545361
f=262.76; logf=5.5713 ; lambda=-10; chi=0.1; psi=1;mode=0.004544516
f=263.96; logf=5.5758 ; lambda=-10; chi=0.1; psi=10;mode=0.004536102
f=275.6; logf=5.619 ; lambda=-10; chi=0.1; psi=100;mode=0.004455231
f=26.263; logf=3.2682 ; lambda=-10; chi=1; psi=0;mode=0.04545455
f=26.264; logf=3.2682 ; lambda=-10; chi=1; psi=0.01;mode=0.04545361
f=26.276; logf=3.2687 ; lambda=-10; chi=1; psi=0.1;mode=0.04544516
f=26.396; logf=3.2732 ; lambda=-10; chi=1; psi=1;mode=0.04536102
f=27.56; logf=3.3164 ; lambda=-10; chi=1; psi=10;mode=0.04455231
f=37.565; logf=3.6261 ; lambda=-10; chi=1; psi=100;mode=0.03866069
f=2.6263; logf=0.96558 ; lambda=-10; chi=10; psi=0;mode=0.4545455
f=2.6276; logf=0.96609 ; lambda=-10; chi=10; psi=0.01;mode=0.4544516
f=2.6396; logf=0.97061 ; lambda=-10; chi=10; psi=0.1;mode=0.4536102
f=2.756; logf=1.0138 ; lambda=-10; chi=10; psi=1;mode=0.4455231
f=3.7565; logf=1.3235 ; lambda=-10; chi=10; psi=10;mode=0.3866069
f=10.1; logf=2.3125 ; lambda=-10; chi=10; psi=100;mode=0.2248134
f=0.26263; logf=-1.337 ; lambda=-10; chi=100; psi=0;mode=4.545455
f=0.26396; logf=-1.332 ; lambda=-10; chi=100; psi=0.01;mode=4.536102
f=0.2756; logf=-1.2888 ; lambda=-10; chi=100; psi=0.1;mode=4.455231
f=0.37565; logf=-0.9791 ; lambda=-10; chi=100; psi=1;mode=3.866069
f=1.01; logf=0.0099154 ; lambda=-10; chi=100; psi=10;mode=2.248134
f=4.4467; logf=1.4922 ; lambda=-10; chi=100; psi=100;mode=0.8960318
f=108.27; logf=4.6846 ; lambda=-1; chi=0.01; psi=0;mode=0.0025
f=108.3; logf=4.6849 ; lambda=-1; chi=0.01; psi=0.01;mode=0.002499984
f=108.48; logf=4.6865 ; lambda=-1; chi=0.01; psi=0.1;mode=0.002499844
f=109.74; logf=4.6981 ; lambda=-1; chi=0.01; psi=1;mode=0.002498439
f=117.5; logf=4.7664 ; lambda=-1; chi=0.01; psi=10;mode=0.002484567
f=159.31; logf=5.0709 ; lambda=-1; chi=0.01; psi=100;mode=0.00236068
f=10.827; logf=2.382 ; lambda=-1; chi=0.1; psi=0;mode=0.025
f=10.848; logf=2.3839 ; lambda=-1; chi=0.1; psi=0.01;mode=0.02499844
f=10.974; logf=2.3955 ; lambda=-1; chi=0.1; psi=0.1;mode=0.02498439
f=11.75; logf=2.4638 ; lambda=-1; chi=0.1; psi=1;mode=0.02484567
f=15.931; logf=2.7683 ; lambda=-1; chi=0.1; psi=10;mode=0.0236068
f=37.356; logf=3.6205 ; lambda=-1; chi=0.1; psi=100;mode=0.01741657
f=1.0827; logf=0.079442 ; lambda=-1; chi=1; psi=0;mode=0.25
f=1.0974; logf=0.092915 ; lambda=-1; chi=1; psi=0.01;mode=0.2498439
f=1.175; logf=0.16123 ; lambda=-1; chi=1; psi=0.1;mode=0.2484567
f=1.5931; logf=0.46571 ; lambda=-1; chi=1; psi=1;mode=0.236068
f=3.7356; logf=1.3179 ; lambda=-1; chi=1; psi=10;mode=0.1741657
f=14.858; logf=2.6985 ; lambda=-1; chi=1; psi=100;mode=0.08198039
f=0.10827; logf=-2.2231 ; lambda=-1; chi=10; psi=0;mode=2.5
f=0.1175; logf=-2.1414 ; lambda=-1; chi=10; psi=0.01;mode=2.484567
f=0.15931; logf=-1.8369 ; lambda=-1; chi=10; psi=0.1;mode=2.36068
f=0.37356; logf=-0.98468 ; lambda=-1; chi=10; psi=1;mode=1.741657
f=1.4858; logf=0.39592 ; lambda=-1; chi=10; psi=10;mode=0.8198039
f=7.4696; logf=2.0108 ; lambda=-1; chi=10; psi=100;mode=0.2968596
f=0.010827; logf=-4.5257 ; lambda=-1; chi=100; psi=0;mode=25
f=0.015931; logf=-4.1395 ; lambda=-1; chi=100; psi=0.01;mode=23.6068
f=0.037356; logf=-3.2873 ; lambda=-1; chi=100; psi=0.1;mode=17.41657
f=0.14858; logf=-1.9067 ; lambda=-1; chi=100; psi=1;mode=8.198039
f=0.74696; logf=-0.29174 ; lambda=-1; chi=100; psi=10;mode=2.968596
f=4.0549; logf=1.3999 ; lambda=-1; chi=100; psi=100;mode=0.9802
f=7.7714; logf=2.0504 ; lambda=-0.1; chi=0.01; psi=0;mode=0.004545455
f=12.725; logf=2.5435 ; lambda=-0.1; chi=0.01; psi=0.01;mode=0.004545361
f=15.232; logf=2.7234 ; lambda=-0.1; chi=0.01; psi=0.1;mode=0.004544516
f=20.172; logf=3.0043 ; lambda=-0.1; chi=0.01; psi=1;mode=0.004536102
f=32.548; logf=3.4827 ; lambda=-0.1; chi=0.01; psi=10;mode=0.004455231
f=76.112; logf=4.3322 ; lambda=-0.1; chi=0.01; psi=100;mode=0.003866069
f=0.77714; logf=-0.25214 ; lambda=-0.1; chi=0.1; psi=0;mode=0.04545455
f=1.5232; logf=0.42084 ; lambda=-0.1; chi=0.1; psi=0.01;mode=0.04544516
f=2.0172; logf=0.7017 ; lambda=-0.1; chi=0.1; psi=0.1;mode=0.04536102
f=3.2548; logf=1.1801 ; lambda=-0.1; chi=0.1; psi=1;mode=0.04455231
f=7.6112; logf=2.0296 ; lambda=-0.1; chi=0.1; psi=10;mode=0.03866069
f=28.032; logf=3.3334 ; lambda=-0.1; chi=0.1; psi=100;mode=0.02248134
f=0.077714; logf=-2.5547 ; lambda=-0.1; chi=1; psi=0;mode=0.4545455
f=0.20172; logf=-1.6009 ; lambda=-0.1; chi=1; psi=0.01;mode=0.4536102
f=0.32548; logf=-1.1225 ; lambda=-0.1; chi=1; psi=0.1;mode=0.4455231
f=0.76112; logf=-0.27296 ; lambda=-0.1; chi=1; psi=1;mode=0.3866069
f=2.8032; logf=1.0308 ; lambda=-0.1; chi=1; psi=10;mode=0.2248134
f=13.556; logf=2.6068 ; lambda=-0.1; chi=1; psi=100;mode=0.08960318
f=0.0077714; logf=-4.8573 ; lambda=-0.1; chi=10; psi=0;mode=4.545455
f=0.032548; logf=-3.425 ; lambda=-0.1; chi=10; psi=0.01;mode=4.455231
f=0.076112; logf=-2.5755 ; lambda=-0.1; chi=10; psi=0.1;mode=3.866069
f=0.28032; logf=-1.2718 ; lambda=-0.1; chi=10; psi=1;mode=2.248134
f=1.3556; logf=0.30425 ; lambda=-0.1; chi=10; psi=10;mode=0.8960318
f=7.2584; logf=1.9822 ; lambda=-0.1; chi=10; psi=100;mode=0.305419
f=0.00077714; logf=-7.1599 ; lambda=-0.1; chi=100; psi=0;mode=45.45455
f=0.0076112; logf=-4.8781 ; lambda=-0.1; chi=100; psi=0.01;mode=38.66069
f=0.028032; logf=-3.5744 ; lambda=-0.1; chi=100; psi=0.1;mode=22.48134
f=0.13556; logf=-1.9983 ; lambda=-0.1; chi=100; psi=1;mode=8.960318
f=0.72584; logf=-0.32043 ; lambda=-0.1; chi=100; psi=10;mode=3.05419
f=4.0184; logf=1.3909 ; lambda=-0.1; chi=100; psi=100;mode=0.9890605
f=7.7918; logf=2.0531 ; lambda=0; chi=0.01; psi=0.01;mode=0.004999875
f=10.299; logf=2.3321 ; lambda=0; chi=0.01; psi=0.1;mode=0.004998751
f=15.12; logf=2.716 ; lambda=0; chi=0.01; psi=1;mode=0.004987562
f=27.101; logf=3.2996 ; lambda=0; chi=0.01; psi=10;mode=0.004880885
f=69.703; logf=4.2442 ; lambda=0; chi=0.01; psi=100;mode=0.004142136
f=1.0299; logf=0.029503 ; lambda=0; chi=0.1; psi=0.01;mode=0.04998751
f=1.512; logf=0.4134 ; lambda=0; chi=0.1; psi=0.1;mode=0.04987562
f=2.7101; logf=0.99697 ; lambda=0; chi=0.1; psi=1;mode=0.04880885
f=6.9703; logf=1.9417 ; lambda=0; chi=0.1; psi=10;mode=0.04142136
f=27.17; logf=3.3021 ; lambda=0; chi=0.1; psi=100;mode=0.02316625
f=0.1512; logf=-1.8892 ; lambda=0; chi=1; psi=0.01;mode=0.4987562
f=0.27101; logf=-1.3056 ; lambda=0; chi=1; psi=0.1;mode=0.4880885
f=0.69703; logf=-0.36092 ; lambda=0; chi=1; psi=1;mode=0.4142136
f=2.717; logf=0.99954 ; lambda=0; chi=1; psi=10;mode=0.2316625
f=13.421; logf=2.5968 ; lambda=0; chi=1; psi=100;mode=0.09049876
f=0.027101; logf=-3.6082 ; lambda=0; chi=10; psi=0.01;mode=4.880885
f=0.069703; logf=-2.6635 ; lambda=0; chi=10; psi=0.1;mode=4.142136
f=0.2717; logf=-1.303 ; lambda=0; chi=10; psi=1;mode=2.316625
f=1.3421; logf=0.29424 ; lambda=0; chi=10; psi=10;mode=0.9049876
f=7.2355; logf=1.979 ; lambda=0; chi=10; psi=100;mode=0.3063858
f=0.0069703; logf=-4.9661 ; lambda=0; chi=100; psi=0.01;mode=41.42136
f=0.02717; logf=-3.6056 ; lambda=0; chi=100; psi=0.1;mode=23.16625
f=0.13421; logf=-2.0083 ; lambda=0; chi=100; psi=1;mode=9.049876
f=0.72355; logf=-0.32359 ; lambda=0; chi=100; psi=10;mode=3.063858
f=4.0144; logf=1.3899 ; lambda=0; chi=100; psi=100;mode=0.99005
f=0.061572; logf=-2.7875 ; lambda=0.1; chi=0; psi=0.01;mode=1
f=0.074104; logf=-2.6023 ; lambda=0.1; chi=0; psi=0.1;mode=1
f=0.059485; logf=-2.822 ; lambda=0.1; chi=0; psi=1;mode=1
f=0.00083192; logf=-7.0918 ; lambda=0.1; chi=0; psi=10;mode=1
f=2.998e-23; logf=-51.862 ; lambda=0.1; chi=0; psi=100;mode=1
f=4.4114; logf=1.4842 ; lambda=0.1; chi=0.01; psi=0.01;mode=0.005555384
f=6.6479; logf=1.8943 ; lambda=0.1; chi=0.01; psi=0.1;mode=0.005553842
f=11.078; logf=2.405 ; lambda=0.1; chi=0.01; psi=1;mode=0.005538514
f=22.405; logf=3.1093 ; lambda=0.1; chi=0.01; psi=10;mode=0.00539392
f=63.819; logf=4.1561 ; lambda=0.1; chi=0.01; psi=100;mode=0.004453624
f=0.66479; logf=-0.40828 ; lambda=0.1; chi=0.1; psi=0.01;mode=0.05553842
f=1.1078; logf=0.10239 ; lambda=0.1; chi=0.1; psi=0.1;mode=0.05538514
f=2.2405; logf=0.80671 ; lambda=0.1; chi=0.1; psi=1;mode=0.0539392
f=6.3819; logf=1.8535 ; lambda=0.1; chi=0.1; psi=10;mode=0.04453624
f=26.341; logf=3.2711 ; lambda=0.1; chi=0.1; psi=100;mode=0.02387856
f=0.11078; logf=-2.2002 ; lambda=0.1; chi=1; psi=0.01;mode=0.5538514
f=0.22405; logf=-1.4959 ; lambda=0.1; chi=1; psi=0.1;mode=0.539392
f=0.63819; logf=-0.44912 ; lambda=0.1; chi=1; psi=1;mode=0.4453624
f=2.6341; logf=0.96854 ; lambda=0.1; chi=1; psi=10;mode=0.2387856
f=13.288; logf=2.5869 ; lambda=0.1; chi=1; psi=100;mode=0.09140418
f=0.022405; logf=-3.7985 ; lambda=0.1; chi=10; psi=0.01;mode=5.39392
f=0.063819; logf=-2.7517 ; lambda=0.1; chi=10; psi=0.1;mode=4.453624
f=0.26341; logf=-1.334 ; lambda=0.1; chi=10; psi=1;mode=2.387856
f=1.3288; logf=0.28428 ; lambda=0.1; chi=10; psi=10;mode=0.9140418
f=7.2126; logf=1.9758 ; lambda=0.1; chi=10; psi=100;mode=0.3073558
f=0.0063819; logf=-5.0543 ; lambda=0.1; chi=100; psi=0.01;mode=44.53624
f=0.026341; logf=-3.6366 ; lambda=0.1; chi=100; psi=0.1;mode=23.87856
f=0.13288; logf=-2.0183 ; lambda=0.1; chi=100; psi=1;mode=9.140418
f=0.72126; logf=-0.32675 ; lambda=0.1; chi=100; psi=10;mode=3.073558
f=4.0104; logf=1.3889 ; lambda=0.1; chi=100; psi=100;mode=0.9910405
f=0.0049751; logf=-5.3033 ; lambda=1; chi=0; psi=0.01;mode=1
f=0.047561; logf=-3.0457 ; lambda=1; chi=0; psi=0.1;mode=1
f=0.30327; logf=-1.1931 ; lambda=1; chi=0; psi=1;mode=1
f=0.03369; logf=-3.3906 ; lambda=1; chi=0; psi=10;mode=1
f=9.6437e-21; logf=-46.088 ; lambda=1; chi=0; psi=100;mode=1
f=0.0049515; logf=-5.3081 ; lambda=1; chi=0.01; psi=0.01;mode=1
f=0.048542; logf=-3.0253 ; lambda=1; chi=0.01; psi=0.1;mode=0.3162278
f=0.45913; logf=-0.77842 ; lambda=1; chi=0.01; psi=1;mode=0.1
f=4.0047; logf=1.3875 ; lambda=1; chi=0.01; psi=10;mode=0.03162278
f=30.559; logf=3.4197 ; lambda=1; chi=0.01; psi=100;mode=0.01
f=0.0048542; logf=-5.3279 ; lambda=1; chi=0.1; psi=0.01;mode=3.162278
f=0.045913; logf=-3.081 ; lambda=1; chi=0.1; psi=0.1;mode=1
f=0.40047; logf=-0.91513 ; lambda=1; chi=0.1; psi=1;mode=0.3162278
f=3.0559; logf=1.1171 ; lambda=1; chi=0.1; psi=10;mode=0.1
f=20.226; logf=3.007 ; lambda=1; chi=0.1; psi=100;mode=0.03162278
f=0.0045913; logf=-5.3836 ; lambda=1; chi=1; psi=0.01;mode=10
f=0.040047; logf=-3.2177 ; lambda=1; chi=1; psi=0.1;mode=3.162278
f=0.30559; logf=-1.1855 ; lambda=1; chi=1; psi=1;mode=1
f=2.0226; logf=0.70438 ; lambda=1; chi=1; psi=10;mode=0.3162278
f=12.172; logf=2.4992 ; lambda=1; chi=1; psi=100;mode=0.1
f=0.0040047; logf=-5.5203 ; lambda=1; chi=10; psi=0.01;mode=31.62278
f=0.030559; logf=-3.4881 ; lambda=1; chi=10; psi=0.1;mode=10
f=0.20226; logf=-1.5982 ; lambda=1; chi=10; psi=1;mode=3.162278
f=1.2172; logf=0.19658 ; lambda=1; chi=10; psi=10;mode=1
f=7.012; logf=1.9476 ; lambda=1; chi=10; psi=100;mode=0.3162278
f=0.0030559; logf=-5.7907 ; lambda=1; chi=100; psi=0.01;mode=100
f=0.020226; logf=-3.9008 ; lambda=1; chi=100; psi=0.1;mode=31.62278
f=0.12172; logf=-2.106 ; lambda=1; chi=100; psi=1;mode=10
f=0.7012; logf=-0.35497 ; lambda=1; chi=100; psi=10;mode=3.162278
f=3.9746; logf=1.3799 ; lambda=1; chi=100; psi=100;mode=1
f=0.00065878; logf=-7.3251 ; lambda=10; chi=0; psi=0.01;mode=1800
f=0.0065878; logf=-5.0225 ; lambda=10; chi=0; psi=0.1;mode=180
f=0.065878; logf=-2.72 ; lambda=10; chi=0; psi=1;mode=18
f=0.65878; logf=-0.41737 ; lambda=10; chi=0; psi=10;mode=1.8
f=6.5878; logf=1.8852 ; lambda=10; chi=0; psi=100;mode=0.18
f=0.00065878; logf=-7.3251 ; lambda=10; chi=0.01; psi=0.01;mode=1800.001
f=0.0065878; logf=-5.0225 ; lambda=10; chi=0.01; psi=0.1;mode=180.0006
f=0.065878; logf=-2.72 ; lambda=10; chi=0.01; psi=1;mode=18.00056
f=0.65878; logf=-0.41737 ; lambda=10; chi=0.01; psi=10;mode=1.800555
f=6.5877; logf=1.8852 ; lambda=10; chi=0.01; psi=100;mode=0.1805539
f=0.00065878; logf=-7.3251 ; lambda=10; chi=0.1; psi=0.01;mode=1800.006
f=0.0065878; logf=-5.0225 ; lambda=10; chi=0.1; psi=0.1;mode=180.0056
f=0.065878; logf=-2.72 ; lambda=10; chi=0.1; psi=1;mode=18.00555
f=0.65877; logf=-0.41737 ; lambda=10; chi=0.1; psi=10;mode=1.805539
f=6.5847; logf=1.8847 ; lambda=10; chi=0.1; psi=100;mode=0.1853939
f=0.00065878; logf=-7.3251 ; lambda=10; chi=1; psi=0.01;mode=1800.056
f=0.0065878; logf=-5.0225 ; lambda=10; chi=1; psi=0.1;mode=180.0555
f=0.065877; logf=-2.72 ; lambda=10; chi=1; psi=1;mode=18.05539
f=0.65847; logf=-0.41784 ; lambda=10; chi=1; psi=10;mode=1.853939
f=6.4527; logf=1.8645 ; lambda=10; chi=1; psi=100;mode=0.2245362
f=0.00065878; logf=-7.3251 ; lambda=10; chi=10; psi=0.01;mode=1800.555
f=0.0065877; logf=-5.0225 ; lambda=10; chi=10; psi=0.1;mode=180.5539
f=0.065847; logf=-2.7204 ; lambda=10; chi=10; psi=1;mode=18.53939
f=0.64527; logf=-0.43809 ; lambda=10; chi=10; psi=10;mode=2.245362
f=5.4207; logf=1.6902 ; lambda=10; chi=10; psi=100;mode=0.4187856
f=0.00065877; logf=-7.3251 ; lambda=10; chi=100; psi=0.01;mode=1805.539
f=0.0065847; logf=-5.023 ; lambda=10; chi=100; psi=0.1;mode=185.3939
f=0.064527; logf=-2.7407 ; lambda=10; chi=100; psi=1;mode=22.45362
f=0.54207; logf=-0.61236 ; lambda=10; chi=100; psi=10;mode=4.187856
f=3.6419; logf=1.2925 ; lambda=10; chi=100; psi=100;mode=1.094042
f=0.00020031; logf=-8.5157 ; lambda=100; chi=0; psi=0.01;mode=19800
f=0.0020031; logf=-6.2131 ; lambda=100; chi=0; psi=0.1;mode=1980
f=0.020031; logf=-3.9105 ; lambda=100; chi=0; psi=1;mode=198
f=0.20031; logf=-1.6079 ; lambda=100; chi=0; psi=10;mode=19.8
f=2.0031; logf=0.69468 ; lambda=100; chi=0; psi=100;mode=1.98
f=0.00020031; logf=-8.5157 ; lambda=100; chi=0.01; psi=0.01;mode=19800
f=0.0020031; logf=-6.2131 ; lambda=100; chi=0.01; psi=0.1;mode=1980
f=0.020031; logf=-3.9105 ; lambda=100; chi=0.01; psi=1;mode=198.0001
f=0.20031; logf=-1.6079 ; lambda=100; chi=0.01; psi=10;mode=19.80005
f=2.0031; logf=0.69468 ; lambda=100; chi=0.01; psi=100;mode=1.980051
f=0.00020031; logf=-8.5157 ; lambda=100; chi=0.1; psi=0.01;mode=19800
f=0.0020031; logf=-6.2131 ; lambda=100; chi=0.1; psi=0.1;mode=1980.001
f=0.020031; logf=-3.9105 ; lambda=100; chi=0.1; psi=1;mode=198.0005
f=0.20031; logf=-1.6079 ; lambda=100; chi=0.1; psi=10;mode=19.80051
f=2.0031; logf=0.69468 ; lambda=100; chi=0.1; psi=100;mode=1.980505
f=0.00020031; logf=-8.5157 ; lambda=100; chi=1; psi=0.01;mode=19800.01
f=0.0020031; logf=-6.2131 ; lambda=100; chi=1; psi=0.1;mode=1980.005
f=0.020031; logf=-3.9105 ; lambda=100; chi=1; psi=1;mode=198.0051
f=0.20031; logf=-1.6079 ; lambda=100; chi=1; psi=10;mode=19.80505
f=2.0031; logf=0.69468 ; lambda=100; chi=1; psi=100;mode=1.985038
f=0.00020031; logf=-8.5157 ; lambda=100; chi=10; psi=0.01;mode=19800.05
f=0.0020031; logf=-6.2131 ; lambda=100; chi=10; psi=0.1;mode=1980.051
f=0.020031; logf=-3.9105 ; lambda=100; chi=10; psi=1;mode=198.0505
f=0.20031; logf=-1.6079 ; lambda=100; chi=10; psi=10;mode=19.85038
f=2.0025; logf=0.69439 ; lambda=100; chi=10; psi=100;mode=2.029279
f=0.00020031; logf=-8.5157 ; lambda=100; chi=100; psi=0.01;mode=19800.51
f=0.0020031; logf=-6.2131 ; lambda=100; chi=100; psi=0.1;mode=1980.505
f=0.020031; logf=-3.9105 ; lambda=100; chi=100; psi=1;mode=198.5038
f=0.20025; logf=-1.6082 ; lambda=100; chi=100; psi=10;mode=20.29279
f=1.9724; logf=0.67925 ; lambda=100; chi=100; psi=100;mode=2.39716
>
> ## --- verify results for special values -------------------------------------
>
> xvals <- c(-Inf,-1,0,Inf)
> lambda <- c(-100, -10, -1, -0.1, 0, 0.1, 1, 10, 100)
> chi <- c(0, 0.01, 0.1, 1, 10, 100)
> psi <- c(0, 0.01, 0.1, 1, 10, 100)
>
> for (x in xvals) { for (l in lambda) { for (c in chi) { for (p in psi) {
+ if (!check.params(l,c,p)) next
+ f <- dgig(x,l,c,p)
+ lf <- dgig(x,l,c,p,TRUE)
+ if (! isTRUE(all(c(f,lf) == c(0,-Inf))) ) {
+ cat ("f=",f,"; logf=",lf,"\t",
+ "; x=",x,"; lambda=",l,"; chi=",c,"; psi=",p,"\n", sep="")
+ stop ("0 and -Inf expected")
+ }
+ }}}}
>
> ## --- verify normalization constant: integrate density ----------------------
>
> lambda <- c(-10, -1, -0.1, 0, 0.1, 1, 10)
> chi <- c(0.01, 0.1, 1, 10, 100)
> psi <- c(0.01, 0.1, 1, 10, 100)
>
> for (l in lambda) { for (c in chi) { for (p in psi) {
+ if (!check.params(l,c,p)) next
+ if ( (c < 0.009 && p < 0.05) || (c < 0.05 && p < 0.009) ) {
+ ## integrate() does not work for these parameters
+ next
+ }
+ a <- integrate(dgig,0,Inf,lambda=l,chi=c,psi=p,rel.tol=1.e-10)
+ if (all.equal(a$value,1, tolerance=1e-4) != TRUE) {
+ cat("lambda =",l,"; chi =",c,"; psi =",p,"\n\t")
+ print(a)
+ stop("normalization constant might be wrong")
+ }
+ }}}
>
> ## special cases
> lambda <- c(-10, -1, -0.5, 0, 0.5,1,10)
> chi <- c(0, 0.1, 1, 10, 100)
> psi <- c(0, 0.1, 1, 10, 100)
>
> for (l in lambda) { for (c in chi) { for (p in psi) {
+ if (!check.params(l,c,p)) next
+ if ( (c < 0.009 && p < 0.05) || (c < 0.05 && p < 0.009) ) {
+ ## integrate does not work for these parameters
+ next
+ }
+ a <- integrate(dgig,0,Inf,lambda=l,chi=c,psi=p,rel.tol=1.e-10)
+ if (all.equal(a$value,1, tolerance=1e-4) != TRUE) {
+ cat("lambda =",l,"; chi =",c,"; psi =",p,"\n\t")
+ print(a)
+ stop("normalization constant might be wrong")
+ }
+ }}}
>
>
> ##############################################################################
> ##
> ## rgig()
> ##
> ##############################################################################
>
> ## --- check handling of invalid parameters ----------------------------------
>
> stopifnot (iserror(rgig(-1, 1,1,1)), TRUE)
Error in rgig(-1, 1, 1, 1) :
sample size 'n' must be non-negative integer.
>
> stopifnot (iserror(rgig(1,-1,-1, 1)), TRUE)
Error in rgig(1, -1, -1, 1) :
invalid parameters for GIG distribution: lambda=-1, chi=-1, psi=1
> stopifnot (iserror(rgig(1,-1, 0, 1)), TRUE)
Error in rgig(1, -1, 0, 1) :
invalid parameters for GIG distribution: lambda=-1, chi=0, psi=1
> stopifnot (iserror(rgig(1,-1, 1,-1)), TRUE)
Error in rgig(1, -1, 1, -1) :
invalid parameters for GIG distribution: lambda=-1, chi=1, psi=-1
> stopifnot (iserror(rgig(1, 0,-1, 1)), TRUE)
Error in rgig(1, 0, -1, 1) :
invalid parameters for GIG distribution: lambda=0, chi=-1, psi=1
> stopifnot (iserror(rgig(1, 0, 1, 0)), TRUE)
Error in rgig(1, 0, 1, 0) :
invalid parameters for GIG distribution: lambda=0, chi=1, psi=0
> stopifnot (iserror(rgig(1, 0, 0, 1)), TRUE)
Error in rgig(1, 0, 0, 1) :
invalid parameters for GIG distribution: lambda=0, chi=0, psi=1
> stopifnot (iserror(rgig(1, 0, 1,-1)), TRUE)
Error in rgig(1, 0, 1, -1) :
invalid parameters for GIG distribution: lambda=0, chi=1, psi=-1
> stopifnot (iserror(rgig(1, 1,-1, 1)), TRUE)
Error in rgig(1, 1, -1, 1) :
invalid parameters for GIG distribution: lambda=1, chi=-1, psi=1
> stopifnot (iserror(rgig(1, 1, 1, 0)), TRUE)
Error in rgig(1, 1, 1, 0) :
invalid parameters for GIG distribution: lambda=1, chi=1, psi=0
> stopifnot (iserror(rgig(1, 1, 1,-1)), TRUE)
Error in rgig(1, 1, 1, -1) :
invalid parameters for GIG distribution: lambda=1, chi=1, psi=-1
>
> ## --- just run generator ----------------------------------------------------
>
> if(!interactive())
+ set.seed(123)
>
> lambda <- c(-100, -10, -1, -0.1, 0, 0.1, 1, 10, 100)
> chi <- c(0, 1e-12, 0.01, 0.1, 1, 10, 100)
> psi <- c(0, 1e-12, 0.01, 0.1, 1, 10, 100)
>
> for (l in lambda) { for (c in chi) { for (p in psi) {
+ if (!check.params(l,c,p)) next
+ x <- rgig(1,l,c,p)
+ cat ("x=",signif(x,5),";\tlambda=",l,"; chi=",c,"; psi=",p,"\n", sep="")
+ if (! isTRUE(is.finite(x))) {
+ cat ("x=",x,";\tlambda=",l,"; chi=",c,"; psi=",p,"\n", sep="")
+ stop ("result not finite")
+ }
+ }}}
x=5.3198e-15; lambda=-100; chi=1e-12; psi=0
x=5.8782e-15; lambda=-100; chi=1e-12; psi=1e-12
x=4.4685e-15; lambda=-100; chi=1e-12; psi=0.01
x=5.0685e-15; lambda=-100; chi=1e-12; psi=0.1
x=5.0803e-15; lambda=-100; chi=1e-12; psi=1
x=4.3483e-15; lambda=-100; chi=1e-12; psi=10
x=4.7721e-15; lambda=-100; chi=1e-12; psi=100
x=3.9583e-05; lambda=-100; chi=0.01; psi=0
x=4.704e-05; lambda=-100; chi=0.01; psi=1e-12
x=4.7187e-05; lambda=-100; chi=0.01; psi=0.01
x=4.5646e-05; lambda=-100; chi=0.01; psi=0.1
x=5.0488e-05; lambda=-100; chi=0.01; psi=1
x=5.8796e-05; lambda=-100; chi=0.01; psi=10
x=5.184e-05; lambda=-100; chi=0.01; psi=100
x=0.00051982; lambda=-100; chi=0.1; psi=0
x=0.00055508; lambda=-100; chi=0.1; psi=1e-12
x=0.00052737; lambda=-100; chi=0.1; psi=0.01
x=0.00060707; lambda=-100; chi=0.1; psi=0.1
x=0.0004717; lambda=-100; chi=0.1; psi=1
x=0.00054891; lambda=-100; chi=0.1; psi=10
x=0.00042325; lambda=-100; chi=0.1; psi=100
x=0.004817; lambda=-100; chi=1; psi=0
x=0.0053621; lambda=-100; chi=1; psi=1e-12
x=0.0044687; lambda=-100; chi=1; psi=0.01
x=0.0047161; lambda=-100; chi=1; psi=0.1
x=0.00449; lambda=-100; chi=1; psi=1
x=0.0046914; lambda=-100; chi=1; psi=10
x=0.0050417; lambda=-100; chi=1; psi=100
x=0.051829; lambda=-100; chi=10; psi=0
x=0.046797; lambda=-100; chi=10; psi=1e-12
x=0.050874; lambda=-100; chi=10; psi=0.01
x=0.046384; lambda=-100; chi=10; psi=0.1
x=0.04602; lambda=-100; chi=10; psi=1
x=0.053794; lambda=-100; chi=10; psi=10
x=0.057177; lambda=-100; chi=10; psi=100
x=0.50106; lambda=-100; chi=100; psi=0
x=0.53202; lambda=-100; chi=100; psi=1e-12
x=0.43132; lambda=-100; chi=100; psi=0.01
x=0.47831; lambda=-100; chi=100; psi=0.1
x=0.45277; lambda=-100; chi=100; psi=1
x=0.54904; lambda=-100; chi=100; psi=10
x=0.37316; lambda=-100; chi=100; psi=100
x=5.7015e-14; lambda=-10; chi=1e-12; psi=0
x=7.9626e-14; lambda=-10; chi=1e-12; psi=1e-12
x=5.7343e-14; lambda=-10; chi=1e-12; psi=0.01
x=7.2658e-14; lambda=-10; chi=1e-12; psi=0.1
x=4.9261e-14; lambda=-10; chi=1e-12; psi=1
x=4.6423e-14; lambda=-10; chi=1e-12; psi=10
x=5.1131e-14; lambda=-10; chi=1e-12; psi=100
x=0.00039824; lambda=-10; chi=0.01; psi=0
x=0.00065872; lambda=-10; chi=0.01; psi=1e-12
x=0.00064581; lambda=-10; chi=0.01; psi=0.01
x=0.00037175; lambda=-10; chi=0.01; psi=0.1
x=0.00060465; lambda=-10; chi=0.01; psi=1
x=0.0010324; lambda=-10; chi=0.01; psi=10
x=0.00039628; lambda=-10; chi=0.01; psi=100
x=0.0055086; lambda=-10; chi=0.1; psi=0
x=0.005917; lambda=-10; chi=0.1; psi=1e-12
x=0.0060243; lambda=-10; chi=0.1; psi=0.01
x=0.005111; lambda=-10; chi=0.1; psi=0.1
x=0.0077388; lambda=-10; chi=0.1; psi=1
x=0.0064675; lambda=-10; chi=0.1; psi=10
x=0.0063866; lambda=-10; chi=0.1; psi=100
x=0.042897; lambda=-10; chi=1; psi=0
x=0.042602; lambda=-10; chi=1; psi=1e-12
x=0.051869; lambda=-10; chi=1; psi=0.01
x=0.050046; lambda=-10; chi=1; psi=0.1
x=0.06084; lambda=-10; chi=1; psi=1
x=0.063388; lambda=-10; chi=1; psi=10
x=0.040537; lambda=-10; chi=1; psi=100
x=0.49961; lambda=-10; chi=10; psi=0
x=0.59059; lambda=-10; chi=10; psi=1e-12
x=0.49503; lambda=-10; chi=10; psi=0.01
x=0.52129; lambda=-10; chi=10; psi=0.1
x=0.28793; lambda=-10; chi=10; psi=1
x=0.34495; lambda=-10; chi=10; psi=10
x=0.18954; lambda=-10; chi=10; psi=100
x=5.1994; lambda=-10; chi=100; psi=0
x=3.2155; lambda=-10; chi=100; psi=1e-12
x=10.144; lambda=-10; chi=100; psi=0.01
x=7.6194; lambda=-10; chi=100; psi=0.1
x=7.1706; lambda=-10; chi=100; psi=1
x=1.9574; lambda=-10; chi=100; psi=10
x=0.92506; lambda=-10; chi=100; psi=100
x=2.5299e-13; lambda=-1; chi=1e-12; psi=0
x=8.2164e-13; lambda=-1; chi=1e-12; psi=1e-12
x=8.4697e-13; lambda=-1; chi=1e-12; psi=0.01
x=8.5986e-13; lambda=-1; chi=1e-12; psi=0.1
x=3.0314e-12; lambda=-1; chi=1e-12; psi=1
x=3.3836e-13; lambda=-1; chi=1e-12; psi=10
x=4.6253e-13; lambda=-1; chi=1e-12; psi=100
x=0.0059525; lambda=-1; chi=0.01; psi=0
x=0.0033479; lambda=-1; chi=0.01; psi=1e-12
x=0.0037231; lambda=-1; chi=0.01; psi=0.01
x=0.0023678; lambda=-1; chi=0.01; psi=0.1
x=0.017296; lambda=-1; chi=0.01; psi=1
x=0.01159; lambda=-1; chi=0.01; psi=10
x=0.0051073; lambda=-1; chi=0.01; psi=100
x=0.031204; lambda=-1; chi=0.1; psi=0
x=0.039645; lambda=-1; chi=0.1; psi=1e-12
x=0.76503; lambda=-1; chi=0.1; psi=0.01
x=0.30344; lambda=-1; chi=0.1; psi=0.1
x=0.0093096; lambda=-1; chi=0.1; psi=1
x=0.051207; lambda=-1; chi=0.1; psi=10
x=0.060959; lambda=-1; chi=0.1; psi=100
x=4.6088; lambda=-1; chi=1; psi=0
x=2.3642; lambda=-1; chi=1; psi=1e-12
x=7.9364; lambda=-1; chi=1; psi=0.01
x=0.19749; lambda=-1; chi=1; psi=0.1
x=0.47418; lambda=-1; chi=1; psi=1
x=0.44547; lambda=-1; chi=1; psi=10
x=0.07695; lambda=-1; chi=1; psi=100
x=6.5574; lambda=-1; chi=10; psi=0
x=7.394; lambda=-1; chi=10; psi=1e-12
x=6.9632; lambda=-1; chi=10; psi=0.01
x=9.9105; lambda=-1; chi=10; psi=0.1
x=0.94265; lambda=-1; chi=10; psi=1
x=1.0247; lambda=-1; chi=10; psi=10
x=0.2205; lambda=-1; chi=10; psi=100
x=58.55; lambda=-1; chi=100; psi=0
x=96.213; lambda=-1; chi=100; psi=1e-12
x=43.965; lambda=-1; chi=100; psi=0.01
x=11.634; lambda=-1; chi=100; psi=0.1
x=14.088; lambda=-1; chi=100; psi=1
x=3.6552; lambda=-1; chi=100; psi=10
x=1.0168; lambda=-1; chi=100; psi=100
x=4.4687e-09; lambda=-0.1; chi=1e-12; psi=0
x=4.5862e-13; lambda=-0.1; chi=1e-12; psi=1e-12
x=5.3712e-06; lambda=-0.1; chi=1e-12; psi=0.01
x=1.1863e-09; lambda=-0.1; chi=1e-12; psi=0.1
x=1.8997e-12; lambda=-0.1; chi=1e-12; psi=1
x=4.1717e-07; lambda=-0.1; chi=1e-12; psi=10
x=8.018e-09; lambda=-0.1; chi=1e-12; psi=100
x=0.039738; lambda=-0.1; chi=0.01; psi=0
x=3.1473e+10; lambda=-0.1; chi=0.01; psi=1e-12
x=0.23561; lambda=-0.1; chi=0.01; psi=0.01
x=0.4198; lambda=-0.1; chi=0.01; psi=0.1
x=0.067642; lambda=-0.1; chi=0.01; psi=1
x=0.038787; lambda=-0.1; chi=0.01; psi=10
x=0.0044857; lambda=-0.1; chi=0.01; psi=100
x=0.085631; lambda=-0.1; chi=0.1; psi=0
x=3928.4; lambda=-0.1; chi=0.1; psi=1e-12
x=6.5582; lambda=-0.1; chi=0.1; psi=0.01
x=0.09847; lambda=-0.1; chi=0.1; psi=0.1
x=0.29626; lambda=-0.1; chi=0.1; psi=1
x=0.10684; lambda=-0.1; chi=0.1; psi=10
x=0.052127; lambda=-0.1; chi=0.1; psi=100
x=86.213; lambda=-0.1; chi=1; psi=0
x=548.16; lambda=-0.1; chi=1; psi=1e-12
x=109.53; lambda=-0.1; chi=1; psi=0.01
x=1.2215; lambda=-0.1; chi=1; psi=0.1
x=0.23694; lambda=-0.1; chi=1; psi=1
x=0.2352; lambda=-0.1; chi=1; psi=10
x=0.12375; lambda=-0.1; chi=1; psi=100
x=293230; lambda=-0.1; chi=10; psi=0
x=196.95; lambda=-0.1; chi=10; psi=1e-12
x=12.69; lambda=-0.1; chi=10; psi=0.01
x=3.7853; lambda=-0.1; chi=10; psi=0.1
x=2.2381; lambda=-0.1; chi=10; psi=1
x=0.50362; lambda=-0.1; chi=10; psi=10
x=0.28161; lambda=-0.1; chi=10; psi=100
x=2946700000; lambda=-0.1; chi=100; psi=0
x=4229; lambda=-0.1; chi=100; psi=1e-12
x=98.266; lambda=-0.1; chi=100; psi=0.01
x=78.511; lambda=-0.1; chi=100; psi=0.1
x=8.3683; lambda=-0.1; chi=100; psi=1
x=3.6636; lambda=-0.1; chi=100; psi=10
x=0.92552; lambda=-0.1; chi=100; psi=100
x=1600200; lambda=0; chi=1e-12; psi=1e-12
x=3.3434e-11; lambda=0; chi=1e-12; psi=0.01
x=4.9032e-07; lambda=0; chi=1e-12; psi=0.1
x=6.2091e-11; lambda=0; chi=1e-12; psi=1
x=7.4185e-10; lambda=0; chi=1e-12; psi=10
x=1.3549e-05; lambda=0; chi=1e-12; psi=100
x=4554300000; lambda=0; chi=0.01; psi=1e-12
x=40.812; lambda=0; chi=0.01; psi=0.01
x=3.4623; lambda=0; chi=0.01; psi=0.1
x=0.08048; lambda=0; chi=0.01; psi=1
x=0.10359; lambda=0; chi=0.01; psi=10
x=0.026447; lambda=0; chi=0.01; psi=100
x=143810000; lambda=0; chi=0.1; psi=1e-12
x=0.12875; lambda=0; chi=0.1; psi=0.01
x=3.6899; lambda=0; chi=0.1; psi=0.1
x=0.04848; lambda=0; chi=0.1; psi=1
x=0.22254; lambda=0; chi=0.1; psi=10
x=0.029077; lambda=0; chi=0.1; psi=100
x=5.0413; lambda=0; chi=1; psi=1e-12
x=1.4174; lambda=0; chi=1; psi=0.01
x=2.3226; lambda=0; chi=1; psi=0.1
x=0.53766; lambda=0; chi=1; psi=1
x=0.16095; lambda=0; chi=1; psi=10
x=0.16294; lambda=0; chi=1; psi=100
x=94713; lambda=0; chi=10; psi=1e-12
x=2.2778; lambda=0; chi=10; psi=0.01
x=3.2252; lambda=0; chi=10; psi=0.1
x=4.509; lambda=0; chi=10; psi=1
x=1.4193; lambda=0; chi=10; psi=10
x=0.37568; lambda=0; chi=10; psi=100
x=43667; lambda=0; chi=100; psi=1e-12
x=36.762; lambda=0; chi=100; psi=0.01
x=29.789; lambda=0; chi=100; psi=0.1
x=7.4166; lambda=0; chi=100; psi=1
x=3.8694; lambda=0; chi=100; psi=10
x=1.0214; lambda=0; chi=100; psi=100
x=5009600000; lambda=0.1; chi=0; psi=1e-12
x=3.4872e-16; lambda=0.1; chi=0; psi=0.01
x=3.1019e-06; lambda=0.1; chi=0; psi=0.1
x=0.040941; lambda=0.1; chi=0; psi=1
x=5.3222e-05; lambda=0.1; chi=0; psi=10
x=1.2521e-19; lambda=0.1; chi=0; psi=100
x=2374400000; lambda=0.1; chi=1e-12; psi=1e-12
x=0.0228; lambda=0.1; chi=1e-12; psi=0.01
x=1.5078e-06; lambda=0.1; chi=1e-12; psi=0.1
x=0.018816; lambda=0.1; chi=1e-12; psi=1
x=2.3662e-10; lambda=0.1; chi=1e-12; psi=10
x=2.1874e-09; lambda=0.1; chi=1e-12; psi=100
x=8.4363e+11; lambda=0.1; chi=0.01; psi=1e-12
x=0.85673; lambda=0.1; chi=0.01; psi=0.01
x=0.0057843; lambda=0.1; chi=0.01; psi=0.1
x=0.59164; lambda=0.1; chi=0.01; psi=1
x=0.069772; lambda=0.1; chi=0.01; psi=10
x=0.027144; lambda=0.1; chi=0.01; psi=100
x=304990000; lambda=0.1; chi=0.1; psi=1e-12
x=0.17238; lambda=0.1; chi=0.1; psi=0.01
x=0.048855; lambda=0.1; chi=0.1; psi=0.1
x=3.4108; lambda=0.1; chi=0.1; psi=1
x=0.14537; lambda=0.1; chi=0.1; psi=10
x=0.063399; lambda=0.1; chi=0.1; psi=100
x=61669000; lambda=0.1; chi=1; psi=1e-12
x=138.29; lambda=0.1; chi=1; psi=0.01
x=26.72; lambda=0.1; chi=1; psi=0.1
x=1.0028; lambda=0.1; chi=1; psi=1
x=0.29028; lambda=0.1; chi=1; psi=10
x=0.090021; lambda=0.1; chi=1; psi=100
x=3.9078e+10; lambda=0.1; chi=10; psi=1e-12
x=10.828; lambda=0.1; chi=10; psi=0.01
x=8.9138; lambda=0.1; chi=10; psi=0.1
x=5.7627; lambda=0.1; chi=10; psi=1
x=1.0396; lambda=0.1; chi=10; psi=10
x=0.23554; lambda=0.1; chi=10; psi=100
x=341870000; lambda=0.1; chi=100; psi=1e-12
x=330.82; lambda=0.1; chi=100; psi=0.01
x=18.119; lambda=0.1; chi=100; psi=0.1
x=17.122; lambda=0.1; chi=100; psi=1
x=2.726; lambda=0.1; chi=100; psi=10
x=1.017; lambda=0.1; chi=100; psi=100
x=1.6802e+11; lambda=1; chi=0; psi=1e-12
x=207.75; lambda=1; chi=0; psi=0.01
x=28.161; lambda=1; chi=0; psi=0.1
x=4.7594; lambda=1; chi=0; psi=1
x=0.10808; lambda=1; chi=0; psi=10
x=0.0092784; lambda=1; chi=0; psi=100
x=1.7313e+12; lambda=1; chi=1e-12; psi=1e-12
x=72.667; lambda=1; chi=1e-12; psi=0.01
x=51.988; lambda=1; chi=1e-12; psi=0.1
x=0.33062; lambda=1; chi=1e-12; psi=1
x=0.061959; lambda=1; chi=1e-12; psi=10
x=0.037929; lambda=1; chi=1e-12; psi=100
x=1.7014e+12; lambda=1; chi=0.01; psi=1e-12
x=299.52; lambda=1; chi=0.01; psi=0.01
x=10.166; lambda=1; chi=0.01; psi=0.1
x=1.2207; lambda=1; chi=0.01; psi=1
x=0.25663; lambda=1; chi=0.01; psi=10
x=0.0058837; lambda=1; chi=0.01; psi=100
x=8.0667e+10; lambda=1; chi=0.1; psi=1e-12
x=58.7; lambda=1; chi=0.1; psi=0.01
x=46.461; lambda=1; chi=0.1; psi=0.1
x=3.1843; lambda=1; chi=0.1; psi=1
x=0.15154; lambda=1; chi=0.1; psi=10
x=0.019606; lambda=1; chi=0.1; psi=100
x=5.0255e+12; lambda=1; chi=1; psi=1e-12
x=241.84; lambda=1; chi=1; psi=0.01
x=2.201; lambda=1; chi=1; psi=0.1
x=2.2627; lambda=1; chi=1; psi=1
x=0.24283; lambda=1; chi=1; psi=10
x=0.084907; lambda=1; chi=1; psi=100
x=1.2464e+12; lambda=1; chi=10; psi=1e-12
x=169.46; lambda=1; chi=10; psi=0.01
x=13.782; lambda=1; chi=10; psi=0.1
x=5.6679; lambda=1; chi=10; psi=1
x=0.76255; lambda=1; chi=10; psi=10
x=0.31157; lambda=1; chi=10; psi=100
x=4.5977e+12; lambda=1; chi=100; psi=1e-12
x=184.07; lambda=1; chi=100; psi=0.01
x=56.119; lambda=1; chi=100; psi=0.1
x=13.893; lambda=1; chi=100; psi=1
x=3.2254; lambda=1; chi=100; psi=10
x=0.9052; lambda=1; chi=100; psi=100
x=1.5586e+13; lambda=10; chi=0; psi=1e-12
x=1854.9; lambda=10; chi=0; psi=0.01
x=100.79; lambda=10; chi=0; psi=0.1
x=18.823; lambda=10; chi=0; psi=1
x=2.1766; lambda=10; chi=0; psi=10
x=0.22276; lambda=10; chi=0; psi=100
x=2.4204e+13; lambda=10; chi=1e-12; psi=1e-12
x=799.49; lambda=10; chi=1e-12; psi=0.01
x=127.65; lambda=10; chi=1e-12; psi=0.1
x=20.106; lambda=10; chi=1e-12; psi=1
x=1.6086; lambda=10; chi=1e-12; psi=10
x=0.09995; lambda=10; chi=1e-12; psi=100
x=2.5864e+13; lambda=10; chi=0.01; psi=1e-12
x=2266.4; lambda=10; chi=0.01; psi=0.01
x=202.14; lambda=10; chi=0.01; psi=0.1
x=26.453; lambda=10; chi=0.01; psi=1
x=2.3342; lambda=10; chi=0.01; psi=10
x=0.12326; lambda=10; chi=0.01; psi=100
x=1.6147e+13; lambda=10; chi=0.1; psi=1e-12
x=866; lambda=10; chi=0.1; psi=0.01
x=240.45; lambda=10; chi=0.1; psi=0.1
x=18.154; lambda=10; chi=0.1; psi=1
x=2.9507; lambda=10; chi=0.1; psi=10
x=0.12418; lambda=10; chi=0.1; psi=100
x=1.8778e+13; lambda=10; chi=1; psi=1e-12
x=795.37; lambda=10; chi=1; psi=0.01
x=116.8; lambda=10; chi=1; psi=0.1
x=27.355; lambda=10; chi=1; psi=1
x=1.5081; lambda=10; chi=1; psi=10
x=0.23705; lambda=10; chi=1; psi=100
x=2.1163e+13; lambda=10; chi=10; psi=1e-12
x=3833.6; lambda=10; chi=10; psi=0.01
x=257.47; lambda=10; chi=10; psi=0.1
x=14.708; lambda=10; chi=10; psi=1
x=3.5502; lambda=10; chi=10; psi=10
x=0.34739; lambda=10; chi=10; psi=100
x=1.4323e+13; lambda=10; chi=100; psi=1e-12
x=2325.5; lambda=10; chi=100; psi=0.01
x=201.31; lambda=10; chi=100; psi=0.1
x=13.02; lambda=10; chi=100; psi=1
x=4.3335; lambda=10; chi=100; psi=10
x=1.1271; lambda=10; chi=100; psi=100
x=1.9957e+14; lambda=100; chi=0; psi=1e-12
x=22580; lambda=100; chi=0; psi=0.01
x=1708.4; lambda=100; chi=0; psi=0.1
x=221.77; lambda=100; chi=0; psi=1
x=18.944; lambda=100; chi=0; psi=10
x=2.0058; lambda=100; chi=0; psi=100
x=1.9137e+14; lambda=100; chi=1e-12; psi=1e-12
x=21662; lambda=100; chi=1e-12; psi=0.01
x=2184.4; lambda=100; chi=1e-12; psi=0.1
x=217.51; lambda=100; chi=1e-12; psi=1
x=20.989; lambda=100; chi=1e-12; psi=10
x=1.8593; lambda=100; chi=1e-12; psi=100
x=2.0228e+14; lambda=100; chi=0.01; psi=1e-12
x=17736; lambda=100; chi=0.01; psi=0.01
x=1992.2; lambda=100; chi=0.01; psi=0.1
x=176.38; lambda=100; chi=0.01; psi=1
x=17.149; lambda=100; chi=0.01; psi=10
x=1.8693; lambda=100; chi=0.01; psi=100
x=1.8749e+14; lambda=100; chi=0.1; psi=1e-12
x=22375; lambda=100; chi=0.1; psi=0.01
x=1634; lambda=100; chi=0.1; psi=0.1
x=194.61; lambda=100; chi=0.1; psi=1
x=19.654; lambda=100; chi=0.1; psi=10
x=2.0385; lambda=100; chi=0.1; psi=100
x=1.8188e+14; lambda=100; chi=1; psi=1e-12
x=19809; lambda=100; chi=1; psi=0.01
x=1730; lambda=100; chi=1; psi=0.1
x=204.02; lambda=100; chi=1; psi=1
x=19.377; lambda=100; chi=1; psi=10
x=2.0725; lambda=100; chi=1; psi=100
x=2.0371e+14; lambda=100; chi=10; psi=1e-12
x=18109; lambda=100; chi=10; psi=0.01
x=1829.2; lambda=100; chi=10; psi=0.1
x=193.45; lambda=100; chi=10; psi=1
x=21.731; lambda=100; chi=10; psi=10
x=1.9884; lambda=100; chi=10; psi=100
x=2.2308e+14; lambda=100; chi=100; psi=1e-12
x=21321; lambda=100; chi=100; psi=0.01
x=1981.2; lambda=100; chi=100; psi=0.1
x=195.86; lambda=100; chi=100; psi=1
x=17.792; lambda=100; chi=100; psi=10
x=2.9358; lambda=100; chi=100; psi=100
>
>
> ## --- end -------------------------------------------------------------------
>
>
> #################################################################################
>
> ## compareHistPdf <- function(n=10^5,lambda=0.5, chi=0.1, psi=2){
> ## # function to compare the histogram of a sample with the theoretic pdf
> ## # demonstrates the correctness of the generator
> ## y <- rgig(n=n,lambda=lambda,chi=chi,psi=psi)
> ## hist(y,breaks=100,freq=FALSE,
> ## main=paste("GIG ( lambda = ",lambda,", chi = ",chi,", psi = ",psi,")"))
> ## xval <- seq(0,max(y),length.out=5000) # to plot the corresponding density
> ## lines(xval,dgig(xval,lambda=lambda,chi=chi,psi=psi))
> ## }
> ## compareHistPdf(n=10^5,1,10,0.1)
> ## compareHistPdf(n=10^5,-5,10,0.5)
> ## compareHistPdf(n=10^5,-5,1,0.2)
> ## compareHistPdf(n=10^5,-0.01,0.1,0.2)
> ##
> ## #########################
> ## # Chi square test for 4 different parameter settings
> ## #
> ## # is slow as the pgig function is slow
> ## # requires the ghyp package
> ## # library(ghyp)
> ## # yu<- pgig(rgig(10000,0.5,0.3,4),0.5,0.3,4)
> ## # res<-hist(yu)
> ## # chisq.test(res$counts)
> ## # yu<- pgig(rgig(10000,2.5,0.3,4),2.5,0.3,4)
> ## # res<-hist(yu)
> ## # chisq.test(res$counts)
> ## # yu<- pgig(rgig(10000,2.5,1.e-17,4),2.5,1.e-17,4)
> ## # res<-hist(yu)
> ## # chisq.test(res$counts)
> ## # yu<- pgig(rgig(10000,-2.5,13,4),-2.5,13,4)
> ## # res<-hist(yu)
> ## # chisq.test(res$counts)
>
>
>
> proc.time()
user system elapsed
0.43 0.12 0.54
Flavor: r-devel-windows-x86_64