Note that there is NOT an intercept term in my pspline basis. One of the features of psplines is that \sum beta_i f_i(x) (where f_i are the spline basis functions) is linear if and only if the coefficients beta are a linear sequence. This makes it easy to decompose the fit into linear and non-linear effects. Terry T