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