Pourquoi est-ce que j'obtiens des résultats différents d'analyse de composant de variance entre SAS et JMP ??Pourquoi est-ce que j'obtiens des résultats différents entre SAS et JMP pour l'analyse des composantes de la variance ?
Cas 1 - Croisé puis croisé
proc varcomp data = data method=TYPE1;
classe Facteur1 Facteur2 Facteur3 ;
variable de modèle = Facteur1 Facteur2 Facteur3 ;
Cours;
= JMP - Qualité et processus - Variabilité / Graphique de jauge d'attribut - Type de modèle - Franchi
(Si "Main Effect" est sélectionné, le résultat est le même avec SAS.)
Cas 2 - Imbriqué puis croisé
proc varcomp data = data method=TYPE1;
classe Facteur1 Facteur2 Facteur3 ;
variable de modèle = Facteur1 Facteur2(Facteur1) Facteur3 ;
Cours;
= JMP - Qualité et processus - Variabilité / Graphique de jauge d'attribut - Type de modèle - Imbriqué puis croisé
Cas 3 - Croisé puis Imbriqué
proc varcomp data = data method=TYPE1;
classe Facteur1 Facteur2 Facteur3 ;
variable de modèle = Facteur1 Facteur2 Facteur3(Facteur2) ;
Cours;
= JMP - Qualité et processus - Variabilité / Graphique de jauge d'attribut - Type de modèle - Croisé puis Imbriqué
Il a été confirmé que les résultats SAS et JMP de Variance Component sont les mêmes lorsque le nombre de facteurs est 1 ou 2.
Cependant, les résultats sont différents pour les cas ci-dessus avec 3 facteurs.\
J'apprécierais vraiment ... si vous pouviez me dire pourquoi les résultats des deux programmes sont différents !
AnalytiqueTableau de donnéesBICHEmacOSAutreles fenêtresRe : Pourquoi est-ce que j'obtiens des résultats différents d'analyse de composant de variance entre SAS et JMP ??Sans voir vos données, je ne peux pas savoir pourquoi elles sont différentes, mais je peux faire une supposition éclairée.
Dans le code SAS, vous utilisez Method=Type I. Ceci utilise des sommes de carrés de type I pour l'estimation. Dans la boîte de dialogue Composants de variance JMP, il y a un bouton pour les paramètres d'analyse. Si vous cliquez dessus, vous verrez que la valeur par défaut de JMP est "Choisir la meilleure analyse". JMP pouvait donc librement basculer entre EMS, REML ou bayésien. JMP utilise donc probablement une méthode d'estimation différente de celle de SAS. Notez que JMP ne propose même pas d'option pour forcer la méthode de type I (ou EMS) car il s'agit d'une approche plus ancienne qui n'est généralement pas recommandée.
Si vous avez besoin qu'ils soient d'accord, dans SAS, définissez la méthode sur REML. Dans JMP, choisissez l'analyse Utiliser REML. Cela devrait les mettre d'accord. Mais gardez à l'esprit que même dans cette situation, en fonction de vos données, il peut y avoir de légères différences en raison de l'approche de la façon dont les grandes matrices peuvent être traitées (séparées et recombinées) à des fins d'efficacité ainsi que la nature itérative de la résolution de telles équations .
Re : Pourquoi est-ce que j'obtiens des résultats différents d'analyse de composant de variance entre SAS et JMP ??Cher Dan,
Merci beaucoup pour votre aimable réponse.
Comme vous l'avez dit, j'ai les mêmes résultats en utilisant REML.
J'ai compris que le même résultat ne peut pas être obtenu dans le cas du non-REML car la méthode d'estimation de la méthode d'erreur de type1 et la méthode EMS sont différentes.
Vous avez dit que la méthode Type1 est une ancienne méthode.
L'utilisation de REML peut-elle obtenir des résultats ANOVA plus précis même si la variacne n'est pas négative ???
Re : Pourquoi est-ce que j'obtiens des résultats différents d'analyse de composant de variance entre SAS et JMP ??Qu'entend-on exactement par précis? L'approche de type I ou EMS et le REML ne sont que des méthodes d'estimation différentes. REML est une méthode d'estimation très bonne et bien respectée qui existe depuis des années. En fait, pour la plupart des situations avec des données équilibrées, il fournira les mêmes résultats que EMS. REML a l'avantage de travailler dans des situations déséquilibrées et dans des situations avec des composantes de variance négatives. De plus, REML vous permettrait même d'exiger que les composantes de la variance soient supérieures ou égales à zéro pour éviter les estimations négatives.
Vous pouvez vous sentir à l'aise avec l'approche REML.
Re : Pourquoi est-ce que j'obtiens des résultats différents d'analyse de composant de variance entre SAS et JMP ??Merci pour votre réponse vraiment gentille et précise !!
