More documentation.  Added inference method that performs all inference steps - Gibbs sampling, calculating theta and phi