Dear List, I want to use the predict() function from the raster-package to apply a GAM to my gridded predictors. Unfortunately I get an error message:
When using a linear model everything works fine: > mod.lm <- lm(y ~ Slope + Boen_By50m + Hoehe_abs, data=dat_nwt) > p <- predict(object = r, model = mod.lm, filename='gebiet_1/sturmwurf.img', progress='text', na.rm = FALSE, overwrite = TRUE, format = 'HFA', type='response') |======================================================================| 100% The same, when using a GLM > mod.glm <- glm(y ~ Slope + Boen_By50m + Hoehe_abs, data=dat_nwt, > family=binomial()) > p <- predict(object = r, model = mod.glm, filename='gebiet_1/sturmwurf.img', progress='text', na.rm = FALSE, overwrite = TRUE, format = 'HFA', type='response') |======================================================================| 100% But, when using a GAM I get an Error. > mod.gam <- gam(y ~ Slope + Boen_By50m + Hoehe_abs, data=dat_nwt, > family=binomial()) > p <- predict(object = r, model = mod.gam, filename='gebiet_1/sturmwurf.img', progress='text', na.rm = FALSE, overwrite = TRUE, format = 'HFA', type='response') | | 0% Error in object[-omit, , drop = FALSE] : incorrect number of dimensions Does anybody know what 'incorrect number of dimensions' means in this context? Thanks TIM *** > sessionInfo() R version 2.13.0 (2011-04-13) Platform: i386-pc-mingw32/i386 (32-bit) locale: [1] LC_COLLATE=German_Germany.1252 LC_CTYPE=German_Germany.1252 [3] LC_MONETARY=German_Germany.1252 LC_NUMERIC=C [5] LC_TIME=German_Germany.1252 attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] rgdal_0.6-33 raster_1.8-12 sp_0.9-81 mgcv_1.7-6 loaded via a namespace (and not attached): [1] grid_2.13.0 lattice_0.19-23 Matrix_0.999375-50 nlme_3.1-100 ---------------- Tim Häring Bavarian State Institute of Forest Research Department of Soil and Climate Hans-Carl-von-Carlowitz-Platz 1 D-85354 Freising E-Mail: [email protected] http://www.lwf.bayern.de _______________________________________________ R-sig-Geo mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-sig-geo
