Gibbs sampling algorithm is introduced for Bayesian inference for multiplicative seasonal autoregressive moving average (SARMA) models.