樓主: wubuyun2009
4107 33

[學習分享] 如何利用SAS做RCS(Restricted cubic spline)曲线? [推廣有獎]

  • 1關注
  • 1粉絲

碩士生

53%

還不是VIP/貴賓

-

威望
0
論壇币
1935 个
通用積分
15.8047
學術水平
11 点
熱心指數
12 点
信用等級
10 点
經驗
1928 点
帖子
49
精華
1
在線時間
110 小时
注冊时间
2011-5-19
最后登錄
2019-8-22

wubuyun2009 发表于 2019-6-1 10:54:30 |顯示全部樓層
本帖最后由 wubuyun2009 于 2019-6-2 22:34 编辑

首先,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學習|主题: 255, 订阅: 59
  • · MyLib|主题: 381, 订阅: 36
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曲線
各類資料下載https://www.docin.com/nowlife
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-8-25 08:49