[R] Gaussian frailty leads to segmentation fault
tlumley at u.washington.edu
Wed Jul 28 20:03:08 CEST 2004
We really need a reproducible example to find segmentation faults. Can
you make one?
On Wed, 28 Jul 2004, Christian Lederer wrote:
> Dear R gurus,
> for a simulation concerning study effects and historical controls
> in survival analysis, i would like to experiment with a gaussian
> frailty model.
> The simulated scenario consists of a randomized trial
> (treatment and placebo) and historical controls (only placebo).
> So the simulated data frames consist of four columns
> $time, $cens, $study, $treat.
> $time, $cens are the usual survival data.
> For the binary thretment indicator we have
> $treat == 0 or 1, if $study == 1,
> $treat == 1 if $study > 1
> Typical parameters for my simulations are:
> sample sizes (per arm): between 100 and 200
> number of historical studies: between 7 and 15
> hazard ratio treatment/placebo: between 0.7 and 1
> variance of the study effekt: between 0 and 0.3
> Depending on the sample sizes, the following call sometimes leads to
> a segmentation fault:
> coxph(Surv(time,cens) ~
> as.factor(treatment) + frailty(study, distribution="gaussian"),
> I noticed, that this segmentation fault occures most frequently, if the
> number of randomized treatment patients is higher than the number of
> randomized placebo patients, and the number of historical studies is
> There seems to be no problem, if there are at least as many randomized
> placebo patients as treated patients. Unfortunately, this is not the
> situation i want to investigate (historical controls should be used
> to decrease the number of treated patients).
> Is there a way to circumwent this problem?
> Is it allowed, to attach gzipped sample data sets in this mailing list?
> R-help at stat.math.ethz.ch mailing list
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Thomas Lumley Assoc. Professor, Biostatistics
tlumley at u.washington.edu University of Washington, Seattle
More information about the R-help