1

我试图从点数来估计一个迁移期的丰度峰值日期。之前的一项研究使用了以下方法:

将截断的正态分布拟合到每 5 天间隔的捕获密度。使用迭代算法来获得分布平均值的最大似然估计(对应于估计的峰值通道)及其标准误差。

我可以弄清楚如何做大部分事情,但我被困在也许是最简单的问题上。我的数据采用以下格式:日期列和该日期的计数列。

在其他程序中,我可以使用我的分布日期和计数作为分布的频率。这将使我能够将分布拟合到计数中。

但是,我无法弄清楚如何在 R 中执行此操作。fitdistr 或 dnorm 等函数只允许使用向量字符串。

是否有任何函数可以让我拥有一个频率向量和一个感兴趣变量的向量并使用它来拟合分布?

数据如下所示:

头(FR.Counts)

Day.of.Year 计数

     112     13780

     122     76727

     123     55458

     125     36746

     126     30079

     127     20715

112-127 是一年中的一天,其余的是计数。

谢谢

4

0 回答 0