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]]

