[R] Survival Plots by Strata
Bret Collier
bret at tamu.edu
Wed Mar 8 14:10:56 CET 2006
All,
I am struggling to create a survival plot using LTRC data for each year
of a 10 year period.
I have a set of individuals (birds) where 'entry' is the day of the
year (1-365) they are released (let out of pens) into the wild (2 year
data snip below). 'Entry' (e.g., day of year the first bird is
released for each year) is highly variable, ranging from 48 to >250.
When I create a plot using the below code I would like to remove the
'solid' line which represent S_hat=1 out to the LC point for each year
and instead have the survival curves formatted in more of a 'hanging'
style, where the LC day (e.g., min(Entry for each year)) is where each
curve begins at S_hat=1 for each year (and does not extend back to the
y-axis)? I could not find anything on this in the archives or MASS or
Survival Analysis using S? Anyone have a suggestion on where to look?
TIA, Bret
#Code snip for R email.
apc.coxfit1<-coxph(Surv(Entry, Exit, Fate)~Sex + Agerelease +
Dayrelease + strata(Year), data=mydat)
coxfit.apc<-survfit(apc.coxfit1)
coxfit.apc
plot(survfit(apc.coxfit1), conf.int=F, log=T, lty=c(1:2), col=c(1:2),
xlim=c(205, 800)) #not run--first entry for this example is day 205 for
1996, 259 for 1997
>mydat
ID Year Dayrelease
Agerelease Survivorship Entry Exit Fate Sex
16240 1996 205 95 164 205 369 1 0
16319 1996 205 88 140 205 345 1 0
16378 1996 248 108 100 248 348 1 0
20383 1996 241 98 204 241 445 1 0
16324 1996 219 90 227 219 446 1 0
16327 1996 219 90 497 219 716 1 0
20373 1996 241 114 413 241 654 1 0
20374 1996 241 111 211 241 452 1 0
16241 1996 205 95 234 205 439 1 1
16321 1996 219 90 118 219 337 1 1
16323 1996 219 90 180 219 399 1 1
20375 1996 241 103 268 241 509 1 1
20384 1996 241 98 299 241 540 1 1
20390 1996 241 93 204 241 445 1 1
20393 1996 241 88 208 241 449 1 1
16313 1996 248 122 512 248 760 0 1
20378 1996 241 103 236 241 477 0 1
20381 1996 241 101 329 241 570 0 1
16328 1996 219 90 224 219 443 0 1
16827 1997 259 127 52 259 311 1 0
16828 1997 259 127 216 259 475 1 0
16831 1997 303 171 19 303 322 1 0
20466 1997 289 149 31 289 320 1 0
20469 1997 289 149 199 289 488 1 0
20483 1997 289 134 18 289 307 1 0
16807 1997 259 137 223 259 482 1 0
16809 1997 259 137 1 259 260 1 0
16819 1997 259 131 237 259 496 1 0
16829 1997 303 171 7 303 310 1 0
20440 1997 289 161 7 289 296 1 0
20470 1997 289 148 257 289 546 1 0
20478 1997 289 143 12 289 301 1 0
16817 1997 259 130 85 259 344 1 0
20454 1997 289 154 4 289 293 1 1
20459 1997 289 153 335 289 624 1 1
20460 1997 289 153 118 289 407 1 1
20465 1997 289 150 31 289 320 1 1
20473 1997 289 147 65 289 354 1 1
20484 1997 289 133 58 289 347 1 1
20489 1997 289 130 56 289 345 1 1
16808 1997 259 137 137 259 396 0 1
16810 1997 303 181 64 303 367 0 1
16816 1997 259 130 1 259 260 0 1
20471 1997 289 147 334 289 623 0 1
16826 1997 259 127 20 259 279 0 1
More information about the R-help
mailing list