if (FALSE) {
# load sample data together with sub-population list
data(RandomVA1)
# extract InterVA style input data
data <- RandomVA1$data
# extract sub-population information.
# The groups are "HIV Positive", "HIV Negative" and "HIV status unknown".
subpop <- RandomVA1$subpop
# run without subpopulation
fit1<- insilico( data, subpop = NULL,
Nsim = 400, burnin = 200, thin = 10 , seed = 1,
external.sep = TRUE, keepProbbase.level = TRUE)
fit1
}
Run the code above in your browser using DataLab