It is traditionally assumed that the random effects in mixed models follow a multivariate normal distribution, making likelihood-based inferences more feasible theoretically and computationally.

For functional regression models with functional responses, we propose a nonparametric random-effects model using Gaussian process priors. The proposed model captures the heterogeneity nonlinearly and ...