==2589555== Memcheck, a memory error detector
==2589555== Copyright (C) 2002-2024, and GNU GPL'd, by Julian Seward et al.
==2589555== Using Valgrind-3.24.0 and LibVEX; rerun with -h for copyright info
==2589555== Command: /data/blackswan/ripley/R/R-devel-vg/bin/exec/R --vanilla
==2589555== 

R Under development (unstable) (2026-03-12 r89611) -- "Unsuffered Consequences"
Copyright (C) 2026 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu

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.

  Natural language support but running in an English locale

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.

> pkgname <- "gctsc"
> source(file.path(R.home("share"), "R", "examples-header.R"))
> options(warn = 1)
> library('gctsc')
> 
> base::assign(".oldSearch", base::search(), pos = 'CheckExEnv')
> base::assign(".old_wd", base::getwd(), pos = 'CheckExEnv')
> cleanEx()
> nameEx("gctsc")
> ### * gctsc
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: gctsc
> ### Title: Fit a Copula-Based Count Time Series Model
> ### Aliases: gctsc
> 
> ### ** Examples
> 
> ## Example 1: Gaussian copula, Poisson marginal, AR(1)
> set.seed(42)
> n <- 500
> sim_dat <- sim_poisson(mu = 10, tau = 0.3, arma_order = c(1, 0),
+                        nsim = n, family = "gaussian")
> 
> dat <- data.frame(y = sim_dat$y)
> 
> fit_gauss <- gctsc(
+   y ~ 1,
+   data = dat,
+   marginal = poisson.marg(lambda.lower = 0),
+   cormat = arma.cormat(p = 1, q = 0), family = "gaussian",
+   method = "CE",
+   options = gctsc.opts(M = 1000, seed = 42)
+ )
> summary(fit_gauss)

--- Summary of Copula Time Series Model ---
Call:
 gctsc(formula = y ~ 1, data = dat, marginal = poisson.marg(lambda.lower = 0),      cormat = arma.cormat(p = 1, q = 0), method = "CE", family = "gaussian",      options = gctsc.opts(M = 1000, seed = 42)) 

Marginal Model Coefficients:
               Estimate Std. Error z value Pr(>|z|)    
mu_(Intercept)    9.777      0.185   52.85   <2e-16 ***

Copula (Dependence) Coefficients:
    Estimate Std. Error z value Pr(>|z|)    
ar1  0.28110    0.04364   6.442 1.18e-10 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Model Fit Statistics:
  Log-likelihood: -1235.8374 
  AIC:             2475.6748 
  BIC:             2484.1040 
> 
> ## Example 2: Student--t copula
> sim_dat_t <- sim_poisson(mu = 10, tau = 0.3, arma_order = c(1, 0),
+                          nsim = 500, family = "t", df = 10)
> 
> dat_t <- data.frame(y = sim_dat_t$y)
> 
> fit_t <- gctsc(
+   y ~ 1,
+   data = dat_t,
+   marginal = poisson.marg(lambda.lower = 0),
+   cormat = arma.cormat(p = 1, q = 0), family ="t",
+   df= 10, method = "CE",
+   options = gctsc.opts(M = 1000, seed = 42)
+ )
> summary(fit_t)

--- Summary of Copula Time Series Model ---
Call:
 gctsc(formula = y ~ 1, data = dat_t, marginal = poisson.marg(lambda.lower = 0),      cormat = arma.cormat(p = 1, q = 0), method = "CE", family = "t",      df = 10, options = gctsc.opts(M = 1000, seed = 42)) 

Marginal Model Coefficients:
               Estimate Std. Error z value Pr(>|z|)    
mu_(Intercept)   9.8087     0.1524   64.36   <2e-16 ***

Copula (Dependence) Coefficients:
    Estimate Std. Error z value Pr(>|z|)    
ar1  0.28304    0.04712   6.006  1.9e-09 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Model Fit Statistics:
  Log-likelihood: -1140.7499 
  AIC:             2285.4999 
  BIC:             2293.9291 
> 
> 
> 
> 
> cleanEx()
> nameEx("marginal.gctsc")
> ### * marginal.gctsc
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: marginal.gctsc
> ### Title: Marginal Models for Copula Time Series
> ### Aliases: marginal.gctsc poisson.marg negbin.marg binom.marg zip.marg
> ###   zib.marg bbinom.marg zibb.marg
> 
> ### ** Examples
> 
> poisson.marg(link = "identity")
$start
function (y, x) 
{
    fit <- glm.fit(x = x, y = y, family = poisson(link))
    lambda <- fit$coefficients
    if (anyNA(lambda)) 
        stop("NA detected in coefficient estimates. Check design matrix.")
    names(lambda) <- prefixed_names(x, "mu_")
    lambda <- attach_bounds(lambda, lambda.lower, lambda.upper, 
        "poisson.marg()")
    lambda
}
<bytecode: 0x1226f0c8>
<environment: 0x176d4670>

$npar
function (x) 
NCOL(x)
<bytecode: 0x122726a8>
<environment: 0x176d4670>

$bounds
function (y, x, lambda, family = "gaussian", df = NULL) 
{
    mu <- invlink(x %*% lambda)
    if (any(mu < 0)) 
        stop("Negative mean detected. Use 'log' link or check predictors.")
    pdf <- dpois(y, mu)
    cdf <- ppois(y, mu)
    bounds <- safe_cdf_bounds(pdf, cdf, family, df)
    cbind(bounds$lower, bounds$upper)
}
<bytecode: 0x12272398>
<environment: 0x176d4670>

attr(,"class")
[1] "marginal.gctsc"
> zibb.marg(link = "logit", size = 24)
$size
[1] 24

