正如 Andy W 建议的那样,最简单的方法可能是使用季节性单变量时间序列模型。如果您使用 R,请尝试auto.arima()
或ets()
从预测包中尝试。
两者都应该可以正常工作,但一般的时间序列方法不会使用提供的所有信息。特别是,您似乎知道每年曲线的形状,因此最好通过相应地对每年的数据进行建模来使用该信息。以下是试图合并此信息的建议。
听起来某种 S 形曲线可以解决问题。例如,一个移位的逻辑:
对于第年和第周,其中、和是要估计的参数。是渐近最大值,控制增长率,时的中点。(将需要另一个参数来允许您描述的不对称性,即到时间的增长率比
ft,j=rteat(j−bt)1+eat(j−bt)
tjatbtrtrtatbtft,j=rt/2btbt. 最简单的方法是允许之前和之后取不同的值。)atbt
可以使用每年的最小二乘法估计参数。每个形成时间序列的参数:,和。这些可以使用标准时间序列方法进行预测,尽管在的情况下,除了使用每个序列的平均值进行预测之外,您可能无能为力。周的值的估计只是,其中使用了、和的预测。a1,…,anb1,…,bnr1,…,rnn=5jf^(6,j)a6b6r6
一旦开始观察第 6 年的数据,您将需要更新此估计值。获得每个新的观测值后,估计从第 6 年开始的数据的 S 形曲线(您至少需要三个观测值才能开始,因为有三个参数)。然后对使用截至第 5 年的数据获得的预测和仅使用第 6 年的数据获得的预测进行加权平均,其中权重分别等于和 . 这是非常临时的,我相信通过将其置于更大的随机模型的上下文中可以使其更加客观。尽管如此,它可能会为您的目的正常工作。(40−t)/36(t−4)/36