%RCS_reg 画限制性立方样图

amei219 发表于 2018-8-8 21:43
就是像这样的

这个图可以从数据中提取，然后采用proc template来绘制

amei219 发表于 2018-7-19 20:58
出现错误？这是怎么回事呢

knots之间用空格隔开，最后只有一个逗号

RCS_Reg的   specif_val设置的含义是啥呢，看说明书也没理解。。。

2286116865 发表于 2018-6-14 19:22
SAS 宏程序 关于画限制性立方样图的，希望对大家有帮助，就不一一给转邮箱了，谢谢大家的热心帮助！

我的币不够，好惨

请问大神们，我的代码运行之后也没报错也没结果显示，请问是啥问题哈？？

libname data "C:\Users\tangminhua\Desktop\1、HP and stroke\data";/*libname语句表示在SAS软件中建立逻辑库data*/
%RCS_Reg( infile=gaofengBPreport5, outfile=out_fin,/*调用限制性立方样条Cox比例风险回归模型宏文件，infile语句选用所需要分析的数据库*/
Dir_data= "C:\Users\tangminhua\Desktop\1、HP and stroke\data\renew strokedata.xlsx",/*Dir_data语句需要包含infile参数中指定的原始SAS数据文件输入目录的名称*/
Main_spline_var=VAR1120, ref_val=146, AVK_msv=0， knots_msv=5 50 95 , /*Main_spline_var语句表示需要分析的自变量，并且必须为连续性变量，ref_val语句表示设置自变量的参考值，在默认情况下，以自变量的中位数作为参考值，AVK_msv语句规定是否自定义节点的选择，AVK_msv=0表示节点是根据自变量的百分位数进行选择，AVK_msv=1表示节点是根据自变量的数值进行选择，并利用knots_msv语句选择节点*/
Oth_spline_var1=age,                     /*Oth_spline_var1语句用来添加需要调整的连续性变量*/
typ_reg=cox,                          /*typ_reg语句选定回归的类型，本研究选取Cox比例风险回归模型*/
dep_var=VAR1107, surv_time_var=VAR1110, /*dep_var和surv_time_var语句分别用来添加因变量和生存时间*/
adjust_var=VAR3 VAR968 VAR1052 VAR1026 VAR1027 activity321 VAR1132 VAR924  VAR32 VAR1068 VAR1069 ,  /*adjust_var语句用来添加需要调整的分类变量*/
Y_ref_line=1,min_Xaxis=0，max_Xaxis=10, /*Y_ref_line=1 表示在Y轴添加一条数值为1的参考线min_Xaxis语句和max_Xaxis语句用来规定x轴的最大值和最小值*/
exp_beta=1, histogram=1,                  /*exp_beta语句用来定义Y轴数值的输出形式，exp_beta=0表示输出的Y轴数值为Ln（HRs），exp_beta=1则表示输出的Y轴数值为HRs*/
no_delete_files=1);

amei219 发表于 2018-8-8 21:43
就是像这样的

请问亲解决这个问题了嘛？

糖豆豆33 发表于 2022-2-27 14:15
请问大神们，我的代码运行之后也没报错也没结果显示，请问是啥问题哈？？

libname data "C:\Users\tangm ...

你好 请问你的问题解决了吗？   我的日志窗口一直在报这个错，要怎麽修改呀？
3618  /*%include "D:\air pollution\RCS_Reg.sas";
3619
3620  %RCS_Reg （infile=oc_6, outfile= rcs_out /* infile=数集(total 为临时数据集)*/
3621
3622  Main_spline_var= P_PM25_lagm, /*主要连续变量*/
      ---------------
      180
ERROR 180-322: 语句无效或未按正确顺序使用。
3623
3624  typ_reg=cox, /*回归类型(log,cox或lin)*/
3625
3626  dep_var=A_death(2),  /*结局分类变量*/
3627
3628  surv_time_var=A_surtime,  /*生存时间*/
3629
3630  adjust_var=A_diagage, /*校正变量*/
3631
3632  avk_msv=0, knots_msv=10 50 90, /*设置节点,按照分数来设置,一般3-5个,本例为3个节点*/
3633
3634  x_ref_line=0,
3635
3636  y_ref_line=1,
3637
3638  no_title=1,
3639
3640  no_label x=1,
3641
3642  no_label y=1,
3643
3644  no_legend=1,
3645
3646  exp_beta=0,
3647
3648  histogram=1,
3649
3650  display_knots=1 /*exp_beta=0显示对数,exp_beta=1显示HR或OR*/);