Added some helper functions for running the gibbs sampler, including with fixed hyperparameters. Added it to script for comparing inference approaches