gam.control(epsilon = 0.001, bf.epsilon = 0.001, maxit = 10, bf.maxit = 10, trace = F)
gam(Kyphosis ~ poly(Age, 2) + s(Start), data = kyphosis, subset = Number>5, control = gam.control(trace = T))gam(formula, family, control = gam.control(bf.maxit = 15)) gam(formula, family, bf.maxit = 15) # these are equivalent