Bonjour,
J'ai developpé un modèle de cox sous sas avec des covariables dépendantes du temps.
Voici le programmela variable LG_ARR_CM est binaire et prend les valeurs 0/1)
PROC PHREG DATA=SERO_EXPO;
MODEL TPsS_HCV*Inf_HCV(0)=LG_ARR_CM ;
ARRAY SUIVI_COL [6] SUIVI_1-SUIVI_6;
ARRAY LG_ARR_COL [6] LG_ARR_1-LG_ARR_6;
DO i=1 TO 6;
IF SUIVI_COL[i]<TPsS_HCV AND SUIVI_COL[i] NE . THEN LG_ARR_CM=LG_ARR_COL[i];
END;
RUN;
J'aimerai maintenant tester l'hypothès de proportionnalité des risques.
Pour cela j'ai essayé de testé l'interaction entre le temps et la covariable mais cela ne marche pas. Voici le modèle que j'ai fait:
PROC PHREG DATA=SERO_EXPO;
MODEL TPsS_HCV*Inf_HCV(0)=LG_ARR_CM LG_ARR_CM_time;
LG_ARR_CM_time=LG_ARR_CM*TPsS_HCV;
ARRAY SUIVI_COL [6] SUIVI_1-SUIVI_6;
ARRAY LG_ARR_COL [6] LG_ARR_1-LG_ARR_6;
DO i=1 TO 6;
IF SUIVI_COL[i]<TPsS_HCV AND SUIVI_COL[i] NE . THEN LG_ARR_CM=LG_ARR_COL[i];
END;
RUN;
Est-ce qu'il y a une erreur dans ma programmation?
Voici le log:
NOTE: 18 event observations have missing values in the
time-dependent covariates. They are not counted as event
observations, but they may still be used in the risk sets
of other event times.
NOTE: Convergence criterion (GCONV=1E-8) satisfied.
NOTE: PROCEDURE PHREG used (Total process time):
real time 0.10 seconds
cpu time 0.01 seconds
Merci de m'avoir lu, bonne journée
Rose
Partager