問題タブ [stochastic]

答えを探してみましたが、私の問題と完全に一致する答えが見つかりませんでした。

私は生物学的システムの確率的シミュレーターをやっています。その結果は、ランダムな時点での濃度レベルを持つ「散布図」時系列です。ここで、複数のシミュレーション実行の平均時系列を取得できるようにしたいと考えていますが、それぞれ数千回の測定を伴う最大 500 回のシミュレーション実行が予想されるため、どのように進めるかについて疑問があります。

当然、間隔を「バケツ」して、おそらく精度がいくらか失われるか、欠落している測定値を補間しようとする可能性があります。しかし、私の場合、好ましい方法は何ですか?

これはJavaで実装する必要があり、メソッドを説明する論文よりも引用を希望します。

ありがとう！

私は小さな遺伝的アルゴリズム フレームワークを実装しています。これは、合理的なものを作成してオープン ソースとして投稿しない限り、主に個人的な使用を目的としています。今は選抜技術に力を入れています。これまでのところ、ルーレット ホイールの選択、確率的ユニバーサル サンプリング、およびトーナメントの選択を実装しました。私のリストの次はランクベースの選択です。すでに実装した他の手法よりも、それに関する情報を見つけるのに少し苦労しましたが、これまでの私の理解は次のとおりです。

1. 次のラウンドで適切な親を取得したい母集団がある場合、まずそれを調べて、各個体の適応度を母集団の合計適応度で割ります。

2. 次に、他の選択手法 (ルーレット ホイールなど) を使用して、実際に誰を繁殖用に選択するかを決定します。

これは正しいです？もしそうなら、ランク調整は一種の前処理ステップであり、その後に候補を選択する実際の選択手順が続く必要があると考えるのは正しいですか? これを誤解している場合は、修正してください。追加のポインタに感謝します。

最小限のコーディングで確率動的システムのシミュレーションに使用できるソフトウェア パッケージ (無料の方がよい) を探しています。たとえば、次のように指定してシステムをモデル化できるようにする必要があります。

ここXで、 は状態、Y観測、A,B行列であり、VnおよびWnは法線、T などの異なる pdf を持つノイズ ソースです。

私はすでに試しましscilabたR。これらは非常に強力に見えますが、上記の方程式を指定してモデルを作成することを直接サポートしているようには見えません。

RのGillespie SSAの構文を知っている人はいますか? 私が見つけたコードは次のとおりです。

SSA のこのレイアウトまたは新しいシンプルなレイアウトについてサポートが必要です

Rで確率的下降勾配を使用してロジスティック回帰をプログラムしようとしています。たとえば、「ex2data1.txt」という名前のAndrew Ngの例に従いました。

要点は、アルゴリズムが適切に機能するということですが、thetas の推定は私が期待したものとまったく同じではありません。そこで、この問題を解決するためにアルゴリズム全体を変更しようとしました。しかし、それは私にはほとんど不可能でした。この問題の原因となっているエラーを検出できませんでした。したがって、誰かが例をチェックして、シータが正しく計算されない理由を教えてくれると非常に便利です。ほんとうにありがとう。

プログラミングに関しては、Rや行列計算で実装されている関数は使用していません。Hadoopでコードを使用したいため、ループで合計と減算を使用しているだけで、行列計算や、「sum」、「sqrt」などのRですでにプログラムされている関数さえ使用できません

確率的勾配降下法は次のとおりです。

そしてロジスティック回帰：

私のコードは次のとおりです。

私のシータ

確率的 GD 中のシータの更新について質問があります。theta を更新するには 2 つの方法があります。

1) 前のシータを使用して、すべてのサンプルのすべての仮説を取得し、サンプルごとにシータを更新します。お気に入り：

2) 別の方法: サンプルのスキャン中に、最新のシータを使用して仮説 [i] を更新します。お気に入り：

SGD コードを確認したところ、2 番目の方法が正しいようです。しかし、私のコーディングでは、最初の方が収束が速く、結果は 2 番目よりも優れています。間違った方法が正しい方法よりも優れたパフォーマンスを発揮するのはなぜですか?

また、完成したコードを次のように添付しました。