r - 在时间序列中使用傅里叶项来表征季节性

Question

我有一个名为“bicoal.tons”的时间序列，其中包含以下信息：

Time Series:
Start = 1920 
End = 1968 
Frequency = 1 
 [1] 569 416 422 565 484 520 573 518 501 505 468 382 310 334 359 372 439 446 349 395
[21] 461 511 583 590 620 578 534 631 600 438 516 534 467 457 392 467 500 493 410 412
[41] 416 403 422 459 467 512 534 552 545

我想从这个年度时间序列中观察季节性影响，我从这段代码开始fourier(time_series,2,h=NULL)

但随后我收到此错误消息：

Error in ...fourier(x, K, 1:NROW(x)) : 
  K must be not be greater than period/2

我做错了什么，我将如何解决这个问题以便我可以获得一个情节？

标签： r