I have similar data set as in Convert multiple irregular time series into regular time series
station Time WaterTemp
1 01-01-1974 5.0000000
1 01-02-1974 5.0000000
1 01-03-1974 8.6000004
1 01-05-1974 8.1333332
1 01-07-1974 12.7999999
2 01-01-1974 5.0000000
2 01-02-1974 5.0000000
2 01-04-1974 8.6000004
2 01-06-1974 8.1333332
2 01-08-1974 12.7999999
If I ran single arima, I can use ACF and PACF plots to determine one station of p,d,q. But what if I have 50 stations univariable time series. How can I get the best p,d,q for each station and do the forecast and validation or RMSE?