Dear Friends, The calculation of penalty criterion in smooth.spline is " pen.crit <- sum(wbar * (ybar - fit$ty) * ybar) " , but I think it should be " pen.crit <- sum(wbar * (ybar - fit$ty) * (ybar-fit$ty)) " right ? Why the integral of the fit$ty is not included ? Thanks a lot . Chuanjun