Running examples in ‘forward-Ex.R’ failed
The error most likely occurred in:

> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: fwdglm
> ### Title: Forward Search in Generalized Linear Models
> ### Aliases: fwdglm print.fwdglm
> ### Keywords: robust regression models
> 
> ### ** Examples
>  
> data(cellular)
> cellular$TNF <- as.factor(cellular$TNF)
> cellular$IFN <- as.factor(cellular$IFN)
> mod <- fwdglm(y ~ TNF + IFN, data=cellular, family=poisson(log), nsamp=200)
lmsglm: found 4 good subsets after trying 10.
lmsglm: found 8 good subsets after trying 20.
lmsglm: found 13 good subsets after trying 30.
lmsglm: found 16 good subsets after trying 40.
lmsglm: found 22 good subsets after trying 50.
lmsglm: found 27 good subsets after trying 60.
lmsglm: found 30 good subsets after trying 70.
lmsglm: found 37 good subsets after trying 80.
lmsglm: found 41 good subsets after trying 90.
lmsglm: found 44 good subsets after trying 100.
lmsglm: found 47 good subsets after trying 110.
lmsglm: found 49 good subsets after trying 120.
lmsglm: found 50 good subsets after trying 130.
lmsglm: found 51 good subsets after trying 140.
lmsglm: found 55 good subsets after trying 150.
lmsglm: found 58 good subsets after trying 160.
lmsglm: found 62 good subsets after trying 170.
lmsglm: found 66 good subsets after trying 180.
lmsglm: found 71 good subsets after trying 190.
lmsglm: found 74 good subsets after trying 200.
lmsglm: found 80 good subsets after trying 210.
lmsglm: found 86 good subsets after trying 220.
lmsglm: found 90 good subsets after trying 230.
lmsglm: found 93 good subsets after trying 240.
lmsglm: found 96 good subsets after trying 250.
lmsglm: found 101 good subsets after trying 260.
lmsglm: found 105 good subsets after trying 270.
lmsglm: found 107 good subsets after trying 280.
lmsglm: found 112 good subsets after trying 290.
lmsglm: found 116 good subsets after trying 300.
lmsglm: found 118 good subsets after trying 310.
lmsglm: found 122 good subsets after trying 320.
lmsglm: found 127 good subsets after trying 330.
lmsglm: found 129 good subsets after trying 340.
lmsglm: found 132 good subsets after trying 350.
lmsglm: found 137 good subsets after trying 360.
lmsglm: found 138 good subsets after trying 370.
lmsglm: found 140 good subsets after trying 380.
lmsglm: found 144 good subsets after trying 390.
lmsglm: found 145 good subsets after trying 400.

	*** Starting Forward Search ***

 ----------- FAILURE REPORT -------------- 
 --- failure: the condition has length > 1 ---
 --- srcref --- 
: 
 --- package (from environment) --- 
forward
 --- call from context --- 
scglm(xbsb, ybsb, family = family, weights = zz$weights, offset = obsb, 
    beta = zz$coefficients, phi = dispersion)
 --- call from argument --- 
if (!(class(Cb) != "try-error")) return("singular")
 --- R stacktrace ---
where 1: scglm(xbsb, ybsb, family = family, weights = zz$weights, offset = obsb, 
    beta = zz$coefficients, phi = dispersion)
where 2: fwdglm(y ~ TNF + IFN, data = cellular, family = poisson(log), 
    nsamp = 200)

 --- value of length: 2 type: logical ---
[1] FALSE FALSE
 --- function from context --- 
function (x, y, family, weights, beta, phi = 1, offset) 
{
    n <- dim(x)[1]
    p <- dim(x)[2]
    if (missing(weights)) 
        weights <- rep(1, n)
    weights <- sqrt(weights)
    if (missing(offset)) 
        offset <- 0
    if (missing(beta)) 
        stop("coefficients estimates must be provided")
    user.error <- options()$show.error.messages
    on.exit(options(show.error.messages = user.error))
    options(show.error.messages = FALSE)
    inverse.fn <- family$linkinv
    variance.fn <- family$variance
    dm.fn <- function(eta) 1/family$mu.eta(eta)
    eta <- offset + x %*% matrix(beta, ncol = 1)
    zz <- eta + (y - inverse.fn(eta)) * dm.fn(eta) - offset
    Xb <- weights * cbind(x, eta^2)
    CXb <- t(Xb) %*% Xb
    Cb <- try(solve(CXb, tol = 1e-20))
    if (!(class(Cb) != "try-error")) 
        return("singular")
    b.add <- Cb %*% t(Xb) %*% (zz * weights)
    se <- sqrt(phi * Cb[p + 1, p + 1])
    b.add[p + 1]/se
}
<bytecode: 0x5572978194c0>
<environment: namespace:forward>
 --- function search by body ---
Function scglm in namespace forward has this body.
 ----------- END OF FAILURE REPORT -------------- 
Fatal error: the condition has length > 1