Extract residuals for each gene from model fit with fitVarPartModel()
Usage
# S4 method for class 'VarParFitList'
residuals(object, ...)
Details
If model is fit with missing data, residuals returns NA for entries that were missing in the original data
Examples
# load library
# library(variancePartition)
library(BiocParallel)
# load simulated data:
# geneExpr: matrix of gene expression values
# info: information/metadata about each sample
data(varPartData)
# Specify variables to consider
# Age is continuous so we model it as a fixed effect
# Individual and Tissue are both categorical, so we model them as random effects
form <- ~ Age + (1 | Individual) + (1 | Tissue)
# Fit model
modelFit <- fitVarPartModel(geneExpr, form, info)
# Extract residuals of model fit
res <- residuals(modelFit)