# Example with galaxy data ================================================
set.seed(123)
# retrieve galaxy data
y <- galaxy
# estimation
bayesmix <- bayes_fit(
data = y,
K = 5, # not many to run the example rapidly
dist = "normal",
nb_iter = 500, # not many to run the example rapidly
burnin = 100
)
# trace plot
# bayes_trace(bayesmix)
Run the code above in your browser using DataLab