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 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@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.
