こんにちは-私はjmpに、予測モデルとして使用している線形回帰直線の周りの標準偏差を計算させようとしています。具体的には、新しい(将来の)サンプルポイントごとに、そのポイントの回帰結果に±1標準偏差を加算したいと思います。信頼区間の計算方法は知っていますが、標準偏差はわかりません。この計算の方法を理解したら、それをスクリプトに配置して、結果を自動化し、時間に対してプロットしたいと思います。前もって感謝します!
これを行う1つの方法の例を次に示します

names default to here(1);
// Open Data Table: big class.jmp
// → Data Table( "big class" )
dt = Open( "/C:/Program Files/SAS/JMPPRO/14/Samples/Data/big class.jmp" );
biv = dt << Bivariate( Y( :Weight ), X( :Height ), invisible );
biv << fit line( {Save Predicteds} );
biv << close window;
// New column: Column 7
Data Table( "big class" ) << New Column( "Upper",
Numeric,
"Continuous",
Format( "Best", 12 ),
Formula(
:Predicted weight + Col Std Dev( :weight )
));
// New column: Column 8
Data Table( "big class" ) << New Column( "lower",
Numeric,
"Continuous",
Format( "Best", 12 ),Formula(
:Predicted weight - Col Std Dev( :weight )
)
);
Graph Builder(
Size( 528, 454 ),
Show Control Panel( 0 ),
Variables(
X( :height ),
Y( :Predicted weight ),
Y( :Upper, Position( 1 ) ),
Y( :lower, Position( 1 ) )
),
Elements(
Points( X, Y( 1 ), Y( 2 ), Y( 3 ), Legend( 12 ) ),
Smoother( X, Y( 1 ), Y( 2 ), Y( 3 ), Legend( 13 ) )
)
);
txnelsonありがとうございます。あなたの例は直線プロットを示していることに気づきました。スクリプトは、回帰式によって予測される将来の値に対して機能しますか?私の回帰式は、濁度(独立)に基づいて大腸菌(依存)の値を予測し、将来のリアルタイムの濁度値を使用して大腸菌を予測するモデルを作成しようとしています.....また、+/1つの標準を示します偏差。以下のプロットが私の目標です。
私が不明確なことの1つは、あなたの例の「col std dev(:weight)」の部分です。stddevがどのように計算されているかを考える必要があります。私の例では、濁りでも大腸菌でもありません。代わりに、それは回帰当てはめの標準偏差になりますよね?(それがスロープの標準偏差になると思いますか?)アドバイスをいただければ幸いです。

描かれた線は予測値から導き出されたものであるため、適合した回帰ライムに従います。 colstddev確率で使用する列の選択が予測列になります
線に関する標準偏差を計算するためにColStdDev()関数を使用しないでください。 この関数は、応答変動の無条件の推定SDを返します。 条件付きSDが欲しいと思います。 この量は、回帰レポートの[適合の概要]テーブルにRMSEとして報告されます。 予測子を条件とします。
マークありがとう。上記のプロット例(スクリーンショット)では、+/-1標準偏差範囲サイズが変動する予測された(すなわち、退行した)大腸菌値の周り。ただし、RMSEは修繕回帰当てはめに基づく値。明らかに、私が見逃しているものがあります。プロット例のように、SDのサイズは時間の経過とともにどのような条件下で変動しますか?(プロットはUSGS Webページからのものであるため、間違いなく完全に精査されています。) https://www2.usgs.gov/water/southatlantic/ga/bacteria/#:~:text=coli)%20concentrations .-、E.、three%20stations%20within%20the%20CRNRA 。