Abstract
The main objective of this article is to develop an approximate Bayesian procedure to identify the orders of seasonal multivariate autoregressive processes. Using either a matrix normal-Wishart prior density or a non informative prior, which is combined with an approximate conditional likelihood function, the foundation of the proposed technique is to derive the joint posterior mass function of the model orders in a convenient form. Then one may easily evaluate the joint posterior probabilities of all possible values of the model orders and choose the orders at which the joint posterior mass function attains its maximum to be the identified orders. A simulation study, with different prior mass functions, is carried out to determine the effectiveness of the proposed technique in handling the identification problem of the seasonal bivariate autoregressive processes.