Additive modeling for location, scale, and shape parameters of the skew normal distribution