[R] HELP

Nekeisha Spencer
Wed, 04 May 2011 04:30:01 -0000

Accelerated Failure Time Model to Proportional Hazard Form

Greetings R users:

I have been working on a problem for  a while and can't seem to get any
result.
I am trying to convert  accelerated failure time estimates to proportional form.
I keep getting an error that I can't understand and don't know if it can
be
debugged or not. My code is as follows:


data<-read.csv("repdata.csv")
da<-na.omit(data)
attach(da)
library(survival)
library(rms)
weib1<-psm(Surv(time,ratified)~island+pop+gdppc+fh+nngo+industryemplimp+icrgfromper+ngoicrgfromper+industryemplimpicrgfromper,
 data=da, dist="weibull",init=NULL, scale=0,
control=survreg.control(maxiter=1000))
print(weib1)
f.ph<-pphsm(weib1)

However, after running the last line of code, I get  the following:

Warning message:
In pphsm(weib1) :
  at present, pphsm does not return the correct covariance matrix


 Any suggestions will be greatly appreciated!! My data is attached.

Thanks in Advance


 NSPENCER