gammod <- gam(Kyp ~ s(Age,3) + s(Start,3) + s(Number), family=binomial) summary(gammod)