Energy intake

This function runs a standard Gibbs sampling algorithm on a posterior distribution associated with a mixture model and 500 datapoints.

1 | ```
gibbsmix(Niter = 10^4, v = 1)
``` |

`Niter` |
Number of MCMC iterations |

`v` |
Scale of the normal prior |

The function returns a plot of the log-posterior surface, along with the MCMC sample represented
both by points and lines linking one value to the next. Evaluating the log-posterior surface on a
*250x250* grid takes some time.

Christian P. Robert and George Casella

Chapter 6 of **EnteR Monte Carlo Statistical Methods**

1 | ```
## Not run: gibbsmix(Nit=10^3)
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.