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

编程问答 2022-05-16

问题描述

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

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

这是我目前的代码：

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")

scatterplot

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

数据链接^

任何帮助将不胜感激！

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

rrspline