nomogram临床预测模型

我想通过两个变量预测肿瘤恶性阳性和阴性,我的数据大概是这样的,样本量67个,阳性的有20个,redscore是连续自变量,PNEP是分类变量,PNEP不显著。这是我的R代码

library(survival)

 library(rms)

 ddist <- datadist(a)

 options(datadist='ddist')

 logistic.lrm <- lrm(outcome~Radscore+PNEP,data=a,x=T,y=T,maxit=1000)

 nomogram <- nomogram(logistic.lrm,fun=plogis,fun.at = c(0,0.05,seq(0.1,0.9,by=0.1),0.95,1))

 plot(nomogram)

call <- calibrate(logistic.lrm,method = "boot",B=1000)

 plot(call,xlim = c(0,1.0),ylim = c(0,1.0))