Constrain first nonempty theta to zero after burnin
Add Beta parameter for population SAMC efficiency
Add population based SAMC sampler