in method update_concentration, DPGMM
i think it should be:
self.gamma[i, 2] = self.gamma_[i + 1, 2] + sz[i+1]
self.gamma_[i, 2] = self.gamma_[i + 1, 2] + sz[i]
I agree. Maybe @alextp can confirm. Otherwise sz[self.n_components] is not included in the sum.
Maybe we should write
self.gamma_[:-1, 2] = np.cumsum(sz[::-1])[::-1][1:]
but that is also not super readable.
Closing: the new Dirichlet process GMM re-write has been merged in master.