v1.6.1

• add likelihood = ‘beta’ (for beta regression)
• faster computations for crossed grouped random effects by enabling matrix_inversion_method = “iterative” ’
• add Vecchia-inducing points full-scale (VIF) approximation (gp_approx = "full_scale_vecchia")
• add metric 'crps_gaussian'
• add option to not estimate / hold fix some covariance parameter via the argument ‘estimate_cov_par_index’ in the ‘params’ parameter (of the ‘fit’ or ‘set_optim_params’ functions)
• add ‘space_time_gneiting’ covariance function
• add objective = ‘mean_scale_regression’ (for plain boosting without a GPModel)
• make iterative methods as default wherever possible