[R] HELP

[Nekeisha Spencer]

 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