Re: [R] survreg 3-way interaction

2011-01-28 Thread tdenes
The fire devoured my laptop, but the PC is still working... :) Sorry for the lazyness, you are totally right. Here goes a reproducible example, which resembles the main features of our dataset: # three experimental factors, 30 subjects expgrid <- expand.grid(F1=0:1,F2=0:1,F3=0:1,id=1:30) # set

Re: [R] survreg 3-way interaction

2011-01-28 Thread Terry Therneau
> I was wondering why survreg (in survival package) can not handle > three-way interactions. I have an AFT . You have given us no data to diagnose your problem. What do you mean by "cannot handle" -- does the package print a message "no 3 way interactions", gives wrong answers, your laptop

[R] survreg 3-way interaction

2011-01-28 Thread tdenes
Dear All, I was wondering why survreg (in survival package) can not handle three-way interactions. I have an AFT (accelerated failure time) model with a frailty term. If the model is restricted to two-way interactions, survreg and gamlss.cens (with random() term) give very similar results, except