問題タブ [poisson]
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.
r - R のポアソン テーブル
Rで 2 つのイベントのポアソン テーブルを生成しようとしています。1 つは平均 1.5 (ラムダ 1) で、もう 1 つは平均 1.25 (ラムダ 2) です。x=0 から x=7+ (7 以上) の両方の場合の確率を生成したいと思います。これはおそらく非常に簡単ですが、その方法を理解できないようです! テーブルのデータ フレームを作成できましたが、これまで関数を作成したことがないため、パラメーターの入力方法がよくわかりません。
私はいくつかの「For」ループが必要になりdpois(x,lambda1)
、ある時点で関与すると推測しています。誰か助けてくれませんか?
r - Conway-Maxwell ポアソン回帰の R パッケージ COMPoissonReg の問題
cmp
パッケージの関数を使用して、COMPoissonReg
非常に分散が不十分な保険死亡率データで Conway-Maxwell Poisson 回帰を実行しようとしていますが、次のエラーが発生します。
関数は機能し始めたようで、係数の推定値を生成しますが、失敗してエラーが発生します。
私のコードは次の形式です
同じデータと予測子に対して標準のポアソン回帰を実行しましたが、すべて正常に実行されました。このエラーの原因は何ですか?
random - ポアソン分布から乱数を生成する
M = 4
平均がで範囲がのポアソン分布から 60 個の乱数を生成したいと考えていますmin = 2
max = 9
。これを解決する方法を知っている人はいますか(できればMatlabで?)。
一番!
c - 一様分布からポアソン分布への変換
乱数の一様分布 (0 から 1 など) をポアソン分布に変換する C プログラムを作成する必要があります。誰でも助けることができますか?
r - ランダムに生成されたサンプルのヒストグラムに真のポアソン分布の pmf を重ねる
私は統計とRの両方に不慣れです.次を使用して、ラムダ= 2.5のポアソン分布から100個のランダムサンプルを生成しました:
サンプルの相対頻度ヒストグラムを正常に作成しました。
ここで、真のポアソン分布の pmf をヒストグラムに重ねる必要がありますが、真の関数を生成する方法がわかりません。これはおそらく非常に単純だと思いますが、何をすべきかわかりません。どんな助けでも本当に感謝しています。
ありがとう!
r - サッカースコアのモデリング (Dixon & Cole)
Dixon, Coles ( 1997 ) では、(4.3) の 2 つの修正独立ポアソン モデルの最尤推定を使用して、サッカーのスコアをモデル化しました。
パッケージを使用せずに、アルファとベータ、およびホーム効果パラメーター (pg. 274、表 4) を「再現」するために R を使用しようとしています (通常の独立したポアソン モデルを使用しても問題ありません)。bivpois パッケージを使用してみましたが、そのパラメーターを変更する方法がわかりません。
誰かがデータをモデル化するための R コードを手伝ってくれたら大変ありがたいです - イングランド プレミア リーグの 2012/13 シーズンのホーム アンド アウェイ チームのスコア。したがって、基本的に、R の最適関数を使用して式 4.3 または 4.5 をコーディングする際に助けが必要です。
ポアソン分布の密度 (1 つの独立したポアソン モデルのみ)
通常のポアソン尤度関数
アルファとベータの適切な見積もりを見つけるための最適化。
アルファ、ベータ、ホーム効果の 3 つのパラメーターの推定値を見つける際にエラーが発生しました。コードを改善する方法を教えてもらえますか?
データは次の形式です。
python - 画像へのポアソン ノイズの追加
より徹底的に分析するために、ポアソン ノイズの増分量を追加する必要がある画像がいくつかあります。MATLAB でこれを行うことができることは知っていますが、Python でそれを行うにはどうすればよいでしょうか? 検索しても、これまでのところ何も得られていません。
r - family="poisson" の lambda=0 の場合の glmnet の収束の欠如
glmnet と glm を比較しているときに、lambda=0 と family="poisson" の収束の問題に遭遇しました。私の理解では、ラムダ= 0(およびデフォルトのアルファ= 1)では、答えは本質的に同じであるはずです。
以下は、glmnet ヘルプ ページ (?glmnet) の poisson の例からわずかに変更されたコードです。唯一の変更点は、nzc = p であるため、すべての変数が真のモデルにあることです。
ここにエラーメッセージがあります
更新: 問題は、family="poisson" の場合に glmnet によって推定されるインターセプトにあり、ラムダ自体の設定には関係ないようです。