PROC MIXED : slope difference between groups over time [Résolu]

[statasma]

Bonjour,

Bonne année et meilleurs vœux à tous.

Je veux calculer la pente entre deux groupes de traitement au fil du temps
j'hésite entre utiliser un modèle avec interaction ou pas, mais je pense qu'une "slope difference between groups over time" c'est le terme d'intéraction temps*groups, pour ceci j'utilise ce code:

Code :

1
2
3
4
5
6
7
 
proc mixed DATA = TABLE ORDER=internal;
	class GROUP;
	model y = PATID temps GROUP temps*GROUP/ solution;
	random intercept / subject=PATID;
	lsmeans GROUP/diff;
	run;

j'obtiens :

Solution for Fixed Effects

Standard
Effect Groupe Estimate Error DF t Value Pr > |t|

Intercept 46.5393 3.7802 112 12.31 <.0001
patid 0.03313 0.03791 333 0.87 0.3828
temps -0.06742 0.02647 333 -2.55 0.0113
GROUP a 2.3455 3.5193 333 0.67 0.5056
GROUP b 0 . . . .
temps*GROUP a 0.1114 0.03790 333 2.94 0.0035
temps*GROUP b 0 . . . .

Est ce que ma slope est bien 0.1114?

[Manoutz]

Je ne suis pas sur ce que ce que tu veux savoir. Ca veut dire quoi

Je veux calculer la pente entre deux groupes de traitement au fil du temps

Tu veux surement calculer les differences entre les pentes deux groupes?

manoutz

[statasma]

Je veux calculer la différence entre les pentes de deux groupes.

[olivier.decourt]

Dans le modèle

Code :

temps GROUP temps*GROUP

tu vas avoir :

• intercept 46.5393 --> ordonnée à l'origine du groupe b
• temps -0.06742 --> pente du temps pour le groupe b
• GROUP a 2.3455 --> différence d'ordonnée à l'origine entre groupes a et b
• GROUP b 0
• temps*GROUP a 0.1114 --> différence de pentes entre groupes a et b
• temps*GROUP b 0

[statasma]

Merci beaucoup pour cette réponse.

J'ai une autre question, pour l'interprétation:

Je veux dire que la différence de pente entre les deux groupes est de 0.1114 par "unités" de temps.



Sauf que ma variable temps comprend (1 ans 2 ans 3ans 5 ans) , j'ai des intervalles de temps irréguliers !

Comment puis je l'exprimer?

[olivier.decourt]

Elle est codée comment ta variable TEMPS ? 1, 2, 3 et 5 ???
Si oui, alors c'est exprimé par années.
Si c'est codé 1, 2, 3 et 4, c'est n'importe quoi.

[Manoutz]

Tu consideres le FACTEUR temps comme variable continuee ou categorielle (simple verif)? Les valeurs sont elles identiques pour tous les patients? Je rejoins la remarque d'Olivier

[olivier.decourt]

Je crois qu'elle est dans le modèle comme variable continue :

class GROUP;
model y = PATID temps GROUP temps*GROUP

Et d'ailleurs dans les sorties copiées dans le 1er message on a clairement une variable TEMPS quanti dans le modèle. Reste à connaître son codage.

J'ai oublié de préciser qu'elle pourrait tout autant être codée 12, 24, 36 et 60 et à ce moment-là les unités seraient des mois, ou idem avec des jours...

[statasma]

la variable temps est continue dans mon modèle et elle est codée en 12, 24, 36 et 60 mois.

Donc l'unité sera mois ?
Et la différence de pente entre les deux groupes est de 0.1114 par MOIS ?.

[olivier.decourt]

Je crois qu'on peut dire que les moyennes de Y entre les traitements a et b divergent à la vitesse de 0.1114 [unités de Y] par mois.

[statasma]

merci encore.

une dernière petite question pour bien comprendre ce modèle une fois pour toute,

si jamais javais la variable temps qui contenait (12 12.30 12.05 24 24.12 36 60 60.23 ...)

A ce moment la je ne peux pas considérer l'unité en mois?

[olivier.decourt]

Si bien sûr. La pente de TEMPS (et la différence entre pentes de TEMPS s'il y a une interaction) est exprimée en unités de Y / unités de temps.
Le tout étant que les unités de la variable TEMPS soient homogènes (ce qui est indispensable si elle est déclarée comme quanti, mais pas obligatoire si elle figurait aussi dans CLASS). Donc si tu as des valeurs de TEMPS qui ont des décimales, pas de souci.
Sauf évidemment si "60.23" signifie "60 mois et 23 jours" et pas "60 mois et 23/100e de mois". Alors là l'unité de temps est bien mystérieuse et je pense qu'un recodage de la variable s'impose pour mettre tout ça en mois et fractions de mois.

[statasma]

merci pour toutes ces réponses et précisions.