robust(family = gaussian, scale=0, k=1.345, maxit=10)
# Fit a linear model robustly glm(formula, family=robust) # Fit a generalized additive model robustly gam(formula, family=robust(binomial))