Кто-нибудь работал с моделью регрессии Кокса в R?
Я пытаюсь использовать ее, но получаю только одну прямую с практически непрерывными значениями.
year <- 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)
d1<-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)
d2<-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)
d3<-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)
dd<-data.frame(cbind(year,d1,d2,d3))
require(survival)
curve <- coxph(Surv(d1, d3) ~ d2, dd)
a<-summary(curve)
coef<-as.numeric(data.frame(a$coef[1]))
ex1<-survfit(curve)$surv
ex2<-survfit(curve)$time
cox<-ex1^exp(coef*dd[6,3])
plot(ex2,ex1,type='S',col="red",ylim=c(0,1.1))
lines(ex2,cox,col="black",type='l',lwd=2)
Я бы хотеложидайте, что график с линиями, расположенными ближе друг к другу, не такими уж далекими и не имеющими результата «кривая Кокса», практически постоянные значения (небольшое отклонение).
Спасибо за помощь!