My model is of the form:
mod-gam(y~s(x0)+s(x1)+s(x2),family=poisson).
But I want to get
estimate and standard error of the difference of two fitted values.
The following code shows you how to do this (a) for differences on the
scale of the linear predictor, and (b) for differences on the
I fitted a GAM model with Poisson distribution
using the function gam() in the mgcv package.
My model is of the form:
mod-gam(y~s(x0)+s(x1)+s(x2),family=poisson).
To extract estimates at a specified set of covariate
values I used the gam `predict' method.
But I want to get
estimate and