我试图从点数来估计一个迁移期的丰度峰值日期。之前的一项研究使用了以下方法:
将截断的正态分布拟合到每 5 天间隔的捕获密度。使用迭代算法来获得分布平均值的最大似然估计(对应于估计的峰值通道)及其标准误差。
我可以弄清楚如何做大部分事情,但我被困在也许是最简单的问题上。我的数据采用以下格式:日期列和该日期的计数列。
在其他程序中,我可以使用我的分布日期和计数作为分布的频率。这将使我能够将分布拟合到计数中。
但是,我无法弄清楚如何在 R 中执行此操作。fitdistr 或 dnorm 等函数只允许使用向量字符串。
是否有任何函数可以让我拥有一个频率向量和一个感兴趣变量的向量并使用它来拟合分布?
数据如下所示:
头(FR.Counts)
Day.of.Year 计数
112 13780 122 76727 123 55458 125 36746 126 30079 127 20715
112-127 是一年中的一天,其余的是计数。
谢谢