[R] Question about GAMs

[pigpigmeow]

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.