Added initialization of parameter matrices, still need to fill them with their values after the Gibbs sampling step.