$start
function (y, x) 
{
    X_mu <- x$mu
    X_pi0 <- x$pi0
    has_intercept <- has_intercept(X_mu)
    y_nonzero <- y[y > 0]
    if (has_only_intercept(X_mu)) {
        df_nonzero <- data.frame(y = y_nonzero, trials = size)
        fit <- VGAM::vglm(cbind(y, trials - y) ~ 1, family = VGAM::betabinomial, 
            data = df_nonzero)
        cf <- VGAM::Coef(fit)
        beta <- qlogis(cf[1])
        rho_logit <- qlogis(cf[2])
        coef_names <- c("mu_Intercept", "logit_rho")
    }
    else {
        if (has_intercept) {
            intercept_idx <- which(is_intercept_col(X_mu))
            X_mu_sub <- X_mu[, -intercept_idx, drop = FALSE]
            coef_names <- c("mu_Intercept", prefixed_names(X_mu_sub, 
                "mu_"), "logit_rho")
        }
        else {
            X_mu_sub <- X_mu
            coef_names <- c(prefixed_names(X_mu_sub, "mu_"), 
                "logit_rho")
        }
        df_nonzero <- data.frame(y = y_nonzero, trials = size, 
            X_mu[y > 0, , drop = FALSE])
        formula <- as.formula(paste("cbind(y, trials - y) ~", 
            paste(colnames(X_mu_sub), collapse = " + ")))
        fit <- VGAM::vglm(formula, family = VGAM::betabinomial, 
            data = df_nonzero)
        cf <- VGAM::Coef(fit)
        n_coef <- length(cf)
        beta <- c(cf[1], cf[3:(n_coef)])
        rho_logit <- (cf[2])
    }
    fit_pi0 <- glm(I(y == 0) ~ X_pi0 - 1, family = binomial(link = "logit"))
    alpha <- coef(fit_pi0)
    lambda <- c(beta, rho_logit, alpha)
    names(lambda) <- c(coef_names, prefixed_names(X_pi0, "pi0_"))
    lambda <- attach_bounds(lambda, lambda.lower, lambda.upper, 
        "zibb.marg()")
    lambda
}
<bytecode: 0x176e3ba8>
<environment: 0x176e9098>

$npar
function (x) 
{
    if (!is.list(x) || is.null(x$mu) || is.null(x$pi0)) {
        stop("x must be a list with elements 'mu' and 'pi0'")
    }
    ncol(x$mu) + 1 + ncol(x$pi0)
}
<bytecode: 0x1091e378>
<environment: 0x176e9098>

$bounds
function (y, x, lambda, family = "gaussian", df = NULL) 
{
    if (!is.list(x)) 
        stop("x must be a list with 'mu' and 'pi0'")
    X_mu <- x$mu
    X_pi0 <- x$pi0
    p_mu <- ncol(X_mu)
    np <- ncol(X_pi0)
    beta <- lambda[1:p_mu]
    rho_logit <- lambda[p_mu + 1]
    rho <- plogis(rho_logit)
    alpha <- lambda[(p_mu + 2):(p_mu + 1 + np)]
    mu <- invlink(X_mu %*% beta)
    pi0 <- plogis(X_pi0 %*% alpha)
    alpha1 <- mu * (1 - rho)/rho
    alpha2 <- (1 - mu) * (1 - rho)/rho
    f0 <- VGAM::dbetabinom.ab(0, size = size, shape1 = alpha1, 
        shape2 = alpha2)
    pmf <- VGAM::dbetabinom.ab(y, size = size, shape1 = alpha1, 
        shape2 = alpha2)
    cdf <- VGAM::pbetabinom.ab(y, size = size, shape1 = alpha1, 
        shape2 = alpha2)
    pdf <- ifelse(y == 0, pi0 + (1 - pi0) * f0, (1 - pi0) * pmf)
    cdf <- pi0 + (1 - pi0) * cdf
    bounds <- safe_cdf_bounds(pdf, cdf, family, df)
    cbind(bounds$lower, bounds$upper)
}
<bytecode: 0x1091d658>
<environment: 0x176e9098>

