[R] Survival: pyears and ratetable: expected events

[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 different variation of the pyears syntax but i  
> get
> either 0.05 expected deaths or something like 25000. I have used  
> different
> age and fu formats but I get no plausible result
>
> Could anyone help me with the syntax?
>
> 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

Have you constructed a ratetable by that name? I do not see one in the  
datasets in the survival package.

>
> sex is a factor (male, female)
> ev=1 is death ev=0 censored at last fu
> fu is in days
> age is in days
> year of surgery is as.Date
>
>
> Any ideas?
> JT
>
> --
> View this message in context: http://r.789695.n4.nabble.com/Survival-pyears-and-ratetable-expected-events-tp3553208p3556443.html
> Sent from the R help mailing list archive at Nabble.com.
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.

David Winsemius, MD
West Hartford, CT