Calculate seasonal ACF
seasonalACF(TS, season, lag.max = 50)
TS
time series
season
name of the season
lag.max
max lag for acf
data('precip') seasonalACF(precip, 'month')