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

Jan de Leeuw deleeuw@stat.ucla.edu
Thu, 14 Nov 2002 08:09:27 -0800


I take that back. I now get the "correct" result on Darwin R-1.6.1, the  
"wrong" result
with Carbon R-1.6.0. On a dual G4 with 10.2.2.

On Thursday, November 14, 2002, at 08:04 AM, Thomas Lumley wrote:

> 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
>
>
===
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._