Notes:

Exploring different ways of fitting random effects GAMs

Why is stochastic lambda > deterministic?

Optimizing the number of meshpoints

Optimizing number of bootstraps