楼主: wubuyun2009
5573 39

[学习分享] 如何利用SAS做RCS(Restricted cubic spline)曲线? [分享]

  • 1关注
  • 3粉丝

硕士生

55%

还不是VIP/贵宾

-

威望
0
论坛币
2016 个
通用积分
16.8047
学术水平
11 点
热心指数
12 点
信用等级
10 点
经验
1928 点
帖子
49
精华
1
在线时间
113 小时
注册时间
2011-5-19
最后登录
2019-9-28

wubuyun2009 发表于 2019-6-1 10:54:30 |显示全部楼层
首先,RCS主要是反应连续变量与结局的关系,目前有很多SAS宏可以帮助做,个人认为最好的还是Desquilbet在2010编写%RCS_Reg宏,另外SAS9.4高级版本也可以绘制。
下面简要讲解下利用%RCS_Reg宏绘制RCS曲线:
1、下载宏程序(见下面链接);
2、将宏程序放在电脑硬盘的一个文件夹内;
3、打开sas,有相应的线性回归、Logistic回归、Cox回归的数据集;
4、输入%include "E:\SASMacro\RCS_Reg.sas"   导入宏程序,具体地址根据自己的喜好来选择;
5、编写宏语句,比如下面的cox回归
%RCS_Reg(infile=total, Main_spline_var=MAP,typ_reg=cox, dep_var=death, surv_time_var=deathtime,
        adjust_var=age sex charlson HDL ALB infection group DVC dialysis_age DVC*dialysis_age  smoking_history logHB logLDL logTG  Hypertension,
         avk_msv=0, knots_msv=5 25 50 75 95,x_ref_line=0,y_ref_line=1,no_title=1,no_label_x=1,no_label_y=1,no_legend=1,
         exp_beta=0,histogram=1,display_knots=1    /*exp_beta=0显示对数,exp_beta=1显示HR或OR*/
          );

说明:1)infile=数据集(total 为临时数据集);
          2)Main_spline_var=主要连续变量;typ_reg=回归类型(log,cox或lin),dep_var=结局分类变量, surv_time_var=生存时间;
          3)adjust_var=校正变量,若校正变量为连续变量,也可以放在Oth_spline_var1=需校正的连续变量1,Oth_spline_var2=需校正的连续变量2中。
          4)avk_msv=0, knots_msv=5 25 50 75 95,设置节点,按照分数来设置,本例为5个节点。
          5)exp_beta=0显示对数(如lnHR或lnOR),exp_beta=1显示HR或OR。
          6)其他的设置就是显示标题,图例,X轴标题,Y轴标题、参考线的设置了。

6、结果图形见附件gplot4

7、RCS一般5个节点足以发现大多数的非线性关系,当然一般需要绘制RCS图形的同志也需要专门学习RCS相关的知识了,有一个ppt和宏程序作者Stat Med文章的链接,专门介绍了相关知识。
8、如果证实不是线性,一般需要根据临床意义(如血压可采用低血压、正常血压、高血压来分类)来分组,或者采用四分位数、五分位数、十分位数(样本很大)来观察不同区间内关注的连续变量与结局的关系。

绘制图形

绘制图形

%RCS_reg_立方样图.zip

28.87 KB

售价: 2 个论坛币

RCS_Reg宏和变量解释

RestrictedCubicSplineforLinearityTest&Continuous:限制三次样条的线性测试和连续.pdf

607.57 KB

售价: 2 个论坛币

SAS中RCS介绍ppt

desquilbet2010.pdf

451.02 KB

售价: 1 个论坛币

引用文献,也有宏的具体介绍

已有 5 人评分经验 论坛币 学术水平 热心指数 收起 理由
np84 + 100 精彩帖子
crystal8832 + 5 精彩帖子
aclyang + 20 精彩帖子
zhdefei + 2 鼓励积极发帖讨论
eijuhz + 1 精彩帖子

总评分: 经验 + 105  论坛币 + 20  学术水平 + 1  热心指数 + 2   查看全部评分

本帖被以下文库推荐

  • · SAS学习|主题: 258, 订阅: 59
  • · MyLib|主题: 388, 订阅: 38
stata SPSS
nkunku 发表于 2019-6-2 09:53:45 |显示全部楼层
利用SAS做RCS(Restricted cubic spline)曲线
回复

使用道具 举报

tianwk 发表于 2019-6-4 00:20:05 |显示全部楼层
thanks for sharing
回复

使用道具 举报

vxiangriv 发表于 2019-6-4 08:26:07 |显示全部楼层
谢谢分享
回复

使用道具 举报

hyq2003 发表于 2019-6-4 08:26:09 |显示全部楼层
谢谢分享
回复

使用道具 举报

xujingjun 发表于 2019-6-4 09:12:49 |显示全部楼层
回复

使用道具 举报

gx666666gx888 发表于 2019-6-4 14:08:08 |显示全部楼层
谢谢楼主。
回复

使用道具 举报

wangyong8935 发表于 2019-6-4 14:49:20 |显示全部楼层
回复

使用道具 举报

panjiquan 发表于 2019-6-5 17:44:57 |显示全部楼层
Reg宏绘制RCS曲线
回复

使用道具 举报

lunaque 发表于 2019-6-6 00:45:25 来自手机 |显示全部楼层
wubuyun2009 发表于 2019-6-1 10:54
首先,RCS主要是反应连续变量与结局的关系,目前有很多SAS宏可以帮助做,个人认为最好的还是Desquilbet在20 ...
请问一下小伙伴们,你们有Desquilbet文章里面appendixB中的sas代码吗?很需要这个,谢谢大家啦
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 我要注册

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2019-12-12 17:45