attr(,"class")
[1] "marginal.gctsc"
> 
> 
> 
> 
> cleanEx()
> nameEx("plot.gctsc")
> ### * plot.gctsc
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: plot.gctsc
> ### Title: Diagnostic Plots for Fitted Copula Count Time Series Models
> ### Aliases: plot.gctsc
> 
> ### ** Examples
> 
> # Simulate data from a Poisson AR(1) model
> set.seed(123)
> n <- 2000
> mu <- 5
> phi <- 0.5
> arma_order <- c(1, 0)
> y <- sim_poisson(mu = mu, tau = phi, arma_order = arma_order, nsim = n)$y
> 
> # Fit the model using the CE method
> fit <- gctsc(y~1,
+   marginal = poisson.marg(link = "identity", lambda.lower = 0),
+   cormat = arma.cormat(p = 1, q = 0), family ="gaussian",
+   method = "CE",
+   options = gctsc.opts(seed = 1, M = 1000),
+   c = 0.5
+ )
> 
> # Produce diagnostic plots
> par(mfrow = c(2, 3), ask = FALSE)
> plot(fit)
> 
> 
> 
> graphics::par(get("par.postscript", pos = 'CheckExEnv'))
> cleanEx()
> nameEx("pmv_ce")
> ### * pmv_ce
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: pmvn_ce
> ### Title: Approximate Log-Likelihood via Continuous Extension (CE)
> ### Aliases: pmvn_ce pmvt_ce
> 
> ### ** Examples
> 
> ## Gaussian copula example
> mu <- 10
> tau <- 0.2
> arma_order <- c(1, 0)
> 
> sim_data <- sim_poisson(mu = mu, tau = tau, arma_order = arma_order, 
+                         nsim = 500, family = "gaussian", seed = 1)
> 
> y <- sim_data$y
> a <- qnorm(ppois(y - 1, lambda = mu))
> b <- qnorm(ppois(y, lambda = mu))
> 
> llk_gauss <- pmvn_ce(lower = a, upper = b,
+                      tau = tau, od = arma_order, c = 0.5)
> 
> 
> ## t copula example
> df <- 8
> 
> sim_data_t <- sim_poisson(mu = mu, tau = tau, arma_order = arma_order,
+                           nsim = 500, family = "t", df = df, seed = 1)
> 
> y_t <- sim_data_t$y
> a_t <- qt(ppois(y_t - 1, lambda = mu), df = df)
> b_t <- qt(ppois(y_t, lambda = mu), df = df)
> 
> llk_t <- pmvt_ce(lower = a_t, upper = b_t, tau = tau, od = arma_order,
+                  c = 0.5, df = df)
> 
> 
> 
> cleanEx()
> nameEx("pmv_ghk")
> ### * pmv_ghk
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: pmvn_ghk
> ### Title: Approximate Log-Likelihood via GHK Simulation
> ### Aliases: pmvn_ghk pmvt_ghk
> 
> ### ** Examples
> 
> ## Gaussian copula example
> mu <- 10
> tau <- 0.2
> arma_order <- c(1, 0)
> 
> sim_data <- sim_poisson(mu = mu, tau = tau, arma_order = arma_order,
+                         nsim = 500, family = "gaussian", seed = 1)
> 
> y <- sim_data$y
> a <- qnorm(ppois(y - 1, lambda = mu))
> b <- qnorm(ppois(y, lambda = mu))
> 
> llk_gauss <- pmvn_ghk(lower = a, upper = b, tau = tau, od = arma_order,
+                       M = 1000)
> 
> 
> ## Student--t copula example
> df <- 8
> 
> sim_data_t <- sim_poisson(mu = mu, tau = tau, arma_order = arma_order,
+                           nsim = 500, family = "t", df = df, seed = 1)
> 
> y_t <- sim_data_t$y
> a_t <- qt(ppois(y_t - 1, lambda = mu), df = df)
> b_t <- qt(ppois(y_t, lambda = mu), df = df)
> 
> llk_t <- pmvt_ghk(lower = a_t, upper = b_t, tau = tau, od = arma_order,
+                   M = 1000, df = df)
> 
> 
> 
> cleanEx()
> nameEx("pmv_tmet")
> ### * pmv_tmet
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: pmvn_tmet
> ### Title: Approximate Log-Likelihood via TMET
> ### Aliases: pmvn_tmet pmvt_tmet
> 
> ### ** Examples
> 
> ## Gaussian copula example
> mu <- 10
> tau <- 0.2
> arma_order <- c(1, 0)
> 
> sim_data <- sim_poisson(mu = mu, tau = tau, arma_order = arma_order,
+                         nsim = 1000, seed = 1)
> 
> y <- sim_data$y
> a <- qnorm(ppois(y - 1, lambda = mu))
> b <- qnorm(ppois(y, lambda = mu))
> 
> # Approximate log-likelihood using TMET
> llk_tmet <- pmvn_tmet(lower = a, upper = b,
+                       tau = tau, od = arma_order)
> llk_tmet
[1] -2578.77
> 
> ## Student--t copula example
> df <- 8
> 
> sim_data_t <- sim_poisson(mu = mu, tau = tau, arma_order = arma_order,
+                           nsim = 500, family = "t", df = df, seed = 1)
> 
> y_t <- sim_data_t$y
> a_t <- qt(ppois(y_t - 1, lambda = mu), df = df)
> b_t <- qt(ppois(y_t, lambda = mu), df = df)
> 
> llk_t <- pmvt_tmet(lower = a_t, upper = b_t, tau = tau, od = arma_order,
+                   M = 1000, df = df)
>                   
> 
> 
> 
> cleanEx()
> nameEx("predict.gctsc")
> ### * predict.gctsc
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: predict.gctsc
> ### Title: One-Step-Ahead Predictive Distribution for Copula Count Time
> ###   Series Models
> ### Aliases: predict.gctsc
> 
> ### ** Examples
> 
> # Simulate Poisson AR(1) data
> set.seed(1)
> y_sim <- sim_poisson(mu = 10, tau = 0.2,
+                      arma_order = c(1, 0),
+                      nsim = 1000,
+                      family = "gaussian")$y
> 
> # Fit Gaussian copula model
> fit <- gctsc(
+   y ~ 1,
+   data = data.frame(y = y_sim),
+   marginal = poisson.marg(lambda.lower = 0),
+   cormat = arma.cormat(p = 1, q = 0),
+   method = "CE",
+   family = "gaussian",
+   options = gctsc.opts(M = 1000, seed = 42)
+ )
> 
> # One-step-ahead prediction
> predict(fit)
==2589555== Conditional jump or move depends on uninitialised value(s)
==2589555==    at 0x1843C06C: mvphi_ (packages/tests-vg/gctsc/src/mvphi.f:34)
==2589555==    by 0x184292E6: gauss_utils::norm_cdf_vec(int, double const*, double*) (packages/tests-vg/gctsc/src/gauss_utils.cpp:34)
==2589555==    by 0x18442208: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:157)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==    by 0x53D1EB: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==2589555==  Uninitialised value was created by a heap allocation
==2589555==    at 0x484577A: operator new[](unsigned long) (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:729)
==2589555==    by 0x18441143: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:44)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==    by 0x53D1EB: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==2589555==    by 0x53D850: Rf_usemethod (svn/R-devel/src/main/objects.c:513)
==2589555== 
==2589555== Conditional jump or move depends on uninitialised value(s)
==2589555==    at 0x4C4109F: __ieee754_exp_avx (/usr/src/debug/glibc-2.39-38.fc40.x86_64/math/../sysdeps/ieee754/dbl-64/e_exp.c:105)
==2589555==    by 0x4C07DD6: exp@@GLIBC_2.29 (/usr/src/debug/glibc-2.39-38.fc40.x86_64/math/./w_exp_template.c:32)
==2589555==    by 0x1843C10A: mvphi_ (packages/tests-vg/gctsc/src/mvphi.f:45)
==2589555==    by 0x184292E6: gauss_utils::norm_cdf_vec(int, double const*, double*) (packages/tests-vg/gctsc/src/gauss_utils.cpp:34)
==2589555==    by 0x18442208: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:157)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==  Uninitialised value was created by a heap allocation
==2589555==    at 0x484577A: operator new[](unsigned long) (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:729)
==2589555==    by 0x18441143: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:44)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==    by 0x53D1EB: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==2589555==    by 0x53D850: Rf_usemethod (svn/R-devel/src/main/objects.c:513)
==2589555== 
==2589555== Use of uninitialised value of size 8
==2589555==    at 0x4C410F1: __ieee754_exp_avx (/usr/src/debug/glibc-2.39-38.fc40.x86_64/math/../sysdeps/ieee754/dbl-64/e_exp.c:145)
==2589555==    by 0x4C07DD6: exp@@GLIBC_2.29 (/usr/src/debug/glibc-2.39-38.fc40.x86_64/math/./w_exp_template.c:32)
==2589555==    by 0x1843C10A: mvphi_ (packages/tests-vg/gctsc/src/mvphi.f:45)
==2589555==    by 0x184292E6: gauss_utils::norm_cdf_vec(int, double const*, double*) (packages/tests-vg/gctsc/src/gauss_utils.cpp:34)
==2589555==    by 0x18442208: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:157)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==  Uninitialised value was created by a heap allocation
==2589555==    at 0x484577A: operator new[](unsigned long) (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:729)
==2589555==    by 0x18441143: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:44)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==    by 0x53D1EB: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==2589555==    by 0x53D850: Rf_usemethod (svn/R-devel/src/main/objects.c:513)
==2589555== 
==2589555== Conditional jump or move depends on uninitialised value(s)
==2589555==    at 0x4C41150: __ieee754_exp_avx (/usr/src/debug/glibc-2.39-38.fc40.x86_64/math/../sysdeps/ieee754/dbl-64/e_exp.c:152)
==2589555==    by 0x4C07DD6: exp@@GLIBC_2.29 (/usr/src/debug/glibc-2.39-38.fc40.x86_64/math/./w_exp_template.c:32)
==2589555==    by 0x1843C10A: mvphi_ (packages/tests-vg/gctsc/src/mvphi.f:45)
==2589555==    by 0x184292E6: gauss_utils::norm_cdf_vec(int, double const*, double*) (packages/tests-vg/gctsc/src/gauss_utils.cpp:34)
==2589555==    by 0x18442208: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:157)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==  Uninitialised value was created by a heap allocation
==2589555==    at 0x484577A: operator new[](unsigned long) (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:729)
==2589555==    by 0x18441143: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:44)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==    by 0x53D1EB: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==2589555==    by 0x53D850: Rf_usemethod (svn/R-devel/src/main/objects.c:513)
==2589555== 
==2589555== Conditional jump or move depends on uninitialised value(s)
==2589555==    at 0x4C07DF3: exp@@GLIBC_2.29 (/usr/src/debug/glibc-2.39-38.fc40.x86_64/math/./w_exp_template.c:33)
==2589555==    by 0x1843C10A: mvphi_ (packages/tests-vg/gctsc/src/mvphi.f:45)
==2589555==    by 0x184292E6: gauss_utils::norm_cdf_vec(int, double const*, double*) (packages/tests-vg/gctsc/src/gauss_utils.cpp:34)
==2589555==    by 0x18442208: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:157)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555== 
==2589555== Conditional jump or move depends on uninitialised value(s)
==2589555==    at 0x4C07DFD: exp@@GLIBC_2.29 (/usr/src/debug/glibc-2.39-38.fc40.x86_64/math/./w_exp_template.c:33)
==2589555==    by 0x1843C10A: mvphi_ (packages/tests-vg/gctsc/src/mvphi.f:45)
==2589555==    by 0x184292E6: gauss_utils::norm_cdf_vec(int, double const*, double*) (packages/tests-vg/gctsc/src/gauss_utils.cpp:34)
==2589555==    by 0x18442208: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:157)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555== 
==2589555== Conditional jump or move depends on uninitialised value(s)
==2589555==    at 0x4C07DFF: exp@@GLIBC_2.29 (/usr/src/debug/glibc-2.39-38.fc40.x86_64/math/./w_exp_template.c:33)
==2589555==    by 0x1843C10A: mvphi_ (packages/tests-vg/gctsc/src/mvphi.f:45)
==2589555==    by 0x184292E6: gauss_utils::norm_cdf_vec(int, double const*, double*) (packages/tests-vg/gctsc/src/gauss_utils.cpp:34)
==2589555==    by 0x18442208: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:157)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555== 
==2589555== Conditional jump or move depends on uninitialised value(s)
==2589555==    at 0x1843C135: mvphi_ (packages/tests-vg/gctsc/src/mvphi.f:47)
==2589555==    by 0x184292E6: gauss_utils::norm_cdf_vec(int, double const*, double*) (packages/tests-vg/gctsc/src/gauss_utils.cpp:34)
==2589555==    by 0x18442208: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:157)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==    by 0x53D1EB: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==2589555==  Uninitialised value was created by a heap allocation
==2589555==    at 0x484577A: operator new[](unsigned long) (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:729)
==2589555==    by 0x18441143: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:44)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==    by 0x53D1EB: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==2589555==    by 0x53D850: Rf_usemethod (svn/R-devel/src/main/objects.c:513)
==2589555== 
==2589555== Conditional jump or move depends on uninitialised value(s)
==2589555==    at 0x1843C06C: mvphi_ (packages/tests-vg/gctsc/src/mvphi.f:34)
==2589555==    by 0x184292E6: gauss_utils::norm_cdf_vec(int, double const*, double*) (packages/tests-vg/gctsc/src/gauss_utils.cpp:34)
==2589555==    by 0x18442228: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:158)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==    by 0x53D1EB: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==2589555==  Uninitialised value was created by a heap allocation
==2589555==    at 0x484577A: operator new[](unsigned long) (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:729)
==2589555==    by 0x18441143: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:44)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==    by 0x53D1EB: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==2589555==    by 0x53D850: Rf_usemethod (svn/R-devel/src/main/objects.c:513)
==2589555== 
==2589555== Conditional jump or move depends on uninitialised value(s)
==2589555==    at 0x4C07DF3: exp@@GLIBC_2.29 (/usr/src/debug/glibc-2.39-38.fc40.x86_64/math/./w_exp_template.c:33)
==2589555==    by 0x1843C10A: mvphi_ (packages/tests-vg/gctsc/src/mvphi.f:45)
==2589555==    by 0x184292E6: gauss_utils::norm_cdf_vec(int, double const*, double*) (packages/tests-vg/gctsc/src/gauss_utils.cpp:34)
==2589555==    by 0x18442228: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:158)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555== 
==2589555== Conditional jump or move depends on uninitialised value(s)
==2589555==    at 0x4C07DFD: exp@@GLIBC_2.29 (/usr/src/debug/glibc-2.39-38.fc40.x86_64/math/./w_exp_template.c:33)
==2589555==    by 0x1843C10A: mvphi_ (packages/tests-vg/gctsc/src/mvphi.f:45)
==2589555==    by 0x184292E6: gauss_utils::norm_cdf_vec(int, double const*, double*) (packages/tests-vg/gctsc/src/gauss_utils.cpp:34)
==2589555==    by 0x18442228: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:158)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555== 
==2589555== Conditional jump or move depends on uninitialised value(s)
==2589555==    at 0x4C07DFF: exp@@GLIBC_2.29 (/usr/src/debug/glibc-2.39-38.fc40.x86_64/math/./w_exp_template.c:33)
==2589555==    by 0x1843C10A: mvphi_ (packages/tests-vg/gctsc/src/mvphi.f:45)
==2589555==    by 0x184292E6: gauss_utils::norm_cdf_vec(int, double const*, double*) (packages/tests-vg/gctsc/src/gauss_utils.cpp:34)
==2589555==    by 0x18442228: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:158)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555== 
==2589555== Conditional jump or move depends on uninitialised value(s)
==2589555==    at 0x1843C135: mvphi_ (packages/tests-vg/gctsc/src/mvphi.f:47)
==2589555==    by 0x184292E6: gauss_utils::norm_cdf_vec(int, double const*, double*) (packages/tests-vg/gctsc/src/gauss_utils.cpp:34)
==2589555==    by 0x18442228: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:158)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==    by 0x53D1EB: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==2589555==  Uninitialised value was created by a heap allocation
==2589555==    at 0x484577A: operator new[](unsigned long) (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:729)
==2589555==    by 0x18441143: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:44)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==    by 0x53D1EB: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==2589555==    by 0x53D850: Rf_usemethod (svn/R-devel/src/main/objects.c:513)
==2589555== 
==2589555== Conditional jump or move depends on uninitialised value(s)
==2589555==    at 0x1843C1A0: mvphnv_ (packages/tests-vg/gctsc/src/mvphi.f:131)
==2589555==    by 0x18429326: gauss_utils::norm_inv_vec(int, double const*, double*) (packages/tests-vg/gctsc/src/gauss_utils.cpp:39)
==2589555==    by 0x184422FD: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:169)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==    by 0x53D1EB: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==2589555==  Uninitialised value was created by a heap allocation
==2589555==    at 0x484577A: operator new[](unsigned long) (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:729)
==2589555==    by 0x18441143: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:44)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==    by 0x53D1EB: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==2589555==    by 0x53D850: Rf_usemethod (svn/R-devel/src/main/objects.c:513)
==2589555== 
==2589555== Conditional jump or move depends on uninitialised value(s)
==2589555==    at 0x4C07CBC: log@@GLIBC_2.29 (/usr/src/debug/glibc-2.39-38.fc40.x86_64/math/./w_log_template.c:32)
==2589555==    by 0x18442350: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:174)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==    by 0x53D1EB: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==2589555==    by 0x53D850: Rf_usemethod (svn/R-devel/src/main/objects.c:513)
==2589555==  Uninitialised value was created by a heap allocation
==2589555==    at 0x484577A: operator new[](unsigned long) (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:729)
==2589555==    by 0x18441143: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:44)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==    by 0x53D1EB: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==2589555==    by 0x53D850: Rf_usemethod (svn/R-devel/src/main/objects.c:513)
==2589555== 
==2589555== Conditional jump or move depends on uninitialised value(s)
==2589555==    at 0x4C418DA: __ieee754_log_avx (/usr/src/debug/glibc-2.39-38.fc40.x86_64/math/../sysdeps/ieee754/dbl-64/e_log.c:61)
==2589555==    by 0x18442350: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:174)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==    by 0x53D1EB: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==2589555==    by 0x53D850: Rf_usemethod (svn/R-devel/src/main/objects.c:513)
==2589555==  Uninitialised value was created by a heap allocation
==2589555==    at 0x484577A: operator new[](unsigned long) (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:729)
==2589555==    by 0x18441143: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:44)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==    by 0x53D1EB: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==2589555==    by 0x53D850: Rf_usemethod (svn/R-devel/src/main/objects.c:513)
==2589555== 
==2589555== Conditional jump or move depends on uninitialised value(s)
==2589555==    at 0x4C418E9: __ieee754_log_avx (/usr/src/debug/glibc-2.39-38.fc40.x86_64/math/../sysdeps/ieee754/dbl-64/e_log.c:85)
==2589555==    by 0x18442350: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:174)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==    by 0x53D1EB: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==2589555==    by 0x53D850: Rf_usemethod (svn/R-devel/src/main/objects.c:513)
==2589555==  Uninitialised value was created by a heap allocation
==2589555==    at 0x484577A: operator new[](unsigned long) (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:729)
==2589555==    by 0x18441143: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:44)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==    by 0x53D1EB: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==2589555==    by 0x53D850: Rf_usemethod (svn/R-devel/src/main/objects.c:513)
==2589555== 
==2589555== Use of uninitialised value of size 8
==2589555==    at 0x4C41966: __ieee754_log_avx (/usr/src/debug/glibc-2.39-38.fc40.x86_64/math/../sysdeps/ieee754/dbl-64/e_log.c:117)
==2589555==    by 0x18442350: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:174)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==    by 0x53D1EB: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==2589555==    by 0x53D850: Rf_usemethod (svn/R-devel/src/main/objects.c:513)
==2589555==  Uninitialised value was created by a heap allocation
==2589555==    at 0x484577A: operator new[](unsigned long) (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:729)
==2589555==    by 0x18441143: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:44)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==    by 0x53D1EB: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==2589555==    by 0x53D850: Rf_usemethod (svn/R-devel/src/main/objects.c:513)
==2589555== 
==2589555== Use of uninitialised value of size 8
==2589555==    at 0x4C4196F: __ieee754_log_avx (/usr/src/debug/glibc-2.39-38.fc40.x86_64/math/../sysdeps/ieee754/dbl-64/e_log.c:117)
==2589555==    by 0x18442350: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:174)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==    by 0x53D1EB: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==2589555==    by 0x53D850: Rf_usemethod (svn/R-devel/src/main/objects.c:513)
==2589555==  Uninitialised value was created by a heap allocation
==2589555==    at 0x484577A: operator new[](unsigned long) (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:729)
==2589555==    by 0x18441143: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:44)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==    by 0x53D1EB: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==2589555==    by 0x53D850: Rf_usemethod (svn/R-devel/src/main/objects.c:513)
==2589555== 
==2589555== Conditional jump or move depends on uninitialised value(s)
==2589555==    at 0x1843C294: mvphnv_ (packages/tests-vg/gctsc/src/mvphi.f:139)
==2589555==    by 0x18429326: gauss_utils::norm_inv_vec(int, double const*, double*) (packages/tests-vg/gctsc/src/gauss_utils.cpp:39)
==2589555==    by 0x184422FD: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:169)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==    by 0x53D1EB: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==2589555==  Uninitialised value was created by a heap allocation
==2589555==    at 0x484577A: operator new[](unsigned long) (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:729)
==2589555==    by 0x18441143: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:44)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==    by 0x53D1EB: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==2589555==    by 0x53D850: Rf_usemethod (svn/R-devel/src/main/objects.c:513)
==2589555== 
==2589555== Conditional jump or move depends on uninitialised value(s)
==2589555==    at 0x4C07CBC: log@@GLIBC_2.29 (/usr/src/debug/glibc-2.39-38.fc40.x86_64/math/./w_log_template.c:32)
==2589555==    by 0x1843C2A4: mvphnv_ (packages/tests-vg/gctsc/src/mvphi.f:140)
==2589555==    by 0x18429326: gauss_utils::norm_inv_vec(int, double const*, double*) (packages/tests-vg/gctsc/src/gauss_utils.cpp:39)
==2589555==    by 0x184422FD: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:169)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==  Uninitialised value was created by a heap allocation
==2589555==    at 0x484577A: operator new[](unsigned long) (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:729)
==2589555==    by 0x18441143: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:44)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==    by 0x53D1EB: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==2589555==    by 0x53D850: Rf_usemethod (svn/R-devel/src/main/objects.c:513)
==2589555== 
==2589555== Conditional jump or move depends on uninitialised value(s)
==2589555==    at 0x4C418DA: __ieee754_log_avx (/usr/src/debug/glibc-2.39-38.fc40.x86_64/math/../sysdeps/ieee754/dbl-64/e_log.c:61)
==2589555==    by 0x1843C2A4: mvphnv_ (packages/tests-vg/gctsc/src/mvphi.f:140)
==2589555==    by 0x18429326: gauss_utils::norm_inv_vec(int, double const*, double*) (packages/tests-vg/gctsc/src/gauss_utils.cpp:39)
==2589555==    by 0x184422FD: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:169)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==  Uninitialised value was created by a heap allocation
==2589555==    at 0x484577A: operator new[](unsigned long) (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:729)
==2589555==    by 0x18441143: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:44)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==    by 0x53D1EB: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==2589555==    by 0x53D850: Rf_usemethod (svn/R-devel/src/main/objects.c:513)
==2589555== 
==2589555== Conditional jump or move depends on uninitialised value(s)
==2589555==    at 0x4C418E9: __ieee754_log_avx (/usr/src/debug/glibc-2.39-38.fc40.x86_64/math/../sysdeps/ieee754/dbl-64/e_log.c:85)
==2589555==    by 0x1843C2A4: mvphnv_ (packages/tests-vg/gctsc/src/mvphi.f:140)
==2589555==    by 0x18429326: gauss_utils::norm_inv_vec(int, double const*, double*) (packages/tests-vg/gctsc/src/gauss_utils.cpp:39)
==2589555==    by 0x184422FD: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:169)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==  Uninitialised value was created by a heap allocation
==2589555==    at 0x484577A: operator new[](unsigned long) (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:729)
==2589555==    by 0x18441143: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:44)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==    by 0x53D1EB: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==2589555==    by 0x53D850: Rf_usemethod (svn/R-devel/src/main/objects.c:513)
==2589555== 
==2589555== Use of uninitialised value of size 8
==2589555==    at 0x4C41966: __ieee754_log_avx (/usr/src/debug/glibc-2.39-38.fc40.x86_64/math/../sysdeps/ieee754/dbl-64/e_log.c:117)
==2589555==    by 0x1843C2A4: mvphnv_ (packages/tests-vg/gctsc/src/mvphi.f:140)
==2589555==    by 0x18429326: gauss_utils::norm_inv_vec(int, double const*, double*) (packages/tests-vg/gctsc/src/gauss_utils.cpp:39)
==2589555==    by 0x184422FD: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:169)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==  Uninitialised value was created by a heap allocation
==2589555==    at 0x484577A: operator new[](unsigned long) (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:729)
==2589555==    by 0x18441143: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:44)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==    by 0x53D1EB: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==2589555==    by 0x53D850: Rf_usemethod (svn/R-devel/src/main/objects.c:513)
==2589555== 
==2589555== Use of uninitialised value of size 8
==2589555==    at 0x4C4196F: __ieee754_log_avx (/usr/src/debug/glibc-2.39-38.fc40.x86_64/math/../sysdeps/ieee754/dbl-64/e_log.c:117)
==2589555==    by 0x1843C2A4: mvphnv_ (packages/tests-vg/gctsc/src/mvphi.f:140)
==2589555==    by 0x18429326: gauss_utils::norm_inv_vec(int, double const*, double*) (packages/tests-vg/gctsc/src/gauss_utils.cpp:39)
==2589555==    by 0x184422FD: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:169)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==  Uninitialised value was created by a heap allocation
==2589555==    at 0x484577A: operator new[](unsigned long) (/builddir/build/BUILD/valgrind-3.24.0/coregrind/m_replacemalloc/vg_replace_malloc.c:729)
==2589555==    by 0x18441143: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:44)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==    by 0x53D1EB: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==2589555==    by 0x53D850: Rf_usemethod (svn/R-devel/src/main/objects.c:513)
==2589555== 
==2589555== Conditional jump or move depends on uninitialised value(s)
==2589555==    at 0x1843C2CF: mvphnv_ (packages/tests-vg/gctsc/src/mvphi.f:141)
==2589555==    by 0x18429326: gauss_utils::norm_inv_vec(int, double const*, double*) (packages/tests-vg/gctsc/src/gauss_utils.cpp:39)
==2589555==    by 0x184422FD: predmvn_ghk(Rcpp::Vector<19, Rcpp::PreserveStorage>, Rcpp::Vector<19, Rcpp::PreserveStorage>) (packages/tests-vg/gctsc/src/pred_mvn.cpp:169)
==2589555==    by 0x1841BF95: _gctsc_predmvn_ghk (packages/tests-vg/gctsc/src/RcppExports.cpp:113)
==2589555==    by 0x4A8029: R_doDotCall (svn/R-devel/src/main/dotcode.c:757)
==2589555==    by 0x4E2103: bcEval_loop (svn/R-devel/src/main/eval.c:8682)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7515)
==2589555==    by 0x4F2257: bcEval (svn/R-devel/src/main/eval.c:7500)
==2589555==    by 0x4F258A: Rf_eval (svn/R-devel/src/main/eval.c:1167)
==2589555==    by 0x4F430D: R_execClosure (svn/R-devel/src/main/eval.c:2389)
==2589555==    by 0x4F4FC6: applyClosure_core (svn/R-devel/src/main/eval.c:2302)
==2589555==    by 0x4F59B8: Rf_applyClosure (svn/R-devel/src/main/eval.c:2324)
==2589555==    by 0x53D1EB: dispatchMethod (svn/R-devel/src/main/objects.c:473)
==2589555== 
==2589555== 
==2589555== More than 10000000 total errors detected.  I'm not reporting any more.
==2589555== Final error counts will be inaccurate.  Go fix your program!
==2589555== Rerun with --error-limit=no to disable this cutoff.  Note
==2589555== that errors may occur in your program without prior warning from
==2589555== Valgrind, because errors are no longer being displayed.
==2589555== 
$mean
[1] 9.969086

