Le test de sphéricité ne s'affiche pas dans les résultats du modèleJ'ai deux questions sur l'exécution d'un test de sphéricité pour une ANOVA à mesures répétées si les gens sont prêts à m'écouter ?
Voici une description de mon expérience et de mes données : mon expérience consistait à tester différents médias filtrants pour éliminer le phosphore de l'eau. J'ai un total de 12 cases (mes "sujets") qui représentent 4 médias filtrants en triple exemplaire. Dans ma table de données JMP, chaque ligne représente une case, et pour mes colonnes, j'ai une colonne "ID" (essentiellement 1-12 pour mes cases de filtre), une colonne qui est le média filtrant appelé "Treatment" (étiqueté "BF" , « NF », « NS » et « WS » et répétés 3 fois pour chaque répétition), et les colonnes suivantes indiquent les concentrations de phosphore pour chaque point dans le temps (total de 11 points dans le temps, espacés de manière inégale). J'ai également inclus mon fichier de données avec ce message.
Alors, mes questions :
1.) J'ai effectué un test de modèle mixte en tant qu'ANOVA à mesures répétées sur mes données, mais un article de JMP indique que je dois tester la sphéricité pour voir si je dois ajuster ma valeur F au cas où je violerais le test de sphéricité. Pour tester la sphéricité, je dois effectuer le test multivarié. Cependant, j'ai des données manquantes pour certains moments, et JMP jette le sujet entier si des données manquent avec l'approche multivariée, c'est pourquoi j'ai opté pour l'approche de modèle mixte (je n'ai que 12 sujets). Comment puis-je alors effectuer un test de sphéricité, qui semble nécessiter l'approche multivariée ?
2.) J'essaie également d'effectuer une ANOVA à mesures répétées en utilisant l'approche multivariée sur un sous-ensemble de mes données juste au moins à ce sous-ensemble pour la sphéricité. J'ai essayé d'exécuter un modèle MANOVA avec mes données, mais le test de sphéricité n'apparaissait pas dans mes résultats. Si je réduis le nombre de points temporels que j'ai (j'ai 11 points temporels), qui sont mes valeurs "y", le test de sphéricité apparaît, mais, selon le nombre de points temporels que j'inclus, le critère de Mauchly est significatif ou non significatif. Je ne suis pas très familier avec les statistiques en général, donc je me demandais si quelqu'un sait ce qui se passe ?
Quelqu'un a-t-il une idée de la raison pour laquelle le test de sphéricité se comporte comme je l'ai décrit ? Existe-t-il d'autres options pour tester la sphéricité qui ne nécessitent pas une approche multivariée ?
Merci beaucoup pour toute contribution ! Je ne suis pas fort en stats du tout, donc toute aide serait appréciée !
Stéphanie
Re : Le test de sphéricité ne s'affiche pas dans les résultats du modèleJ'ai un problème similaire avec le test de sphéricité qui ne s'affiche pas dans JMP. J'ai l'impression que c'est parce que mes données ne disent pas grand-chose, ou à cause de problèmes avec la taille de mes échantillons, mais je suis malheureusement trop vert pour en être sûr. Puis-je cocher la case de corrélation partielle pour voir que l'hypothèse de sphéricité est de toute façon violée (sans le test) ? J'ai joint mon jeu de données. Sera reconnaissant pour toute aide.
Re : Le test de sphéricité ne s'affiche pas dans les résultats du modèleMerci beaucoup pour votre réponse détaillée! J'apprécie vraiment que vous preniez le temps d'expliquer les concepts et de fournir des exemples concrets basés sur mes données !
Stéphanie
Re : Le test de sphéricité ne s'affiche pas dans les résultats du modèleStéphanie,
J'ai vu des zéros dans vos données, mais aucune valeur manquante. L'erreur que j'ai reçue était "degrés de liberté d'erreur insuffisants".Avec 2 df pour chaque groupe de traitement qui vous donne 8 df pour chaque période (répéter).Je m'attendais donc à ce que seulement 4 périodes de temps soient autorisées pour le test. Pour être juste ou pour montrer ce qui est vu dans les données, utilisez 4, 12, 72 et 120, cela montrera le problème de sphéricité.
En fait, votre jeu de données est assez intéressant.Je ne pense pas que vous ayez besoin d'un test statistique pour constater une différence importante dans les résultats des matériaux. Les données montrent un modèle de croissance.
À l'aide du vaste ensemble de données que vous aviez joint, créez la matrice de covariance à l'aide de l'analyse multivariée et sélectionnez Covariance. Vous verrez que la structure de covariance (et la corrélation) sont assez différentes pour les premières périodes par rapport aux périodes ultérieures où les valeurs mesurées ne changent pas beaucoup.Encore une fois, vous devriez regarder les données et la covariance, vous pouvez voir que la règle de sphéricité est violée sans test.À mon avis, une valeur p n'a aucune valeur ici.
La prochaine étape consiste donc à choisir un test différent. Vous devez au moins modéliser un composant de série chronologique pour vos mesures répétées.(Voir suggestion ci-dessous.)
Measured ValuesMeasured ValuesValeurs mesurées
Weibull Growth ModelsWeibull Growth Models Modèles de croissance de Weibull
Predicted vs. ActualPredicted vs. Actual Prédit vs Réel
Si vous voulez quantifier la différence, je vous suggère de modéliser la croissance pour chaque groupe de traitement.Les images ci-dessus affichent les valeurs mesurées, les résultats des courbes d'ajustement utilisant une courbe de croissance de Weibull et la performance. Notez que le modèle Weibull Growth n'est pas dans JMP13, mais est disponible dans JMP14. Les valeurs prédites ont été enregistrées dans le tableau de données et tracées par rapport aux valeurs réelles avec la ligne y=x.Les scripts pour les deux sont enregistrés dans la version empilée de votre table (jointe).Le script Fit Curves peut ne pas fonctionner, selon votre version de JMP.
Vous devez rechercher l'instabilité du traitement BF. Y a-t-il quelque chose dans le matériau, l'environnement comme l'emplacement de la boîte ou la température qui pourrait affecter la précision de vos lectures pour les concentrations de phosphore.
