<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />
Est-ce que quelqu'un connaît un script JMP pour effectuer la procédure Hodges-Lehmann sur des données? Certains documents de normalisation européens demandent que l'analyse de Hodges-Lehmann soit effectuée sur des types de données spécifiques. Les documents montrent comment faire Hodges-Lehmann à la main. Je peux obtenir la sortie désirée en utilisant le test Wilcoxon à 1 échantillon dans Minitab, même si cela nécessite de nombreuses explications sur les parties de la sortie Minitab correspondant aux informations HL souhaitées. Ces documents recommandent spécifiquement à StatXact ou à SAS de faire Hodges-Lehmann. Je n'ai pas de SAS, mais JMP. J'espérais trouver un script JMP pour faire Hodges-Lehmann. Le test que j'effectue est un test unilatéral selon lequel la médiane d'un ensemble de données est une valeur spécifique. Et je dois utiliser Hodges-Lehmann.
Je vous remercie!
Bruce
La méthode manuelle que j'ai trouvée provient d'un document de normes européennes. Ici est un lien vers le document. Je ne suis pas sûr des informations de copyright, ou je posterais une image de la page appropriée ou deux. J'ai pu trouver un algorithme généralisé ici . J'ai aussi trouvé quelques étapes SAS dans ce document et ce un.
Je pense pouvoir citer un paragraphe du document européen. Une fois la médiane trouvée, la procédure ci-dessous est décrite pour un ensemble de données de taille n = 20.
"Les différences moyennes par paires qui ne dépassent pas la médiane ... sont calculés. Le tableau E.5 des valeurs critiques pour les rangs signés par paires appariées de Wilcoxon teste l'entrée pour n = 20 et un niveau de signification unilatéral de 0,025 indique la valeur critique de 52. D'où c = 52 + 1 = 53. Les différences par paires sont triées par ordre décroissant. La 53ème valeur est 0,11 [NOTE: cela provient des exemples de données dans le document]. par conséquent, la limite de confiance unilatérale supérieure de 97,5% de Hodges-Lehmann pour la différence de lg Rs entre RP et PP est de 0,11, ce qui est inférieure à la marge d'infériorité convenue de 0,6. Par conséquent, l'hypothèse d'infériorité de PP est rejetée et on peut en conclure que la préparation du test PP n'est pas inférieure à RP. "
Le document européen contient un petit tableau des valeurs critiques de Wilcoxon, mais je n'ai pas été en mesure de trouver un tableau plus complet nulle part.
Merci!
Bruce
bwwhite,
Avez-vous un lien pour faire le HL à la main? De plus, avez-vous installé R sur le même ordinateur que JMP?
