gistfile1.txt

PHassump <- cox.zph(fitSURV)
PHassump
par(mfrow=c(5,2))
plot(PHassump[1], xlab="Time", ylab="Residual", main="Schoenfeld residuals (scaled) for Gender"); abline(h=0, lty=2, lwd=1.5, col="red")
plot(PHassump[2], xlab="Time", ylab="Residual", main="Schoenfeld residuals (scaled) for NRS"); abline(h=0, lty=2, lwd=1.5, col="red")
plot(PHassump[3], xlab="Time", ylab="Residual", main="Schoenfeld residuals (scaled) for Disease Stage"); abline(h=0, lty=2, lwd=1.5, col="red")
plot(PHassump[4], xlab="Time", ylab="Residual", main="Schoenfeld residuals (scaled) for Age"); abline(h=0, lty=2, lwd=1.5, col="red")
plot(PHassump[5], xlab="Time", ylab="Residual", main="Schoenfeld residuals (scaled) for Tumorsite"); abline(h=0, lty=2, lwd=1.5, col="red")
plot(PHassump[6], xlab="Time", ylab="Residual", main="Schoenfeld residuals (scaled) for Non-Opioid"); abline(h=0, lty=2, lwd=1.5, col="red")
plot(PHassump[7], xlab="Time", ylab="Residual", main="Schoenfeld residuals (scaled) for Opioid"); abline(h=0, lty=2, lwd=1.5, col="red")
plot(PHassump[8], xlab="Time", ylab="Residual", main="Schoenfeld residuals (scaled) for Surgery"); abline(h=0, lty=2, lwd=1.5, col="red")
plot(PHassump[9], xlab="Time", ylab="Residual", main="Schoenfeld residuals (scaled) for RT"); abline(h=0, lty=2, lwd=1.5, col="red")
plot(PHassump[10], xlab="Time", ylab="Residual", main="Schoenfeld residuals (scaled) for Opioid_OME"); abline(h=0, lty=2, lwd=1.5, col="red")
plot(PHassump[11], xlab="Time", ylab="Residual", main="Schoenfeld residuals (scaled) for meds"); abline(h=0, lty=2, lwd=1.5, col="red")