Hi all,
I am using the coxme function to fit random effects survival models. The
base survival models show a typical "J" shaped curve where survival drops
sharply then plateaus over time interval. I want to include a time
covariate which could account for this non-linear decrease in survival over
time. However, I have not found a clear way to do this in the coxme
models. Any suggestions would be greatly appreciated. Below is some
sample code and data.
-Jared
selection=coxme(Surv(start,stop,status)~prob+(1|year);
fawn
status
start
stop
prob
year
3
0
0
86
0.83
2009
4
1
0
9
0.6
2009
5
0
0
355
0.63
2009
6
1
0
6
0.52
2010
7
0
0
68
0.79
2010
8
0
0
354
0.69
2010
9
1
0
3
0.72
2010
10
1
0
4
0.77
2010
11
1
0
12
0.62
2010
13
0
0
82
0.7
2011
14
1
0
7
0.67
2011
15
1
0
74
0.49
2011
16
0
0
351
0.64
2011
17
0
0
351
0.74
2011
18
0
0
82
0.72
2011
19
1
0
342
0.69
2011
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
