Displaying 1 result from an estimated 1 matches for "m1tp".
Did you mean:
1tp
2012 Feb 13
3
mgcv: increasing basis dimension
...sis dimension k=10
m2<-gam(y~x1+x2+te(d,bs="ts",k=10)+te(x3,bs="ts",k=10)+te(d,x3,bs="ts",k=10),family="poisson")
#gcv increased
m1$gcv
m2$gcv
summary(m1)
summary(m2)
gam.check(m1)
gam.check(m2)
#is this due to bs="ts"?
#basis dimension k=5
m1tp<-gam(y~x1+x2+te(d,bs="tp")+te(x3,bs="tp")+te(d,x3,bs="tp"),family="poisson")
#basis dimension k=10
m2tp<-gam(y~x1+x2+te(d,bs="tp",k=10)+te(x3,bs="tp",k=10)+te(d,x3,bs="tp",k=10),family="poisson")
m1tp$gcv
m2tp$g...