$median
[1] 10

$mode
[1] 9

$variance
[1] 9.581621

$p_y
 [1] 3.314302e-05 3.669702e-04 1.966024e-03 6.896520e-03 1.791651e-02
 [6] 3.687656e-02 6.275380e-02 9.092820e-02 1.146240e-01 1.277940e-01
[11] 1.276539e-01 1.154520e-01 9.536170e-02 7.246176e-02 5.096773e-02
[16] 3.336189e-02 2.041705e-02 1.172992e-02 6.349320e-03 3.248549e-03
[21] 1.575572e-03 7.262849e-04 3.189511e-04 1.337291e-04 5.363762e-05
[26] 2.061776e-05 7.607941e-06 2.699085e-06 9.219585e-07 3.036193e-07
[31] 9.651819e-08 2.965201e-08 9.058665e-09 9.058665e-09

$lower
[1] 4

$upper
[1] 16

> 
> 
> 
> 
> cleanEx()
> nameEx("residuals.gctsc")
> ### * residuals.gctsc
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: residuals.gctsc
> ### Title: Randomized Quantile Residuals for Copula Count Time Series
> ###   Models
> ### Aliases: residuals.gctsc
> 
> ### ** Examples
> 
> # Simulate Poisson AR(1) data under a Gaussian copula
> set.seed(1)
> y <- sim_poisson(mu = 5, tau = 0.7,
+                  arma_order = c(1, 0),
+                  nsim = 500,
+                  family = "gaussian")$y
> 
> fit <- gctsc(
+   y ~ 1,
+   data = data.frame(y = y),
+   marginal = poisson.marg(),
+   cormat = arma.cormat(1, 0),
+   family = "gaussian",
+   method = "CE",
+   options = gctsc.opts(seed = 1, M = 1000)
+ )
> 
> res <- residuals(fit)
> hist(res$residuals, main = "Randomized Quantile Residuals")
> hist(res$pit, main = "PIT Histogram")
> 
> 
> 
> 
> cleanEx()
> nameEx("sim_gctsc")
> ### * sim_gctsc
> 
> flush(stderr()); flush(stdout())
> 
> ### Name: sim_gctsc
> ### Title: Simulate from Gaussian and t Copula Time Series Models
> ### Aliases: sim_gctsc sim_poisson sim_negbin sim_zip sim_binom sim_bbinom
> ###   sim_zib sim_zibb
> 
> ### ** Examples
> 
> # Poisson example
> sim_poisson(mu = 10, tau = c(0.2, 0.2), 
+   arma_order = c(1, 1), nsim = 100, 
+   family = "gaussian", seed = 42)
$y
  [1]  9  3  5 10 10  9 12  9 11  7  5  4  8  9 14 11 14 18 11 11  8 11  7  6  9
 [26] 10 11 10  8  9 15 12  7  6  6  6  5  2  8 12 10  8  7 10 11 11 10 12 12 12
 [51] 12 16  9 14 10 10 10  9  3  6  9  5  9  8 14 12  8 10 11 13 11  5  4  6 10
 [76]  9  8  7  9  7  8  7  5 11  7  9  8  6  9 11 14 12 10 18 11  7  7  9 11 12

