我无法弄清楚我无法分解下面的时间序列。有什么想法吗?
data(co2,package = "datasets")
library(lubridate)
library(xts)
plot(co2)#ok
decompose(co2)# ok
my_co2 <- xts(co2,order.by = seq(from=dmy(01011959),to=dmy(01121997),length.out = 468))
plot(my_co2)
my_co2
和co2
似乎相同......
decompose(my_co2) #
分解时出错(my_co2):时间序列没有或少于2个句号
我不明白为什么会这样。有什么想法吗?
编辑:
my_co2_v2 <- xts(frequency = 12,co2,order.by = seq(from=dmy(01011959),to=dmy(01121997),length.out = 468))
frequency(my_co2_v2) # still 1 !
attr(my_co2_v2, 'frequency')<- 12
frequency(my_co2_v2) # ok, it's 12
plot(decompose(my_co2_v2)) # still not good
你是对的plot(decompose(ts(my_co2_v2,frequency = 12)))
工作,但似乎有点扭曲。是否存在关于频率的错误?