mean_cond_var: Mean conditional variance

In gmonette/Tcells2: Fits Tcell dynamic models with mixed models and STAN

Description

Mean conditional variance

Usage

mean_cond_var(lambda = parms$lambda, birth = parms$birth,
  death = parms$death, tran = parms$tran, parms, verbose = FALSE)

Arguments

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

Value

E ( Var( Y_{t+1} | Y_t)) = T E(D) T' where D is the diagonal matrix of variances of new cell births, deaths and transitions and T is the coding matrix of 0s and \pm 1s to linearly transform the birth, death and transition components into the cell counts

Examples

parms <- list(lambda = 10, birth = rep(.01,3),
              death = .03, tran = c(.01,.02))
mean_cond_var(parms = parms)

gmonette/Tcells2 documentation built on May 17, 2019, 7:25 a.m.