glmmtmb

问题 I'm trying to predict to a raster stack in R, using a GLMM I fit with the glmmTMB package, at the population level (i.e. setting random effects to 0). I followed Ben Bolker's solution in this thread which works for models fit with lme4 , but even though the function argument re.form=~0 appears to be applicable to predict.glmmTMB in addition to predict.merMod , it's not working for me when I predict using a glmmTMB model. Here is an example using the same example code provided by Robert

问题 As a reproducible example, let's use the next no-sense example: > library(glmmTMB) > summary(glmmTMB(am ~ disp + hp + (1|carb), data = mtcars)) Family: gaussian ( identity ) Formula: am ~ disp + hp + (1 | carb) Data: mtcars AIC BIC logLik deviance df.resid 34.1 41.5 -12.1 24.1 27 Random effects: Conditional model: Groups Name Variance Std.Dev. carb (Intercept) 2.011e-11 4.485e-06 Residual 1.244e-01 3.528e-01 Number of obs: 32, groups: carb, 6 Dispersion estimate for gaussian family (sigma^2):