Re: [R] Use of survreg.distributions

From: Abdus Sattar <upsattar_at_yahoo.com>
Date: Fri, 25 Apr 2008 10:02:21 -0700 (PDT)


Hello Dr. Therneau:
 

Thank you for your response. Let me explain to you want I want. My y variable(time) is normal and would like to fit the model on logarithmic transformation of y (log(y)). I tried to run codes according to your suggestion: tfit=survreg(Surv(y, y>=-5, type="left")~x + cluster(id), dist="lognormal", data=y.data, scale=0, weights=w)

The following error message is giving:
 

Error in survreg(Surv(y, y >= -5, type = "left") ~ x +  :
  Invalid survival times for this distribution
In addition: Warning messages:

1: In log(dlist$dtrans(Y[exactsurv, 1])) : NaNs produced 2: In log(y) : NaNs produced

Note, the data file y.data does not contain any missing data! Do you know why it is giving me such an error message please? Thank you again for your helpful comment/suggestion.
 

Best Regards,
 

Abdus Sattar
upsattar_at_yahoo.com

--begin included message ---

I am using survreg(Surv()) for fitting a Tobit model of left-censored longitudinal data. For logarithmic transformation of y data, I am trying use survreg.distributions in the following way: tfit=survreg(Surv(y, y>=-5, type="left")~x + cluster(id), dist="gaussian", data=y.data, scale=0, weights=w)

my.gaussian<-survreg.distributions$gaussian
my.gaussian$name="lognormal"
my.gaussian$dist<-my.gaussian

tfit=survreg(Surv(y, y>=-5, type="left")~x + cluster(id), dist=my.gaussian, data=y.data, scale=0, weights=w)

If I run these codes then I got the following error message,

Error in survreg(Surv(y, y >= -5, type = "left") ~ x +  : Invalid distribution object

Does anybody can help me in identifying the error(s) in these code please?

  Can you tell us what you are trying to do? 
  Your first model was a fit of y ~ x + eps, eps ~ Gaussian.  If what you want
is log(y) ~ x + eps, then all that you need do is use dist="loggaussian" in the survreg call.  (Or 'lognormal'; which is the same distribution.)
 
      Terry Therneau


[[elided Yahoo spam]]
=Ahu06i62sR8HDtDypao8Wcj9tAcJ

        [[alternative HTML version deleted]]



R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Fri 25 Apr 2008 - 17:31:27 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Fri 25 Apr 2008 - 18:30:33 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive