用SAS做Logistic回归和COx回归限制性立方样图

如何使用这个宏，先运行宏代码，然后运行一下命令：感谢@wangyiyisheng 网友提供宏代码和数据

COX回归

ods rtf file="C:\Users\ALIENWARE\Desktop\Figure\export_curves.rtf";/*结果输出路径，输出到word*/

            ods rtf select where=(_NAME_='export_curves');/*输出的文件名为export_curves*/

%RCS_Reg( infile= Initial_datafile, outfile= out_fin, Dir_data= D:\RCS\Data,

Main_spline_var= HDL_CHOL, ref_val= 1.0, AVK_msv= 1, knots_msv= 0.8 1.0 1.5 2.0,

typ_reg= cox,dep_var= CVD, surv_time_var= SURVIVAL_T,

adjust_var= SEX BLACK_NH HISPANIC ETHN_OTHER CURR_SMOKER INF_HS COLLEGE,

Y_ref_line= 1, max_Xaxis= 5,

exp_beta= 1, round= 0.01);

ods rtf close;

Logistic回归：仍以原数据为例，因变量仍未CVD。

ods rtf file="C:\Users\ALIENWARE\Desktop\Figure\export_curves.rtf";/*j结果输出路径，输出到word*/

ods rtf select where=(_NAME_='export_curves');/*输出的文件名为export_curves*/

%RCS_Reg( infile= Initial_datafile, outfile= out_fin,Dir_data= D:\RCS\Data,

Main_spline_var= HDL_CHOL, ref_val= 0.21, AVK_msv= 1, knots_msv= 0.78 1.03 1.27 1.53 2.04,

typ_reg= log,dep_var= CVD,

adjust_var= SEX BLACK_NH HISPANIC ETHN_OTHER CURR_SMOKER INF_HS COLLEGE,

Y_ref_line= 1, max_Xaxis= 5,

exp_beta= 1, round= 0.01);

ods rtf close;

说明

infile=数据文件名，格式为.sas7bdat，

nfile= Initial_datafile, outfile= out_fin, Dir_data= D:\RCS\Data,表示在D盘文件夹下的文件夹下的名字为Initial_datafile的数据集

ref_val= 1.0, 一般变量的中位数，也可以设置其他

AVK_msv= 1，默认

knots_msv= 0.8 1.0 1.5 2.0,设置节点，数字可改，建议以百分位数

typ_reg= cox,表示做cox回归

dep_var= CVD, 因变量

surv_time_var= SURVIVAL_T,生存时间

adjust_var= SEX BLACK_NH HISPANIC ETHN_OTHER CURR_SMOKER

INF_HS COLLEGE,调整的变量

Y_ref_line= 1, 参考线默认

max_Xaxis= 5,横坐标刻度，可以设置为条样变量的最大值的整数

exp_beta= 1, round= 0.01，默认

结果输出到word里面，打开word里面有图，再word里把图缩小一下，把整个word用转换成pdf，然后用PS打开pdf，单独选择图这一页，稍微剪切修改之后就可以另存为tif格式了，300dpi。

COX回归

Logistic回归