$z
  [1] -0.13300515 -2.70724545 -1.86871373  0.20859964 -0.08672751 -0.33518193
  [7]  0.68546744 -0.20587789  0.44018808 -1.11235648 -1.62229233 -2.10167018
 [13] -0.65218033 -0.30564379  1.23561793  0.47542362  1.34068129  2.31973154
 [19]  0.36961155  0.34264292 -0.64176564  0.39844936 -0.78332468 -1.30692785
 [25] -0.18408090  0.10073290  0.28010339  0.16884070 -0.44282260 -0.30699412
 [31]  1.41446795  0.56632103 -1.11073331 -1.24819423 -1.42610336 -1.32255869
 [37] -1.59871289 -2.98431364 -0.46364095  0.64223841  0.01865270 -0.71277775
 [43] -1.02851085  0.01306728  0.28620182  0.48955789  0.18046386  0.61005634
 [49]  0.63577729  0.72388277  0.70732684  1.77190677 -0.31466945  1.13814897
 [55]  0.19750016 -0.08254386  0.04652605 -0.39064767 -2.52548859 -1.26769011
 [61] -0.27316851 -1.65579460 -0.14377844 -0.47989464  1.13377594  0.74004510
 [67] -0.53321791 -0.09740792  0.25005201  0.94209731  0.21206221 -1.86917478
 [73] -1.95359972 -1.42503118 -0.04909324 -0.19459679 -0.66045553 -1.02806697
 [79] -0.21976290 -0.80376692 -0.68681133 -1.08129274 -1.67558173  0.29215810
 [85] -0.98446635 -0.37323418 -0.58716234 -1.20257495 -0.19642428  0.35418610
 [91]  1.31078222  0.79742707  0.09150955  2.40208076  0.36025336 -0.97712580
 [97] -0.87659782 -0.37662743  0.46012003  0.58819978

