library(data.table)
data('precip')
# \dontshow{
precip <- precip[between(date, as.POSIXct('1990-1-01', format('%Y-%m-%d'), tz = 'America/Regina'),
as.POSIXct('1990-1-10', format('%Y-%m-%d'), tz = 'America/Regina'))]
# }
x <- seasonalACF(precip, 'month')
seasonalAR(precip$date, x)
Run the code above in your browser using DataLab