select_model: select_model

In Malindrie/scShapes: A Statistical Framework for Modeling and Identifying Differential Distributions in Single-cell RNA-sequencing Data

select_model

Description

This function is used to select the distribution of best fit for scRNA-seq count data

Usage

select_model(lrt.value)

Arguments

lrt.value

A list of genes with the p-values from performing the GOF tests from gof_model

Value

A list of selected model distributions for genes scShapes selects.

Examples

data(scData)

# apply the select_model function to the best fit model from the results of
# the gof_model function

library(BiocParallel)
scData_models <- fit_models(counts=scData$counts, cexpr=scData$covariates,
                            lib.size=scData$lib_size, BPPARAM=bpparam())
scData_bicvals <- model_bic(scData_models)
scData_least.bic <- lbic_model(scData_bicvals, scData$counts)
scData_gof <- gof_model(scData_least.bic, cexpr=scData$covariates, lib.size=scData$lib_size,
BPPARAM=bpparam())

scData_fit <- select_model(scData_gof)

Malindrie/scShapes documentation built on Nov. 21, 2022, 8:58 a.m.