使用matlab 2023a计算cox模型，导入csv文件后，这两个分级变量在求取cox单因素模型时，均以第一个变量等级为依据求取相对的P值与HR值 文件中有列"NAFLD", "sex", "age", "CACSgrades", "CADRADS", "Hypertension", "Diabetes", "Dyslipidemia", "smoking", "RCA","PCATgrade", "SIS", "SSS",MACE，time。其中CACSgrades，CADRADS为等级变量，CACSgrades分位1、2、3、4，共4个等级，CADRADS分为0、1、2、3、4、5，共6个等级。time为生存时间，事件结果为MACE（定义为1）。分别求其中每个单因素cox回归模型，建立4个多因素模型，第1组age + Diabetes + Hypertension + CACSgrades + CADRADS + SIS + SSS，第2组age + Diabetes + Hypertension + CACSgrades + CADRADS + SIS + SSS + NAFLD第3组age + Diabetes + Hypertension + CACSgrades + CADRADS + SIS + SSS + PCATgrade第4组age + Diabetes +Hypertension + CACSgrades + CADRADS +SIS + SSS + PCATgrade + NAFLD求取每个多因素模型的C-INDEX，对每组两两之间的C-INDEX进行检验。

您可以按照以下步骤进行操作：

1. 首先，您需要导入csv文件，您可以使用MATLAB中的readmatrix函数读取csv文件，例如：

   data = readmatrix('your_file_path.csv');

2. 对于等级变量，您需要将其转换为虚拟变量，您可以使用MATLAB中的dummyvar函数将其转换为虚拟变量，例如：

   CACSgrades_dummy = dummyvar(data(:,4));
   CADRADS_dummy = dummyvar(data(:,5));

3. 然后，您可以使用MATLAB中的coxphfit函数计算单因素cox回归模型，例如：

   [b_age,logl_age,H,stats_age] = coxphfit(data(:,3),data(:,14),'censoring',~data(:,15));
   [b_diabetes,logl_diabetes,H,stats_diabetes] = coxphfit(data(:,7),data(:,14),'censoring',~data(:,15));
   [b_hypertension,logl_hypertension,H,stats_hypertension] = coxphfit(data(:,6),data(:,14),'censoring',~data(:,15));
   [b_CACSgrades,logl_CACSgrades,H,stats_CACSgrades] = coxphfit(CACSgrades_dummy,data(:,14),'censoring',~data(:,15));
   [b_CADRADS,logl_CADRADS,H,stats_CADRADS] = coxphfit(CADRADS_dummy,data(:,14),'censoring',~data(:,15));
   [b_SIS,logl_SIS,H,stats_SIS] = coxphfit(data(:,12),data(:,14),'censoring',~data(:,15));
   [b_SSS,logl_SSS,H,stats_SSS] = coxphfit(data(:,13),data(:,14),'censoring',~data(:,15));

4. 接着，您可以使用MATLAB中的coxphfit函数计算多因素cox回归模型，例如：

   X1 = [data(:,3),data(:,7),data(:,6),CACSgrades_dummy,CADRADS_dummy,data(:,12),data(:,13)];
   [b1,logl1,H,stats1] = coxphfit(X1,data(:,14),'censoring',~data(:,15));
   
   X2 = [data(:,3),data(:,7),data(:,6),CACSgrades_dummy,CADRADS_dummy,data(:,12),data(:,13),data(:,1)];
   [b2,logl2,H,stats2] = coxphfit(X2,data(:,14),'censoring',~data(:,15));
   
   X3 = [data(:,3),data(:,7),data(:,6),CACSgrades_dummy,CADRADS_dummy,data(:,12),data(:,13),data(:,11)];
   [b3,logl3,H,stats3] = coxphfit(X3,data(:,14),'censoring',~data(:,15));
   
   X4 = [data(:,3),data(:,7),data(:,6),CACSgrades_dummy,CADRADS_dummy,data(:,12),data(:,13),data(:,11),data(:,1)];
   [b4,logl4,H,stats4] = coxphfit(X4,data(:,14),'censoring',~data(:,15));

5. 最后，您可以使用MATLAB中的coxphint函数计算C-INDEX，例如：

   [cindex1,~,~,~,~,~,~,~,~,~] = coxphint(b1,X1,data(:,14),'censoring',~data(:,15));
   [cindex2,~,~,~,~,~,~,~,~,~] = coxphint(b2,X2,data(:,14),'censoring',~data(:,15));
   [cindex3,~,~,~,~,~,~,~,~,~] = coxphint(b3,X3,data(:,14),'censoring',~data(:,15));
   [cindex4,~,~,~,~,~,~,~,~,~] = coxphint(b4,X4,data(:,14),'censoring',~data(:,15));

6. 最后，您可以使用MATLAB中的anova函数对两两之间的C-INDEX进行检验，例如：

   p12 = anova(cindex1,cindex2);
   p13 = anova(cindex1,cindex3);
   p14 = anova(cindex1,cindex4);
   p23 = anova(cindex2,cindex3);
   p24 = anova(cindex2,cindex4);
   p34 = anova(cindex3,cindex4);