ポアソン混合モデルの gibbs sampling を書いたので供養のため以下に載せておきます。 for i in range(MAXITER): for n in range(N): # (4.37), (4.38) for k in range(K): etas[n][k] = exp(Xs[n] * ln(lambdas[k]) - lambdas[k] + ln(pis[k])) # normaliza…