[R] plot.gam help

[Pollyanna Fisher]

 Hello, I've created a list of gams and for my last step I'd like to plot them but I keep getting an error method. Here is an example of one of themI appreciate your help

Gam1s=gam(E1~chlag +chla + waveMax + SST + IRR + site + day +chlag:site, family=gaussian)
summary(Gam1s)

Family: gaussian 
Link function: identity 

Formula:
E1 ~ chlag + chla + waveMax + SST + IRR + site + day + chlag:site

Parametric coefficients:
                Estimate Std. Error t value Pr(>|t|)    
(Intercept)    1.059e+02  6.153e-01 172.154  < 2e-16 ***
chlag         -3.112e+01  4.356e+00  -7.144 1.25e-12 ***
chla          -1.380e+01  1.650e+00  -8.359  < 2e-16 ***
waveMax        6.718e-03  8.450e-04   7.950 3.02e-15 ***
SST            2.627e-02  9.210e-03   2.853  0.00438 ** 
IRR           -4.702e-02  6.665e-03  -7.055 2.33e-12 ***
siteJOH       -6.356e+00  7.508e-01  -8.465  < 2e-16 ***
siteKUR       -4.291e+00  4.497e-01  -9.540  < 2e-16 ***
day            1.208e-03  4.348e-04   2.777  0.00554 ** 
chlag:siteJOH  3.582e+01  7.436e+00   4.817 1.56e-06 ***
chlag:siteKUR  3.083e+01  4.602e+00   6.699 2.68e-11 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1


R-sq.(adj) =  0.294   Deviance explained = 29.8%
GCV = 3.9675  Scale est. = 3.9469    n = 2126

plot(Gam1s)
Error in plot.gam(Gam1s) : 
  No terms to plot - nothing for plot.gam() to do.
> plot.gam(Gam1s)
Error in plot.gam(Gam1s) : 
  No terms to plot - nothing for plot.gam() to do.
> plot.gam(E1~chlag +chla + waveMax + SST + IRR + site + day +chlag:site, family=gaussian)
Error in if (se && x$Vp[1, 1] < 0) { : 
  missing value where TRUE/FALSE needed
	[[alternative HTML version deleted]]