VariantSelection_Group — VariantSelection

We get variants of interest by selecting variants with a high VAF difference between two groups. This function is adapted from the Peter van Galen. Source: https://github.com/petervangalen/MAESTER-2021

Usage

VariantSelection_Group(
  SE,
  min_coverage = 2,
  quantiles = c(0.1, 0.9),
  thresholds = c(0.1, 0.9),
  min_quality = NULL,
  mean_allele_frequency = 0,
  group_of_interest = NULL,
  group1 = NULL,
  group2 = NULL,
  group_factor = 5,
  remove_nocall = TRUE,
  verbose = TRUE
)

Arguments

SE: SummarizedExperiment object.
min_coverage: Minimum coverage needed.
quantiles: The lower and upper quantile you want to use.
thresholds: The VAF thresholds you want to use for the quantiles.
min_quality: The minimum quality you want for the Variants of Interest. Can be ignored by setting it to NULL.
mean_allele_frequency: The minimum mean allele frequency. Default = 0
group_of_interest: The column data that divides the cells.
group1: The first group of interest. If set, the quantiles are only calculated for this group.
group2: The second group of interest.
group_factor: How much higher has the mean allele frequency to be in group 1 when compared to group 2?
remove_nocall: Should NoCall cells (consensus = 0) be disregarded during the analysis?
verbose: Should the function be verbose? Default = TRUE