ウィルコクソンの署名付き順位テストJMP で Wilcoxon Signed-Rank オプションを実行すると、SAS で得られるものとは大きく異なる答えが得られます。 各オプション (SAS と JMP) で使用されるコード、または結果が異なる理由を解明するために使用される式を取得する方法はありますか?
事前にご協力いただきありがとうございます。Re: ウィルコクソン符号順位検定上記の通り、
こんにちは。以前の参考文献の説明に従って分析を行うために、以下の SAS コードを設定し、同じ分析を JMP で実行しました。
JMP で同じことを実行するには、 www.personal.psu.edu/ jda188/JMP/Signedrank.doc にある James Abbey の次のノートをダウンロードし、データ セットの DIFF 列にある指示に従います。
つまり、WILCライブラリからFERTILIZER1.SAS7BDATまたはFERTILIZER3.SAS7BDATのいずれかをJMPにインポートし、「分析」→「分布」でDIFF列を分析します。DIFF列を「Y, 列」プレースホルダーにコピーし、「OK」をクリックします。次に、2番目の赤い三角をクリックし、「検定平均値」を選択します。「仮説平均値の指定」に0を入力し、「Wilcoxon符号順位検定」にチェックを入れて、「OK」をクリックします。出力結果の「検定平均値=値」セクションには、検定統計量=74、P値=0.0017が表示されます。これはSASが生成する結果と同じです。
それでいいですか？
Re: ウィルコクソン符号順位検定どちらかの設定が間違っている可能性が高いので、まずは教科書的な例題を入手し、そのデータをSASとJMPの両方の分析で実行して、どちらが期待通りの結果を生んでいないかを確認するのが最善です。分かりやすい例題として個人的に気に入っているのは、シドニー・シーゲル著の『行動科学のためのノンパラメトリック統計』ですが、他にもたくさんあります。
2つ目の可能性として、SAS検定かJMP検定のいずれかが、実際には対応のあるサンプルの符号順位検定ではなく、独立サンプルの順位和検定（Mann-Whitney U検定とも呼ばれる）である可能性が考えられます。これら2つの検定は、それぞれ対応のあるt検定と対応のないt検定のノンパラメトリックな類似物であり、どちらもWilcoxonの名で知られています。これらの検定の詳細は、以下のリンク先のWikipediaで、実例とともにご覧いただけます。
http://en.wikipedia.org/wiki/ウィルコクソン符号順位検定
http://en.wikipedia.org/wiki/マン・ホイットニー・ウィルコクソンテストRe: ウィルコクソン符号順位検定以下に、これに一致する UCLA Stat Academic Tech Services のサンプル データの例を示します。
Re: ウィルコクソン符号順位検定フィードバックをいただきありがとうございました。
David Q/RI は本の例を使用していたので、異なる回答が返ってきたことにとても驚きました。
ソース コードを取得して手動で問題を処理し、コードがどのステップで異なる動作をするかを確認する方法はありますか?
再度、感謝します。Re: ウィルコクソン符号順位検定実際のソースコードを入手する方法があるとすれば、それがどのようなものかは残念ながら分かりません。ただ、MattFさんの上記の投稿で言及されているWilcoxon検定をSASとJMPの両方でざっと調べてみたところ、どちらもデータペア間の差に対してMann-Whitney-Wilcoxon検定を実行しているように見えます。つまり、差の集合に対して一元配置比較を実行する前に、実際にはペア間の差を自分で計算する必要があるということです。実際にそうしているのですか？
以下は、SAS で Wilcoxon の Matched-Pairs Signed-Ranks テストを実行する方法を示した別の参考資料です。これは、基本的なテストが実際には簡単な手順として提供されるものではないことを示しています (または、なぜこのように少し回りくどい方法で行う手間をかけるのでしょうか)。
http://udel.edu/~mcdonald/statsignedrank.html
実行中のコードのコピーを投稿していただければ、私が正しいかどうか喜んで確認させていただきます。
