問題タブ [poisson]

たとえば、X が発生する可能性を表す P(x1...n) 個の離散的な独立した確率値を指定しました。

質問の普遍的なコードが必要です: X が同時に 0 回から n 回発生する確率はどれくらいですか?

例: 与えられた: 各車 (A、B、C) が駐車する 3 つの確率 P(A)、P(B)、P(C)。質問は次のとおりです: 車が 1 台も駐車しない、1 台が駐車する、2 台が駐車する、3 台が駐車する確率はどれくらいですか?

たとえば、同時に 2 台の車を駐車する場合の答えは次のようになります。

私はすべての可能性についてコードを書きましたが、より多くの値を取得すると、可能な組み合わせが増えるため、もちろん遅くなります。

この関数を高速化する方法を知っている人はいますか? それとも、Matlab にその機能が実装されているのでしょうか。

計算の名前を見つけました: ポアソン二項分布

点群から視覚的に魅力的な表面再構成を生成したいと思います。

点群ライブラリを使用しています。ポアソン再構成法を使用してメッシュを作成しようとしましたが、後で完全に再構成できることがわかりました。

例: 私の場合、部屋の点群があります。

http://justpaste.it/code1のコードを使用して、このような再構成を取得できました

_{(出典: pcl-users.org )}

上の写真は、平面図を覆っている表面を持っています。これは、MeshLab を使用して視覚化されました。その後、MeshLab GUI でポイントを押すと、このように表示されます。 _{(出典: pcl-users.org )}

しかし、2番目の写真では、表面にも点があります（添付の写真でははっきりと見えません）.

上部にポイントがなく、内部の構造だけを持つモデルの作成を手伝ってもらえますか?

再構成の品質を改善するための他の提案はありますか?

部屋の点群と生成されたプライ ファイルは、https://dl.dropboxusercontent.com/u/95042389/temp_pcd_ply_files.tar.bz2からダウンロードできます。

時間t ₀ ... t _Nに多くのプロセスx ₀ ... x _Nによって取得されたxの一連の測定値があるとします。時間tで、私が知っている長期的な傾向はなく、xは指数平滑法などのアルゴリズムから予測できるという仮定に基づいて、 xの現在の値を推定したいとします。多くのプロセスがあり、Nが非常に大きくなる可能性があるため、いくつかの値 (以前の状態など) しか保存できません。

ここでの 1 つのアプローチは、通常の指数平滑化アルゴリズムを適応させることです。サンプルが定期的に取得される場合、次のような推定量y _{nを維持します。}

y _n = α . y _n-1 + ( 1 - α ). × _n

サンプリングが不規則な場合、このアプローチは適切ではありません。多くのサンプルを一緒に使用すると、不均衡な影響が生じるからです。したがって、この式は次のように適用できます。

y _n = α _n . y _n-1 + (1 - α _n )。× _n

どこ

α _n = e ^{-k.(t _n - t _n-1 )}

前の 2 つのサンプル間の間隔に応じて平滑化定数を動的に調整する IE。私はこのアプローチに満足しており、うまくいっているようです。これはここで与えられた最初の回答であり、この種の手法の優れた要約は、この 2012 年の論文 ( PDF ) で Eckner によって与えられています。

さて、私の質問は次のとおりです。上記を適応させて、発生率を推定したいと思います。ときどきイベントが発生します。同様の指数手法を使用して、イベントの発生率を推定したいと考えています。

2 つの明白な戦略は次のとおりです。

• データ系列x _nとして最後の 2 つのイベント間の遅延を使用して、最初または 2 番目の手法を使用するには。
• 最後の 2 つのイベント間の遅延の逆数 (つまり、速度) をデータ系列x _nとして使用して、最初または 2 番目の手法を使用します。

私が知る限り、これらはどちらも良い戦略ではありません。まず、500 ミリ秒ごとに発生するイベント (一方で) と、他方で 200 ミリ秒の遅延と 800 ミリ秒の遅延で発生するイベントを取り上げます。明らかに、これらは両方とも 1 秒間に 2 回発生するため、指定されたレート推定値は同じになるはずです。最後のサンプルの時間を無視するのは無謀に思えるので、2 番目の戦略に集中します。200 ミリ秒/800 ミリ秒のサンプル ストリームをシミュレートすると、約 1.5 の推定値が生成されるため (逆数ではなく) 遅延を使用しても、適切な予測因子にはなりません (逆数の平均は平均の逆数ではないことに基づく)。

しかし、はるかに重要なことは、どちらの戦略も、実際に実際に起こっていること、つまり突然すべてのイベントが長時間停止することに対応していないことです。したがって、 yの「最新」の値は最後のイベントの値であり、レートの推定値は計算されません。したがって、レートは一定に見えます。もちろん、データを遡及的に分析していれば問題ありませんが、リアルタイムで分析しています。

これを行う別の方法は、スレッドを定期的に (たとえば 10 秒ごとに) 実行し、この 10 秒間隔で発生回数をカウントすることです。統計は頻繁に必要とされないため、これは非常にリソースが多く、ミューテックスの問題のためにすべてをポーリングするスレッドを実行するのは嫌いです。したがって、（どういうわけか）最後のサンプルが取得されてからの時間（たとえば）によって読み取られた状態を調整するアルゴリズムを使用したいと思います。これは、パフォーマンスがサンプルとは無関係に選択された時間に測定される場合、合理的なアプローチのように思われます。測定時間は平均してサンプル間の期間の半分になるため、レートの非常に粗い平滑化されていない推定値は、最後のサンプルからの時間。さらに複雑なことに、測定時間はサンプルとは無関係ではありません。

これには簡単な答えがあると感じていますが、私にはわかりません。イベントがポアソン分布であると仮定し、最後のサンプルからの間隔と何らかの形の移動平均に基づいてλの推定値を導き出すのが正しいルートだと思いますが、私の統計はさびすぎてこれを機能させることができません。

ここにこの質問のほとんどだましがありますが、答えはあまり満足のいくものではないようです（理由を説明したいと思います）。推定する変数が 1 つあり、それについて何も知らないことを考えると、カルマン フィルターは重く見えると付け加えておきます。他にも多くの似たようなものがありますが、そのほとんどは値の大きなビンを保持することを提案するか (ここではメモリの観点からは現実的ではありません)、上記の 2 つの問題に対処していません。

私は 500 の平均値で poisson_distribution を使用します。問題は、1200 回の反復の後、分布が 423-571 の間にあることです。たとえば、結果が 1 から 4.000.000 になるように最大制限を変更するにはどうすればよいですか?ただし、平均値は500?

これが私のコードです：

ポアソン回帰を行い、モデルを視覚化しました。

次に、これらのシミュレートされたポアソン分布をプロットに追加します。

プロットは次のようになります (赤い点はフォトショップ処理されています)。

最適化中にのみ問題が発生するネストされたモデルのより大きなセットをデバッグしています。rpois()エラーの原因と思われるものをゼロにするプロセス中に、関数で異常な動作に遭遇しました。

平均値が非常に大きい場合、値の代わりにrpois()返されるようです。NAこの問題は警告を生成しません。再現可能なコード セットについては、以下を参照してください。

私の質問は 2 つあります:
1 - なぜこの動作を示しているのですか (rpois 関数の整数のサイズに最大制限はありますか?)、
2 - NA の生成を防ぐための回避策はありますか (それが平均入力のサイズをより小さい値に制限するには)?

64x Windows 7 で 32x R バージョン 3.0.2 を実行しています。