in method update_concentration, DPGMM
i think it should be:
self.gamma[i, 2] = self.gamma_[i + 1, 2] + sz[i+1]
instead of:
self.gamma_[i, 2] = self.gamma_[i + 1, 2] + sz[i]