rule: which rule to use in approx() of the response before calculating the global mean, rule=1 means no extrapolation, rule=2 means to extrapolate the closest non-missing value, see approx
silent: print error messages of model fit?
cyclic: defaults to FALSE, if TRUE cyclic splines are used