fitting models to clustered correlated data