如何解决RMS包中的Function是否打印出正确的拟合模型?
当我探索和研究rms
包以拟合回归样条时,我注意到了一些奇怪的东西,但不确定是不是错误。
在下面的示例中,我拟合了线性回归,其中一个组件是具有3节的自然三次样条。
我使用Function
包中的rms
打印出拟合的模型。我还使用rcsspline.eval
打印出结。但是,在拟合模型的打印输出中,我看不到预测值 bp 的2级项。另外,括号中的结与rcsspline.eval
中的结不匹配(我假设括号中的术语表示截断的幂基)。
向我指出正确方向的任何帮助将不胜感激。
> library(rms)
>
> # Generate the response and the predictors
> set.seed(123)
> age <- rnorm(500,50,15)
> bp <- rnorm(500,120,10)
> y <- abs(age-50)*.08 + pmax(bp,100)*.09 + rnorm(500)
>
>
> f <- ols( y ~ age + rcs(bp,nk =3))
>
# print out the cofficients
> coef(f)
Intercept age bp bp' bp'' bp'''
2.10418 0.00821 0.07531 0.08501 -0.66785 1.44522
>
> # print out the knots
> rcspline.eval(bp,nk = 3,knots.only = TRUE)
[1] 106 120 132
>
> # print out the fitted model
>
> Function(f)
function(age = NA,bp = NA) {2.1041812+0.0082086026*age+0.075309422* bp+7.6118728e-05*pmax(bp-103.13774,0)^3-0.00059797264*pmax(bp-113.95708,0)^3+0.0012940095*pmax(bp-119.98812,0)^3-0.00097515678*pmax(bp-125.82031,0)^3+0.00020300115*pmax(bp-136.55717,0)^3 }
<environment: 0x0000019cd3ff1c68>
>
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。