fit <- ets(log(tseries1[1:32]), model="AAN")
fit
# 预测结果准确性度量
accuracy(fit)
# 预测时序未来的K步。深色为置信区间为95%的预测范围,浅色为80%。
fit1 <- forecast(fit,4)
fit1
plot(fit1) #这里画的图是取对数之后的预测值,如何还原预测值并画出图?
fit1$mean <- exp(fit1$mean)
fit1$upper <- exp(fit1$upper)
fit1$lower <- exp(fit1$lower)
p <- cbind(fit1$mean, fit1$upper, fit1$lower)
p