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

> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: cAIC
> ### Title: Conditional Akaike Information for 'lme4'
> ### Aliases: cAIC
> ### Keywords: regression
> 
> ### ** Examples
> 
> 
> ### Three application examples
> b <- lmer(Reaction ~ Days + (Days | Subject), sleepstudy)
> cAIC(b)
 ----------- FAILURE REPORT -------------- 
 --- failure: the condition has length > 1 ---
 --- srcref --- 
: 
 --- package (from environment) --- 
cAIC4
 --- call from context --- 
calculateGaussianBc(model, sigma.estimated, analytic)
 --- call from argument --- 
if (class(Lambday) == "try-error") {
    Rchol <- try(chol(B))
    L1 <- backsolve(Rchol, C, transpose = TRUE)
    Lambday <- backsolve(Rchol, L1)
}
 --- R stacktrace ---
where 1: calculateGaussianBc(model, sigma.estimated, analytic)
where 2: biasCorrectionGaussian(mer, sigma.estimated, analytic)
where 3: bcMer(object, method = method, B = B, sigma.estimated = sigma.estimated, 
    analytic = analytic)
where 4: cAIC(b)

 --- value of length: 2 type: logical ---
[1] FALSE FALSE
 --- function from context --- 
function (model, sigma.estimated, analytic) 
{
    C <- model$C
    B <- model$B
    e <- model$e
    A <- model$A
    tye <- model$tye
    V0inv <- model$V0inv
    if (analytic) {
        WAlist <- lapply(model$Wlist, function(w) {
            if (!model$isREML) 
                return(w %*% V0inv)
            else return(w %*% A)
        })
        for (j in 1:length(model$theta)) {
            Wj <- model$Wlist[[j]]
            eWje <- model$eWelist[[j]]
            C[j, ] <- as.vector((e %*% Wj) %*% A - eWje * e/(2 * 
                tye))
            for (k in j:length(model$theta)) {
                Wk <- model$Wlist[[k]]
                WkAWjA <- sum(t(WAlist[[j]]) * WAlist[[k]])
                eWke <- model$eWelist[[k]]
                if (!model$isREML) {
                  B[j, k] <- B[k, j] <- -tye * WkAWjA/(2 * model$n) - 
                    eWje * eWke/(2 * tye) + as.numeric(e %*% 
                    Wk %*% (A %*% (Wj %*% e)))
                }
                else {
                  B[j, k] <- B[k, j] <- -tye * WkAWjA/(2 * (model$n - 
                    ncol(model$X))) - eWje * eWke/(2 * tye) + 
                    as.numeric(e %*% Wk %*% (A %*% (Wj %*% e)))
                }
            }
        }
    }
    else {
        if (!model$isREML) {
            np <- model$n
        }
        else {
            np <- model$n - ncol(model$X)
        }
        for (j in 1:length(model$theta)) {
            Wj <- model$Wlist[[j]]
            eWje <- model$eWelist[[j]]
            C[j, ] <- 2 * np/tye * as.vector((e %*% Wj) %*% A - 
                eWje * e/tye)
        }
    }
    Lambday <- try(solve(B) %*% C)
    if (class(Lambday) == "try-error") {
        Rchol <- try(chol(B))
        L1 <- backsolve(Rchol, C, transpose = TRUE)
        Lambday <- backsolve(Rchol, L1)
    }
    df <- model$n - sum(diag(A))
    for (j in 1:length(model$theta)) {
        df <- df + sum(Lambday[j, ] %*% (A %*% (model$Wlist[[j]] %*% 
            e)))
    }
    if (sigma.estimated) {
        df <- df + 1
    }
    return(df)
}
<bytecode: 0xaed79c0>
<environment: namespace:cAIC4>
 --- function search by body ---
Function calculateGaussianBc in namespace cAIC4 has this body.
 ----------- END OF FAILURE REPORT -------------- 
Error in if (class(Lambday) == "try-error") { : 
  the condition has length > 1
Calls: cAIC ... bcMer -> biasCorrectionGaussian -> calculateGaussianBc
Execution halted