$marginal
[1] "poisson"

$parameters
$parameters$mu
[1] 10


$cormat
$cormat$arma_order
[1] 1 1

$cormat$tau
[1] 0.2 0.2


> 
> # Negative Binomial example
> sim_negbin(mu = 10, dispersion = 2, tau = c(0.5, 0.5),
+   arma_order = c(1, 1),family = "gaussian",
+   nsim = 100, seed =1)
$y
  [1]  1  0  0  0  1  1  0  0  5 10 13 22 28 22  3  0  2  0  1  2  1  0  0  5 12
 [26]  8 13 41 30 20 16  3  0 10 32 34 22 18 14  9  9  2 21 46 20 26 12 26 18  0
 [51]  1 19  8  0  0  0  0  0  0  4  2  1  3 24 11  0  0  0  1  7  7  8 10  8  3
 [76] 10  9 10 12 34 45 18  8  5  0  3 23  9  6  9  1  8 11  7 10  2  0  0  3  6

$z
  [1] -0.52920915 -1.39723738 -2.17940118 -1.33491569 -0.66686346 -0.77774644
  [7] -1.45578947 -1.88290815  0.03231338  0.46889043  0.63428176  1.08440088
 [13]  1.29881595  1.09223654 -0.17724232 -1.15942620 -0.41824976 -0.81969281
 [19] -0.49274956 -0.42298039 -0.64693033 -0.97075887 -0.89448944  0.05167619
 [25]  0.61000232  0.30133670  0.64412891  1.69094953  1.38335978  0.99760282
 [31]  0.83494209 -0.13164827 -0.99251103  0.45723276  1.43394051  1.49124268
 [37]  1.08513409  0.90863319  0.71261041  0.43350640  0.37676476 -0.39999018
 [43]  1.04122143  1.83160238  1.01385059  1.21872067  0.59573477  1.21276518
 [49]  0.90831259 -1.19470526 -0.46421325  0.95562415  0.34151203 -1.00600564
 [55] -2.16563661 -1.85941911 -2.21835298 -2.24266501 -1.28084665 -0.02681298
 [61] -0.27016396 -0.48439492 -0.11575316  1.15087626  0.51057481 -1.10873722
 [67] -2.28773739 -1.38417080 -0.54953495  0.21100303  0.26602506  0.34425240
 [73]  0.44758190  0.31259556 -0.24382109  0.50035613  0.39762614  0.50249169
 [79]  0.57727714  1.50881092  1.80314506  0.91330119  0.34012582  0.08503408
 [85] -1.36728010 -0.25682494  1.12366380  0.40820036  0.20324863  0.38502368
 [91] -0.46987424  0.36965368  0.50999764  0.29392463  0.47210448 -0.28222985
 [97] -1.17709883 -0.97561813 -0.12656179  0.17362081

