détermination de groupes de trajectoires avec proc traj [Résolu]

[mathis86]

Bonjour,

dans le cadre de mon stage je souhaite utiliser proc traj, un outil tout nouveau pour moi et que je ne maitrise pas encore.
Ma problématique est la suivante : j'ai à ma disposition 912 individus pour lesquels on observe un taux A (taux de créatinine) à plusieurs dates. Je souhaite avec proc traj pouvoir identifier deux groupes
Voici mon programme:

Code :

1
2
3
4
5
6
7
8
9
10
11
proc traj DATA=dn out=out outstat=os outplot=op;
var v1-v9;
indep years01-years09;
model cnorm;
min 20;
max 818;
ngroups 2;
ORDER 1 1;
id ID;
run;
%trajplot (OP, OS, “Title of graph”, “Subtitle”, “Y-axis label”, “X-axis label”);

Je ne suis pas sure de choix du model (je ne comprend vraiment l'option cnorm) et l'option order aussi.
Ce qui me pose le plus problème c'est que sas me retourne ce message d'erreur :

ERROR: Floating Point Zero Divide.
ERROR: Termination due to Floating Point Exception

Par ailleurs les tableau os et op sont vide.
Merci pour votre aide

[Manoutz]

première chose: identifier d'où vient l'erreur.

Tu utilises la proc traj puis une macro. A quel moment ca plante?

si ca vient de la macro nous ne pouvons pas t'aider car nous n'avons pas le contenu de la macro.

Bon courage,

Manoutz

[mathis86]

l'erreur arrive avant l'utilisation de la macro.
J'obtient

*** SAS PROC TRAJ *** Bobby L. Jones, Version Date: Mar 11 2009

NOTE: 912 observations read.
NOTE: 912 observations used in trajectory model.


ERROR: False convergence

Parameter estimates for adding risk factors

61.743610 0.000000 126.469013 0.000000 48.544052 0.000000
. . . . . .


Parameter estimates

61.743610 0.000000 126.469013 0.000000 48.544052 50.000000
50.000000 . . . . .

ERROR: Floating Point Zero Divide.
ERROR: Termination due to Floating Point Exception
seconds

etc.

Merci pour votre aide

[Manoutz]

Je ne connais pas plus que cela la proc traj, mais de toute évidence tu as un problème de convergence, et donc les résultats ne sont pas fiables.

On peut chercher, d'une manière générale, la cause de la non convergence part diverses pistes:

Mais avant de commencer: Es tu sure que cette procédure va te permettre de répondre à ton objectif sur le taux de créatinine? Pourquoi utilises tu une procédure si avancée soit la plus recommandée pour analyser des tests de laboratoire(je ne connais pas le contexte)?

Si tu veux nécéssairement utiliser la proc traj, Es tu sur de ton modèle? Si oui essaye de manipuler tes paramètres entrée (si possible) afin de voir si cela te permet de converger.

Manoutz

[mathis86]

je ne suis pas sure de mon modèle et je ne comprend pas d'ou vient la non convergence. Ce dont je suis sure c'est que proc traj peut répondre à ma question. La problèmatique est la suivante :
une cohorte de patients pour lesquels on mesure la créatine à plusieurs reprises. (plusieurs observations par patients, le nombre d'observation varie d'un individu à l'autre et les intervalles, les dates d'observations ne sont pas les memes).
On souhaite predire un événement (de type cardiovasculaire et autre). L'idée dans un premiers temps c'est de déterminer des groupes de trajectoires du paramètre biologique (la créatinine) avec proc traj (une trajectoire moyenne pour les patients qui font un événement et une autre pour les patients qui ne font pas d'événement) puis de modéliser le taux de créatinine avec une proc mixed.


Je te remercie pour ton aide

[mathis86]

ça fonctionne, j'ai tout simplement modifier quelques paramètres et réorganiser mes données

Code :

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
proc traj  DATA=b.don out=b.b outstat=b.os outplot=b.op ;
   id id;
   var v1-v8;
 
   indep tmp1-tmp8;
 
   model cnorm;
 
/*min 0; actually unnecessary because default is 0 */ 
min 10;
max 520;  /* always necessary but only pertinent to cnorm option */
ngroups 2;
ORDER 1 1;
run;
 
 proc print DATA=b.op;
 run;
 
  proc print DATA=b.os;
 run;
 
%trajplot(b.op,b.os,'titre,'sous titre','Creat (micromol/l)','Time until end of follow-up (years)')

Pour ceux que ça interesse voici un doc permettant d'expliquer l'utilisation de proc traj
https://circle.ubc.ca/handle/2429/834

et ici le lien pour la macro
http://www.andrew.cmu.edu/user/bjones/download.htm