Calcul de la taille de l'échantillon pour TOSTSalut tout le monde!
Contexte: J'ai un processus A qui est ma référence, et j'ai développé un processus B dont je souhaite évaluer la comparabilité avec le processus A. Pour ce faire, je prévois d'utiliser un TOST (test d'équivalence) en utilisant le processus A's 3xSTD (écart type) comme seuil. Maintenant, je veux savoir combien de fois dois-je répéter le processus B pour être sûr de pouvoir détecter une telle variation (sachant qu'une course représente une somme d'argent importante).
Option 1 : Je voulais utiliser l'outil Sample size & Power, avec 2 comparaisons de moyennes d'échantillons, avec un alpha de 0,1 (2x0,05 pour chaque processus), un STD = 1, pas de paramètre supplémentaire, Différence à détecter = 3 (3x STD de 1), Puissance = 0,95, & j'obtiens 7, donc 4 exécutions du processus B à comparer à 4 exécutions du processus A pour assurer la comparabilité.
Option 2 : Pour utiliser la distribution t (comme indiqué ci-dessous) et donc suggérer 6 exécutions :
Que pensez-vous des deux options, laquelle est la meilleure pour ce cas? Avez-vous un autre outil auquel je n'aurais pas pensé ?
Merci beaucoup ! 
@martindemel peut être ?
Modélisation statistique avancéeAnalyse et modélisation de données de baseExploration et visualisation des donnéesPersonnalisation de masseModélisation prédictive et apprentissage automatiqueIngénierie de la qualité et des procédésRe : Calcul de la taille de l'échantillon pour TOSTOui, c'est dans JMP 16. Utilisez le menu DOE --> Sample Size Explorers --> Power pour trouver des tests d'équivalence pour les cas à 1 et 2 échantillons.
Re : Calcul de la taille de l'échantillon pour TOSTChère équipe JMP,
cette fonction a-t-elle été ajoutée dans JMP16 ? Si oui, où ? Je n'ai pas pu le trouver.
Merci d'avance
Anne
Re : Calcul de la taille de l'échantillon pour TOSTSalut John,
Veuillez contacter votre représentant commercial JMP SE/JMP pour vous envoyer une invitation au programme d'adoption précoce (EA) pour JMP 16 (tant que vous êtes un utilisateur et que vous êtes prêt à faire part de vos commentaires). Dans l'EA, notre développement partage de nouvelles fonctionnalités potentielles et demande des commentaires avant qu'il ne soit livré dans la nouvelle version. C'est une excellente occasion de fournir des commentaires sur la "nouvelle" façon de faire des calculs de puissance et des tests d'équivalence, ainsi que s'il y a encore quelque chose de différent de ce à quoi vous vous attendez. Et du milieu à la fin du mois de mai, la nouvelle plate-forme de puissance et de taille d'échantillon sera probablement présentée dans l'EA, y compris les tests d'équivalence pour autant que je sache. Vos remarques sont très appréciées.
Meilleur,
Martin
Re : Calcul de la taille de l'échantillon pour TOST@calking - Y a-t-il une chance que des tests d'équivalence soient ajoutés ?
Actuellement, cela peut être fait dans un logiciel similaire, tel que Minitab -- https://support.minitab.com/en-us/minitab/18/help-and-how-to/statistics/equivalence-tests/supporting-topics/power-for-equivalence-tests/
*canards*
Re : Calcul de la taille de l'échantillon pour TOSTCertainement!
Pour répondre à vos questions de suivi, les valeurs que j'ai choisies pour l'alpha et la puissance à utiliser dans le calculateur de puissance étaient en fait basées sur les valeurs de double alpha et de puissance que vous recherchiez. Pour clarifier, étant donné que JMP ne dispose pas actuellement d'un calculateur pour un test d'équivalence, nous utilisons une astuce (voir l'article) pour approximer le calcul de puissance à l'aide du calculateur de puissance du test t à deux échantillons. L'astuce implique le fait que le test d'équivalence a changé les hypothèses que vous utilisez généralement dans un test t à deux échantillons, nous devrons donc changer certaines valeurs. Donc pour votre cas, nous avons procédé comme suit :
1. Vous voulez un alpha de 0,1 (2x0,05). Pour l'astuce, cela signifie que vous réglerez la puissance de la calculatrice sur 0,9 ; c'est-à-dire 1-(alpha destiné au test d'équivalence).
2. Vous voulez une puissance de 0,95. Pour l'astuce, cela signifie que vous réglerez le niveau alpha de la calculatrice sur 0,05 ; c'est-à-dire 1-(puissance destinée au test d'équivalence).
Notez que nous avons dû changer les valeurs que nous avons utilisées pour la puissance et l'alpha. J'espère que cela aide à clarifier d'où viennent ces valeurs que j'ai utilisées.
Je me suis assuré d'utiliser la valeur alpha totale plutôt que les valeurs individuelles pour correspondre spécifiquement à vos valeurs particulières. Cependant, si 0,05 est effectivement le niveau de confiance global que vous souhaitez atteindre, alors c'est probablement la valeur que vous devriez utiliser. Une valeur de 0,1 se traduira par une taille d'échantillon totale plus petite, mais uniquement parce que vous assouplissez la restriction sur la confiance globale (confiance à 90 % contre confiance à 95 %). De plus, le test t à deux échantillons tient déjà compte de l'existence de deux populations et s'ajuste donc déjà pour garantir le maintien du niveau de confiance global.
J'ai couru quelques chiffres rapides et, dans tous les cas (alpha = 0,1 contre alpha = 0,05), la différence globale dans la taille totale de l'échantillon n'est que de 8 contre 9. Donc, en fin de compte, cela dépendra de la façon dont votre direction décide de peser les risques et les rendements (8 échantillons au total pour une confiance globale de 90 % contre 9 échantillons au total pour une confiance globale de 95 %).
Re : Calcul de la taille de l'échantillon pour TOSTD'accord super merci beaucoup !
Je n'ai pas R donc je vais essayer ça avec le script !Re : Calcul de la taille de l'échantillon pour TOSTSalut @calking ,
Ça sonne bien, merci ! Quelques questions :
1. À propos de l'alpha, on nous a dit que puisque nous avons 2 populations, nous devons appliquer l'alpha de 1-0,95 aux deux, donc pour utiliser 0,05x2=0,1 dans l'alpha global de l'outil, qu'en pensez-vous ?
2. A propos de la puissance, est-ce habituel d'utiliser la 0.9 ?
Merci beaucoup pour votre aide!
Re : Calcul de la taille de l'échantillon pour TOSTPour le cas à un échantillon, vous spécifiez les valeurs de ces variables :
n = 15; //total sample size
conf = 0.95; //confidence
mu = 505; //hypothesized mean
muu = 510; //upper equivalence limit
mul = 490; //lower equivalence limit
sigma = 4; //estimated standard deviation
Pour le cas à deux échantillons, vous spécifiez les valeurs de ces variables :
N = 36; //total sample size
W = 1; //weighting for sample size (e.g. use 2 for 1 of the groups to have 2x as many samples as the other
conf = 0.95; //confidence level
mudiff_l = -2; //lower equivalence limit on mean difference
mudiff_u = 2; //upper equivalence limit on mean difference
mudiff = 0; //estimated mean difference
sigma = 2; //estimated standard deviation
Pour trouver une taille d'échantillon, il vous suffirait d'implémenter une boucle pour augmenter N jusqu'à ce que la puissance souhaitée soit atteinte.
Je comprends votre inquiétude quant à la nécessité de soumettre cela à une autorité pour examen.Je travaille moi-même dans une entreprise réglementée et nous avons un protocole de validation pour le complément JMP avec les calculatrices que nous avons construites sur la base de ce script. La documentation SAS pour PROC POWER montre les formules pour le calcul de puissance pour 1 et 2 échantillons TOST.Vous pouvez faire référence à ce document pour montrer que les calculs sont exacts, le reproduire dans un autre logiciel (par exemple R, j'ai un script pour cela si vous le souhaitez) ou confirmer les estimations de puissance par simulation.Si vous avez SAS n'importe où dans votre organisation, vous pouvez simplement obtenir une impression d'alimentation proc et avoir terminé.
Re : Calcul de la taille de l'échantillon pour TOSTMerci pour cette réponse, je ne suis pas sûr d'avoir compris tout le script, qu'est-ce que je suis censé changer ?
Mon souci est que cette expérience est "officielle" (sera soumise aux autorités), donc je ne suis pas sûr d'utiliser un tel script ...Re : Calcul de la taille de l'échantillon pour TOSTNous avons développé une application qui effectue des calculs de taille d'échantillon TOST.Je ne peux pas envoyer cela, mais je peux partager le script que j'ai utilisé pour répliquer les calculs de puissance pour TOST dans SAS en utilisant le Q d'Owen.Il contient le code des scénarios à 1 et 2 côtés avec un exemple pour chacun.Notez que cela ne gérera pas bien les petits scénarios signal sur bruit.D'après mon expérience avec cela jusqu'à présent, cela ne fonctionne que jusqu'à ce que nu (paramètre de taille d'échantillon) atteigne environ 200-250 car Gamma( nu/2 ) rencontre un problème de débordement numérique et renvoie une valeur manquante.
//Owen's Q
OwensQ = Function( {t, delta, a, b, nu},
c = Sqrt( 2 * Pi() ) / (Gamma( nu / 2 ) * Power( 2, (nu - 2) / 2 ));
integrand = Expr(
Normal Distribution( ((t * x / Sqrt( nu )) - delta) ) * Power( x, nu - 1 ) * Normal Density( x )
);
c * Integrate( integrand, x, a, b, <<STARTING value=""></STARTING>
Re : Calcul de la taille de l'échantillon pour TOSTBonjour Elofar,
j'ai trouvé un article ce qui semblait indiquer que vous pouviez approximer le calcul de puissance en utilisant un test t où alpha_test = 1- power_equiv, power_test = 1-alpha_equiv et la différence à détecter correspond à la valeur seuil du test d'équivalence. À l'aide du calculateur de puissance et de taille d'échantillon pour deux moyennes d'échantillons, vous spécifiez alpha=0,05 (1-0,95), StdDev=1, Différence à détecter = 3 et Puissance=0,9 (1-0,1). Cela vous donne une taille d'échantillon totale de 8, ce qui semble être en accord avec votre option 1.
J'espère que cela pourra aider!
