Multivariate smoothing using iterative bias reduction with kernel, thin plate splines, Duchon splines or low rank splines.
