I'd be inclined to use predict(ozonea,type="terms") to extract the
estimates of s(ratio,bs = "cr")
that you need. But do you really want newozone - s(ratio,bs="cr") when
you've used a log link?
best,
Simon
On 10/13/2011 09:05 AM, pigpigmeow wrote:
hi! I hope all of you can help me this question
for example GAMs:
ozonea<-gam(newozone~
pressure+maxtemp+s(avetemp,bs="cr")+s(ratio,bs="cr"),family=gaussian
(link=log),groupA,methods=REML)
formula(ozonea)
newozone ~ pressure + maxtemp + s(avetemp, bs = "cr") + s(ratio,bs = "cr")
#formula of gams
coef(ozonea) # extract the coefficient of GAMs
(Intercept) pressure maxtemp s(avetemp).1 s(avetemp).2
s(avetemp).3 s(avetemp).4 s(avetemp).5 s(avetemp).6 s(avetemp).7
s(avetemp).8 s(avetemp).9 s(ratio).1
2.848204561 -0.024662520 -0.095528572 -0.052613005 -0.135309283
0.066243642 0.030937179 0.176267227 0.243405931 0.357027920 0.600628869
0.758581430 0.004653696
s(ratio).2 s(ratio).3 s(ratio).4 s(ratio).5 s(ratio).6
s(ratio).7 s(ratio).8 s(ratio).9
0.156547789 -0.132527837 -0.179947367 -0.208918785 -0.468014803
-0.195498523 -0.813667830 -1.505844690
then, I want to calculate
newozone - s(ratio,bs="cr")
we define this term (newozone - s(ratio,bs="cr")) is X
can i use the coefficient to calculate?
--
View this message in context:
http://r.789695.n4.nabble.com/Question-about-GAMs-tp3900848p3900848.html
Sent from the R help mailing list archive at Nabble.com.
______________________________________________
[email protected] mailing list
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.
______________________________________________
[email protected] mailing list
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.
