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]] ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.