Je ne rentre pas plus dans les explications, voici mon code:
> Li<-c(17,19,23,24) # age à l'entrée dans létude
> ti<-c(25,27,26,28) # age à la sortie de létude
> deltai<-c(0,0,1,1) # 1=décédé ; 0=censuré
> Zi<- c(3.18,2.74,1.76,2.11) # variable explicative
> donnees<-data.frame(Li=Li,ti=ti,deltai=deltai,Zi=Zi)
> donnees
Li ti deltai Zi
1 17 25 0 3.18
2 19 27 0 2.74
3 23 26 1 1.76
4 24 28 1 2.11
> Zplus=c(2.13) # valeur de la variable explicative pour l'individu que l'on veut prédire
> covar<-data.frame(Zplus)
> regress1.cph <- coxph(Surv(Li,ti,deltai,type="counting")~Zi, data=donnees,method=c("efron"))
Warning message:
Ran out of iterations and did not converge in: fitter(X, Y, strats, offset, init, control, weights = weights,
> regress1.cph
Call:
coxph(formula = Surv(Li, ti, deltai, type = "counting") ~ Zi,
data = donnees, method = c("efron"))
coef exp(coef) se(coef) z p
Zi -56.6 2.67e-25 57062 -0.000992 1
Likelihood ratio test=2.2 on 1 df, p=0.138 n= 4
> s0<-survfit(regress1.cph,newdata=covar)
Warning message:
'newdata' had 1 rows but variable(s) found have 4 rows
Donc d'après ce que je comprends, ça ne veut pas me donner d'estimation pace qu'il n'y a pas le même nombre d'individus à prédire que d'individus qui ont servi à faire le modèle. Ou alors je ne pars pas dans la bonne direction?
Merci de votre aide.