Discussion :

[craow87000]

Bonjour,

j'ai un fichier dans lequel pour chaque individu (1 ligne par individu) j'ai une variable oui/non qui se repete à differents temps : mes patients prennent t'il trop de médicaments (oui/non) 1 mois, 2 mois, 3 mois..; aprés leur sortie d'hospitalisation. La question est : il y a t'il une variation au niveau du nombre de personnes pernant trop de médicament dans le temps?
Apres une recherhce sur internet j'ai vu qu'il fallait faire un test du Q de cochran (option agree de proc freq) mais en pratique comment organiser ma commande. comme cela?

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
proc freq data=X;
tables variableT1*variableT2*...*variable Tn/agree;
run;

quand je fait cette commande je ne trouve pas le test voulu dans les sorties et j'ai l'impression qu'il teste les variables 2à 2.

Sinon, j'ai pour certaines variables des valeurs manquantes car certains individus sont DCD durant l'étude ou autre.

Merci par avance pour votre aide

[frajen]

Bonjour,

Vu le test que tu veux faire, regarde la procédure TTEST dans l'aide qui est très bien faite, avec l'option cochran et tu auras les résultats souhaités.

[olivier.decourt]

Bonjour.
A priori la syntaxe que tu proposes (hors valeurs manquantes) est la bonne. Elle rejoint celle de la doc de la proc Freq.
Ce sont les valeurs manquantes qui compliquent la donne. Tu peux essayer en te restreignant aux individus complets ? Est-ce que ça a un sens dans ton étude ? Est-ce qu'ils sont assez nombreux ?

@Frajen : le Q de Cochran et son approximation pour les tests de comparaison de moyennes à variances inégales ne sont pas du tout les mêmes tests. Ici c'est bien un Q dont on a besoin, donc la proc Freq.

Bon courage.
Olivier

[craow87000]

Bonjour,

Si je me restreint aux seuls patients qui ont des données sur les 12 mois de suivi il ne me restera pas grand monde. Si c'est la bonne commande pourquoi je n'ai pas le Q à la fin comme dans les exemple de commande pour l'option Agree?

Merci pour votre aide

[craow87000]

du coup j'ai fait mon test par un autre moyen :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
 
proc genmod data=medinap.k descending;
class nir _name_(param=ref ref='polymedicsortiessrg');
model polymedication1=_NAME_/dist=bin link=log type3;
repeated subject=nir /corrw type=unstr;
run;

nIir étant l'identifiant patient, _name_ la période