非線形スクリプトプラットフォームを介してデータテーブルに対してNLRを実行すると予想通りの結果が得られますが(Croucher http://www.walkingrandomly.com/?p=5215を参照)、対応するスクリプト(添付)を実行するとJMP 7でクラッシュ(アクセス違反)が発生します。 [テーブルの付け方は見ませんでしたが、必要に応じて送ることができます。]
私は問題なく同様のスクリプト[データと方程式の変更]を使用しました。
何かご意見は?
スクリプトの方程式セクションは次のとおりです。Parameter({b0 = 1、b1 = 0.2}、
b0 *余弦(b1 *:X)+ b1 *正弦(b0 *:X))
Clear Log(); Clear Globals();
xv = [ -2, -1.64, -1.33, -.07, 0, 0.45, 1.2, 1.64, 2.32, 2.9];
yv = [.699369,.700462,.695354,1.030905,1.97389,2.411430,1.910910,0.919576,-0.770975,-1.42001];
dt = New Table( "Croucher Data",
New Column( "x", SetValues( xv ) ),
New Column( "y", SetValues( yv ) ),
); //close New Table
dt <<
Nonlinear(
Y( :y ),
X( :x ),
Model( Parameter( {b0 = 1., b1 = 0.2},
b0 * Cosine( b1 * :x ) + b1 * Sine( b0 * :x ) ) ),
QuasiNewton SR1,Finish,plot(1),
Confidence Limits
);//close NonLinear
非線形スクリプトありがとう、タイピングの問題!私は自分のファイル用にそれらを修正しました。
非線形スクリプトスクリプトをどうもありがとう。それは私にとってはうまくいきます。私はその機能をあまり使ったことはありませんが、明らかに多くのモデルにとって理想的です。
非線形スクリプトxの値の1つは-0.07ではなく-0.7になるはずです。
非線形スクリプトあなたのスクリプトはうまくいきます。リンクされた例と比較すると、入力データにわずかな違いがあります(xv [4]&yv [9]を参照)。
プラットフォームを使用せずにこの非線形フィットを実行することも可能です。
Model = Expr(
b0 * Cosine(b1 * xv) + b1 * Sine(b0 * xv)
);
xv = [-2, -1.64, -1.33, -0.7, 0, 0.45, 1.2, 1.64, 2.32, 2.9];
yv = [0.699369, 0.700462, 0.695354, 1.03905, 1.97389, 2.41143, 1.91091, 0.919576, -0.730975, -1.42001];
b0 = 1;
b1 = 0.2;
SSE = Minimize(Sum((Model - yv) ^ 2), {b0, b1});
Eval List({b0, b1, SSE});