Olá,
Estou tentando desenvolver uma Regressão Cox no R, mas consigo apenas uma reta com valores praticamente contínuos.
ano<- c(1982,1983,1984,1985,1986,1987,1988,1989,1990,1991,1992,1993,1994,1995,1996,1997,1998,1999,2000,2001,2002,2003,2004,2005,2006,2007,2008,2009,2010)
dado1<-c(117.6423,116.3174,120.5568,140.6938,172.2240,143.6083,140.9587,121.3517,155.5315,145.1981,127.4458,126.6509,133.8048,155.2666,174.8736,116.5824,130.3603,125.0611,124.0013,121.6166,130.8902,157.6512,119.2320,111.2832,144.9331,160.5658,125.3261,166.3949,145.1981)
dado2<-c(237.2314,226.8339,237.7374,246.8556,245.0277,226.8549,240.7147,242.8530,235.3038,243.4697,228.0853,237.0662,234.8285,233.6033,245.6090,237.1481,234.6894,239.9852,237.6996,234.6507,229.7693,239.0660,236.2122,243.6228,233.9454,242.9659,239.3584,242.5270,227.0022)
dado3<-c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1)
dados<-data.frame(cbind(ano,dado1,dado2,dado3))
require(survival)
curva <- coxph(Surv(dado1, dado3) ~ dado2, dados)
a<-summary(curva)
coef<-as.numeric(data.frame(a$coef[1]))
eixo1<-survfit(curva)$surv
eixo2<-survfit(curva)$time
cox<-eixo1^exp(coef*dados[6,3])
plot(eixo2,eixo1,type='S',col="red",ylim=c(0,1.1))
lines(eixo2,cox,col="black",type='l',lwd=2)
Eu esperaria um gráfico com as linhas mais próximas uma da outra, não assim tão distante e não também com a 'curva cox' resultado valores praticamente constantes (pouca variação).
Como corrigir o modelo de regressão Cox que fiz acima?
Agradeço a ajuda!
Surv(dado1, dado3)
, não seriaano
ao invés dedado1
? Não sei o que significadado1
mas a função funçãoSurv()
pede tempo e evento como entradas.coxph(Surv(dado1, dado3) ~ dado2, dados)
? A funçãocoxph
já ajusta o modelo de riscos proporcionais de Cox, então não consegui entender o que você está chamando de "curva Cox".