fit_nzi_model {switchde}R Documentation

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

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)

[Package switchde version 1.18.0 Index]