library(rms)
library(rmda)

setwd("")

non_tumor<-read.table("input.txt",header=T,sep="\t")

modul<- decision_curve(MDR~CREA+Fever+	FIB+Invasiveoperation+	
                         LYMPH+	Age	+	DD	+APTT,data= non_tumor, 
family = binomial(link ='logit'),
thresholds= seq(0,1, by = 0.01),
confidence.intervals = 0.95)

pdf("DCA.pdf")
plot_decision_curve(modul,
curve.names="Nonadherence prediction nomogram",xlab="Threshold probability",
cost.benefit.axis =FALSE,col= "blue",
confidence.intervals=FALSE,
standardize = FALSE)
dev.off()