Re: [R] Survival: pyears and ratetable: expected events

2011-05-31 Thread Terry Therneau
> Should the following syntax work getting the expected number of > deaths? I am > having no luck getting the thing to work. > pyears(Surv(fu,ev)~ratetable(sex=sex, age=age,year=yearofsurgery), > ratetable=expsurv.us, scale=365.25)$expected Your syntax is fine. It looks like you built your own r

Re: [R] Survival: pyears and ratetable: expected events

2011-05-28 Thread Jim Trabas
Hi thanks, The conventional method is solid (cross multiplication of patient years at each age and year and gender, with the corresponding risk for death at each age year and gender). I get about 26 expected deaths (verified by many different sources) Now what I am trying to do is get the same a

Re: [R] Survival: pyears and ratetable: expected events

2011-05-27 Thread David Winsemius
On May 27, 2011, at 5:35 PM, Jim Trabas wrote: I tried again but I cannot get a valid expected number of death directly from the pyear function. Using cross multiplication from Lifetables I am expecting in my cohort about 26 expected deaths (to match the 25 observed deaths). I have tried di

Re: [R] Survival: pyears and ratetable: expected events

2011-05-27 Thread Jim Trabas
I tried again but I cannot get a valid expected number of death directly from the pyear function. Using cross multiplication from Lifetables I am expecting in my cohort about 26 expected deaths (to match the 25 observed deaths). I have tried different variation of the pyears syntax but i get either

Re: [R] Survival: pyears and ratetable: expected events

2011-05-27 Thread Jim Trabas
Thanks, of course I am familiar with your work. I've read 63 but I think my answer may be in 81 (reading now). I've tried running the examples I've found in the help of the survival package, they seem to work fine, but when I parse my own variables everything gets mixed up. pyears wotks fine but wh

Re: [R] Survival: pyears and ratetable: expected events

2011-05-27 Thread Terry Therneau
>I am having a (really) hard time getting pyears to work together with a > ratetable to give me the number of expected events (deaths). The most complete discussion of pyears and rate tables is in a technical report from the Biostatistics group at Mayo http://mayoresearch.mayo.edu/mayo/researc

[R] Survival: pyears and ratetable: expected events

2011-05-26 Thread Jim Trabas
Dear all, I am having a (really) hard time getting pyears to work together with a ratetable to give me the number of expected events (deaths). I have the following data: dos, date of surgery, as.Date dof, date of last follow-up, as.Date dos, date of surgery, as.Date sex, gender, as.factor (femal