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]]
_______________________________________________
R-sig-ecology mailing list
R-sig-ecology@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-ecology
