> 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
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
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
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
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
>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
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
7 matches
Mail list logo