JMP 14.1.0 で二項分布分析を実行する方法こんにちは、
誰かがこのクエリを手伝ってくれませんか (サンプルデータセットを添付)
二項分布を実行したいと思います。治療（4回の治療）および日（3時点）ごとの死亡率データ数に関する分析（JMP）。
ありがとう
消費者および市場調査Re: JMP 14.1.0 で二項分布分析を実行する方法こんにちは、
私は生存分析に関するあなたの提案に従いました。
添付の Excel には、(すべてのペンの) 変換されたデータが表示されます。これは、ペン 5 について説明したもの (「Time Int Left」、「Time Int Right」、および「Count」の列) に従っていると思います。
次に、治療別の生存分析、生存分布分析を実行しました (統計出力は添付されています)。
これはあなたが提案していたことですか？
私は生存出力について詳しくないので、意味についてのガイダンスが必要です。
Re: JMP 14.1.0 で二項分布分析を実行する方法確かに、データは列見出し (元のファイル) に適切に記述されていませんでした。 列には生きた動物に関するデータが含まれていました。
代わりに死亡率に関するデータを添付しました。
あなたが提案した生存テストを検討してみます。
詳細な説明をありがとうございます。
Re: JMP 14.1.0 で二項分布分析を実行する方法説明により多くのことが明確になりました。 そして、表にした観察が独立していないことを確認します。 Pen=1 の極端なケースを見てください。 Day=0 での死亡動物の合計が 0 であると仮定するのはエラーであり、0 であるべきです。 4 つの観察は、22 匹の動物からなる同じグループからの反応を記録しています。
この問題は別の方法で処理できるにもかかわらず、このタイプの研究は生存分析の対象となります。 以下は、研究において興味のある最も単純な質問です。
- Pen = 1 から Day = 22 までの生存率はどれくらいですか? 答えは100%です。 そして補体死亡率は0%です。
- Pen = 5 から Day = 22 までの生存率はどれくらいですか? Day = 0、10、22、または 35 日に正確に死亡した動物がないと仮定すると、Day = 10 と Day = 22 の間に 1 頭の動物が死亡したため、答えは 21/22 = 95.5% になります。補体死亡率は4.5%です。
- Pen = 5 から Day = 35 までの生存率はどれくらいですか? Day = 22 と Day = 35 の間に別の動物が死亡したため、答えは 20/22 = 90.9% となります。 補体死亡率は9.1%です。
これは Pen=5 データのスクリーンショットです。
あなたの研究では、どの変数が二項分布を持っていますか? これは、指定された時間間隔内での囲い内の死亡 (または生存) の数です。
檻内の時間の経過に伴う死亡率も分布です。 しかし、これは二項式ではありません。 これが先ほどコメントした方法の理由です。
ここで、このデータに対して上記で行ったような単純な生存分析を実行する方法を説明します。 まず、データを別の形式に変換する必要があります。 以下は、Pen = 5 の変換されたデータです。
ステップ 1、「Time Interval Left」と「Time Interval Right」の 2 つの列を作成し、次のスクリーンショットのように値を入力します。 これはすべてのペンで同じである必要があります。 2 つの列の最初の行は、Day = 0 から Day = 10 までの時間間隔を定義します。 2 行目と 3 行目も同様に後続の時間間隔を定義します。 最後の行は、Day = 35 以降の期間を定義します。
ステップ 2、「カウント」を作成し、スクリーンショットのようにペン = 5 の値を入力します。 値は、個々の囲いにおける死亡数と生存数に基づく必要があります。 Pen = 5 の場合、意味は次のとおりです。第 1 時間間隔で死亡者がいない場所。 2 番目の時間間隔で 1 人の死亡者が発生します。 3 番目の時間間隔では 1 人の死亡者が発生します。残りの 20 人は全員、Day = 35 の時点で生存していました。
この形式の詳細については、Life Distribution プラットフォームのドキュメントを参照してください。 https://www.jmp.com/support/help/en/15.2/index.shtml#page/jmp/event-plot.shtml#
Pen = 5 のテーブルのみを作成していることに注意してください。 そして、次のイラストはそのペンのみにあります。
次に、「分析」メニューに移動し、「信頼性と生存」、「寿命分布」を見つけて、次のような起動ダイアログを指定して、「OK」をクリックします。
デフォルトのレポートには「死亡率」の結果が表示されます。 レポートの赤い三角ボタンのメニューに移動すると、「生存曲線の表示」が表示され、「死亡率」と「生存率」の結果が切り替わります。
レポートの下部近くに「ノンパラメトリック推定」アウトライン ノードがあるので、それを開きます。 先ほど述べた死亡率の推計値が載っています。
上記はペン1本のみの場合です。 これらの時間間隔ですべての檻の死亡 (生存) 率を収集した後、この道をたどって、かなりの距離を進むことができるかもしれません。 ペンやトリートメントなどで比較できます。
私が説明したようにすべてのデータを変換すると、より包括的な分析を実行して、グループ化変数を使用して死亡率を同時に推定することもできます。 変換されたデータをフォローアップしていただければ、データ内で何が見つかったかに応じて、追加の提案をさせていただく可能性があります。
Re: JMP 14.1.0 で二項分布分析を実行する方法それは正しい。研究の開始時 (0 日目)、64 のすべての囲いは各囲いに 20 匹の動物から始まります。
各時点 (0、10、22、35 日目) での死亡動物の数が列 6 (シート 1) に記録されます。
4 つの処理があり、各処理が 16 回繰り返されます。
死亡動物の数に対する治療効果の経時変化が関心のある変数です。
データは次のように表示されます (Sheet1)。
- 列 1: ペン ID
- 列 2: 複製
- 列 3: 治療 ID
- 列 4: 時点
- 列 5: 生きている動物の数
- 列 6: 死亡した動物の数
Re: JMP 14.1.0 で二項分布分析を実行する方法データには 64 個の固有のペンがあります。 生きている動物の合計はすべて22以下です。 ただし、Day=0 では、Sum Live Animals はすべて 22 です。 生きている動物の合計は、個々の囲いの時間の経過とともに減少するだけです。 私の質問: あなたのすべての囲いは、day=0 に 22 匹の生きた動物から始まり、時間の経過とともにいくつかの動物が死亡します。また、個別の囲いで生きている動物と死んだ動物を追跡していますか?
Re: JMP 14.1.0 で二項分布分析を実行する方法ありがとう。
それが私が探していたものですが、関数 logit または probit についてはわかりませんでした。
Re: JMP 14.1.0 で二項分布分析を実行する方法ありがとう、リュー
分析の目的は、治療 (4 つのレベル) が死亡率に影響を与えるかどうか、また治療効果が時間の経過とともに変化するかどうかを確認することです。
Re: JMP 14.1.0 で二項分布分析を実行する方法この場合に役立つ可能性があるもう 1 つの分析は、プロビット分析です。 [分析] > [モデルの適合] を使用して、これらの効果を仮定できます。
GLM プラットフォームを使用すると、重大な効果、適合度の欠如、過剰分散を調査し、応答 Prob( Dead ) をプロファイリングできます。
Re: JMP 14.1.0 で二項分布分析を実行する方法- データを JMP にインポートします。
- データ、またはデータを収集した方法について説明します。 変数と観測値を説明します。 必要に応じてデータをクリーンアップします。 時点は 4 つあり、time=0 での観察は他のものと矛盾すると思います。
- 二項当てはめの目的を決めたら、試行回数と死亡または生存の数はいくらかという質問をしてください。 「どの時点で」ということも聞かなければいけないと思います。
- データを変更せずにロジスティック回帰または同様の分析に投入する場合は注意してください。観察はそのままでは独立していないと考えられるからです。 確認するには説明を参照する必要があります。
- 目的が経時的な死亡率をモデル化することである場合、二項分布は関連性のある主題ではない可能性があります。 時間の経過に伴う死亡率は曲線です。