$marginal
[1] "negbin"

$parameters
$parameters$mu
[1] 10

$parameters$dispersion
[1] 2


$cormat
$cormat$arma_order
[1] 1 1

$cormat$tau
[1] 0.5 0.5


> 
> # Zero Inflated Beta-Binomial example with seasonal covariates
> n <- 100
> xi <- numeric(n)
> zeta <- rnorm(n)
> for (j in 3:n) {
+   xi[j] <- 0.6 * xi[j - 1] - 0.4 * xi[j - 2] + zeta[j]
+ }
> prob <- plogis(0.2 + 0.3 * sin(2 * pi * (1:n) / 12) +
+              0.5 * cos(2 * pi * (1:n) / 12) + 0.3 * xi)
> sim_zibb(prob, rho = 1/6, pi0 = 0.2, size = 24, tau = 0.5,
+  arma_order = c(1, 0),family = "t", df = 10, nsim = 100)
$y
  [1] 20 10 21 16  9  7 15 11 10 17 12  0  0 21 14 16  3  0  6  8 15 18 10 15 23
 [26] 24 16 11  6 17  6  0  3  0  0  0  8 19 23 19  3 19 18 13 18 20 18 12 20  0
 [51]  0  0  0  0  0  0  0  7 15 17 21 21 22 17  0  6 14 22 15 19  0  0 16 15 22
 [76] 21  5  8  8  7 10  3  0 11 21 20 16  7 16 17 23 19 13 18 21  0 23 13 22  0

