# NOT RUN {
data(Mtb80)
# simulate 44% of genes to be essential
essential <- rep(0,nrow(Mtb80))
essential[sample(1:nrow(Mtb80),ceiling(nrow(Mtb80)*0.44))] <- 1
# simulate 759 mutants
counts <- sim.mutants(Mtb80[,1], essential, Mtb80[,2], 759)
# run the Gibbs sampler without returning detailed output
# }
# NOT RUN {
output <- negenes(Mtb80[,1], counts[,1], Mtb80[,2], counts[,2])
# }
# NOT RUN {
# run the Gibbs sampler, returning the detailed output
# }
# NOT RUN {
output2 <- negenes(Mtb80[,1], counts[,1], Mtb80[,2], counts[,2], return=TRUE)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab