if (FALSE) {
data(RandomVA1)
fit1a<- insilico(RandomVA1, subpop = NULL,
Nsim = 1000, burnin = 500, thin = 10 , seed = 1,
auto.length = FALSE)
summary(fit1a, id = "d199")
# The following script updates credible interval for individual
fit1b <- updateIndiv(fit1a, CI = 0.95)
summary(fit1b, id = "d199")
}
Run the code above in your browser using DataLab