Лагирование данных о безработице
Поскольку экономические тенденции могут влиять на туризм с некоторой задержкой, перед анализом может быть полезно применить лаг к данным о безработице.
В пакете xts вычислить лаг несложно с помощью команды lag(). Нужно указать данные для сдвига (аргумент x) и значение k, которое задаёт направление и величину лага.
Обратите внимание на разницу в форматировании. В базовом R и пакете zoo лаг задаётся отрицательным значением: лаг 1 записывается как "-1" (а опережение 1 — как "1", что может показаться нелогичным). В пакете xts всё наоборот: лаг задаётся положительным значением, поэтому лаг 1 записывается как "1" (а опережение 1 — как "-1").
Это упражнение является частью курса
Практический пример: анализ городских временных рядов в R
Инструкции к упражнению
- Используйте
lag(), чтобы вычислить одномесячный лаг безработицы в США. Для месячных данных с лагом в один месяц задайте аргументk = 1. Помните, что объектunemploymentсодержит временные ряды как по безработице в США (us), так и по безработице в Массачусетсе (ma). Укажите нужный столбец явно. Сохраните полученный объект xts какus_monthlag. - Ещё раз вызовите
lag(), чтобы вычислить однолетний лаг безработицы в США. Снова укажите нужный столбец вunemploymentи подберите подходящее значениеkдля лага в один год. Сохраните результат какus_yearlag. - С помощью
merge()объедините исходные данные о безработице (unemployment) с новыми лагами (us_monthlagиus_yearlag). Сохраните объединённые данные какunemployment_lags. - Используйте
head(), чтобы вывести первые15строкunemployment_lags.
Интерактивное практическое упражнение
Попробуйте выполнить это упражнение, дополнив этот пример кода.
# Create a one month lag of US unemployment
us_monthlag <- lag(___$___, k = ___)
# Create a one year lag of US unemployment
us_yearlag <-
# Merge your original data with your new lags
unemployment_lags <- merge(unemployment, ___, ___)
# View the first 15 rows of unemployment_lags