search for: x_t0

Due to missing data with some of my predictor variables I first do multiple imputation as follows: library(foreign) library(Amelia) library(norm) set.seed(666) M=10 impdat <- NVP[,c("X_t0","X_t","nvp","adstatus","t0rwfa","ageatran","whostage","t0rhfa","vlsupp","t0rwfh","t0rvl","t0rcd4pc","postrantb","resistance","id")] round(apply(apply(imp...

...Markovian series of poisson processes, which generate (or kill) some objects. At t0, say I have initially x = 3 objects. In the next time step, t1, the number of total objects is Poisson-distributed, subject to a function of x at t0, covariates, and parameters. So x_t1 ~ Pois, with E(x_t1) = f(x_t0, covariates, parameters). Let's choose a very simple function for f, say just f = x * par1 * Covariate1. Now let this process be repeated for say 6 times, always with the number of objects obtained in a previous step as input (x) for the next step. The problem is, at all time steps the tota...