Re: [R] parametric proportional hazard regression

From: Thomas Lumley <tlumley_at_u.washington.edu>
Date: Sat 08 Jul 2006 - 01:47:59 EST

On Fri, 7 Jul 2006, Valentin Dimitrov wrote:
> I am trying to find a suitable R-function for
> parametric proportional hazard regressions. The
> package survival contains the coxph() function which
> performs a Cox regression which leaves the base hazard
> unspecified, i.e. it is a semi-parametric method. The
> package Design contains the function pphsm() which is
> good for parametric proportional hazard regressions
> when the underlying base distribution is "weibull" or
> "exponential". But what if I need a parametric
> proportional hazard model with the other "usual"
> distributions like '"gaussian"', '"logistic"',
> '"lognormal"'and '"loglogistic"?

Those are not proportional hazards families of distributions. That is, if the distribution is gaussian for one value of the hazard ratio parameters it will not be gaussian for any other value.

You can get accelerated failure models with these distributions using survreg() in the survival package.

         -thomas

Thomas Lumley			Assoc. Professor, Biostatistics
tlumley@u.washington.edu	University of Washington, Seattle

______________________________________________
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Sat Jul 08 02:06:30 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Sat 08 Jul 2006 - 12:15:41 EST.

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