多重仮説検定<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />こんにちは、みんな、
私はすぐに多数の仮説検定を実行する必要があるでしょう。そのような場合、タイプIのエラーのため、p値に頼るのは間違ったことであることはよく知られているので、誤検出率、より正確には正の誤検出率（Storey）を使用したいと思います。私の "調整されたp値"としてq値を計算します。
私の質問は2つに分かれています。
1。JMP 9にはpFDRとそれに対応するq値を計算するためのコードがありますか。私はそのような機能がJMP Clinicalに既に存在していると私が思うところから私はインターネット上で何かを見ました。私はJMP 9のライセンスしか持っていませんが、JMPの臨床費用がいくらなのか、または予算があるのかどうかはわかりません。単純なJSLコードがどこかにあるとしたら、それは私にとって非常に便利です。
2。そのようなコードがある場合、たぶんそれは出力としてp値のベクトルを必要とするでしょう。1つのカテゴリカル変数と500の連続変数を持つJMPファイルがあり、500の各変数に対してt検定を実行し、次の内容を含むp値の変数を新しいファイルに保存するとします。それらの500 JSLでそれを行うにはどうすればいいですか？それは可能ですか？そのJMPは自動的に500のテストを実行し、p値のみを保存しますか？
前もって感謝します....
多重仮説検定みんな助けてくれてありがとう。
週末が終わったら、コードを使って作業し、スライドを読みます。
リンクをありがとうDave、それはそれが私が必要としているものだけであるように見えます、そしてアシスタントをみんなにありがとう、私はそれを感謝します！
多重仮説検定それらのスライドにはたくさんの良い情報があります。Daveを共有してくれてありがとう。
多重仮説検定以下はその方法を説明したスライドです。 JSLを使ってRに接続する
多重仮説検定行列用||水平連結演算子です。Concat（m、owr ...）は等価です。
私はFDRの経験がありませんが、私はJMP-Rで少し遊んだことがあります。以下は一例です。前半は上記の例と同じです。出力はjmpログに表示されます。出力を文字列として取得するには、Log Capture（）関数内にR Submitを入れます。
dt = 開いた （ "$ ENGLISH_SAMPLE_DATA / Big Class.jmp" ） ;
メートル = [] ;
にとって （ 私 = 4 、私 <= Nコル （ dt ） 、私 ++、
わ = 一方通行 （ Y （ カラム （ 私 ） ） 、バツ （ ： 性別 ） 、 tテスト （ 1 ） 、見えない ） ;
フクロウ = 報告書 （ わ ） ;
メートル = メートル || フクロウ [ "tテスト" ] [ ナンバーコルボックス （ 4 [] 1 ] ;
） ;
DPTP = 新しいテーブル （ "t検定" 、新しいコラム （ "p" 、数値 、設定値 （ メートル ） ） ） ;
varname = "p" ;
R Init （） ; // R接続を初期化する
R送信 （ DPTP ） ; //データテーブルをRに送信
// Rコマンドで文字列を送信する
R送信する （ "データ< - dtp
名前（データ）
添付（データ）
平均（" || varname || 「）」 ） ;
Rterm （） ;
多重仮説検定ありがとう、それは良い考えです…
もう1つやらなければならないことがあります（私は忙しくてJSLを正しく学ぶ時間がないため、謝罪を受け入れます）。あなたのコードから私はp値の、1列を含むJMPテーブルになるでしょうね？Rに接続するJSLコードを作成する必要がある場合は、少し手引きが必要になるでしょうか。それとも、テーブルをエクスポートしてRから直接作業する必要があるのでしょうか。
新しいテーブルを再検索した後、Rに接続してmean（p）のような単純なものでも実行する方法の簡単な簡単な例を誰かに書いてもらえたら、あとでpFDR用に変更することができます。
そして最後の質問（約束）の最後の1行：
メートル = メートル || フクロウ [ "tテスト" ] [ ナンバーコルボックス （ 4 [] 1 ] ;
なぜm または owr [...]つまり、なぜ||必要ですか？
たくさんの人に感謝します、あなたは私に大きなデータセットをRにエクスポートしてそこでコーディングさせるのに多くの時間を節約しています....
多重仮説検定1。私は自由に利用できるJMPコードを知りませんが、Rにはいくつかの選択肢があります。 http://strimmerlab.org/notes/fdr.html いくつかのFDR関連のパッケージを一覧表示します。JMP 9はRに接続してRから結果を取得することができます。これは、JSLのアルゴリズムを最初からコーディングするよりもはるかに簡単な作業です。
2。それは可能なはずです。以下はあなたにいくつかのアイデアを与えるかもしれない例です。これは、一方向近似で固定のx列の右側にあるすべての数値y列をループ処理し、レポートからp値（両側t検定）を取得してそれを行列に挿入することによって機能します。行列は、最後に新しいデータテーブルの列に追加されます。
dt = 開いた （ "$ ENGLISH_SAMPLE_DATA / Big Class.jmp" ） ;
メートル = [] ;
にとって （ 私 = 4 、私 <= Nコル （ dt ） 、私 ++、
わ = 一方通行 （ Y （ カラム （ 私 ） ） 、バツ （ ： 性別 ） 、 tテスト （ 1 ） 、見えない ） ;
フクロウ = 報告書 （ わ ） ;
メートル = メートル || フクロウ [ "tテスト" ] [ ナンバーコルボックス （ 4 [] 1 ] ;
） ;
新しいテーブル （ "t検定" 、新しいコラム （ "p" 、数値 、設定値 （ メートル ） ） ） ;
多重仮説検定pFDRはJMP Genomicsには存在しますが、標準（またはPro）JMPには存在しないと思います
