The inference on shape parameters is a challenging issue for the family of skew normal distributions, especially under multivariate skew normal settings. In this paper, the plausibility regions for the shape parameter of multivariate skew normal family are constructed based on inferential models (IMs) and its generalization, and on generalized IMs when location parameters and scale parameters are known. For illustrations of our results, simulation studies are conducted as well.