Modeling the joint effect of intensity and duration of alcohol drinking with bivariate spline models