RCS曲线,除了可以反映自变量与因变量间的非线性关系,还可以用作阈值效应分析,即研究自变量达到某个阈值时对因变量的影响,帮助研究发现连续数据的转折点或临界点。
像是这篇Lancet子刊,IF=24.1的文章中,就分别绘制了4个RCS曲线,并在图像中重点标注了拐点。
但是R语言寻找拐点的方法,根据网络上的方法来说,还比较质朴,通过查表,确定OR/HR=1或者β=0时,自变量对应的值。
得出拐点大约在age=44.5附近,此时HR≈1
但是这个方法不够便捷,并且要求使用者对代码有一定基础进行实操!这对于科研压力重的医务人员来说属实不够友好!
因此,这里这里向大家介绍一个2步走绘制RCS曲线,并给出拐点值的工具,无需代码、无需统计基础,1分钟直接出结果!
下面我们就结合一份实操数据来为大家详细介绍一下浙中大郑老师研发的风暴统计平台绘制RCS曲线并计算出拐点的具体步骤!
实操具体网址:https://www.zstats.cn/software2/trend1/
或者百度、必应Bing搜索“风暴统计”
本平台上线的所有工具都是免费的
1.进入风暴统计平台
首先,浏览器搜索风暴统计,依次点击"风暴智能统计"——"临床研究趋势分析"——"趋势性分析与RCS"!进入分析界面后,根据提示,完成数据的导入与整理。
这里我们不再赘述数据的导入与整理过程,详细教程大家可以点击下方链接:
我们绘制RCS曲线就可以直接来到"非线性趋势RCS图"模块啦!
风暴统计平台提供了“线性回归”、“logistic回归”、“COX回归”三大常见的回归的RCS绘制!一站式满足不同研究者的统计需求。
2.选择变量
这里我们以logistic回归为例,进行下面的实操!
第一步:
绘制单因素回归的RCS曲线,只需要选入”因变量“,”连续型自变量“(因变量只能是二分类以0、1赋值哦!)
第二步:
如果想要绘制调整协变量后的RCS曲线,只需要在”潜在混杂变量筛选设置“里模块的下拉框中,依次勾选想要调整的协变量!
平台就会实时出现调整协变量后的RCS曲线啦!超快速!
接着,在"模型拟合信息"这里就可以分别找到单因素、多因素RCS曲线的拐点对应坐标啦!
3.下载结果
最后在曲线下载设置中,选择需要下载单因素、多因素图像,以及图片的保存类型等,选择完毕后,就可以下载啦!
最后下载的pdf也是十分的清晰,如果需要将变量名修改为中文,也可以在pdf中直接编辑修改,非常便捷!
风暴Zstats统计
www.medsta.cn/software
郑老师团队致力于让忙碌的医务人员或科研人员不擅长统计,也能轻松绘制发表级图表!
关于郑老师团队及公众号
全国最大的医学统计公众号平台,专注于医学生、医护工作者学术研究统计支持,我们是你们统计助理
我们开展对临床预测模型、轨迹增长模型、医学免费数据库NHANES、GBD数据库、孟德尔随机化方法1对1R语言指导(一年内不限时间,周末、晚上均统计师一对一指导)。
①指导学习R语言基本技巧
②全程指导课程学习
③课程R语言代码运行bug修复
④支持学员一篇SCI论文的数据分析
详情联系助教小董咨询(微信号aq566665)