## Not run:# simulating data from mixture of gamma with two componentsn =500# number of observationsweight = c(0.6,0.4)alpha = c(12,1)beta = c(3,2)data <- rmixgamma( n = n, weight = weight, alpha = alpha, beta = beta )# plot for simulation data hist( data, prob =TRUE, nclass =50, col ="gray")x = seq(0,10,0.05)truth = dmixgamma( x, weight, alpha, beta )lines( x, truth, lwd =2)# Runing bdmcmc algorithm for the above simulation data set bmixgamma.obj <- bmixgamma( data )plot( bmixgamma.obj )## End(Not run)