GEE en SAS et R: résultats différents

Bonjour

En utilisant Proc GENMOD pour SAS et geeglm pour R dans le cadre de données répétées, j'obtiens des résultats différents si j'utilise une matrice AR1 (les résultats sont les mêmes avec la matrice exchangeable) et je me demande pourquoi et quel résultat est le bon.

Voici mes syntaxes:

geeglm(y ~ x , id=sujet,
waves=jour, data = pol, family = binomial, corstr = "ar1", na.action=na.omit)

proc genmod data=pol descend;
class num_pat ;
model y = x / dist=binomial link=logit ;
repeated subject=sujet / corr = AR WITHIN=jour ;
run;

je précise que y est une variable binaire, x est une variable quantitative, jour (1,2,3, etc.) correspond aux jours des mesures répétées chez un même sujet, et qu'il manque des jours pour certains sujets.

Merci de votre aide

Attention avec SAS qui se base sur des formules américaines.
Je connais pas la PROC que tu utilises ni les algorithmes qu'elle fait intervenir mais par exemple si tu prends l'analyse discriminante l'aglorithme qui fait intervenir les matrices de covariances en les divisant par n dans la littérature or SAS les divise par n-1 ou encore n'utilise pas la matrice de la formule d'Huygens mais celle de la dispersion inter-groupe pour calculer les composantes canoniques discriminantes.

Ensuite voir également comment SAS et R gèrent les données manquantes que tu as, as-tu consulté la doc du support SAS? tu devrais avoir une section consacrée à la gestion des données manquantes par cette procédure.