Re: [R] relsurv package

2008-09-11 Thread Giulia Barbati
ls ? Giulia - Original Message From: Prof Brian Ripley <[EMAIL PROTECTED]> To: Yihui Xie <[EMAIL PROTECTED]> Cc: Giulia Barbati <[EMAIL PROTECTED]>; r-help@r-project.org Sent: Wednesday, September 10, 2008 11:44:43 AM Subject: Re: [R] relsurv package On Wed, 10 Sep 2008,

Re: [R] relsurv package

2008-09-10 Thread Giulia Barbati
Dear R-users, I have a couple of questions about the relsurv package: 1) when I try to run the example: fit <- rsmul(Surv(time,cens)~sex+as.factor(agegr)+ratetable(age=age*365.24,sex=sex,year=year),ratetable=slopop,data=rdata) with the datasets in the package (rdata and slopop) it gives me

[R] hazard ratio of interaction Cox model

2007-12-18 Thread Giulia Barbati
Dear Forum, I have a question about interaction estimate in the Cox model: why the hazard ratio of the interaction is not produced in the summary of the model? (Instead, the estimate of the coefficient is given in the print of the model.) # Example: modINT <-cph( Surv(T_BASE, T_FIN,STATUS)~ N

[R] datadist options error, DESIGN library

2007-10-31 Thread Giulia Barbati
Hello, using the Design library, and the following command (from the Harrell's book, example at Ch.20): dd <- datadist(rx, age, wt, pf, pf.coded, heart, map, hg, sz, sg, ap, bm) options(datadist=='dd') I get the following error: Error in datadist == "dd" : comparison (1) possible only for atomi