[R] Question about GAMs

pigpigmeow glorykwok at hotmail.com
Thu Oct 13 17:20:10 CEST 2011


I'm confused...
I type ..
predict.gam(ozonea,type=s(ratio,bs="cr"))

  pressure      maxtemp   s(avetemp)     s(ratio)
1   -0.0459102290 -0.185178463  0.263358446 -0.164558673
2   -0.0286464652 -0.194731320  0.199315027  0.727823293
3    0.0478073459 -0.013227033  0.002228896  0.342373202
4   -0.0089164494  0.082301539 -0.037331159 -0.067260889
5    0.0675373617  0.024984396 -0.047067786 -0.357569389
6    0.0823348735  0.101407254 -0.075884852 -1.485036738
7   -0.0977015204  0.177830112 -0.094755158  0.236575309
8   -0.0903027645  0.225594398 -0.113346667  0.435141242
9    0.0206785742  0.187382969 -0.066346157 -0.256133513
10  -0.1371615520  0.101407254 -0.131656887  0.145057584
11  -0.0477674066 -0.181001505  0.260279546  0.180513043
12  -0.0921599421 -0.009050075 -0.020511366  0.281470433
13   0.0681464361 -0.219212934  0.335348247  0.270813178
......

I want to show s(ratio,bs="cr") term, and show the warning message
Warning message:
In predict.gam(ozonea, type = s(ratio, bs = "cr")) :
  Unknown type, reset to terms.

I don't understand what does it mean.

By the way, i use log-link function, 
should I convert log-link function of newozone to fitted value of newozone?
1. log(newozone) - s(ratio,bs="cr") = x
then   exp(x)

2. exp(newozone) - s(ratio,bs="cr") =X
then x

which one is correct?
so confused

--
View this message in context: http://r.789695.n4.nabble.com/Question-about-GAMs-tp3900848p3901842.html
Sent from the R help mailing list archive at Nabble.com.



More information about the R-help mailing list