simulate_new | R Documentation |
See vignette("sim", package = "glmmTMB")
for more details and examples,
and vignette("covstruct", package = "glmmTMB")
for more information on the parameterization of different covariance structures.
simulate_new(
object,
nsim = 1,
seed = NULL,
newdata,
newparams,
...,
show_pars = FALSE
)
object |
a one-sided model formula (e.g. |
nsim |
number of simulations |
seed |
random-number seed |
newdata |
a data frame containing all variables listed in the formula, including the response variable (which needs to fall within the domain of the conditional distribution, and should probably not be all zeros, but whose value is otherwise irrelevant) |
newparams |
a list of parameters containing sub-vectors
( |
... |
other arguments to |
show_pars |
(logical) print structure of parameter vector and stop without simulating? |
## use Salamanders data for structure/covariates
simulate_new(~ mined + (1|site),
zi = ~ mined,
newdata = Salamanders, show_pars = TRUE)
sim_count <- simulate_new(~ mined + (1|site),
newdata = Salamanders,
zi = ~ mined,
family = nbinom2,
newparams = list(beta = c(2, 1),
betazi = c(-0.5, 0.5), ## logit-linear model for zi
betad = log(2), ## log(NB dispersion)
theta = log(1)) ## log(among-site SD)
)
head(sim_count[[1]])
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.