基于rms包的限制性立方样条回归（RCS）R代码实现

1.原理

��       限制性立方样条(Restricted cubic spline，RCS)是分析非线性关系的最常见的方法之一。RCS用三次函数拟合不同节点之间的曲线并使其平滑连接，从而达到拟合整个曲线并检验其线性的过程。可以想见，RCS的节点数对拟合结果来说非常重要。通常，小于30个样本数的小样本取3个节点，大样本取5个节点。

2.R实现

1.cox回归

#Used for RCS(Restricted Cubic Spline)#我们使用rms包library(ggplot2)library(rms)library(survminer)library(survival)

在这里我们使用survival包中的lung数据

#####基于cox回归#这里用survival包里的lung数据集来做范例分析head(lung)# inst time status age sex ph.ecog ph.karno pat.karno meal.cal wt.loss#    3  306      2  74   1       1       90       100     1175      NA#    3  455      2  68   1       0       90        90     1225      15#    3 1010      1  56   1       0       90        90       NA      15#    5  210      2  57   1       1       90        60     1150      11#    1  883      2  60   1       0      100        90       NA       0#   12 1022      1  74   1       1       50        80      513       0#status:censoring status 1=censored, 2=dead#sex:Male=1 Female=2#ph.ecog：医生对患者的体能状态评级#ph.karno：医生对患者的另一种体能状态评级karnofsky#pat.karno:患者karnofsky自评#meal.cal:摄入卡路里#wt.loss:过去半年体重减轻# 对数据进行打包，整理dd