BBOptimはScalaで同等か代替のものです

Question

Rコードの一部をScalaに変換する予定で、BBOptimという関数がありました。 SPGのラッパーであると思われますが、数学についての知識はほとんどありません。Scalaに相当するコードは何か分かりません。

たとえば、以下のコードをScalaに変換できますか？あるいはこれに代わる方法がありますか？私は、SpectralProjectedGradientまたはBreezeライブラリのNonLinearMinimizerが使用されていると思われます。

# Use a preset seed so test values are reproducable. 
require("setRNG") 
old.seed <- setRNG(list(kind="Mersenne-Twister", normal.kind="Inversion", 
    seed=1234)) 

rosbkext <- function(x){ 
# Extended Rosenbrock function 
n <- length(x) 
j <- 2 * (1:(n/2)) 
jm1 <- j - 1 
sum(100 * (x[j] - x[jm1]^2)^2 + (1 - x[jm1])^2) 
} 

p0 <- rnorm(500) 
BBoptim(par=p0, fn=rosbkext) 

ありがとうございます。

編集

私だけJVMを使用して制限していますのでScalaからRを呼び出しても可能ではありません。

Answer 1

このRライブラリは、spgのラッパーです。これは、スペクトラル投影勾配アルゴリズムです。

TANGO projectによる非線形グラジエント最適化のためのSPG implementationは、Javaを含む多くの言語へのインタフェースを備えたFortran 77で記述されています。すべてのJavaライブラリはScalaで動作します。したがって、これは適切なソリューションでなければなりません。この実装は、Spectral Projected Gradientsにいくつかのacademic articlesの焦点です。

また、Open Optimization LibraryのJava APIをチェックアウトすることもできます。

Answer 2

最も慣れ親しの国のScalaではなく、JVMに限定されている場合は、Apache Commons Mathから適切な最適化ルーチンを選ぶのが最も安全です。あなたが使っている正確なものはありませんが、一般的なソルバは、私がそれらに投げたほとんどのクラスの問題ではうまく機能します。適切なボックスの制約を得るには、いくつかの異なるクラスのソリューションを試してみる必要があります。壁の近くで急なペナルティを加え、壁の近くにいる場合は、検索スペースをそれらのサーフェスに制限します。

それ以外の場合は、rScalaを使用してScalaからRを呼び出すのはかなり簡単です。そうすれば、アルゴリズムをまったく変更する必要はありません。 （しかし、あなたはRをインストールしなければなりません）