問題タブ [stochastic-process]

私はPyMCを初めて使用し、最大事後推定値を使用して、不均一なポアソンプロセスを区分定数レート関数に適合させようとしています。

私のプロセスは、1 日の間にいくつかのイベントを記述します。したがって、私は 1 日を 24 時間に分割しています。つまり、レート関数 (区分定数) 内に 24 の定数があることを意味します。

以下のアイデアを組み合わせる:

• ハッカーのための確率的プログラミングとベイズ法: 第 1 章
• SO の質問: キャップ付きポアソン過程を可変レートでフィッティングする

次のコードを思いつきましたが、満足のいくものではありません(結果的には間違っていると確信しています):

a0、a1... の値が私のデータに適合していないようです (指定されたラムダを使用した非均質ポアソン過程からのサンプリングによって生成されました -> テストケース!)

ラムダをどのように適合/推定できますか? 私は何を間違っていますか？

(私は pyMC 2.3.2 を使用しています!)

jmeter を使用して、確率的パフォーマンス テストの実装を簡素化しようとしています。以前にこれを行ったことがありますが、時間がかかります。

問題は、(合計訪問数と経過時間に基づいて) 同時ユーザーの数があれば、経過時間は変わらないと仮定します。しかし、ストレスの多い状況でサーバーが予想以上に処理を行った場合、一部のユーザーは予想された時間内にサーバーを離れないため、ユーザーは積み重なっていきます。したがって、実際の一致は計算されたものとは異なります。

従来のテストでは、サーバーの応答が遅いために経過時間が増加しても、事前に確立された仮想ユーザーが維持されます (増加するかどうかは問題ではありません)。サーバーの応答が遅い場合、jmeter の反復回数は少なくなります。私はそれをしたくありません。到着率を維持するために、必要に応じてさらにスレッドを作成したいと考えています。

時間の経過とともに一定のスレッドではなく、ユーザーの到着率を設定する方法はありますか?

一連の方程式を解いて、エルゴード マルコフ行列の定常分布を決定しようとしています。

つまり、行列は

方程式のセットは、以下の定理からのものです

上記の方程式を有効な Matlab 構文に変換するにはどうすればよいですか?

私は生物学の大学院生で、特定の動作を R のモデルにコーディングしようとしていますが、「翻訳で失われた」問題がいくつかあります。私が持っているコードは投稿に従います。私はこのシステムをモデル化しようとしています:

大理石が振動するバスタブを想像してみてください。バスタブの表面は関数「bathtub」で与えられます。「温度」パラメータを指定して、浴槽のビー玉の振動をシミュレートします。低温/ゼロ温度では、底に置き、温度としてください。増加すると、側面をより高く探索する必要があります。これらをベクトルに保存し、この確率的プロセスの特定の実現がどのように見えるかを確認できるようにしたいと考えています。

私の問題は、大理石が動いている表面を記述する関数ではなく、開始値、時間、およびその他のパラメーター (曲率パラメーター a と b、温度パラメーター、およびおそらく他のパラメーター)。

私は基本的に正規分布のためにこのコードに類似したものを探しています:

アイデア、コード、または役立つリソースへのリンクをありがとうございます。

Sim.DiffProcStratonovich の確率積分をシミュレートするために Rのパッケージを使用しています。次のコードを使用して、t=0 から t=5 までの確率積分の 5 つのパスをシミュレートできます。

fun=expression(w) strat=st.int(fun, type="str", M=5, lower=0, upper=5)

st.int()関数がさまざまな t の値を出力として与えない場合、t=5 の確率積分の値を取得するにはどうすればよいですか?

隣接する 1 を持たないすべての NxN 行列に均一に分布するマルコフ連鎖の MC シミュレーターを実行しようとしています。私のアルゴは、チェーンを何度も実行して状態空間を埋めることになっています。しかし、私のロジックのどこかにひどく問題があり、状態空間がいっぱいになっていません。どんな助けでも大歓迎です。これが私のコードです。

2 人がテーブルに座り、3 冊の本を共有するシステムを考えてみましょう。いつでも両方とも本を読んでいて、テーブルには 1 本の本が残されています。人は今読んでいる本を読み終えると、テーブルの上にある本と入れ替えて読み始めます。読書時間は指数関数的に分布し、人 i が本 j を読む平均時間を bi,j で表します。

このマルコフ連鎖の状態空間とは何ですか?また、どのようにレート行列 Q を構築できますか?

講義ノートからこの演習を取得しましたが、連続時間マルコフ連鎖であるため、どういうわけか状態空間が混乱しています。

これらは私が考えることができる可能な状態です:

人 i1 と i2、本 A、B、C

しかし、どうすればこれをグラフィカルに表現できますか? 試してみましたが、各ユーザーには個別のマルコフ連鎖（還元可能）があり、これは正しいとは思えません。そこから、マトリックス b のレートに基づいてレート マトリックスを作成しても問題ないと思います。