#Oi segue abaixo um exemplo, espero ter ajudado!
#Alexandro
#modelo
tempo=c(40,30,20,10,5,45,35,25,15,5)
censura=c(1,1,1,1,1,1,1,0,1,1)
sexo=c("m","f","m","f","m","f","m","m","m","f")
doente=factor(c(1,1,0,0,0,0,1,1,1,1))
idade=c(20,30,40,50,60,25,35,45,55,65)
dados=data.frame(tempo,censura,sexo,doente,idade);dados
require(survival)
.Survfit <- survfit(Surv(tempo, censura) ~ doente + sexo, conf.type="log",
conf.int=0.95, type="kaplan-meier", error="greenwood", data=dados)
summary(.Survfit)
plot(.Survfit, col=1:4, lty=1:4, mark.time=TRUE)
legend("bottomleft", legend=c("0.f","1.f","0.m","1.m"), title="doente + sexo",
col=1:4, lty=1:4, bty="n")#INTERACAO
>________________________________
> De: Ana Luiza Cassin <[email protected]>
>Para: [email protected]
>Enviadas: Segunda-feira, 6 de Maio de 2013 15:35
>Assunto: [R-br] Survival
>
>
>
>Prezados,
>
>
>Estou começando a usar or programa R e preciso fazer curva de
sobrevivência (Kaplan-Meier) sobre o acompanhamento de pacientes ao
longo do tempo.
>Gentileza, se alguém tiver conhecimento ou algum tutorial sobre essa
análise, o mais detalhado possível, agradeceria.
>
>Atenciosamente,
>
>--
>
> Ana Luiza Cassin
>Monitora de Pesquisa
>Monitor Research
>Tel.com.: +55(31)3313-1251
>Tel.cel.: +55(31)9308-0072
>
>Labcor
>www.labcor.com.br
>_______________________________________________
>R-br mailing list
>[email protected]
>https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
>mínimo reproduzível.
>
>_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
