我正在使用limma分析R中的甲基化数据集。
我为数据拟合了以下线性模型：
#methylation dataset is a 400000x20 table (400000 CGs and 20 samples)
X <- model.matrix(~variable1+variable2+variable3)
fit <- lmFit(mValues,X)
resids <- residuals(fit,mValues)
#resids returns as expected a 400000x20 table
我需要一种将样本中的残差转换为标准化残差的方法。我知道rstandard，但它仅适用于lm（）。 MarrayLM对象有类似的东西吗？
我是否应该对残渣的每一行都使用standard（）函数？