bs_binomial_model.RdEta is a matrix of individual-level parmaters for a B-spline with equivalent knots over all individuals
bs_binomial_model(observation_model_parameters, MegaLMM_state = list())This function should pre-calculate design matrices for each individual (assuming they are all unique)