search for: pyear

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 (female,male) ev, event(death), 0= censored at time point dof, 1...

Dear List, Is there an automated way to use the survival package to generate survival rate estimates and their standard errors? To be clear, *not *the survivorship estimates (which are cumulative), but the survival *rate * estimates... Thank you in advance for any help. Best, Brian [[alternative HTML version deleted]]

...urray Jorgensen <maj at waikato.ac.nz> To: Prof Brian Ripley <ripley at stats.ox.ac.uk> References: <44C063E5.3020703 at waikato.ac.nz> <Pine.LNX.4.64.0607210716270.12611 at gannet.stats.ox.ac.uk> Apologies for a non-selfcontained example. Here is what I should have sent: pyears <- scan() 18793 52407 10673 43248 5710 28612 2585 12663 1462 5317 deaths <- scan() 2 32 12 104 28 206 28 186 31 102 l <- log(pyears) Smoke <- gl(2,1,10,labels=c("No","Yes")) Age <- gl(5,2,10,labels=c("35-44","45-54","55-64","6...

...ted the offending expression has been changed (mostly for other reasons). The problem case in question no longer is generated. If fixing this issue was on anyone's list of things to do, feel free to drop the priority to zero. (You could still generate this with just exactly the right call to pyears, but I'll be boarding up that final door in 2.36-2.) Terry Therneau

Consider the following example (based on an example in Pat Altham's GLM notes) pyears <- scan() 18793 52407 10673 43248 5710 28612 2585 12663 1462 5317 deaths <- scan() 2 32 12 104 28 206 28 186 31 102 Smoke <- gl(2,1,10,labels=c("No","Yes")) Age <- gl(5,2,10,labels=c("35-44","45-54","55-64","65-74","75-...

...ng a package on CRAN is that it have to pass R CMD check with no errors, warnings *or* notes. So, in that sense notes are no different from warnings. --------------------------------- Getting rid of these notes would be very hard in the survival package. The population survival routines (survexp, pyears) that have to deal with dates try to be very accommodating wrt the date format the user chose, so have constructions like if (user date looks like a chron object) { calculation using chron functions } else if (is from the date package) { use mdy.date } else if (D...

...agexact.o agfit2.o agfit3.o agfit5.o agfit_null.o agmart.o agmart2.o ags core.o agsurv1.o agsurv2.o agsurv3.o char_date.o chinv2.o chinv3.o cholesky2.o c holesky3.o chsolve2.o chsolve3.o cox_Rcallback.o coxdetail.o coxfit2.o coxfit5.o coxmart.o coxph_wtest.o coxscho.o coxscore.o dmatrix.o doloop.o pyears1.o pyear s2.o pyears3.o pystep.o surv_callback.o survdiff2.o survfit2.o survfit3.o survin dex2.o survindex3.o survreg2.o survreg3.o survreg4.o survreg5.o ** R ** data ** moving datasets to lazyload DB ** inst ** preparing package for lazy loading /bi/common/src/R-2.1.1/bin/INSTALL[325]: 3667 M...

...ithout warnings: 1) unbalanced braces: * Rd files with unbalanced braces: * man/Surv.Rd * man/cluster.Rd * man/cox.zph.Rd * man/coxph.Rd * man/coxph.detail.Rd * man/date.ddmmmyy.Rd * man/lines.survfit.Rd * man/plot.cox.zph.Rd * man/predict.survreg.Rd * man/pyears.Rd * man/residuals.coxph.Rd * man/stanford2.Rd * man/strata.Rd * man/summary.survfit.Rd * man/survdiff.Rd * man/survexp.Rd * man/survfit.Rd * man/survreg.control.Rd * man/survsum.Rd * man/untangle.specials.Rd 2) missing docs: * checking for undocumented obj...