A nonparametric regression method for multiple longitudinal phenotypes using multivariate adaptive splines