我想使用 Cox 比例风险回归来调查一些连续风险因素与我的偏好疾病之间的假定 U 形关联。到目前为止，我已经完成了这个R和survival这样的包
library(survival)
coxph(Surv(entry, exit, status) ~ pspline(my.risk.factor, df = df))
然而，虽然不幸的是文档并不太冗长，但在我看来，这些样条曲线不受限制（即，被迫在尾部是线性的），因此在所述尾部中表现出很大的差异，而且我必须这样做也让我很恼火提供所需的拟合复杂度（通过df参数）。我宁愿把它放在mgcv包里，它可以自动调节拟合的复杂性，但 AFAIK 不支持 cox 模型。
我可以用 来做吗R，如果可以，怎么做？