カスタムテストを使用して、ポストホックのペアワイズスロープ比較をジュリーリグしようとしています。 カスタム テストに関するドキュメントを読みましたが、テストの設定に使用する係数をどのように決定すればよいのかまだ理解できません。 3 つの治療法 [PT CS OF] があり、これらの治療法と連続変数 [期間] の間の相互作用をペアごとに検定したいと考えています。 カスタム テストを正しく理解していれば、各テストは独自の列に入力されます。
テストしたい:
PT*期間-CS*期間=0;
PT*期間-OF*期間=0;
CS*期間-OF*期間=0;
以下に、テスト係数と結果の画面キャプチャを添付しました。 何か間違ったことをしていると思われる警告が表示されます。
よろしくお願いいたします。
自由度の問題ではないかと思いますが、MANOVA 手順のテストについてはよくわかりません。 もしかしたら他のメンバーがアドバイスしてくれるかも知れません。 ヘルプが必要な場合は、JMP または SAS ソフトウェアのドキュメントを参照することをお勧めします。 JMP テクニカル サポート ( support@jmp.com ）あなたの質問と一緒に。
さて、カスタムテストとMANOVAの質問です。
ANOVA で交互作用の傾きに対して 3 列のペアごとの検定を行うと、すべての検定統計量と、対比ができないという警告が表示されますが、少なくともすべての推定値と検定統計量は得られます。 MANOVAでも同じことを試してみました。 因子とその相互作用を 4 つの応答変数すべてに対して行う場合と同様に、モデルを当てはめます。 (高次の有意でない相互作用の一部を削除しました) 応答オプションとして「各列もテストする」と「アイデンティティ」を選択しました。 特定の列について、「見出しの合計変更」とします。ANOVA と同じ方法でカスタム テストを試みましたが、結果は 1 つの F 統計と、コントラストをテストできないという同じ警告でした。 単一のペアごとのテストを実行すると、期待どおりの結果が得られます。 関心のある列テストのそれぞれで 10 億回のペアワイズ テストを実行するだけで、ボンフェローニは正しいでしょうか?
あなたは正しい道を進んでいると思います！
効果試験を見ましたか？ データを分析してモデル化する方法は複数ありますが、私はモデル全体から始めて、効果テストとその他の情報を進めてモデルを縮小し、次にカスタム テストを使用して重要な効果内で興味深い比較を行うことを好みます。
JMP サイトには役立つ資料がたくさんあります。 私は利用可能なものすべてについて知っているふりをするつもりはありません。 ただし、検索可能です。 たとえば、次のことがわかりました。彼のページ学術ユーザー向けのチュートリアルとビデオが含まれています (ただし、誰でも使用できます)。
明確にしておきたいのですが、私はすべての独立因子と重要な交互作用項を単一のモデルに含めています。 同じモデル内のすべての依存因子を収集する MANOVA を実行しましたが、モデル全体が重要です。 さらに、MANOVA を個々の標準最小二乗分析に分解して、依存因子の挙動を解析するように言われました。 MANOVA を使用して作品を解釈する方法を説明したビデオまたはスライドデッキへのリンクがある場合は、それをいただければ幸いです。
多重比較の問題は 1 つの研究と 1 つのデータセット内にあることを理解しています。 研究を繰り返す場合、最初の研究で行われた比較を調整に含める必要はありません。 ただし、同じ研究に複数の要素があるように聞こえるため、累積テスト数に基づいて調整することになります。
要因ごとに個別に分析するのではなく、すべての要因を 1 つのモデルに同時に含めない理由はありますか?
追加の質問が 1 つあります。
同じ基本モデルをさまざまな要素に当てはめています。 そして、最終的には同じインタラクションを異なる状況でテストすることになります。
ボンフェローニはモデル内の結合テスト (つまり 3 つのレベル) 内で修正する必要がありますか、ボンフェローニは単一モデル内のすべてのテストを修正する必要がありますか、またはこれらすべてのカスタム テストをすべてのモデルにわたってプールし、1 つの大きなボンフェローニ修正を実行する必要がありますか?
ありがとう
アルファの選択に対してボンフェローニ調整を使用する場合の質問の意味を私は解釈します。 はい、これは複数の比較を使用してタイプ I エラー率を制御する 1 つの方法です。
共同テストは、一部のアナリストが高い偽陽性率を防ぐために使用しています。 多くのテストを実行する場合があります。 テスト数に応じて個別に調整されることはありません。
警告は合同テストに関するものです。 結合テストの数に対して十分な自由度がありません。
ありがとう、これが私が探している解決策であり、結果は理にかなっています。 唯一の問題は、依然として「警告: テストできないコントラスト」が表示されることです。 それは、結合テストがテスト不可能である、または対照の 1 つがテスト不可能であると言っているのでしょうか? ペアワイズコントラストを個別に実行すると、同じ応答が得られません。
混乱を招いてしまい申し訳ございません。 混乱していることはよくわからないので、最初から始めましょう。
あなたは欲しいと言うことから始めましたある事後的なカスタムテストによるペアごとの傾き比較。 したがって、各グループの傾向線が異なって見えるので、その違いをテストしたいとします。 やってみよう。
(注: テクニカル サポートの同僚から素晴らしい支援を受けました。 私の返答は彼の説明に基づいています。)
カスタム テストは 1 つのステップで実行できますが、正しい係数の展開を確認できるように、2 つのステップに分割します。
ステップ 1: 各治療の線の傾きのカスタム テストを使用します。 カスタム テストの最初の列は [ 0, 0, 0, 1, 1, 0, 0 ] である必要があります。 つまり、1 のパラメータはパルス持続時間と治療[PT]*パルス持続時間だけです。 その理由は、線の傾きが共通の傾きであるパルス持続時間に治療[PT]の差分傾きを加えたものであるためです。 それらの合計は、0 に対してテストする傾きです。 パルス持続時間と治療[CS]*パルス持続時間を加算しているため、2 番目の列は [ 0, 0, 0, 1, 0, 1, 0 ] である必要があります。 最後に、3 番目の処理の微分勾配を追加しますが、これは明示的に報告されていません。 しかし、JMP がカテゴリ因子をパラメータ化する方法により、それが他の 2 つのパラメータの合計のマイナスであることがわかります。 したがって、3 列目は [ 0, 0, 0, 1, -1, -1, 0 ] となるはずです。
繰り返しますが、必要な実際のカスタム テストにはこの手順は必要ありません。 実際のテストの説明に役立つだけです。 この時点で、個々の勾配のカスタム テストの正しい係数が決定されました。
ステップ 2: 個々の傾きをゼロではなく、一度に 2 つずつ比較します。 ステップ 1 で使用した列を一度に 2 つずつ減算するだけです。 PT と CS の傾きを比較するには、ステップ 1 の CS の係数の列を PT の係数の列から減算します。 したがって、[ 0, 0, 0, 1, 1, 0, 0 ] - [ 0, 0, 0, 1, 0, 1, 0 ] = [ 0, 0, 0, 0, 1, -1, 0 ] となります。 PT から OF をテストする列は [ 0, 0, 0, 1, 1, 0, 0 ] - [ 0, 0, 0, 1, -1, -1, 0 ] = [ 0, 0, 0, 0 、２、１、０]。 最後に、CS から OF をテストする列は [ 0, 0, 0, 1, 0, 1, 0 ] - [ 0, 0, 0, 1, -1, -1, 0 ] = [ 0, 0, 0] です。 、０、１、２、０］。
以前と同様に、パラメータ推定値の合計と差の手動計算を使用して、JMP がテストに対してレポートする値と比較することで、いつでも係数を確認できます。
まだ、正しく理解できているとは思えません。
まず、依然として非対照的な警告が表示されます (以下を参照)。
次に、CS*Duration -(-OF*Duration)=0 の 2 番目の係数を -1 にすることによって、2 つの値が加算されるようですが、これは私が望むものではありません。 最終的には -0.0447 になります。
3 番目に、PT*Duration を OF*Duration (OF*Duration-PT*Duration=0) と比較しようとしている場合、値は 0.00777 である必要があり、代わりにそれは単なる OF*Duration の値 (0.0174926) になります。 したがって、係数を 0 のままにしても機能しないようです。
最初の推定値の係数を 1 に設定し、残りは 0 のままだと思います。 結果を確認できます。 [値] は、比較する推定値間の差異をゼロの差とする必要があります。
ご説明いただきありがとうございます。 下に別のウィンドウを貼り付けました。 私の理解が正しければ、最初のコントラスト (CS*Duration-OF*Duration=0) が正しく設定されています。 負の係数を使用すると、3 番目のカテゴリ効果 (PT*期間) に対してテストできるはずだと思いました。拡張された推定値ではそれがわかりますが、カスタム テストの設定ではわかりません。 PT*Duration に対してテストするにはどうすればよいですか?
これらのカスタム テストは、他の場所ではコントラストとしても知られています。
指定が間違っているのは、最初のテストで両方のレベルに +1 を入力したためです。 係数の合計がゼロになることを確認してください。 次に、「=」の後の値が帰無仮説であることを確認します。 有意差をテストする場合は通常、ゼロを入力するのが適切ですが、状況によってはそうでない場合もあります。
