Dear all,
I have a data
set<http://yaap.it/paste/c11b9fdcfd68d02b#gIVtLrrme3MaiQd9hHy1zcTjRq7VsVQ8eAZ2fol1lUc=>with
6 clusters, each containing 48 (possibly censored, in which case
"event = 0") survival times. The "x" column contains a binary explanatory
variable. I try to describe that data with a gamma frailty model as follows:
library(survival)
mod <- coxph(Surv(time, event) ~
x + frailty.gamma(cluster, eps=1e-10, method="em", sparse=0),
outer.max=1000, iter.max=10000,
data=data)
Here is the error message:
Error in if (history[2, 3] < (history[1, 3] + 1)) theta <-
mean(history[1:2, :
missing value where TRUE/FALSE needed
Does anyone have an idea on how to debug?
Yours sincerely,
Marco
[[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.
