R/expand.models.R

# expand.models

# all additive combinations

expand.models <- function (model, ...) {
    # model is a list of models
    # LHS, RHS defined in utility.R
    oneformula <- function(model1) {
        coll <- "+"  
        if(any(grepl('\\*', as.character(model1))))
            warning ("ignoring interactions in ", LHS(model1), RHS(model1))
        lhs <- LHS(model1)
        vars <- all.vars(RHS(model1))
        tmp <- lapply(0:length(vars), combn, x = vars)
        tmp[[which(sapply(tmp,nrow)==0)]] <- matrix('1')   ## constant
        rhs <- unlist(lapply(tmp, function(x) apply(x, 2, paste, collapse = coll)))
        lapply(paste(lhs, rhs, sep='~'), as.formula, ...)
    }
    tmp <- lapply(model, oneformula)
    names(tmp) <- sapply(model, LHS)
    apply(do.call(expand.grid, tmp), 1, as.list)
}

## expand.models(list(p~t+b+h2, phi~t))
## expand.models(list(p~1, phi~1))

## could construct input for par.openCR.fit

Try the openCR package in your browser

Any scripts or data that you put into this service are public.

openCR documentation built on Aug. 14, 2021, 9:08 a.m.