Thiago Gentil Ramires - Departamento de Estatística (UEM
Thiago Gentil Ramires - Departamento de Estatística (UEM
Thiago Gentil Ramires - Departamento de Estatística (UEM
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
ekm = survfit(Surv(tempo,censura)~cor,type="kaplan-meier")<br />
amarela=1<br />
branca=1<br />
negra=1<br />
sort(tempo)<br />
mu = 0.9951<br />
gama= 1.0050<br />
beta0= 5.9623<br />
beta1= -1.1069<br />
beta2= -1.3950<br />
beta3= -1.0080<br />
S1 = exp(-(tempo/(exp(beta0+beta1*amarela)))^gama)<br />
S2 = exp(-(tempo/(exp(beta0+beta2*branca)))^gama)<br />
S3 = exp(-(tempo/(exp(beta0+beta3*negra)))^gama)<br />
S4 = exp(-(tempo/(exp(beta0)))^gama)<br />
plot(ekm,lty=c(1,2,3,4),xlab="t",ylab="S(t)")<br />
lines(tempo,S1,type="l",lty=2,col="blue")<br />
lines(tempo,S2,type="l",lty=2,col="red")<br />
lines(tempo,S3,type="l",lty=2,col="black")<br />
lines(tempo,S4,type="l",lty=2,col="green")<br />
legend(250,0.8, lty=c(1,2,3,4),c("Amarela", "Branca","Negra","parda"), bty="n",<br />
cex=0.8)<br />
text(150.5,0.93,c("Ajuste da variavel Cor"),bty="n",cex=1.2)<br />
#Analise <strong>de</strong> residuos<br />
# ajuste geral<br />
Cor<br />
dialisee