Includes random intercpt, option to include random slope, second random intercept, and time-correlation.
aldex.mem.sim(
D,
days,
subjects,
depth = 10000,
location = 1,
random_slope = FALSE,
corr = 0,
rho_ar1 = 0,
sd_resid = 0.1
)number of taxa/genes
num days (i.e., repeated measurements) for each subject
num of subjects to simulate
sum of counts for each multinomial draw
second random intercept, if 0 ignore, else simulate num of locations
If true, simulate a random slope for each subject.
The correlation between slope/random intercept
The ar1 time-correlation, if 0, don't simulate
The residual error time-correlation.
Kyle McGovern