Displaying 3 results from an estimated 3 matches for "agerelease".
2006 Feb 13
1
JRG Console Output
...[13] 77.554024 48.579639 48.111326
[16] 44.994786 65.241722 65.852035
[19] 53.254482 43.217719 30.255150
Reading in some data, the 'mydat' line extends across the console, but
the output is wrapped around again?
> mydat<-data.frame(ID, Year = factor(Year), Dayrelease, Agerelease,
Survivorship, Entry, Exit, Fate, Gender)
> mydat
ID Year Dayrelease Agerelease
1 16240 1996 205 95
2 16319 1996 205 88
3 16378 1996 248 108
.
.
.
576 3094 2005 251 136
577 2667 2005 264 861
578 3035 2005...
2006 Mar 08
2
Survival Plots by Strata
...rve 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 Survivors...
2006 Mar 20
0
Estimating Daily Survival
...0.605 0.524 -0.96 0.34
Likelihood ratio test=0.92 on 1 df, p=0.337 n= 19
> survfit(brettest1)$surv
[1] 0.9502716 0.9010356 0.8508737 0.7996929 0.7473818 0.6957615 0.5872297 0.5324263 0.4756545
[10] 0.4201535 0.3547775 0.2916481 0.2189030 0.1374015
>mydat
ID Year Dayrelease Agerelease Survivorship Entry Exit Fate Sex
1 16240 1996 205 95 164 205 369 1 0
2 16319 1996 205 88 140 205 345 1 0
3 16378 1996 248 108 100 248 348 1 0
4 20383 1996 241 98 204 241...