survreg (survival) reports erroneous results for left-censored data (PR#2287)

Thomas Lumley tlumley@u.washington.edu
Thu, 14 Nov 2002 08:04:37 -0800 (PST)


On Wed, 13 Nov 2002, Jan de Leeuw wrote:

>
> No problemo. And, in fact, I get the same results in
> the R-1.6.0 Carbon version.

I don't. Could there be a G3/G4 issue?

	-thomas


> --- Jan
>
> On Wednesday, November 13, 2002, at 02:05 PM, tim@timcohn.com wrote:
>
> > Full_Name: Tim Cohn
> > Version: 1.6.1
> > OS: Macintosh OS X
> > Submission from: (NULL) (130.11.34.250)
> >
> >
> > The Mac version of survreg does not handle left-censored data
> > correctly (at
> > least the results are not what I get doing it other ways, and they are
> > not the
> > same as I get running R 1.6.1 in Windows 98se;  the Windows 98 results
> > are
> > correct).  On the windows version of R 1.6.1.
> >
> >> survreg(Surv(c(0,1,-1),c(T,T,F),type='left')~1,dist='gaussian')
> > Call:
> > survreg(formula = Surv(c(0, 1, -1), c(T, T, F), type = "left") ~
> >     1, dist = "gaussian")
> >
> > Coefficients:
> > (Intercept)
> >  -0.2321569
> >
> > Scale= 1.161135
> >
> > Loglik(model)= -4.1   Loglik(intercept only)= -4.1
> > n= 3
> >
> >
> > =================
> >
> > The Macintosh OS X result for the same command is:
> >
> >> survreg(Surv(c(0,1,-1),c(T,T,F),type='left')~1,dist='gaussian')
> > Call:
> > survreg(formula = Surv(c(0, 1, -1), c(T, T, F), type = "left") ~
> >     1, dist = "gaussian")
> >
> > Coefficients:
> > (Intercept)
> >  -0.7373613
> >
> > Scale= 1.366282
> >
> > Loglik(model)= -5.6   Loglik(intercept only)= -5.6
> > n= 3
> >
> >
> > -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
> > .-.-.-.-.-
> > r-devel 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-devel-request@stat.math.ethz.ch
> > _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
> > ._._._._
> >
> >
> ===
> Jan de Leeuw; Professor and Chair, UCLA Department of Statistics;
> Editor: Journal of Multivariate Analysis, Journal of Statistical
> Software
> US mail: 9432 Boelter Hall, Box 951554, Los Angeles, CA 90095-1554
> phone (310)-825-9550;  fax (310)-206-5658;  email: deleeuw@stat.ucla.edu
> homepage: http://gifi.stat.ucla.edu
>
> ------------------------------------------------------------------------
> -------------------------
>            No matter where you go, there you are. --- Buckaroo Banzai
>                     http://gifi.stat.ucla.edu/sounds/nomatter.au
>
> ------------------------------------------------------------------------
> -------------------------
>
> -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
> r-devel 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-devel-request@stat.math.ethz.ch
> _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
>

Thomas Lumley			Asst. Professor, Biostatistics
tlumley@u.washington.edu	University of Washington, Seattle
^^^^^^^^^^^^^^^^^^^^^^^^
- NOTE NEW EMAIL ADDRESS

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-devel 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-devel-request@stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._