Discussion :

[xavier_dcf]

Bonjour,

J'ai un doute sur l'interprétation de résultats (odds ratios) pour une proc logistic :

je teste une régression logistique sur une variable d'intérêt Y(rémuneration_ifaq)=0 ou 1 VS une variable d'interaction (statut*esatis)
J'ai pris comme références statut="ex_dgf" et esatis="NA" comme indiqués dans "capture sas 1"

Je veux étudier les odds ratios ' capture sas 2" .
Selon les sorties de "capture sas 2", est ce que je peux dire cela: si esatis=10 on a 0.5 fois moins de chance d'avoir remuneration_ifaq=1 si on est " ex dg" par rapport à "ex oqn".
Est ce la bonne interprétation? Je pense faire une erreur car si j'inverse la référence pour statut , je ne trouve pas
si esatis=10 on a 2 fois plus de chances d'avoir remuneration_ifaq=1 si on est " ex oqn" par rapport à "ex dg".



En gros j'ai un problème d'interprétation pour les odds ratios appliqués à une variable d'interaction.

Merci de votre aide.

[Haache]

Bonjour
C'est normal que le contraire ne donne pas 2.
Pour donner ma réponse je veux m'assurer d'une chose. Quand vous changez de référence pour statut, le odd ratio donne 1.84 ou 1.83 environ et non 2. Est - ce bien ça ?

[xavier_dcf]

Bonjour,

Merci pour votre réponse.

Non , cela ne donne pas 2 . Il y a un soucis dans la cohérence des calculs ou dans ma proc logistic.
Car dans mon message précédent, j'ai bien pris comme référence statut= "ex dgf" et dans les odds ratios SAS me sort des résultats pour le cas suivant "statut ex-DGF vs ex-OQN dans esatis=10" , (cela devrait être l'inverse)
Si je mets comme référence statut="ex OQN" , SAS me sort encore des résultats pour le même cas suivant "statut ex-DGF vs ex-OQN dans esatis=10" avec des chiffres différents ce qui n'est pas logique.

Le problème vient peut être de ma proc logistic :
ODS output globaltests= ods;
proc logistic data=base_ifaq;
class statut (ref="ex-DGF") esatis(ref="NA")/param=ref;
model remuneration_ifaq(event='1')=STATUT*esatis;
oddsratio statut;
run;
ODS output close;

Il faut peut être mettre les variables explicatives seules hors interaction+ la variable interaction dans le modèle pour que la proc fonctionne.
Je teste.

[xavier_dcf]

rebonjour,


En effectuant les codes suivants (c'est à dire je mets les variables hors interaction dans le modèle), j'obtiens des résultats cohérents mais exactement les mêmes résultats. j'ai dans les tables ods et ods2 "statut ex-DGF vs ex-OQN dans esatis=10" avec les mêmes chiffres. Les chiffres sont les mêmes mais je ne comprends pas pourquoi la référence ne change pas.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
 
	ODS output   OddsRatiosWald= ods;	
			proc logistic data=base_ifaq;	
				class statut (ref="ex-DGF") esatis/param=ref;
				model remuneration_ifaq(event='1')=statut esatis STATUT*esatis;
                oddsratio statut;	
			run;	
	ODS output close;	
 
 
	ODS output   OddsRatiosWald= ods2;	
			proc logistic data=base_ifaq;	
				class statut (ref="ex-OQN") esatis/param=ref;
				model remuneration_ifaq(event='1')=statut esatis STATUT*esatis;
                oddsratio statut;	
			run;	
	ODS output close;

Le problème est à moitié résolu, je ne sais pas si je peux prendre les odds ratios obtenus comme des bons résultats. Qqun a une idée?

Merci

[Haache]

Bonjour
Si ça ne donne pas 2 je ne pense pas que le résultat soit incohérent. Ici il s'agit d'une interaction qui est la variable considérée dans le modèle. Au plan théorique cette variables a 2*4=8 modalités. Dans la résolution du modèle, une des 8 modalités sera prise comme référence (c'est à dire le produit des modalités considérées comme référence pour chacune des deux modèles). Du moment où la variable n'est plus binaire, le faire de changer la modalité de référence ne veut pas dire que les odd ratio seront inversés car il y a plus de deux modalités. C'est vrai que la variable pour laquelle vous avez changer la modalité est binaire mais ce n'est pas cette variable qui est utilisée dans le modèle.

[xavier_dcf]

oui merci.

Mais pourquoi dans ma sortie sur les odss ratios , j'ai alors que 4 lignes?

[Haache]

parce que vous avez demander seulement pour Statut