powered by
An internal function
GP.csh(Z_mat, fixed_effects, control)
data frame
formula specifying fixed effects to be included in model
a list