This article introduces a bimodal model based on the epsilon-skew-normal distribution. This extension generates bimodal distributions similar to those produced by the mixture of normal distributions. We study the basic properties of this new family. We apply maximum likelihood estimators, calculate the information matrix and present a simulation study to assess parameter recovery. Finally, we illustrate the results to three real data sets, suggesting this new distribution as a plausible alternative for modelling bimodal data.
