estDisp {airpart}R Documentation

Estimate overdispersion parameter of a beta-binomial

Description

Estimate overdispersion parameter of a beta-binomial

Usage

estDisp(sce, genecluster, type = c("plot", "values"))

Arguments

sce

SingleCellExperiment with a1 matrix and counts

genecluster

the gene cluster for which to estimate the over-dispersion parameter. Default is the cluster with the most cells

type

whether to output the over-dispersion estimates as a plot or a value

Value

A ggplot object of the dispersion estimates over the mean, or a data.frame of the mean and dispersion estimates (theta)

Examples


sce <- makeSimulatedData()
sce <- preprocess(sce)
sce <- geneCluster(sce, G = seq_len(4))
estDisp(sce, genecluster = 1)

[Package airpart version 1.2.0 Index]