ConjComp: Conjugate complement of span(X) in span(Z) with respect to...

In gmonette/spida15: Collection of miscellaneous functions for mixed models etc. prepared for SPIDA 2009+ (development version)

Usage

ConjComp(X, Z = diag(nrow(X)), ip = diag(nrow(X)), tol = 1e-07)

Arguments

X
Z
ip
tol

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (X, Z = diag(nrow(X)), ip = diag(nrow(X)), tol = 1e-07) 
{
    help <- "\n    ConjComp returns a basis for the conjugate complement of the\n    conjugate projection of X into span(Z) with respect to inner product with\n    matrix ip.\n    Note: Z is assumed to be of full column rank but not necessarily X.\n    "
    xq <- qr(t(Z) %*% ip %*% X, tol = tol)
    if (xq$rank == 0) 
        return(Z)
    a <- qr.Q(xq, complete = T)[, -(1:xq$rank)]
    Z %*% a
  }

gmonette/spida15 documentation built on May 17, 2019, 7:26 a.m.