Comment afficher sur le graphique uniquement le test T ayant échouéSalut tout le monde,
J'utilise le script suivant afin d'effectuer un test T entre deux groupes (contrôle abd stressé).
Oneway(
Y( :Shift ),
X( :Group ),
by(:Test),
t Test( 1 ),
Box Plots( 1 ),
X Axis Proportional( 0 ),
Points Jittered( 1 ),
Grand Mean( 0 ),
SendToReport(
Dispatch(
{},
"Oneway Plot",
FrameBox,
{DispatchSeg(
Box Plot Seg( 1 ),
{Box Type( "Outlier" ), Line Color( "Red" )}
), DispatchSeg(
Box Plot Seg( 2 ),
{Box Type( "Outlier" ), Line Color( "Red" )}
)}
)
)
);comme vous pouvez le voir, j'effectue le test T pour chaque "Test" et j'ai des milliers de tests. ma question s'il existe une option pour afficher sur le graphique uniquement le test où le test T a échoué ?
Merci!!
Re : Comment afficher sur le graphique uniquement le test T ayant échoué@Adir52 , si possible, il est recommandé d'avoir le même nombre d'unités dans chaque groupe pour une comparaison de deux échantillons.
La taille de l’échantillon est le facteur utilisé pour réduire le risque d’erreur de type II.
" Dans tests d'hypothèses statistiques , un erreur de type I est le rejet d'un vrai hypothèse nulle (également connu sous le nom de résultat « faux positif »), alors qu'un erreur de type II est l'incapacité à rejeter une fausse hypothèse nulle (également connue sous le nom de résultat « faux négatif »). [1] Plus simplement, une erreur de type I consiste à déduire faussement l'existence de quelque chose qui n'est pas là (confirmant une croyance commune avec de fausses informations), tandis qu'une erreur de type II consiste à déduire faussement l'absence de quelque chose qui est présent (en allant à l'encontre de la croyance commune avec de fausses informations). "
Extrait https://en.wikipedia.org/wiki/Type_I_and_type_II_errors . Même si vous déclarez que votre hypothèse nulle est que le contrôle et le non stressé sont différents, l'hypothèse nulle du test à deux échantillons (ANOVA pour 2 groupes) est que les moyennes du groupe 1 et du groupe 2 sont équivalentes, égales.Vous devez donc concevoir votre expérience de manière à ce que le test ait le pouvoir de détecter une différence spécifique. Par exemple, si vous souhaitez être sûr à 95 % que votre expérience détectera une différence de 0,5 sigma, vous pouvez utiliser JMP pour vous aider à sélectionner le nombre d'unités.
Dans le menu principal de JMP, sélectionnez DOE et selon votre version de JMP, le menu ou l'un des sous-menus contient un Taille et puissance de l'échantillon calculatrice. Dans JMP 14, vous le trouverez en sélectionnant DOE > Design Diagnostics. Sélectionnez ensuite Deux exemples de moyens . Ci-dessous, une capture d'écran. Étant donné que chaque test a probablement un écart type différent, cette étape de conception sera effectuée avec des valeurs relatives. Spécifier Développement standard comme 1 puis le Différence comme la différence à détecter. Ci-dessous, il est spécifié comme 1. Si vous appuyez Continuer , la taille de l'échantillon sera de 54 ou 27 unités pour chaque groupe. Supprimez la taille de l'échantillon et spécifiez la Différence pour qu'il soit 0,5 ou la moitié de Std Dev, puis appuyez sur Continuer, la taille de l'échantillon sera de 210 ou 105 unités pour chaque groupe. Ou supprimez la différence et entrez 160 pour la taille de l'échantillon (80 unités par groupe), la différence est alors signalée comme 0,57 ou votre test sera sûr qu'une différence de 0,57 sigma sera détectée dans 95 % du temps.
Si vous disposez d'un statisticien ou d'un professeur de statistiques universitaire à proximité qui connaît votre secteur, je vous suggère d'en discuter.Votre étude spécifique peut avoir un long historique de contrôle, et le nombre d'unités de contrôle ne constitue qu'une vérification que le contrôle se situe dans les limites de référence.Mais le plus important est la sélection des unités, la manière dont l'expérience est menée pour en faire une comparaison équitable et la population à laquelle ces conclusions s'appliquent.
Re : Comment afficher sur le graphique uniquement le test T ayant échouépValue est la probabilité que votre test statistique soit significatif.La signification est historiquement définie comme 1 sur 20 (pValue <= 0,05) ou 2 sur 20 (pValue <= 0,1).Et oui, dans votre cas, puisque votre Ho est Group1=Group2, vous voudriez trouver où la pValue est > 0,05 ou > 0,1.
Re : Comment afficher sur le graphique uniquement le test T ayant échouéSalut @gzmorgan0
Avez-vous des recommandations concernant la taille de l'échantillon des unités de contrôle.
le groupe d'unités streesed sera de 80 unités.
Re : Comment afficher sur le graphique uniquement le test T ayant échouéSalut Txnelson
tout d'abord merci beaucoup pour votre contribution, elle a été très utile.
j'avais une mauvaise question. je souhaite afficher les tests où les groupes ne sont pas statistiquement égaux. c'est simple en utilisant votre code et supprimez la ligne suivante :
dt << inverser la sélection de ligne ;
pourriez-vous nous indiquer la signification de la valeur P ?
est-ce = Prob>|t|?
je pose cette question parce que mon H0 est que Group1=Group2Re : Comment afficher sur le graphique uniquement le test T ayant échoué@txnelson (Jim), C'est tellement drôle, j'ai aussi commencé une réponse en utilisant un tableau créé à partir de Semiconductor Capability.jmp, réalisé à partir de résumés de lots de résumés de plaquettes et en choisissant au hasard les lots pour être dans le groupe A ou le groupe B. J'ai vu votre réponse, je n'écrirai donc pas les détails. Vous trouverez ci-joint le tableau de données et le script permettant de produire le graphique ci-dessous, par exemple tous les 100 tests, ou un regroupement logique de tests.
En pratique, il est bon d'ajouter le test F et la valeur p pour comparer les écarts types et une comparaison sur n, le nombre d'unités testées pour chaque groupe. Depuis @Adir52 Les groupes de étaient contrôlés et stressés, le stress pouvait produire des résultats non mesurables (vides).De plus, s'il existe des milliers de tests, peut-être comme les tests de semi-conducteurs, certains tests pourraient être corrélés, par exemple tous les tests PNP augmentés, mais tous n'étaient pas significatifs. Cela serait affiché sous forme de ratios t supérieurs à 0.
En tant que statisticien, je dois mettre en garde contre la signification statistique. Notez que pour ces données regroupées de manière aléatoire, 118 tests t (128 tests, 10 avaient un écart type nul, pas de rapport t), 11 ont été signalés avec une valeur p <0,10, ce qui est le taux d'erreur attendu (et 5 tests avec p -valeur < 0,05).Il existe donc de nombreuses mises en garde concernant le fait de fonder le jugement uniquement sur les valeurs p et sur le seul test des moyennes.
Dernière mise en garde à @Adir52 , si les mesures sont rapportées par un équipement de test automatisé, l'ordre peut être important.S'il y a un élément temporel, je recommande quelque chose comme cet affichage pour chaque test.
Re : Comment afficher sur le graphique uniquement le test T ayant échouéIl existe plusieurs façons de résoudre cette question dans JMP.Je préfère utiliser Response Screening pour effectuer un parcours très rapide des données afin d'obtenir uniquement les valeurs p des analyses, puis prendre ces résultats et sélectionner les données souhaitées pour l'analyse.Vous trouverez ci-dessous un script simple qui effectue cette opération sur certains exemples de données que je crée pour correspondre aux noms de colonnes spécifiques de votre table de données.
names default to here(1);
// Modify a sample data table to mimic your t-test code
dtSemi = open("$SAMPLE_DATA/semiconductor capability.jmp", invisible);
dtSemi << select where( :Site > 2 );
dtSemi << delete rows;
dtSemi:Site << set name( "Group" );
colNames = dtSemi << get column names( continuous );
dt = dtSemi << stack( columns(colNames) ,
Source Label Column( "Test" ),
Stacked Data Column( "Shift" ),
Name( "Non-stacked columns" )(Keep( :Group ))
);
close( dtSemi, nosave );
// The above code just creates a sample data table....The code below generates
// the desired output
// Now run the Response Screening to get the failed t-tests
rs = dt << Response Screening( Y( :Shift ), X( :Group ), Grouping( :Test ) );
dtRS = rs << Get PValues;
//dtRS << show window( 0 );
report(rs) << close window;
// select all of the failed t-tests and get their test names
// We are using failue of .1 as the target.....just change the
// target value if you want .05 or something else as the failure value
dtRS << select where( :PValue <= .1 );
If( N Rows(dtRS << get selected rows ) > 0,
dtRS << delete rows;
);
failingTests = dtRS:Test << get values;
// close the dtRS data table that is no longer needed
close( dtRS, nosave );
// Now select the Tests in the original data table that are not significant
dt << select where( contains( failingTests, :Test ) );
If( N Rows( dt << get selected rows ) > 0,
dt << invert row selection;
dt << hide and exclude ;
dt << clear select;
);
// Run your code
Oneway(
Y( :Shift ),
X( :Group ),
by(:Test),
t Test( 1 ),
Box Plots( 1 ),
X Axis Proportional( 0 ),
Points Jittered( 1 ),
Grand Mean( 0 ),
SendToReport(
Dispatch(
{},
"Oneway Plot",
FrameBox,
{DispatchSeg(
Box Plot Seg( 1 ),
{Box Type( "Outlier" ), Line Color( "Red" )}
), DispatchSeg(
Box Plot Seg( 2 ),
{Box Type( "Outlier" ), Line Color( "Red" )}
)}
)
)
);
