residualize | R Documentation |
Compute a regression model for each column in a matrix and return residual matrix
residualize(form, X, design)
form |
the formula defining the model to fit for residuals |
X |
the response matrix |
design |
the |
X <- matrix(rnorm(20*10), 20, 10) des <- data.frame(a=rep(letters[1:4], 5), b=factor(rep(1:5, each=4))) xresid <- residualize(~ a+b, X, design=des) ## design is saturated, residuals should be zero xresid2 <- residualize(~ a*b, X, design=des) sum(xresid2) == 0
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.