dropCoef: Ensure Full Rank Design Matrix

Description Usage Arguments Details Value Author(s) See Also Examples

Description

Coefficients (columns) are dropped from a design matrix to ensure that it has full rank.

Usage

1
drop.coef(X, silent = FALSE)

Arguments

X

a design matrix, e.g., the result of model.matrix possibly of less than full column rank, i.e., with redundant parameters. Works for ncol(X) >= 0 and nrow(X) >= 0.

silent

should a message not be issued if X is column rank deficient?

Details

Redundant columns of the design matrix are identified with the LINPACK implementation of the qr decomposition and removed. The returned design matrix will have qr(X)$rank columns.

Value

The design matrix X without redundant columns.

Author(s)

Rune Haubo B Christensen

See Also

qr and lm

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
X <- model.matrix( ~ PRODID * DAY, data = soup)
ncol(X)
newX <- drop.coef(X)
ncol(newX)

## Essentially this is being computed:
qr.X <- qr(X, tol = 1e-7, LAPACK = FALSE)
newX <- X[, qr.X$pivot[1:qr.X$rank], drop = FALSE]
## is newX of full column rank?
ncol(newX) == qr(newX)$rank
## the number of columns being dropped:
ncol(X) - ncol(newX)

ordinal documentation built on May 2, 2019, 5:47 p.m.