Skip to contents

Plot histogram of allelic ratio distribution across all genes with beta function based on alpha and beta global parameters

Usage

plot_glob_params(a1_counts, tot_counts, glob_params, min_counts = 5)

Arguments

a1_counts

Integer matrix (genes x cells): allele-1 counts.

tot_counts

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

glob_params

Named vector/list with elements alpha, beta, mu, theta as returned by glob_disp().

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

Examples

plot_glob_params(a1_counts,
                 tot_counts,
                 glob_params,
                 min_counts = 5)
#> Error: object 'a1_counts' not found