Comment puis-je déterminer la différence dans les médianes et l'intervalle de confiance pour cette différenceOn m'a demandé de rapporter la différence et les intervalles de confiance à 95 % autour de la différence pour les moyennes et les médianes dans une étude. Il est facile d'obtenir ces valeurs pour les moyennes, mais j'ai eu du mal à comprendre comment procéder pour les médianes et les données non paramétriques.
Je travaille avec deux groupes de sujets, des données entières non appariées.
Toute aide est appréciée.J'utilise JMP 12.
Re : Comment puis-je déterminer la différence entre les médianes et l'intervalle de confiance pour cette différence ?C'est merveilleux - merci beaucoup !
Re : Comment puis-je déterminer la différence entre les médianes et l'intervalle de confiance pour cette différence ?Si vous ne disposez pas de JMP Pro, vous pouvez produire le bootstrap vous-même, mais cela nécessite les étapes supplémentaires d'échantillonnage de votre table avec remplacement afin que vous puissiez calculer la médiane pour chaque groupe à chaque itération. Vous pouvez le faire via des scripts, mais voici comment vous pouvez procéder de manière interactive en utilisant mon Échantillon avec complément de remplacement (téléchargez ceci en premier).
- Lancez l'échantillon avec le complément de remplacement (Menu Complément > Échantillon avec remplacement > Échantillon maintenant)
- Sélectionnez les colonnes dont vous avez besoin pour la comparaison particulière des médianes (votre Y et votre facteur X)
- Spécifiez au moins 5000 pour le nombre d'échantillons, puis cliquez sur OK (laissez la taille de l'échantillon seule, ceci est uniquement pour le sur- ou le sous-échantillonnage). Vous aurez un nouveau tableau contenant autant de lignes que votre tableau d'origine multiplié par le nombre d'échantillons que vous avez dessinés. Pour une grande table initiale, cela pourrait être énorme.
- Tableaux de lancement > Résumé
- Regrouper par numéro d'échantillon
- Sous-groupe selon votre facteur
- Sélectionnez votre Y, puis cliquez sur "Statistiques" et sélectionnez "Médiane". Cliquez sur OK. Vous aurez maintenant le même tableau récapitulatif que nous avons obtenu du bootstrap du clic droit, donc les étapes sont maintenant les mêmes
- Sélectionnez les deux colonnes médianes, faites un clic droit dans l'en-tête > Nouvelle colonne de formule > Combiner > Différence
- Lancez Analyser> Distribution, placez la colonne de différence dans pour Y, puis cliquez sur OK
- Il n'y aura pas de section de limites Bootstrap, mais vous pouvez utiliser les quantiles pour trouver vos limites de couverture. Pour une confiance de 95 %, lisez les valeurs des quantiles 2,5 % et 97,5 %.
Voici une autre vidéo montrant cela rapidement:
J'espère que cela t'aides!
@julien
Re : Comment puis-je déterminer la différence entre les médianes et l'intervalle de confiance pour cette différence ?Salut @strickerp ,
La détermination de l'erreur standard de la médiane en tant qu'estimateur présente des complications, de sorte que les intervalles de confiance exacts ne sont généralement pas définis (bien qu'il existe des formules approximatives qui fonctionnent bien pour les grands échantillons). La meilleure façon de former un intervalle de confiance pour une telle statistique est de l'amorcer. Vous utilisez JMP Pro 12 par hasard ? Si c'est le cas, vous pouvez générer des intervalles de confiance amorcés concernant la différence médiane à l'aide de l'amorçage par clic droit de JMP Pro :
Utilisation de JMP Pro :
- Lancez Fit Y by X avec votre variable
- Allez dans le Triangle Rouge > Quantiles.
- Faites un clic droit sur les médianes rapportées et sélectionnez "Bootstrap". Entrez le nombre d'échantillons bootstrap que vous souhaitez (je suggérerais 5000, ce qui prendra un certain temps). Vous obtiendrez un tableau avec une ligne par échantillon bootstrap et des colonnes pour la médiane trouvée à chaque itération pour chaque niveau de votre variable X.
- Puisque vous souhaitez former un intervalle de confiance sur la différence des médianes, sélectionnez les deux colonnes dans le tableau, cliquez avec le bouton droit sur l'en-tête de ligne > Nouvelle colonne de formule > Combiner > Différence.
- Pour obtenir votre intervalle de confiance bootstrap, accédez à Analyser> Distribution, placez votre nouvelle colonne de score de différence dans pour Y et cliquez sur OK.
- JMP sait qu'il s'agit d'une table bootstrap, il renverra donc une section intitulée "Bootstrap Confidence Limits". Pour des intervalles de 95 %, regardez les lignes de couverture de 0,95.
Voici une vidéo rapide montrant ces étapes :
J'espère que ça aide!
@julien