lin function

Generate orthogonal polynomial base for a numeric covariate without intercept