approximate_gene_mean: Approximate Gene Mean Distribution

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 means.

Usage

approximate_gene_mean(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 grand means

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 'fitdistrplus' package is required for 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_mean_params <- approximate_gene_mean(data_summaries)

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