approximate_gene_drop: Approximate Gene Dropout as a gamma

In kdzimm/hierarchicell: Simulating Cell-Type Specific and Hierarchical Single-Cell RNA-Seq Data

Description

This function estimates the parameters for the gamma distribution of gene dropout.

Usage

approximate_gene_drop(data_summaries, plot = FALSE)

Arguments

data_summaries	an R object that has been output by the package's compute_data_summaries function.
plot	a TRUE/FALSE statement for the output of a plot to observe a histogram of dropout

Details

Prior to estimating the simulation parameters, it is important to run the compute_data_summaries function to build an object that is in the right format for the following estimation functions to properly compute. The 'ggplot2' package is required for the plotting component of this function to work. Please install it.

Value

A vector of length two, where the first number is the shape of the gamma distribution and the second number is the rate of the gamma distribution

Examples

clean_expr_data <- filter_counts()
data_summaries <- compute_data_summaries(clean_expr_data)
gene_drop_params <- approximate_gene_drop(data_summaries)

kdzimm/hierarchicell documentation built on Dec. 21, 2021, 5:23 a.m.