Re: [R] Fitting spline using Pspline

From: guy33 <david.reshef_at_magd.ox.ac.uk>
Date: Sun, 29 May 2011 15:28:18 -0700 (PDT)

Ravi,

Thanks so much! You're right, smooth.spline does work on larger n.

Although, for some reason it's results are different (slightly less good?, but I'm not sure). For example, on the simple doppler function below, sm.spline seems to be closer to the true function than smooth.spline:

x=array(0,1000)
y=array(0,1000)
for (i in 1:1000){

x[i] = i/1000
y[i] = (x[i]*(1-x[i]))^.5 * sin(2*pi*(1.05/(x[i]+.05)))
}
plot(x,y)

fit = sm.spline(x, y, norder=2, cv=FALSE) lines(fit\$x,fit\$y)

fit2 = smooth.spline(x, y, cv=FALSE)
lines(fit2\$x,fit2\$y)

What do you make of that?
-guy33

--
View this message in context: http://r.789695.n4.nabble.com/Fitting-spline-using-Pspline-tp3559202p3559610.html
Sent from the R help mailing list archive at Nabble.com.

______________________________________________
R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Received on Sun 29 May 2011 - 22:30:17 GMT

This quarter's messages: by month, or sorted: [ by date ] [ by thread ] [ by subject ] [ by author ]

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Mon 30 May 2011 - 06:10:10 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.