Dear all,

I have used predict.gam (mgcv) to create lookup tables for further use in ArcGIS. So far so good. However, I would also like to predict the standard error in the same manner, and that I can't figure out how to do.

Previously I have done the predictions directly on point files in R (type="response"), and retrieved the standard error via se.fit=TRUE. The se.fit=TRUE for type="terms" (used for the lookup tables) does return an SE for each predictor variable, but I do not understand how to go from the separate SE:s to the single value retrieved by se.fit for type="response". Adding the SE-terms + model intercept and using the link function (in the same manner as for the smooth components) is apparently wrong, as it does not result in the same value.

Best regards,
Antonia




Antonia Nyström Sandman
Ph.D./Marine Ecology

+46 (0)8-522 302 52 (direct)
+46 (0)70-633 03 66 (mobile)
+46 (0)8-522 302 40 (switchboard)
antonia.sandman@aquabiota.se<mailto:antonia.sandman@aquabiota.se>
AquaBiota Water Research
Löjtnantsgatan 25
SE-115 50 Stockholm
Sweden
www.aquabiota.se<http://www.aquabiota.se/>


	[[alternative HTML version deleted]]

