ポイントカウントから、移動期の個体数がピークに達した日付を推定しようとしています。以前の研究では、次の方法が使用されています。
切り捨てられた正規分布を 5 日間隔ごとのトラップ密度に当てはめました。反復アルゴリズムを使用して、分布の平均の最尤推定値 (推定されたピーク通過に対応) とその標準誤差を取得しました。
これのほとんどを行う方法を考え出すことはできますが、おそらく最も単純な問題に行き詰まっています。私のデータは次の形式です: 日付の列とその日付のカウントの列。
他のプログラムでは、配布の日付とカウントを配布の頻度として使用できます。これにより、分布をカウントに合わせることができます。
ただし、Rでこれを行う方法がわかりません。fitdistr や dnorm などの関数では、ベクトル文字列のみが許可されます。
頻度のベクトルと対象変数のベクトルを取得し、それを使用して分布に適合させる関数はありますか?
データは次のようになります。
head(FR.Counts)
年間通算日
112 13780 122 76727 123 55458 125 36746 126 30079 127 20715
112 ～ 127 は通年日で、残りはカウントです。
ありがとう