CheckStability {mlica} | R Documentation |
Performs a correlation test to see which of the inferred ICA modes are reproducible across multiple runs using different random initialisations. Returns a set of consensus ICA modes and stability scores for each following the algorithm of Chiappetta,...et.al (2004).
CheckStability(a.best.l, corr.th)
a.best.l |
List of a.best objects from mlica runs. |
corr.th |
Correlation threshold to use to decide whether a mode is reproducible. |
A list with the following components
consS |
Consensus source matrix with columns labeling the consensus ICA modes. Has same number of rows as a.best$S . |
consA |
Consensus mixing matrix with rows labeling the consensus ICA modes. |
stabM |
Vector of same length as consM giving the
stability measures of each consensus ICA mode. Stability or
reproducibility measures are given as fractions, that is, the number
of times the ICA mode correlates with one of the other runs at
threshold level corr.th divided by the number of runs (length of a.best.l ). |
Andrew Teschendorff aet21@cam.ac.uk