如何选择RCS曲线最佳节点数？

样条曲线本质是一个分段多项式函数，此函数受限于某些控制点，称为 “节点”，节点放置在数据范围内的多个位置，多项式的类型以及节点的数量和位置决定了样条曲线的类型。

其中，RCS节点的数量比位置更重要。由于节点个数的选择和自由度有关, 所以当样本量比较大的时候可以设置较多的节点。但是节点越多, 自由度越大, 模型越复杂, 越难解释。

那么如何确定RCS曲线的最优节点数呢？

在«Regression Modeling Strategies»这本书中，Harrell建议节点数为4时，模型的拟合效果较好，即同时可以兼顾曲线的平滑程度以及避免过拟合造成的精确度降低。当样本量较大时，5个节点是更好的选择。小样本（n<30）可以选择3个节点。当节点的个数为2时，得到的拟合曲线就是一条直线。因此，大多数研究者推荐的节点为3-5个。但是一个一个试验跑代码，未免有点太耗时耗力。

这里推荐大家使用风暴统计平台进行分析，可以帮助您快速确认最佳节点数！

风暴统计是由浙江中医药大学的郑卫军教授基于R语言开发的，平台将复杂的参数设置转换为了菜单式操作模式，只需在对应选框更改节点数，实时更新曲线变化结果，超快速！还提供AIC与BIC智能确定最佳节点数的功能！

下面我们就结合一份实操数据来为大家详细介绍一下风暴统计平台绘制RCS曲线的具体步骤！

实操具体网址：https://shiny.medsta.cn/trend/

或者百度、必应Bing搜索“风暴统计”

本平台上线的所有工具都是免费的

1.进入风暴统计平台

首先，浏览器搜索风暴统计，依次点击"风暴智能统计"——"临床研究趋势分析"——"趋势性分析与RCS"！进入分析界面后，根据提示，完成数据的导入与整理。

这里我们不再赘述数据的导入与整理过程，详细教程大家可以点击下方链接：

详细指南！风暴统计如何高效导入数据，统计分析快人一步？

详细版！如何利用风暴统计进行数据的整理转换？

我们绘制RCS曲线就可以直接来到"非线性趋势RCS图"模块啦！

风暴统计平台提供了“线性回归”、“logistic回归”、“COX回归”三大常见的回归的RCS绘制！一站式满足不同研究者的统计需求。

2.选择变量

这里我们以logistic回归为例，进行下面的实操！

绘制单因素回归的RCS曲线，只需要选入”因变量“，”连续型自变量“（因变量只能是二分类以0、1赋值哦！）

如果，想要绘制调整协变量后的RCS曲线，只需要在”潜在混杂变量筛选设置“里模块的下拉框中，依次勾选想要调整的协变量！

平台就会实时出现调整协变量后的RCS曲线啦！超快速！这时默认的节点数是4！

最最最重要的节点数自定义设置来啦！！关于RCS曲线的节点数选择，平台给出了3种形式！

• 自定义节点数
  

• 根据AIC
  

• 根据BIC

如果选择自定义节点数，平台支持3-7个节点设置，可以逐一调整，观察曲线变化，不到半分钟，不同节点数对应的图形全部得到！

此外，平台还可以根据AIC值、BIC值确定最佳的节点数，更加智能！

AIC和BIC都是用来衡量模型拟合优良性的指标，一般认为，AIC或BIC越小，则模型的拟合效果越好，根据这一原理，平台会根据不同节点模型的比较，选择最佳节点数制图！

如下图，我们选择根据AIC值来定义节点数，可以发现曲线发生了变化。

另外，在"模型拟合信息"中，可以看到此时定义的最佳节点数是3，这时的AIC值相比4-7个节点数也都是最小的！

3.下载结果

最后在曲线下载设置中，选择下载图片，保存类型，曲线页面，中文字体。

注：如果图形显示不完全，可以在曲线页面设置自定义。自由调整图片的长度和宽度，避免下载图形不完整的情况。

选择完毕后，就可以下载单因素RCS或者调整协变量后的RCS曲线啦！

最后下载的pdf也是十分的清晰，如果需要将变量名修改为中文，也可以在pdf中直接编辑修改，非常便捷！

如果您在风暴统计平台的使用过程中有任何的建议或疑问，欢迎加入我们的讨论群！群里郑老师与助教会在群内解答！

统计机器人交流群