compare_many {seqCAT} | R Documentation |
Overlap and compare genotypes in many SNV profiles.
compare_many(many, one = NULL, a = 1, b = 5)
many |
SNV profiles to be compared (list of dataframes). |
one |
SNV profile to be compared to all others (dataframe). |
a |
Similarity score parameter a (integer). |
b |
Similarity score parameter b (integer). |
This is a function that compares all the combinations of the SNV profiles input to it, either in a one-to-many or many-to-many manner. It returns both a dataframe containing summary statistics for all unique combinations and a list of dataframes with all the performed comparisons, for easy re-use and downstream analyses of said comparisons.
A list of summary statistics and comparisons.
# Load test data data(test_profile_1) data(test_profile_2) # Perform many-to-many comparisons profiles <- list(test_profile_1, test_profile_2) comparisons <- compare_many(profiles) # View aggregate similarities ## Not run: comparisons[[1]]) # View data of first comparison ## Not run: head(comparisons[[2]][[1]])