#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
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,
--
![Labcor]() |
Ana Luiza Cassin
Monitora de
Pesquisa
Monitor
Research
Tel.com.:
+55(31)3313-1251
Tel.cel.:
+55(31)9308-0072
Labcor
|