標準ベータの 95% 信頼区間の計算やあみんな、
標準ベータとその 95% 信頼区間をレポートする必要がある分析を行っています。 回帰モデルを実行すると、JMP によって標準ベータ自体が表示されますが、その横にリストされている信頼区間はベータそのものではなく推定値を指します。
最初に新しい標準化された変数を作成し、代わりにこの変数を使用して分析を実行することで CI を取得できることは理解しています。その後、CI はベータ版を参照します。しかし、特に多くの変数を使用しているため、これは少し面倒です。変数の。
JMP でこれを行う簡単な方法はありますか? あるいは、これを行うために Excel シートを作成することは可能でしょうか? 私はそれを試してみましたが、標準ベータ値をその推定値で割り、その結果に 95% CI を乗算すると、通常は小数点以下数桁以内で非常に近くなるように見えますが、わずかにずれており、なぜなのか疑問に思っています。
ご協力いただければ幸いです。 ありがとう！
PS - 標準ベータ版をデフォルトで表示する方法はありますか? CI では実行できますが、標準ベータ版では実行できないようです。
Re: 標準ベータの 95% 信頼区間の計算@デバサー
私はモデルのフィット インターフェイスを使用し、モデルとウィンドウをキャプチャするスクリプトの作成方法を知っているため、次の 2 つのオプションがあります。
- ブートストラップを実行するコードを追加します (これには、置換を伴うサンプリングとモデルの繰り返し実行が必要です)。 これは非常に遅くなるのではないかと心配しています。
- 参照した記事で提案されているデルタ法を計算するコードを追加します。
もちろん、モデルの仮定が厳密に正しい場合にのみ、参考記事で正確であると証明された計算を追加するという代替案もあります。
デルタ法のアルゴリズムについては勉強したことがありません。 数日以内にこの投稿に返信して、これらの方法のいずれかがパフォーマンスに大きな影響を与えることなく機能するかどうかを確認する予定です。
Re: 標準ベータの 95% 信頼区間の計算ディベーサー、
私はモデルのフィット インターフェイスを使用し、モデルとウィンドウをキャプチャするスクリプトの作成方法を知っているため、次の 2 つのオプションがあります。
- ブートストラップを実行するコードを追加します (これには、置換を伴うサンプリングとモデルの繰り返し実行が必要です)。 これは非常に遅くなる可能性があると懸念しています。
- 参照した記事で提案されているデルタ法を計算するコードを追加します。
もちろん、別の方法としては、モデルの仮定が厳密に真実であることのみが正確であると参考記事で実証された計算を追加することもできます。
デルタ法のアルゴリズムについては勉強したことがありません。 数日以内にこの投稿に返信して、これらの方法のいずれかがパフォーマンスに大きな影響を与えることなく機能するかどうかを確認する予定です。
Re: 標準ベータの 95% 信頼区間の計算ああ、それで私の問題があるんです！ 私はJMPのProエディションではなく、通常バージョンを使用しています。
ブートストラップせずにこれを行う方法はありますか、それとも運が悪いのでしょうか?
再度、感謝します！
Re: 標準ベータの 95% 信頼区間の計算まず第一に、スクリプトを実行すると信頼区間が表示テーブルに表示されますが、インターフェイスを使用して列を確認すると、新しい列は表示されません。これは、使用したインターフェイスが組み込みであり、組み込みのもののみをリストするためです。 -列内。
次に、JMP Pro 14.1 を使用しています。 オンライン ヘルプによると、Bootstrap は JMP Pro でのみ利用できるようです。それが問題なのかもしれない。
JMP PRO をお持ちの場合、エラーを特定するにはログを確認する必要があります。 ログを表示するには、CTRL+SHIFT+L を押します。 右クリックして「ログのクリア」を選択し、スクリプトを再度実行します。 ログにはエラーが報告されるはずです。
スクリプトを実行してブートストラップを実行できます。
- 入力からランダムに選択し、標準ベータ版をキャプチャします
- 分布から間隔を計算します。
今後数日で時間ができたら、脚本を書きます。 しかし、おそらく誰かがそれを行うためのアプリを構築し、それに対応するでしょう。
すみません、JMP PROについて言及するべきでした。
Re: 標準ベータの 95% 信頼区間の計算やあ、
本当にありがとうございました！ これはまさに私たちが探しているもののようです。 JMP 14.1 を使用していますが、何らかの理由でスクリプトが正しく機能しません。 単に正しく使用していないだけなのかどうかはわかりません。
データ ファイルは期待どおりに開きますが、スクリプトを実行すると、何らかの理由で標準ベータ CI の列が表示されません。 これは、自分のテーブルで動作するように何も変更せず、「Big Class」サンプル テーブルを使用するだけです。 ここで何らかの素人の間違いを犯しているのではないかと思います。 取得した出力のスクリーンショットを添付しました。 私が何を間違っているのかわかりますか？
ご協力いただきましたこと、重ねて御礼申し上げます。本当に感謝しています。
Re: 標準ベータの 95% 信頼区間の計算ps この論文には、推奨される計算のための R コードも含まれています。JMP は R を呼び出すことができますが、R ユーザーではない (R をインストールしていない) 場合はそのパスをたどりませんでした。
Re: 標準ベータの 95% 信頼区間の計算@デバサー標準ベータ版の信頼区間を計算して報告するための JSL メッセージが見つかりませんでした。 私の経験はスケール係数と推定値のほうが多いので、このトピックについてインターネットで少し検索しました。 これらの信頼区間の最適な方法については、少し議論があるようです。この紙が気に入りました標準化回帰係数の信頼区間の計算。 彼らのシミュレートされた逸脱モデルからのこの論文の結論は、「デルタ法 (Yuan & Chan、2011) またはノンパラメトリック ブートストラップ (Efron & Tibshirani、1998) を推奨する」です。
このトピックに関する私の検索は広範囲にわたるものではなかったので、JMP の技術専門家からの意見を聞くのは興味深いでしょう。 この文書をサポートとして使用し、以下のスクリプトは JMP が提供する素晴らしい Bootstrap() メソッドを使用して、標準ベータの 95% の下限と上限を計算します。 ブートストラップには時間がかかり、このスクリプトは 1000 回の反復を実行します。 これを 500 に減らすこともできます。
添付のスクリプトを参照してください。 JMP 14 でのみテストしました。 異なるバージョンがある場合は、その旨を明記してください。
