similar to: Survival: pyears and ratetable: expected events

A few weeks ago I reported a problem with model.frame, whose root lay in a formula expression "....+ ratetable(x1=x1, x2=x2, ....x100=x100)" that was really long and caused model.frame to fail. Brian had some indefinite ideas on what might need to change in the base code to handle it. In survival_2.36-1 the bit of code that generated the offending expression has been changed (mostly

I am not sure if this is the right place for that kind of questions, but I wondered that the recommended package survival did not pass R's check procedure without 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 *

Hello, I'm attempting to use the ratetable argument to survexp in the survival package. I use the example from the ?survexp help page below, and then slightly modify it to produce an error. library(survival) data(pbc) #fit a model without any factors pfit1 <- coxph(Surv(time, status > 0) ~ trt + log(bili) + log(protime) + age + platelet, data=pbc) #this works as expected

Dear All, reposting, because I did not find a solution, maybe someone could check the example below. It's taken from the help page of survdiff. Executing it, gives the error "Error in floor(temp) : Non-numeric argument to mathematical function" best regards, Heinz library(survival) ## Example from help page of survdiff ## Expected survival for heart transplant patients based

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-84"),

Hello, I probably have a syntax error in trying to generate an expected survival curve from a weighted cox model, but I can't see it. I used the help sample code to generate a weighted model, with the addition of a "weights=albumin" argument (I only chose albumin because it had no missing values, not because of any real relevance). Below are my code with the resulting error

This is a follow-up on my note of Saturday. Let me start with two important clarifications - I think this would be a nice addition, but I've had exactly one use for it in the 15+ years of developing the survival package. - I have a work around for the current case. Prioritize accordingly. The ideal would be to change survexp as follows: fit <- survexp( ~ gender,

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]]

Dear R-users, I have a couple of questions about the relsurv package: 1) when I try to run the example: fit <- rsmul(Surv(time,cens)~sex+as.factor(agegr)+ratetable(age=age*365.24,sex=sex,year=year),ratetable=slopop,data=rdata) with the datasets in the package (rdata and slopop) it gives me an error: Error in nrow(x) : object "x" not found 2) If I have a date format

Dear All, Does someone have, or know of survexp ratetables for european contries, especially Austria and Germany? I know only about slopop in the package relsurv. Thanks in advance Heinz T??chler

The survexp function can fail when called from another function. The "why" of this has me baffled, however. Here is a simple test case, using a very stripped down version of survexp: survexp.test <- function(formula, data, weights, subset, na.action, rmap, times, cohort=TRUE, conditional=FALSE, ratetable=survexp.us, scale=1, npoints, se.fit,

survexp in survival5 doesn't seem to work for me. see below: > library(survival5) Attaching Package "package:survival5": The following object(s) are masked from package:base : sort.list > library(chron) > data(ratetables) > survexp(~ratetable(year=julian(6,1,1991), + sex=1,age=35*365.24),times=(0:30)/6*365.24) Error in as.character(as.date(c(min(R[, 3]),

Dear R-Users, I've recently used relsurv package for relative survival analysis. In particular I've tried to reproduce the examples proposed in the R-documentation about rsadd, rsmul and rstrans functions in R latest version (R 2.14.1). These examples don't run and the error message is always the following: data(slopop) data(rdata)

Hello, I have a question about the function lifetab in package KMsurv. The description of the output value surv says "the estimated survival function at the start of the intervals". Are these estimates the ones calculated via Kaplan-Meier probability of survival ? Thanks in advance! -- View this message in context:

Bother! This cold has made me accident-prone. I meant to hit Reply-all. Clarification below. -------- Original Message -------- Subject: Re: [R] Parameterization puzzle Date: Fri, 21 Jul 2006 19:10:03 +1200 From: Murray Jorgensen <maj at waikato.ac.nz> To: Prof Brian Ripley <ripley at stats.ox.ac.uk> References: <44C063E5.3020703 at waikato.ac.nz>

> Several changes in print.survfit, plot.survfit and seemingly in the structure > of ratetabels effect some of my syntax files. > Is there somewhere a documentation of these changes, besides the code itself? I agree, the Changelog.09 file is not as comprehensive as one would like. Specific comments: 1. The ratetables were recently changed to accomodate a new option. I thought

Henrik wrote: I think what people are also thinking about is that the policy for publishing 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)

Hi everybody, trying to get the Dovecot antispam_plugin to work and I must be doing something wrong, because as soon as it is enabled with a certain backend, imap logins do not work anymore (the session is immediately closed after a successful login). Interestingly, pipe and spool2dir are working (that is, the session won't be closed), dspam-exec and crm114-exec are not. If this happens,

Dear All, when I try to reproduce an example of survexp, taken from the help page of survdiff, I receive the error message "Error in floor(temp) : Non-numeric argument to mathematical function" . It seems to come from match.ratetable. I think, it has to do with character variables in a ratetable. I would be interested to know, if it works for others. With an older version of

The model.frame function has trouble with a certain type of really long formula. Here is a test: tname <- paste('var', 1:50, sep='') tmat <- matrix(rnorm(500), ncol=50, dimnames=list(NULL, tname)) tdata <- data.frame(tmat) temp1 <- paste( paste(tname, tname, sep='='), collapse=', ') temp2 <- paste("~1 + cbind(", temp1, ")")