Re: R-beta: Survival 4

Gregory R. Warnes (warnes@biostat.washington.edu)
Wed, 24 Sep 1997 09:06:45 -0700 (PDT)


Date: Wed, 24 Sep 1997 09:06:45 -0700 (PDT)
From: "Gregory R. Warnes" <warnes@biostat.washington.edu>
To: Ennapadam Venkatraman <VENKAT@biosta.mskcc.org>
Subject: Re: R-beta: Survival 4
In-Reply-To: <199709241411.QAA14201@aragorn.ci.tuwien.ac.at>




>>>>> Ennapadam Venkatraman writes:

> I installed the latest version of R (Version 0.50 Alpha-4) and it breaks 
> survival4 as follows (it worked fine under Version 0.50 Alpha-3)

I initially had the same problems with the survival 4 code and 0.50-r4,
but I have the large advantage of being classmates with Thomas Lumley who
does the port from Splus.  Working with him we resolved the problems. The
solution involves 2 steps. 
1) get the newest version of the survuval4 port from
     http://www.biostat.washington.edu/survival-4R1.0.tar.gz
2) and make a slight change to the terms.formula function in
     $RHOME/src/library/base/funs/models
   by adding two lines to the beginning to ensure that the data argument
   is a data frame:

        if(!is.null(data) & !is.data.frame(data))  ## Added 9/12/97 GRW
          data _ as.data.frame(data)               ## Added 9/12/97 GRW

   a patch to do this is at the end of this message.  

This corrected the problems I was having.  I hope it helps.

Greg

-------------------------------------------------------------------------------
    Gregory R. Warnes          | It is high time that the ideal of success
warnes@biostat.washington.edu  |  be replaced by the ideal of service.
                               |                       Albert Einstein
-------------------------------------------------------------------------------


====patch startes here-- run in $RHOME/src/library/base/funs =====
--- models.orig Wed Sep 24 08:57:29 1997
+++ models      Wed Sep 24 08:59:05 1997
@@ -58,6 +58,8 @@
 terms.formula <-
 function (x, specials = NULL, abb = NULL, data = NULL, keep.order =
FALSE) 
 {
+        if(!is.null(data) & !is.data.frame(data))  ## Added 9/12/97 GRW 
+          data _ as.data.frame(data)               ## Added 9/12/97 GRW
        new.specials <- unique(c(specials, "offset"))
        terms <- .Internal(terms.formula(x, new.specials, abb, data,
keep.order))
        offsets <- attr(terms,"specials")$offset
=====patch ends here===========



=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request@stat.math.ethz.ch
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=