Description Usage Arguments Value Examples
View source: R/means_variances.R
Marginal (equilibrium) expectation of process
1 2 |
lambda |
autonomous rate of cell formation for first cell type per unit of time |
birth |
vector for each type of cell with proportions of cells that replicate per unit of time |
death |
vector for each type of cell with proportions of cells that die per unit of time |
tran |
vector of length equal to length(Y) -1 with proportions of cells that transform to the next cell type per unit of time |
parms |
option list with parameters lambda, birth, death and tran |
verbose |
(default FALSE) print additional output for debugging |
marginal, i.e. long run, expectation of process
1 2 3 4 5 6 7 8 9 10 11 | parms <- list(lambda = 10, birth = rep(.01,3),
death = .03, tran = c(.01,.02))
marg_ev(parms=parms)
# comparing with asymptotic conditional expectation:
Nsteps <- 2000
mat <- matrix(NA, 3, Nsteps)
mat[,1] <- c(100,100,100)
for ( i in 2:Nsteps) mat[,i] <- cond_ev(mat[,i-1],parms = parms)
# Equilibrium profile with these parameters:
head(t(mat),10)
tail(t(mat),10)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.