Dear r-sig-ecology list users,

I have been using the mgcv package for fitting a GAM model to a large set of
spatial/temporal data with several predictors. Each variable contains >
1,000,000 values...
I have been able to successfully fit a model using the bam() function for
large data sets. It takes a while but the results seem to be favorable. My
setup looks something like this:

require(mgcv)
bs="cr";k=10
res <- bam(CHLA ~ s(SSS, bs = bs, k = k) + s(SST, bs = bs, k = k) +
s(SIarea,
   bs = bs, k = k) + s(ZOOC, bs = bs, k = k) + s(DIN, bs = bs,
   k = k) + s(DSI, bs = bs, k = k) + s(DFE, bs = bs, k = k), data=DB)

I was hoping someone might be able to suggest a way to extract the values of
plot.gam (smooth function and standard error). My reason for wanting to do
this is to be able to create a single graphic of several different gam fits
at a later time. The "res" object above becomes very large, and I would like
to avoid saving the entire object.

Many thanks for your help.

Cheers,
Marc

	[[alternative HTML version deleted]]

