問題タブ [stochastic]
For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.
python - pymc の確率変数
私はそのような用語に出くわしました
多くの場所で。
私はどこにもそれらについて見つけることができませんでした。私は彼らの使用法しか見ることができません。
それらは何を意味し、uniform, binomialそれぞれどのように異なるのですか
r - 浴槽内の確率論的ビー玉をシミュレートします。e
私は生物学の大学院生で、特定の動作を R のモデルにコーディングしようとしていますが、「翻訳で失われた」問題がいくつかあります。私が持っているコードは投稿に従います。私はこのシステムをモデル化しようとしています:
大理石が振動するバスタブを想像してみてください。バスタブの表面は関数「bathtub」で与えられます。「温度」パラメータを指定して、浴槽のビー玉の振動をシミュレートします。低温/ゼロ温度では、底に置き、温度としてください。増加すると、側面をより高く探索する必要があります。これらをベクトルに保存し、この確率的プロセスの特定の実現がどのように見えるかを確認できるようにしたいと考えています。
私の問題は、大理石が動いている表面を記述する関数ではなく、開始値、時間、およびその他のパラメーター (曲率パラメーター a と b、温度パラメーター、およびおそらく他のパラメーター)。
私は基本的に正規分布のためにこのコードに類似したものを探しています:
アイデア、コード、または役立つリソースへのリンクをありがとうございます。
graph - 絶対加重による加重グラフでのページランクの計算
このリンク ( Networkx PageRank - Equal Ranks with Different Weights )で表現されているのと同じ問題に直面しています。
基本的に、グラフのページランクを計算するために networkx を使用しています。したがって、ページランクの計算では、最初にグラフが右確率行列に変換されます (すべての出力エッジが 1 に正規化されます)。
私が必要としているのは、エッジの重みを正規化しない方法です。したがって、1 つのノードが重み 0.1 の 1 つの発信エッジのみであり、別のノードが重み 0.05 の発信エッジを 1 つしか持たない場合、この情報をページランクの計算に使用する必要があります (それぞれ 1 に正規化するのではなく)。
これを達成するためにページランクを変更する正しい方法を知っている人はいますか?
前もってありがとう、アミット
r - %dopar% または逐次確率計算を高速化する代替方法
確率過程シミュレーターを作成しましたが、かなり遅いので高速化したいと思います。
シミュレーターの主要部分はforループで構成されていますが、これforeachを `%dopar%.
単純化されたループでそうしようとしましたが、いくつかの問題が発生しています。私のforループがこのように見えるとします
これは、反復ごとにrandの値を更新しs、2 つの結果のいずれかで vectorを設定することを意味しtます。ループとして書き直すいくつかの異なる方法を試しましたforeachが、繰り返しのたびに値が更新されず、かなり奇妙な結果が得られるようです。使ってみreturnたけど効果なし！
これは私が思いついた例です。
使用foreachできない場合、すべてのコアを使用して高速化できるようにループを書き直すには、他にどのような方法がありますか?
matlab - 確率的時間依存入力を持つ ODE
論文で見つけた例を再現しようとしています。
この ODE を解決する必要があります。
25 a + 15 v + 330000 x = p(t)
ここで、p(t)は 10 ～ 25 Hz の範囲に帯域制限されたホワイト ノイズ シーケンスです。aは加速度、vは速度、xは変位です。
次に、システムはルンゲクッタ手順を使用してシミュレートされていると表示されます。サンプリング周波数は 1000 Hz に設定され、ノイズが信号 rms 値の 5% に寄与するようにガウス ホワイト ノイズがデータに追加されます (この最後の情報をどのように使用すればよいでしょうか?)。
主な問題は、帯域制限されたホワイト ノイズに関連しています。ここで見つけた指示に従いましたhttps://dsp.stackexchange.com/questions/9654/how-to-generate-band-limited-gaussian-white-noiseそして、次のコードを書きました：
ここで、頌歌の関数を定義する必要がありますが、 p (ホワイト ノイズ)を与える方法がわかりません...次の方法を試しました。
実際のところ、入力励起は適切に機能していません。方程式の固有周波数は約 14 Hz であり、ホワイト ノイズが 10 ～ 25 Hz の範囲にあるため、応答に共振が見られると予想されます。
この Q/A も見ましたが、まだ機能させることはできません。