Re: [R] can we include nonparametric component for survival regression?

From: gallon li <>
Date: Tue, 19 Feb 2008 19:19:55 +0800

Thanks a lot, Prof Lumley.

Now I can fit a model like


but I am not sure how to extract the estimated function for x1 alone. I tried to use predict function but couldn't find appropriate option to do this.

if i only have one covariates, then the example in help manule can be used as
plot(x1, predict(coxfit))

but with more than 1 predictor, i am not sure how to select the one i want.

On 2/19/08, Thomas Lumley <> wrote:
> On Mon, 18 Feb 2008, gallon li wrote:
> > i am trying to fit a survival regression model (cox model or parametric
> > model) in R by including the covariate effects as a function m(x)
> instead of
> > just beta*x. is it possible to fit such a model? can someone recommend
> some
> > reference? I searched but only found a package called addreg where
> > the hazard is actually modeled additively. That is not what i want.
> The survival package has pspline() for this purpose (assuming that by
> 'non-parametric' you mean a flexible smooth curve)
> -thomas
> Thomas Lumley Assoc. Professor, Biostatistics
> University of Washington, Seattle

        [[alternative HTML version deleted]] mailing list PLEASE do read the posting guide and provide commented, minimal, self-contained, reproducible code. Received on Tue 19 Feb 2008 - 11:22:13 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 Tue 19 Feb 2008 - 12:30:15 GMT.

Mailing list information is available at Please read the posting guide before posting to the list.

list of date sections of archive