A data frame with the columns resp (the response), id and
cask (random effects), w1, w2, and w3
(nonparametric confounders), and x1 if beta0 is of length
1 and x1 and x2 if beta0 is of length 2.
The random effects are modelled with "(1|id) + (1|cask:id)".