Thiago Gentil Ramires - Departamento de Estatística (UEM

ekm = survfit(Surv(tempo,censura)~sexo,type="kaplan-meier")
masculino=0
feminino=1
sort(tempo)
mu = 0.9897
gama= 1.0104
beta0= 4.892
beta1= -0.4462
S1 = exp(-(tempo/(exp(beta0+beta1*masculino)))^gama)
S2 = exp(-(tempo/(exp(beta0+beta1*feminino)))^gama)
plot(ekm,lty=c(1,4),xlab="t",ylab="S(t)")
lines(tempo,S1,type="l",lty=2,col="blue")
lines(tempo,S2,type="l",lty=2,col="red")
legend(250,0.8, lty=c(1,2),c("Masculino", "Feminino"), bty="n", cex=0.8)
text(150.5,0.93,c("Ajuste da variavel Sexo"),bty="n",cex=1.2)
AntiHBS
ekm = survfit(Surv(tempo,censura)~anthbs,type="kaplan-meier")
false=1
true=0
sort(tempo)
mu = 0.8951
gama= 1.1172
beta0= 5.1288
beta1= -0.9557
S1 = exp(-(tempo/(exp(beta0+beta1*false)))^gama)
S2 = exp(-(tempo/(exp(beta0+beta1*true)))^gama)
plot(ekm,lty=c(1,4),xlab="t",ylab="S(t)")
lines(tempo,S1,type="l",lty=2,col="blue")
lines(tempo,S2,type="l",lty=2,col="red")
legend(250,0.8, lty=c(1,2),c("False", "True"), bty="n", cex=0.8)
text(150.5,0.93,c("Ajuste da variavel AntiHBS"),bty="n",cex=1.2)
Pressão
ekm = survfit(Surv(tempo,censura)~pres,type="kaplan-meier")
false=1
true=0
sort(tempo)
mu = 0.9781
gama= 1.0224
beta0= 4.2038
beta1= 0.6054
S1 = exp(-(tempo/(exp(beta0+beta1*false)))^gama)
S2 = exp(-(tempo/(exp(beta0+beta1*true)))^gama)
plot(ekm,lty=c(1,4),xlab="t",ylab="S(t)")
lines(tempo,S1,type="l",lty=2,col="blue")
lines(tempo,S2,type="l",lty=2,col="red")
legend(250,0.8, lty=c(1,2),c("False", "True"), bty="n", cex=0.8)
text(150.5,0.93,c("Ajuste da variavel Pressão"),bty="n",cex=1.2)
64