Hi, I want to perform a survival analysis using survreg procedure from
survival library in R for a pareto distribution for a time variable, so I
set the new distribution using the following sintax:

     library(foreign)
     library(survival)
     library(VGAM)

     mypareto <- list(name='Pareto',
                  init= function(x, weights,parms){
etc.

The survreg routine fits location-scale distributions such that (t(y) - Xb)/s ~ F, where t is an optional transformation, F is some fixed distribution and X is a matrix of covariates. For any distribution the questions to ask before trying to add the distribution to survreg are
  - can it be written in a location-scale form?
  - if so, how do the parameters of the distribution map to the location (Xb) 
and scale (s).

In fitting data we normally have per-subject location (X b) but an intercept-only model is of course possible.

If y is Weibull then log(y) fits into the framework, which is how survreg fits it. The transformation of parameters location and scale parameters for log(y) back to the usual Weibull parameterization for y often trips people up (see comments in the Examples section of ?survreg).

The log of a Pareto can be written in this form (I think?). The two parameters are the scale a and lower limit b, with survival function of S(x)= (b/x)^a, for x >= b. If y = log(x) the survival function for y is S(y) = (b/exp(y))^a = exp[-(y - log(b))/(1/a)], which has location log(b) and scale 1/a. But even if I am correct the discontinuity at b will cause the underlying Newton-Raphson method to fail.

 Terry Therneau

______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
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.

Reply via email to