$z
  [1]  0.826127205 -0.487578730  1.218164086  0.965308659 -0.001355709
  [6] -0.044213057  1.351040676  0.553323975  0.510791967  1.157962318
 [11] -0.048023130 -2.367758143 -1.208943397  0.809135621  0.203699677
 [16]  0.869801459 -0.198690948 -1.160395290 -0.074671436 -0.482563666
 [21]  0.732951040  1.468668104  0.259425764  0.246300746  1.188032570
 [26]  2.282082915  0.682080896 -0.217660815 -0.474547532  1.162132909
 [31] -0.289778758 -0.895786858 -0.741934937 -1.870258620 -1.426258259
 [36] -1.174764971 -0.518135317  0.528157167  2.079162982  1.977885256
 [41] -0.481317715  1.499544949  1.683077806  0.916224162  2.012551417
 [46]  1.632997101  0.638317606 -0.354325612  0.848204139 -2.133168777
 [51] -1.245717436 -1.926399989 -3.299527300 -1.392288693 -0.844068230
 [56] -2.762904392 -1.704909397 -0.365982440  0.560524926  0.620723560
 [61]  1.142071249  0.793673213  0.821791814  0.279746576 -1.058979019
 [66] -0.237360082  0.948272033  3.105169206  1.042047245  0.534615209
 [71] -2.192625123 -2.012417350  0.484348964 -0.051973513  0.631947044
 [76]  1.087900298 -0.554319468 -0.259223479  0.139319117 -0.290818515
 [81]  0.092545838 -0.710083016 -1.688819415 -0.234350283  1.494935058
 [86]  0.957528838  0.398615081 -0.141250710  1.564744094  2.054553157
 [91]  2.651521251  2.046432215  1.340039469  1.247875417  0.928500211
 [96] -0.923723488  2.026740064 -0.311400692  0.945670395 -0.956548849

$marginal
[1] "zibb"

$parameters
$parameters$prob
  [1] 0.6863291 0.6703586 0.6600723 0.4718064 0.4653147 0.3914002 0.3577349
  [8] 0.3998695 0.3609605 0.4790356 0.5926847 0.6293253 0.6358514 0.7342628
 [15] 0.5978461 0.4901190 0.2823565 0.2865110 0.3633346 0.5722659 0.5107741
 [22] 0.4437483 0.4394952 0.6188499 0.7842371 0.6771907 0.5468435 0.6034885
 [29] 0.4906652 0.4623140 0.4513921 0.3433873 0.4967060 0.6493561 0.6358584
 [36] 0.5751872 0.5981944 0.7071886 0.6428896 0.4089080 0.3756557 0.4769856
 [43] 0.3926909 0.3784495 0.3641832 0.5063529 0.6288524 0.6861921 0.6520420
 [50] 0.5917143 0.6502589 0.6103909 0.4890829 0.2081664 0.1973339 0.4308099
 [57] 0.5510027 0.4995067 0.5284760 0.6089262 0.6515489 0.7041994 0.7528124
 [64] 0.6924097 0.6002045 0.4143521 0.4126447 0.3617568 0.4234765 0.6950612
 [71] 0.7196628 0.6526697 0.5944313 0.6744654 0.7863936 0.6916811 0.4690866
 [78] 0.5145915 0.4048462 0.4712615 0.4728648 0.4503118 0.6498093 0.5979730
 [85] 0.6117509 0.6544484 0.6120517 0.4305603 0.3655316 0.2989266 0.5124615
 [92] 0.4046398 0.2787239 0.4798987 0.7177203 0.7348871 0.6753834 0.6833629
 [99] 0.7262032 0.4985306

$parameters$rho
[1] 0.1666667

$parameters$size
[1] 24

$parameters$pi0
[1] 0.2


$cormat
$cormat$arma_order
[1] 1 0

$cormat$tau
[1] 0.5


> 
> 
> 
> 
> ### * <FOOTER>
> ###
> cleanEx()
> options(digits = 7L)
> base::cat("Time elapsed: ", proc.time() - base::get("ptime", pos = 'CheckExEnv'),"\n")
Time elapsed:  310.78 1.011 313.652 0.006 0.027 
> grDevices::dev.off()
null device 
          1 
> ###
> ### Local variables: ***
> ### mode: outline-minor ***
> ### outline-regexp: "\\(> \\)?### [*]+" ***
> ### End: ***
> quit('no')
==2589555== 
==2589555== HEAP SUMMARY:
==2589555==     in use at exit: 260,806,358 bytes in 56,172 blocks
==2589555==   total heap usage: 132,968 allocs, 76,796 frees, 679,083,547 bytes allocated
==2589555== 
==2589555== LEAK SUMMARY:
==2589555==    definitely lost: 0 bytes in 0 blocks
==2589555==    indirectly lost: 0 bytes in 0 blocks
==2589555==      possibly lost: 0 bytes in 0 blocks
==2589555==    still reachable: 260,806,358 bytes in 56,172 blocks
==2589555==         suppressed: 0 bytes in 0 blocks
==2589555== Reachable blocks (those to which a pointer was found) are not shown.
==2589555== To see them, rerun with: --leak-check=full --show-leak-kinds=all
==2589555== 
==2589555== For lists of detected and suppressed errors, rerun with: -s
==2589555== ERROR SUMMARY: 10000000 errors from 26 contexts (suppressed: 0 from 0)