8.0 DCA&CIC

install.packages("rmda")

library(rmda)

#######################训练集

t.mod <- decision_curve(MACE~Age +Killip +Hypertention+Hhcy +HF +Sent+VD3LM
                           +LDLC+LVEF_F,
                       data = Vdata2 
                       )

# 画图
plot_decision_curve(t.mod, curve.names = "t.mod",
                    cost.benefit.axis = F, # 是否需要损失：获益比 轴
                    confidence.intervals = "none" # 不画可信区间
)
                    
summary(t.mod,measure= "NB")  

summary(t.mod,measure= "sNB")   
                    
#######################验证集
  
  v.mod <- decision_curve(MACE~Age +Killip +Hypertention+Hhcy +HF +Sent+VD3LM
                          +LDLC+LVEF_F,
                          data = Vdata1 
  )  
  

  # 画图
  plot_decision_curve(v.mod, curve.names = "v.mod",
                      cost.benefit.axis = F, # 是否需要损失：获益比 轴
                      confidence.intervals = "none" # 不画可信区间
  )
  
  
  summary(v.mod,measure= "NB")  
  summary(v.mod,measure= "sNB") 
  
  
###############################################################################33                  
                    
                    
######临床影响曲线
  
  # 1次只能绘制1个模型
 
#内部训练
  
  plot_clinical_impact(t.mod,population.size = 1000,cost.benefit.axis = F,
                       n.cost.benefits= 8,col = c('red','blue'),
                       confidence.intervals= T)
  
#外部验证
  plot_clinical_impact(v.mod,population.size = 1000,cost.benefit.axis = F,
                       n.cost.benefits= 8,col = c('red','blue'),
                       confidence.intervals= T)