Use a projection of the given variance-covariance matrix.
proj_matrix(VarRubin.matrix, metadata)
VarRubin.matrix | A variance-covariance matrix. |
---|---|
metadata | Metadata of the experiment. |
A list of variance-covariance matrices.
M. Chion, Ch. Carapito and F. Bertrand (2021). Accounting for multiple imputation-induced variability for differential analysis in mass spectrometry-based label-free quantitative proteomics. arxiv:2108.07086. https://arxiv.org/abs/2108.07086.
library(mi4p) data(datasim) datasim_imp <- multi.impute(data = datasim[,-1], conditions = attr(datasim,"metadata")$Condition, method = "MLE") VarRubin.matrix <- rubin2.all(datasim_imp[1:5,,], attr(datasim,"metadata")$Condition) proj_matrix(VarRubin.matrix, attr(datasim,"metadata"))#> [1] 0.9556735 0.6325980 1.2207511 0.7328072 0.6709007