Estimates global beta-binomial distribution parameters on all genes together to evaluate the degree of skewness towards the reference allele. For the estimation, genes located on the sex chromosomes and imprinted genes are excluded

Estimates global beta-binomial distribution parameters on all genes together to evaluate the degree of skewness towards the reference allele. For the estimation, genes located on the sex chromosomes and imprinted genes are excluded

Usage

glob_disp(a1_counts, tot_counts, genes.excl, min_counts = 5)

Arguments

a1_counts

Reference allele counts

tot_counts

Integer matrix (genes x cells): total counts (same dimenstions and rownames as a1_counts).

genes.excl

Character vector of gene IDs to exclude (eg. sex chromosome or imprinted genes)

min_counts

Integer >= 0. Minimum reads per cell to include (default 0). Cells with a number of mapped reads less than min_counts are excluded from the estimation

Examples

glob_disp(a1_counts,
          tot_counts,
          genes.excl,
          min_counts = 5)
#> Error: object 'a1_counts' not found