install.packages("rmda")

library(rmda)

#######################訓(xùn)練集

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" # 不畫可信區(qū)間
)
                    
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" # 不畫可信區(qū)間
  )
  
  
  summary(v.mod,measure= "NB")  
  summary(v.mod,measure= "sNB") 
  
  
###############################################################################33                  
                    
                    
######臨床影響曲線
  
  # 1次只能繪制1個模型
 
#內(nèi)部訓(xùn)練
  
  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)