Calculate seasonal ACF
seasonalACF(TS, season, lag.max = 50)
TS
season
lag.max
library(data.table) data('precip') seasonalACF(precip, 'month')
Related functions from the same R package