He intentado replicar los resultados de la entrada de ayer con GAM (vía mgcv) así (véase el enlace anterior para la definición de los datos):
library(mgcv)
modelo_gam <- gam(
y ~ x + s(id, bs = "re"),
data = datos,
method = "REML",
family = "poisson")
Y nada:
Error in gam(y ~ x + s(id, bs = "re"), data = datos, method = "REML", :
Model has more coefficients than data
Sí, ya sé que tengo más variables que observaciones. Pero, ¿no es para eso que estoy usando efectos aleatorios?
En fin…