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…