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