get.resids: Estimate Residuals
In Sage-Bionetworks/snm: Supervised Normalization of Microarrays

Description Usage Arguments Author(s) Examples

Get Residuals from data given a model

1	get.resids(Y, mod)

`Y`	Matrix of data
`mod`	Model matrix parameterizing coefficients

Brig Mecham <brig.mecham@sagebase.org> and John D. Storey <jstorey@princeton.edu>

##---- 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 (data.matrix, design.matrix) 
{
    x.qr <- qr(design.matrix)
    x.det <- prod(diag(x.qr$qr))
    x.cov.unscaled <- solve(t(qr.R(x.qr)) %*% qr.R(x.qr))
    x.df <- nrow(design.matrix) - x.qr$rank
    Fit.coef <- qr.coef(x.qr, t(data.matrix))
    Fit.res <- qr.resid(x.qr, t(data.matrix))
  }