We describe parameter estimation for the multivariate sub-Gaussian symmetric stable distribution using Monte Carlo EM algorithm. Two augmented vectors are employed in the construction of the posterior joint density of the stable parameters. Gibbs sampling enables the generation of these vectors from their respective conditional posterior distributions and thus facilitates the expectation step of the algorithm.