On Aug 4, 2012, at 10:45 AM, hafida wrote:
Dear R-community,
I have tried to estimate an EXPONENTIEL accelerated failure time(AFT)
power rule model with time-independent . For that purpose, I have
used
the eha package.
Please, consider this example:
vi Ti
1 26 5.79
2 26 1579.52
3 26 2323.70
snipped
74 38 1.13
75 38 0.09
76 38 2.38
aftexp<-aftreg(Surv(time,status) ~ vi, data=data.frame(dataexp),
dist="exponentiel")
Error in Surv(time, status) : Time variable is not numeric
Your variables being given to the `Surv` function do not match those
in you data object.
aftexp<-aftreg(Surv(Ti,status) ~ vi, data=data.frame(dataexp),
dist="exponentiel")
Error in Surv(Ti, status) : object 'status' not found
status<- rep(1, 76)
status
[1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
1 1 1 1
1 1
[39] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
1 1 1 1
1 1
cbind(dataexp, status)
vi Ti status
1 26 5.79 1
2 26 1579.52 1
snipped
72 38 0.74 1
73 38 0.39 1
74 38 1.13 1
75 38 0.09 1
76 38 2.38 1
aftexp<-aftreg(Surv(Ti,status) ~ vi, data=data.frame(dataexp),
dist="exponentiel")
Error in aftreg.fit(X, Y, dist, strats, offset, init, shape, id,
control, :
exponentiel is not an implemented distribution
You are spelling 'exponential' incorrectly (in English anyway.)
Bon chance;
David.
aftexp<-aftreg(Surv(time,status) ~ vi, data=data.frame(dataexp),
dist="exponentiel")
Error in Surv(time, status) : Time variable is not numeric
pleas help me to find a solution to my problem
--
David Winsemius, MD
Alameda, CA, USA
______________________________________________
R-help@r-project.org 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.