I am currently fitting a Poisson model with multiple-membership random effects.
For the continuous covariates, I fit a smoothing spline, s(...).
I would like to plot the smoothing terms in a similar way as plot.gam() does.
However, if I apply the marginal_smooths(...) function, I get the following error: "Error in if (any(weights < 0)) { : missing value where TRUE/FALSE needed"
Does anyone know how to fix this and/or circumvent the error? I can't really wrap my mind around what is going wrong here.
Kind regards,
Hans