nullspace function

Nullspace

Compute the (right or left) nullspace of matrix using a (semi-complete) Singular Value Decomposition.

nullspace(A, type = c("right", "left"), tol = sqrt(.Machine$double.eps))

Arguments

• A: a numeric matrix.
• type: "right" (default) gives is the standard nullspace, "left" gives left nullspace of A.
• tol: tolerance multiple of the first singular value to determine if subsequent singular values are (sufficiently) positive to be determined greater than zero.

Returns

a matrix with as many rows as there are columns in A. The number of columns (which may be zero) determine the dimensionality of the nullspace of A.

Details

This implementation is fastest on matrices with more rows than columns such as a typical design matrix for a linear model.

Examples

# FIXME: We need some examples here

Author(s)

Rune Haubo B. Christensen