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
______________________________________________ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
