testMRCA - Tests for MRCA annotation enrichment in lineage trees
testMRCA performs a permutation test on a set of lineage trees to determine
the significance of an annotation’s association with the MRCA position of the lineage
testMRCA(graphs, field, root = "Germline", exclude = c("Germline", NA), nperm = 200, progress = FALSE)
- list of igraph object containing annotated lineage trees.
- string defining the annotation field to test.
- name of the root (germline) node.
- vector of strings defining
fieldvalues to exclude from the set of potential founder annotations.
- number of permutations to perform.
TRUEshow a progress bar.
An MRCATest object containing the test results and permutation realizations.
# Define example tree set graphs <- ExampleTrees[1-10] # Perform MRCA test on isotypes x <- testMRCA(graphs, "ISOTYPE", nperm=10)
ANNOTATION COUNT EXPECTED PVALUE 1 IgA 16 13.6 0.0 2 IgA,IgG 1 1.2 0.2 3 IgG 31 33.5 1.0