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
______________________________________________
[email protected] mailing list
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.
