fit models to clustered, correlated data (R package)