Closed-form marginal likelihood in Gamma-Poisson factorization