在 R 中的散点图数据上生成平滑样条的问题

如何解决在 R 中的散点图数据上生成平滑样条的问题

试图在一些骨密度数据上生成平滑样条，但遇到了一些麻烦。基本上是从我的统计学习教科书中复制了这个例子，并试图修改一些代码以适应，但并没有真正理解它到底发生了什么或者为什么它不起作用。

为同一图中的男性和女性数据生成两条拟合线将是理想的选择。

这是我目前的代码：

bone <- read.csv('/Users/myname/Desktop/ML Extra Credit/spnbmd.csv')
library(splines)

plot(spnbmd ~ age,data=bone,col = 
       ifelse(sex=="mal","blue","red"),xlab="Age",ylab="Relative Change in Spinal BMD")
legend("topleft",c("male","female"),fill=c("blue","red")
)


fit=lm(spnbmd~bs(age,knots=c(10,15,20)),data=bone)
pred=predict(fit,newdata=list(age=age.grid),se=T)
plot(age,wage,col="gray")
lines(age.grid,pred$fit,lwd=2)
lines(age.grid,pred$fit+2*pred$se,lty="dashed")
lines(age.grid,pred$fit-2*pred$se,lty="dashed")

https://github.com/khammerberg53/MLEC/blob/main/spnbmd.csv

数据链接^

任何帮助将不胜感激！