#######先调用spline包
library ( splines )
###########用lm拟合，主要注意部分是bs（age,knots=c(...)）这部分把自变量分成不同部分
fit =lm(wage～bs(age ,knots =c(25 ,40 ,60) ),data=Wage )
############进行预测，预测数据也要分区
pred= predict (fit , newdata =list(age =age . grid),se=T)
#############然后画图
plot(age ,wage ,col =" gray ")