fit_zi_model {switchde} | R Documentation |
Fits a zero-inflated sigmoidal model for a single gene vector, returning MLE model parameters and p-value.
fit_zi_model(y, pst, maxiter = 10000, log_lik_tol = 0.001, verbose = FALSE)
y |
Vector of gene expression values |
pst |
Pseudotime vector, of same length as y |
maxiter |
Maximum number of iterations for EM algorithm if zero inflation enabled. Default 100 |
log_lik_tol |
If the change in the log-likelihood falls below this for zero inflated EM the algorithm is assumed to have converged |
verbose |
Print convergence update for EM algorithm |
A vector with 6 entries: maximum likelihood estimates for μ_0, k t0, λ, σ^2 and a p-value
data(synth_gex) data(ex_pseudotime) y <- synth_gex[1, ] fit <- fit_zi_model(y, ex_pseudotime)