“…An approach typically used for fitting multimodal data is the mixture of normal or asymmetric-normal models, which present difficulties such as identifiability problems (see McLachlan and Peel [3]; Marin et al [4]) and complicated numerical implementation. Bimodal distributions generated from skew-symmetric distributions can be found in Azzalini and Capitanio [5], Ma and Genton [6], Arellano-Valle et al [7], Kim [8], Lin et al [9,10], Elal-Olivero et al [11], Arnold et al [12,13], Gómez et al [14], Braga et al [15], Venegas et al [16] and Gómez-Déniz et al [17], among others.…”