モデル プロットの信頼帯用の JSL、より良いアイデアはありますか?モデルの結果を重ねて、信頼帯を陰影付けしたデータのプロットを取得したいと思います。
たとえば、Fit Model からの次のプロットに似ていますが、適合線の周囲に信頼区間が陰影付けされています。
この特定の例では、FitYbyX で探しているものを取得できることはわかっています。
これで私が求めているものが示されましたが、一部の適合では機能しません。たとえば、この適合を実行して、性別ごとに同じ傾斜に制限したい場合などです。
また、他の人が同様のことを要求しているが、フィット曲線や非線形フィットも含まれている投稿もいくつかあります。
解決済み: ロジスティック回帰プロットに信頼区間を追加するにはどうすればいいですか? - JMP ユーザー コミュニティ
JMP の「フィット曲線」で信頼区間をプロットする - JMP ユーザー コミュニティ
これらのスレッドの提案のうち 2 つは、私にとって完全に満足できるものではありません。
1 つは、プロファイラーに信頼帯があることです。これは便利ですが、適合度を重ねて表示したいこの種の例には役立ちません。
(モデルの適合では、回帰プロットを作成するためにも要因が多すぎる可能性が十分にあることも理解しています。これが、機能が充実していない理由かもしれません。プロファイラーで調査する必要があるより複雑なモデルでこれを使用することが期待されています。ただし、単純なプロットや曲線が重ねられたプロットが役立つ場合があります)。
もう 1 つは、保存された予測値を使用してグラフ ビルダーの曲線間に陰影を付ける便利な方法でした。
解決済み: グラフビルダーのプロットに陰影付きの信頼区間を作成するにはどうすればいいですか? - JMP ユーザー コミュニティ
場合によっては良いのですが、予測は実際のデータを超えて滑らかな曲線になるようにしたいのです。
そうなると、プラットフォームから信頼帯の数式を取得し、それをグラフィック スクリプトで使用することになるようです。
その道筋は見えていると思いますが、あまりにも面倒なので、あまり深く考えすぎないうちに、もっと賢い方法を見逃しているのではないかと思っています。それとも、これはすでにどこかで実証されているのでしょうか?
私の例では、式は次のようになります (信頼限界の下限のみ)。
これにいくつかの代入を行って、性別ごとに Y 曲線を作成する必要があります。
下限ラインを描いたところですが、ここまでは順調です...
予測自体と上側信頼線を繰り返し、次にシェーディングすることができます...各曲線のペアに沿った一連のポイントを使用してポリゴンを作成できます...ああ、それはできることはわかりますが、過度に複雑に思えます... Y関数間でシェーディングしたいというこの部分だけでも、もっと簡単な方法を見逃していますか? それとも、これが本当に正しいパスですか?
要するに：
実際に、これらのモデリング プラットフォームからこのタイプのプロットをより直接的に取得する方法はありますか? (ない場合は、価値のあるウィッシュ リスト項目になると思います)。
そうでない場合、私が説明したグラフィック スクリプトよりも、目的のプロットを作成するためのより良いアイデアはありますか?
そうでない場合、グラフィック スクリプトで Y カーブ間のシェーディングを行うためのポリゴンよりも簡単なものが欠けているのでしょうか?
ありがとう！
高度な統計モデリング自動化とスクリプト基本的なデータ分析とモデリングRe: モデル プロットの信頼帯用の JSL、より良いアイデアはありますか?こんにちは@hardner さん
少なくとも私が理解している限りでは、Fit Y by X プラットフォームを通じて、やりたいことを実現できるはずです。非線形プラットフォームについてはよくわかりませんし、同じことができるかどうかもわかりません。
とにかく、二変量フィットの横にある赤い三角形に移動して、フィットラインの代わりにフィットスペシャルを選択すると、データの変換を実行したり、特定の傾きまたは切片を持つように線形フィットを制限したりできます。あなたの例を使用して、0.113 の傾きを使用してそれを実行しました。
また、信頼帯はデータ ポイントを超えて拡張されるため、軸の範囲を調整するだけで済みます。
あなたの投稿で説明されていた内容から判断すると、これがあなたが求めていたものだと思います。あるいは、私があなたが望んでいることを完全に理解していないのかもしれません。
お役に立てれば、
DS
Re: モデル プロットの信頼帯用の JSL、より良いアイデアはありますか?こんにちは@SDF1 さん
ありがとうございます。問題は、傾きを指定したために信頼帯が非常に狭くなることです。その部分には非常に自信があり、切片のみを適合させます。適合モデルで傾きを適合させてから、ここで指定すると、適合線がきれいに重ね合わされ、適合モデルと比較して適合の方程式がきれいに表示されます。また、信頼帯もきれいに描画されます。つまり、まさに私が探していた種類の画像であり、一部の適合では FitYbyX が機能しますが、それらは私が求めている信頼帯ではありません。
Re: モデル プロットの信頼帯用の JSL、より良いアイデアはありますか?こんにちは@hardner さん
あなたが何を求めているのか、まだ理解できていないようです。
元の投稿で、「...たとえば、これを適合させて、性別ごとに同じ傾きに制限したい場合」と述べていました。これが、私が示した例で、性別ごとに傾きを固定するというまさにその操作を実行した理由です。
最新の投稿で、「しかし、それらは私が求めている信頼区間ではありません」と述べています。どの信頼区間をお探しですか? Fit Y by X プラットフォームには、平均の期待値を表す信頼曲線 Fit と、誤差とパラメータ推定値の変動を反映する個々の予測値の信頼限界を表す信頼曲線 Indiv の 2 つの信頼曲線があります (JMP のオンライン ヘルプ ページから要約)。探しているのは Indiv 曲線ですか?
DS
Re: モデル プロットの信頼帯用の JSL、より良いアイデアはありますか?こんにちは@SDF1 さん
私の問題について引き続き検討していただきありがとうございます。私の例をもっと明確に説明しなかったことをお詫びします。私の知る限り、Fit Y by X は信頼区間には機能しません。
問題は、各性別で傾きが同じになるように制約する場合、その傾きを適合させたいということです。Fit Y by X では、傾きを制約できますが、グループごとに同じであるだけでなく、値を指定する必要があります (そう思います)。
Fit Model でフィッティングを行うと、各性別の傾き (または切片) が同じになるように強制しながら、傾きと切片の両方を実際にフィッティングすることができます。
次に、「列を保存」の下で、平均信頼限界の値または数式を保存できます。
以下に、右側のグラフ ビルダーを使用してこれらの制限をプロットしました (これは、テーブル内の特定のポイントに制限されるため、問題に対する解決策としては気に入らなかったのですが、少なくとも限られた範囲では、このケースに必要な信頼帯を示しています)。左側は FitYby X で、Fit Model から取得した傾斜値に制限するように指示しています。これらの信頼帯 (平均値、つまりより暗く狭い帯が欲しいのですが、Indiv も追加しました) は線形です。指定された傾斜に 100% の信頼度があるため、重みが高くても低くても扇状に広がりません。データの中央では、Fit Y by X のものはすべて垂直方向の変動、つまり切片の適合に関するものであるため、非常に接近しています。
個々の曲線は、バンド内に含まれる個々のポイントに関する信頼性を求めるものではありません。平均に関する信頼バンドは必要ですが、傾きと切片の適合に関する信頼バンドも必要です。
FitYbyX は、元のデータ、フィッティング、信頼帯をうまくオーバーレイする唯一のプラットフォームのようですが、すべてのフィッティングをそこで実行できない理由があります。Fit Model、Nonlinear、Fit Curve などの他のプラットフォームでは、FitYbyX では実行できないいくつかのフィッティングを実行できると思います (これはその一例です)。また、X が 1 つしかないためこのようなプロットを作成できる場合もありますが、代わりに、フィッティングの信頼限界を保存するオプションはありますが、オーバーレイ プロット オプションはありません。
Re: モデル プロットの信頼帯用の JSL、より良いアイデアはありますか?こんにちは@hardner さん
あなたのジレンマがよく理解できました。あなたが直面している問題を言い換えるために時間を割いていただき、ありがとうございます。
つまり、Fit Y by X は、プロットを広げて信頼帯がどのように変化するかを確認できる、フィットの優れたオーバーレイを可能にする唯一のプラットフォームであるというあなたの意見は正しいと思います。また、データ テーブルに保存される予測式は、そのデータ テーブルの行数だけになるというのは理にかなっています。
ただし、希望にかなり近い回避策があるかもしれません。完璧ではありませんが、使用する必要があるさまざまなプラットフォームの制限 (Fit Y by X ではない他のモデル フィッティング プラットフォームなど) があるため、おそらくこれが最も近いものになります。
私が思いついたアイデアは、データ テーブルに、値が予想される平均値であるダミー行を持つ「ダミー」列をいくつか作成することでした。そうすれば、信頼区間が混乱することはありません。
以下は、Graph Builder と対応する JSL コードを使用した結果です (示したものと同様です)。
Graph Builder(
Graph Spacing( 4 ),
Variables(
X( :hidden weight ),
Y( :Lower 95% Mean height ),
Y( :Upper 95% Mean height, Position( 1 ) ),
Y( :Pred Formula height, Position( 1 ) ),
Y( :hidden height, Position( 1 ) ),
Overlay( :sex )
),
Elements(
Area( X, Y( 1 ), Y( 2 ), Legend( 2 ), Area Style( "Range" ) ),
Line( X, Y( 3 ), Legend( 1 ) ),
Points( X, Y( 4 ), Legend( 1 ) )
),
SendToReport( Dispatch( {}, "X title", TextEditBox, {Set Text( "weight" )} ) )
);
そこに到達するために私が実行した手順は次のとおりです (もちろん、これはすべて JSL でコーディングできますが、作業するケースごとに時間がかかります)。
1. 以前と同じように適合モデルを実行した後、目的の数式列 (予測平均、および +/-95% 平均 CI) を保存します。
2. ダミー列を 2 つ作成します。名前には :hidden height と :hidden weight を選択しました。
2a. データ テーブルに 4 つの行を追加します。これにより、M と F の両方の性別に対して 2 つのポイント (最低値と最高値) を設定できます。
3. 各列を数式にする
3a. :hidden の高さの場合は次のようになります。
3b. :hidden ウェイトの場合は次のようになります。
注: 値として 0 と 200 を選択したのは、簡単だったためです (切片が 1 つ)。もう 1 つは、単に高い方の極値として選択しましたが、これらは必要に応じて設定されます。
4. 予測平均と +/-95% 信頼区間の数式を変更して、計算に :weight を使用する代わりに、計算に :hidden weight を使用するようにします。予測平均の例を以下に示します。
5. 戻って JSL コードを再実行しますが、:height 列と :weight 列を対応するダミー列に変更して、データ テーブルにこれらの追加の 4 行が含まれるようにします。
6. 最後に、最後の 4 行をデフォルトのマーカー スタイルから「/」マーカー スタイルに変更して、これらの平均値の終了点がデータ自体ほど目立たないようにしました。[その他...] を選択して別の文字を使用することもできます。MS Word の二重右向き引用符 (») を試したところ、問題なく動作しました。
「非表示の重量」が表示されないように、x 軸の名前も「重量」に変更したことに注意してください。
しかし、あなたが説明した内容に基づくと、これがあなたが求めているものに対応できると思います。これは最も簡単でエレガントな回避策ではありませんが、目的は達成できます。たとえば、Graph Builder でポリゴンを操作するよりも簡単だと思います。
考慮し、考えるべきこと。
お役に立てれば、
DS
Re: モデル プロットの信頼帯用の JSL、より良いアイデアはありますか?ありがとう！
グラフ ビルダーでとてもきれいに表示できました。テーブルを「改変」したくはありません (私の使用例は他の人のためのツールになります) が、テーブルのプライベート コピーは可能です。とにかく検討すべきことですが、少なくともグラフィックス スクリプトの方法を試して、ポリゴンを滑らかにしてきれいに表示できるかどうか確認する必要があると思います。
他の人の意見も参考にして、確かにこの状況のようで、私が見逃している明らかな組み込みオプションはないということに特に感謝しています。
