ARIMA中单个实例中的多个产品的预测

Question

我有一些CSV格式的SKU级别的24个月的销售记录（该格式粘贴在下面）。我想在ARIMA中通过CSV上传在单个实例中对这多个产品执行3个月的预测。 编写的代码（摘自Hyndman教授的帖子http://robjhyndman.com/hyndsight/batch-forecasting/。）如下：

library(forecast)
print(system.time({ 
sales <- read.csv("C:\\Users\\Red Eye\\Desktop\\ARIMA\\Sales History1.csv", header = FALSE)
sales <- ts(sales[,-1], f = 12, s = c(2016,4))
  h <- 3
fit <- apply(sales, 2, auto.arima)
  fcast <- lapply(fit, forecast, h = h)
fcast <- sapply(fcast,"[",4)
fcast <- (as.data.frame(fcast))
write.csv(t(fcast), file = "Sales Forecast.csv")  
}))

问题在于预测导致Sales Forecast.csv文件仅具有最后一个SKU。 输入的CSV文件（分别具有24个月销售数据的SKU1，SKU2和SKU3）的快照如下所示。

输出预测CSV文件仅包含SKU3预测值....

任何人都可以帮助我如何通过R中的ARIMA对单个实例中的所有SKU进行预测。在此先感谢