Goals:
(1) Start coding a simple HMMSBM where the Dirichlet hyperparameters are estimated and come from the same distribution whose hyperparameters are also estimated. We will hold \xi fixed to be 1/g for now, and try to do just estimate \gamma_k for each network now.
(2) Make sure things for the single network guy are working out. Maybe keep \xi fixed? See how that works maybe…