fic (version 1.0.0)

all_inds.default: Form indicator matrix describing all submodels of a general linear wide model

Description

Form indicator matrix describing all possible submodels of a general linear wide model, where the submodels are defined by selected covariates.

Usage

# S3 method for default
all_inds(wide, inds0 = NULL, auxpars = NULL, ...)
all_inds(wide, inds0, ...)
# S3 method for lm
all_inds(wide, inds0, ...)
# S3 method for glm
all_inds(wide, inds0, ...)
# S3 method for coxph
all_inds(wide, inds0, ...)

Value

A matrix in the format required by the inds argument of fic(), representing all possible submodels of the wide model.

The number of rows is the number of models, and the number of columns is the number of parameters in the wide model. The \(r,s\) entry of the matrix is a 1 if the \(r\)th submodel includes parameter \(s\), and 0 otherwise.

If a factor is included (excluded) from the submodel, then all corresponding parameters are included (excluded).

See the package vignettes for examples.

Arguments

wide

A fitted model of standard R format, such that terms(wide) returns information about the terms of the model formula. Models outside standard R packages may not support this.

inds0

Narrow model indicators, in format described in fic.

auxpars

Names of parameters in the wide model other than the covariate effects being selected from. By default, for linear and generalised linear models this is c("(Intercept)"), and for Cox regression this is omitted.

...

Other arguments. Currently unused.

Examples

bwt.glm <- glm(low ~ lwtkg + age + smoke, data=birthwt, family="binomial")
all_inds(bwt.glm, inds0=c(1,0,0,0))

# note no intercept term in Cox models, so inds0 has two elements here
library(survival)
wide <- coxph(Surv(years, death==1) ~ sex + thick_centred, data=melanoma)
all_inds(wide, inds0=c(0,0))