失敗したTテストのみをグラフに表示する方法こんにちは、みんな、
2つのグループ間のTテストを実行するために次のスクリプトを使用しています（ストレスアブドコントロール）。
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" )}
)}
)
)
);ご覧のとおり、私は各「テスト」に対して T テストを実行しており、何千ものテストがあります。私の質問は、T テストが不合格となったテストのみをグラフに表示するオプションがあるかどうかです。
ありがとう！！
Re: 失敗したTテストのみをグラフに表示する方法@Adir52 、可能であれば、2 つのサンプルを比較するために各グループのユニット数を同じにすることをお勧めします。
サンプル サイズは、タイプ II エラーのリスクを軽減するために使用される要素です。
」で統計的仮説検定、タイプ I エラーそれは真実の拒絶です帰無仮説(「偽陽性」所見としても知られています)、タイプ II エラー偽の帰無仮説 (「偽陰性」結果とも呼ばれます) を棄却できないことです。 [1]より簡単に言うと、タイプ I の過誤は、存在しない何かが存在すると誤って推論すること (誤った情報で通説を裏付ける) であり、タイプ II の過誤は、存在する何かが存在しないと誤って推論することです (常識に反する)誤った情報を含む一般的な信念）。 」
からの抜粋https://en.wikipedia.org/wiki/Type_I_and_type_II_errors 。 帰無仮説はコントロールとストレスなしでは異なると述べていますが、2 標本検定 (2 つのグループの分散分析) の帰無仮説は、グループ 1 とグループ 2 の平均が同等であり、等しいことになります。したがって、テストが特定の違いを検出できるように実験を設計する必要があります。 たとえば、実験で 0.5 シグマの差が検出されることを 95% の信頼したい場合は、JMP を使用して単位数を選択できます。
JMP メイン メニューから DOE を選択すると、JMP のバージョンに応じて、メニューまたはサブメニューの 1 つにサンプルサイズと検出力電卓。 JMP 14 では、DOE > Design Diagnostics を選択すると見つかります。 次に選択します2 つのサンプルの平均値。 以下はスクリーンショットです。 各テストの標準偏差は異なるため、この設計ステップは相対値を使用して実行されます。 特定標準開発1として、その後、違い検出する違いとして。 以下では 1 として指定します。 押したら続くの場合、サンプル サイズは各グループの 54 単位または 27 単位になります。 サンプルサイズを削除し、違いを 0.5 または 1/2 Std Dev にして、 を押します。続く、サンプル サイズは各グループの 210 単位または 105 単位になります。 または、「差」を削除し、「サンプル サイズ (各グループ 80 単位)」に 160 を入力すると、「差」は 0.57 として報告されます。そうでない場合、テストでは 0.57 シグマの差が 95% の確率で検出されると確信できます。
統計学者がいる場合、または近くにあなたの業界に詳しい大学の統計教授がいる場合は、話し合ってみることをお勧めします。特定の研究には、コントロールの長い履歴が含まれている可能性があり、コントロール ユニットの数は、コントロールがベースライン内にあるかどうかの単なるチェックにすぎません。しかし、さらに重要なのは、単位の選択、これを公平に比較するための実験の実行方法、およびこれらの結論がどの母集団に適用されるかということです。
Re: 失敗したTテストのみをグラフに表示する方法pValue は、統計検定が有意である確率です。有意性は歴史的に 20 分の 1 (pValue <= .05) または 20 分の 2 (pValue <= .1) として定義されています。そしてはい、あなたの場合、Ho は Group1=Group2 であるため、pValue が > .05 または > .1 である場所を見つける必要があります。
Re: 失敗したTテストのみをグラフに表示する方法こんにちは@gzmorgan0
コントロールユニットのサンプルサイズに関する推奨事項はありますか。
並べられたユニット グループは 80 ユニットになります。
Re: 失敗したTテストのみをグラフに表示する方法こんにちは、テクネルソン
まず最初に、ご意見いただきありがとうございました。非常に役に立ちました。
質問が間違っていました。グループが統計的に等しくないテストを表示したいと考えています。コードを使用して次の行を削除するのは簡単です。
dt << 行選択を反転;
PValueの意味を教えていただけますか？
それは=確率>|t|ですか?
私の H0 が Group1=Group2 であるため、この質問をしています。Re: 失敗したTテストのみをグラフに表示する方法@txnelson (ジム)、とても面白いですね。私も、Semiconductor Capability.jmp で作成したテーブルを使用して返信を開始しました。このテーブルは、ウェーハ概要のロット概要から作成され、グループ A またはグループ B に入るロットをランダムに選択しました。 返信拝見しましたので詳細は書きません。 添付されているのは、100 個のテストごと、またはテストの論理グループごとに、以下のグラフを作成するためのデータ テーブルとスクリプトです。
実際には、F 検定と p 値を追加して標準偏差を比較し、n (各グループでテストされたユニットの数) を比較することをお勧めします。 以来@Adir52のグループは対照群とストレス群であり、ストレスを与えると測定不能な（空虚な）結果が生じる可能性がありました。また、おそらく半導体テストのように、何千ものテストがある場合、一部のテストは相関している可能性があります。たとえば、すべての PNP テストが上方にシフトしましたが、すべてが有意だったわけではありません。 これは、0 より大きい t 比として表示されます。
統計学者として、私は統計的有意性について警告しなければなりません。 このランダムにグループ化されたデータ、118 個の t 検定 (128 個のテストのうち 10 個は stdev がゼロで、t 比がありません)、11 個には予想される誤り率である < 0.10 の p 値でフラグが付けられました (および 5 個のテストでは p が設定されました) ことに注意してください。 -値 < 0.05)。したがって、p 値のみと平均値の検定のみに基づいて判断を行う場合には、多くの注意点があります。
最後の注意点@Adir52 、測定値が自動テスト装置によって報告される場合、順序が重要になる可能性があります。時間要素がある場合は、テストごとにこのような表示を推奨します。
Re: 失敗したTテストのみをグラフに表示する方法JMP ではこの問題を解決する方法が複数あります。私は、応答スクリーニングを使用してデータを非常に高速に実行して分析の p 値だけを取得し、その結果を取得して分析に必要なデータを選択することを好みます。以下は、特定のデータ テーブルの列名と一致するように作成したサンプル データに対してこの処理を実行する簡単なスクリプトです。
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" )}
)}
)
)
);
