fit_nzi_model: Fit a (non-zero-inflated) model for a single gene

In switchde: Switch-like differential expression across single-cell trajectories

Description

Fits a sigmoidal expression model for a single gene vector, returning MLE model parameters and p-value.

Usage

fit_nzi_model(y, pst)

Arguments

y	Vector of gene expression values
pst	Pseudotime vector, of same length as y

Value

A vector with 5 entries: maximum likelihood estimates for μ_0, k t0, σ^2 and a p-value

Examples

data(synth_gex)
data(ex_pseudotime)
y <- synth_gex[1, ]
fit <- fit_nzi_model(y, ex_pseudotime)

switchde documentation built on Nov. 8, 2020, 5:19 p.m.