Except for boundary effects, the function will reproduce polynomials up to degree 3.
spence.21(y)
spence.15
,
spencer
,data(spencer)
yy <- spence.21(spencer$mortality)
plot(spencer$age, spencer$mortality)
lines(spencer$age, yy)
Run the code above in your browser using DataLab