Running examples in 'bzinb-Ex.R' failed
The error most likely occurred in:

> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: bzinb.se
> ### Title: The bivariate zero-inflated negative binomial distribution -
> ###   Standard error estimation
> ### Aliases: bzinb.se
> 
> ### ** Examples
> 
> set.seed(1)
> data1 <- rbzinb(n = 20, a0 = 1, a1 = 1, a2 = 1, 
+                 b1 = 1, b2 = 1, p1 = 0.5, p2 = 0.2, 
+                 p3 = 0.2, p4 = 0.1)
> bzinb.se(xvec = data1[,1], yvec = data1[,2], 
+          param = c(5.5, 0.017, 0.017, 0.33, 0.36, 
+                    0.53, 0.30, 0.08, 0.09))
 ----------- FAILURE REPORT -------------- 
 --- failure: the condition has length > 1 ---
 --- srcref --- 
: 
 --- package (from environment) --- 
bzinb
 --- call from context --- 
bzinb.se(xvec = data1[, 1], yvec = data1[, 2], param = c(5.5, 
    0.017, 0.017, 0.33, 0.36, 0.53, 0.3, 0.08, 0.09))
 --- call from argument --- 
if (class(cov.mat) == "try-error") {
    std.param = setNames(rep(NA, 11), c(abp.names, "rho", "logit.rho"))
    cov.mat <- NA
} else {
    var.p4 <- sum(cov.mat[6:8, 6:8])
    d.g <- rho * c(1/param[1] - 1/{
        2 * (param[1] + param[2])
    } - 1/{
        2 * (param[1] + param[3])
    }, -1/{
        2 * (param[1] + param[2])
    }, -1/{
        2 * (param[1] + param[3])
    }, 1/{
        2 * param[4] * (param[4] + 1)
    }, 1/{
        2 * param[5] * (param[5] + 1)
    })
    var.rho <- t(d.g) %*% cov.mat[1:5, 1:5] %*% d.g
    var.logit.rho <- var.rho/rho^2/(1 - rho)^2
    std.param = sqrt(c(diag(cov.mat), p4 = var.p4, rho = var.rho, 
        logit.rho = var.logit.rho))
}
 --- R stacktrace ---
where 1: bzinb.se(xvec = data1[, 1], yvec = data1[, 2], param = c(5.5, 
    0.017, 0.017, 0.33, 0.36, 0.53, 0.3, 0.08, 0.09))

 --- value of length: 2 type: logical ---
[1] FALSE FALSE
 --- function from context --- 
function (xvec, yvec, a0, a1, a2, b1, b2, p1, p2, p3, p4, param = NULL, 
    ...) 
{
    .check.input(xvec, yvec)
    if (!is.null(param)) {
        if (length(param) != 9) {
            stop("length(initial) must be 9.")
        }
        .check.ab(param[1:5])
        .check.p(param[6:9])
    }
    if (any(!is.finite(param))) {
        return(rep(NA, 10))
    }
    xy.reduced <- as.data.frame(table(xvec, yvec))
    names(xy.reduced) <- c("x", "y", "freq")
    xy.reduced <- xy.reduced[xy.reduced$freq != 0, ]
    xy.reduced$x <- as.numeric(as.character(xy.reduced$x))
    xy.reduced$y <- as.numeric(as.character(xy.reduced$y))
    xy.reduced$freq <- as.numeric(as.character(xy.reduced$freq))
    n <- sum(xy.reduced$freq)
    n.reduced <- as.integer(length(xy.reduced$freq))
    if (is.null(param)) {
        warning("param was not provided. It will be estimated.")
        est <- bzinb(xvec, yvec, ...)
        param <- unlist(est$coefficients[, 1])
        iter <- est$iter
    }
    else {
        iter <- NA
    }
    if (any(is.na(param))) {
        warning("params include NA.")
        result <- list(rho = matrix(rep(NA, 4), ncol = 2, dimnames = list(c("rho", 
            "logit.rho"), c("Estimate", "Std.err"))), coefficients = matrix(rep(NA, 
            18), ncol = 2, dimnames = list(abp.names, c("Estimate", 
            "Std.err"))), lik = NA, iter = NA, info = NA, vcov = NA)
        return(result)
    }
    rho <- param[1]/sqrt((param[1] + param[2]) * (param[1] + 
        param[3])) * sqrt(param[4] * param[5]/(param[4] + 1)/(param[5] + 
        1))
    logit.rho <- qlogis(rho)
    expt = setNames(as.double(rep(0, 12)), expt.names)
    s_i = setNames(as.double(rep(0, 8)), abp.names[-9])
    info <- matrix(0, ncol = 8, nrow = 8, dimnames = list(abp.names[-9], 
        abp.names[-9]))
    dBvZINB_Expt_vec(xvec = xy.reduced$x, yvec = xy.reduced$y, 
        freq = xy.reduced$freq, n = n.reduced, a0 = param[1], 
        a1 = param[2], a2 = param[3], b1 = param[4], b2 = param[5], 
        p1 = param[6], p2 = param[7], p3 = param[8], p4 = param[9], 
        expt = expt, s_i = s_i, info = info, se = 1, bnb = 1)
    qr.info <- try(qr(info))
    if (class(qr.info) == "try-error") {
        warning("The information matrix has NA/NaN/Inf and thus the standard error is not properly estimatd.")
        std.param = setNames(rep(NA, 11), c(abp.names, "rho", 
            "logit.rho"))
        cov.mat <- NA
    }
    else if (qr(info)$rank < 8) {
        warning("The information matrix is (essentially) not full rank, and thus the standard error is not reliable.")
        std.param = setNames(rep(NA, 11), c(abp.names, "rho", 
            "logit.rho"))
        cov.mat <- NA
    }
    else {
        cov.mat <- try(solve(info))
        if (class(cov.mat) == "try-error") {
            std.param = setNames(rep(NA, 11), c(abp.names, "rho", 
                "logit.rho"))
            cov.mat <- NA
        }
        else {
            var.p4 <- sum(cov.mat[6:8, 6:8])
            d.g <- rho * c(1/param[1] - 1/{
                2 * (param[1] + param[2])
            } - 1/{
                2 * (param[1] + param[3])
            }, -1/{
                2 * (param[1] + param[2])
            }, -1/{
                2 * (param[1] + param[3])
            }, 1/{
                2 * param[4] * (param[4] + 1)
            }, 1/{
                2 * param[5] * (param[5] + 1)
            })
            var.rho <- t(d.g) %*% cov.mat[1:5, 1:5] %*% d.g
            var.logit.rho <- var.rho/rho^2/(1 - rho)^2
            std.param = sqrt(c(diag(cov.mat), p4 = var.p4, rho = var.rho, 
                logit.rho = var.logit.rho))
        }
    }
    result <- list(rho = matrix(c(rho, logit.rho, std.param[c("rho", 
        "logit.rho")]), ncol = 2, dimnames = list(c("rho", "logit.rho"), 
        c("Estimate", "Std.err"))), coefficients = matrix(c(param, 
        std.param[1:9]), ncol = 2, dimnames = list(abp.names, 
        c("Estimate", "Std.err"))), lik = expt[1], iter = iter, 
        info = info, vcov = cov.mat)
    return(result)
}
<bytecode: 0xbd91640>
<environment: namespace:bzinb>
 --- function search by body ---
Function bzinb.se in namespace bzinb has this body.
 ----------- END OF FAILURE REPORT -------------- 
Error in if (class(cov.mat) == "try-error") { : 
  the condition has length > 1
Calls: bzinb.se
Execution halted