if (FALSE) {
## define model code
code <- nimbleCode({
lambda ~ dgamma(0.001, 0.001)
for(i in 1:N) {
AC[i, 1, 1] ~ dunif(0, 100)
AC[i, 2, 1] ~ dunif(0, 100)
for(t in 2:T) {
AC[i, 1:2, t+1] ~ dDispersal_exp(s = AC[i, 1:2, t], rate = lambda)
}
}
})
constants <- list(N = 10, T = 6)
## create NIMBLE model object
Rmodel <- nimbleModel(code, constants)
## use model object for MCMC, etc.
}
Run the code above in your browser using DataLab