gam.fit(x = X, y = Y, smooth.frame = mf, weights = weights, start = start,
etastart = etastart, mustart = mustart, offset = offset,
family = family, control = control)

gam(Life ~ s(Doc) + s(TV), data = life, control = list(maxit = 100))