Функция R 'STL'. Как сделать сезонный?

Question

Я хочу знать, как сделать «сезонный» в функции STL.

Я думаю, что «тренд» - это временной ряд. и как сделать сезонное ??

в R 'stl' описание равно

which.cycle <- cycle(x)
z$seasonal <- tapply(z$seasonal, which.cycle, mean)[which.cycle]

и R документации

fdrfourier  Calculation of the false discovery rates (FDR) for periodic expression
backgroundData  Generation of background expression set
ar1analysis     Performs AR1 fitting
fourierscore    Calculation of the Fourier score

AR и преобразование Фурье ??

Answer 1

Временной ряд, строго говоря, не может быть сделан сезонным. Я подозреваю, что вы имеете в виду извлечение модели сезонности с использованием функции stl () в R.

Давайте рассмотрим пример временного ряда, измеряющего максимальную зарегистрированную температуру воздуха в Дублине, Ирландия, с 1941 по 2019.

Вот график исходного временного ряда:

weatherarima <- ts(mydata$maxtp, start = c(1941,11), frequency = 12)
plot(weatherarima,type="l",ylab="Temperature in Celsius")
title("Maximum Air Temperature - Dublin")

сезонный , тренд и случайные компоненты могут быть извлечены с помощью stl () следующим образом:

stl_weather = stl(weatherarima, "periodic")
seasonal_stl_weather   <- stl_weather$time.series[,1]
trend_stl_weather     <- stl_weather$time.series[,2]
random_stl_weather  <- stl_weather$time.series[,3]

plot(as.ts(seasonal_stl_weather))
title("Seasonal")
plot(trend_stl_weather)
title("Trend")
plot(random_stl_weather)
title("Random")

Сезонные

Тренд

Случайный

Как можно заметить, в данных о погоде существует явная сезонность - учитывая изменение температуры в разные сезоны - но это было просто извлечено из ряда с помощью stl () - не создан как таковой.

Может оказаться информативным следующее: Извлечение сезонности и тренда из данных: декомпозиция с использованием R * 105 5 *