Your problem is that you've specified a parametric linear model. Your
output is a set of parametric coefficients. You can replace `gam` with
`lm` in your example, and you'll get identical results. If you wish to
specify one of the terms nonparametrically, you'd use s(chlag), for
example (assuming you're using mgcv).
On 11/13/2015 10:59 PM, Pollyanna Fisher via R-help wrote:
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.
______________________________________________
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.
