nomogram

javascript
library(rms)
library(foreign)
library(survival)
tcga<-read.csv("9.csv")
tcga$age<-factor(tcga$age,labels=c("<50","50-59","60-69",">=70"))
tcga$sex<-factor(tcga$sex,labels=c("FEMALE","MALE"))
tcga$visiondecreased<-factor(tcga$visiondecreased,labels=c("NO","YES"))
tcga$movement<-factor(tcga$movement,labels=c("NO","YES"))
tcga$PNI<-factor(tcga$PNI,labels=c("NO","YES"))
tcga$invasiontobone<-factor(tcga$invasiontobone,labels=c("NO","YES"))
tcga$T_group<-factor(tcga$T_group,labels=c("T1","T2","T3+T4"))
tcga$M<-factor(tcga$M,labels=c("M0","M1"))
tcga$N<-factor(tcga$N,labels=c("N0","N1"))
tcga$N<-factor(tcga$freemagrin,labels=c("negative","positive"))
tcga$histopathologicaltype = factor(tcga$histopathologicaltype,

                                    levels = c("NONBASOLID","BASOLID"))
head(tcga$histopathologicaltype)
tcga$metastasistime = as.numeric(tcga$metastasistime)
class(tcga$metastasistime)
metastasis <- coxph(Surv(metastasistime,metastasis) ~age + sex + visiondecreased + movement + PNI + invasiontobone + T_group + N + M + histopathologicaltype,

                    data=tcga)
install.packages("regplot")
library(regplot)
regplot(metastasis,failtime = c(12,36,60),prfail = T,observation = tcga[2,])
library(survminer)
ggforest(metastasis)