This paper proposes to structure the spectrogram with nonnegative matrix factorization (NMF) in a probabilistic framework, and derives two expectation-maximization algorithms for estimating the parameters of the model, depending on whether the problem in the time- or time-frequency domain.