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