Para bucle para pronosticar varios conjuntos de datos a la vez en R

Tengo un conjunto de datos con variables "Hora, Región, Ventas" y quiero pronosticar las ventas para cada región usando ARIMA o ETS (SES) usandolibrary(forecast). Hay un total de 70 regiones y todas tienen 152 observaciones cada una y (3 años de datos). Algo como esto:

  Week      Region    Sales 
01/1/2011      A       129
07/1/2011      A       140
14/1/2011      A       133
21/1/2011      A       189
...           ...      ...
01/12/2013     Z       324
07/12/2013     Z       210
14/12/2013     Z       155
21/12/2013     Z       386
28/12/2013     Z       266 

Entonces, quiero que R trate a cada región como un conjunto de datos diferente y realice unauto.arima. Supongo que un bucle for debería ser un ajuste ideal aquí, pero fracasé miserablemente con él. Lo que idealmente querría que hiciera es un bucle for para ejecutar algo como esto (una arima automática por cada 152 observaciones):

fit.A <- auto.arima(data$Sales[1:152])  
fit.B <- auto.arima(data$Sales[153:304])
....
fit.Z <- auto.arima(data$Sales[10490:10640])

Me encontré conesta pero al convertir el marco de datos en series de tiempo, todo lo que obtuve son NA.

Cualquier ayuda es apreciada! Gracias.