Hello. Using the smooth.spline function one can obtain the 1st derivatives by specifiying < deriv=1 > in the predict function (i.e. predict(fit,fit$x,deriv=1) where fit is the object created by smooth.spline). Can one obtain the first derivatives using gam() and predict.gam? If so, how? Bill Shipley [[alternative HTML version deleted]]