この式をRの非線形ソルバーに実装しようとするのが最も困難です。nleqslv
とBB
の両方のパッケージを試していますが、私は目を覚ますまで文書を検索して読んでいますが、私の脳を包むことはできません。式自体は次のように動作:Rの和を含む非線形方程式R
s2
、
s2.bar
、及び
price
は等しい長さのベクトルを知られている
s2 * sum(price^(2*x+2)) - s2.bar * sum(price^(2*x)) = 0
。私はBB
にしようとした
最後の試みはこれだった:すべてのヘルプ
gamma = function(x){
n = len(x)
f = numeric(n)
f[n] = s2*sum(price^(2*x[n]+2)) - s2.bar*sum(price^(2*x[n]))
f
}
g0 = rnorm(length(price))
results = BBsolve(par=g0, fn=gamma)
大幅にあなたが式を混乱しているように見える機能で使用される各種部品のご説明から
f [n]はベクトルの1つの要素ですが、長さ(s2)に等しい長さのベクトル全体を割り当てようとしています。あなたはs2とs2.barの代わりにs2 [n]とs2.bar [n]を使用しているのですか?あなたが指数を取った後にベクトル全体を合計しているので、価格は問題を引き起こさない。 – psychometriko
私はもう一度それを試してみましょう: 私はおそらく最初から間違っています。私はそれが同じように簡単になってしまった信じることができない http://i1330.photobucket.com/albums/w561/lompoc421/CodeCogsEqn%201_zpsrwocvotk.gif – Lompoc42