In recent years, the skew-normal models introduced by Azzalini (1985)—and their multivariate generalizations from Azzalini and Dalla Valle (1996)—have enjoyed an amazing success, although an important literature has reported that they exhibit, in the vicinity of symmetry, singular Fisher information matrices and stationary points in the profile log-likelihood function for skewness, with the usu...