Блог

Главная - Вопросы по программированию - Как интегрировать несколько независимых переменных для анализа временных рядов и сравнить временные ряды в R?

Как интегрировать несколько независимых переменных для анализа временных рядов и сравнить временные ряды в R?

#r #time-series

Вопрос:

Я пытаюсь определить влияние вмешательства (вызывающего отвращение) на количество сообщений, полученных в конкретных общественных парках. Я подумывал об использовании дизайна типа «до-после», но у меня не всегда есть отличные данные «до», и часто бывает много «вмешательств» за короткий промежуток времени на парк. Поэтому я думал об использовании временных рядов и 1) количестве отчетов на парк за 2 недели и 2) количестве вмешательств на парки за 2 недели. До сих пор я использовал тест причинно-следственной связи Грейнджера и функцию автокорреляции, чтобы изучить корреляцию между моим «вмешательством» (AC) и временными рядами моих отчетов.

Как я могу интегрировать компонент даты, компонент типа данных (вмешательство против отчета) и идентификатор парка во временных рядах? И как я должен сравнивать временные ряды?

Вот как выглядят мои данные:

 PARK                  REPORT_AC Date      
  <chr>                     <dbl> <chr>     
1 Nosehill gravel pit           0 13-11-1930
2 Arbour Lake West              0 31-05-2002
3 Arbour Lake West              0 17-06-2002
4 Nosehill Macewin              1 02-05-2018
5 Winston-Victoria Park         1 28-06-2018
6 Winston-Victoria Park         1 29-06-2018

Где отчет = 1, а вмешательство (AC) = 0.

До сих пор мне удавалось построить два моих временных ряда отдельно, используя отчетную дату, только используя следующий метод:

Отделение докладов от мероприятий

 reports_only <- All_years_park_type_date_only[All_years_park_type_date_only$REPORT_AC==1,]

AC_only <- All_years_park_type_date_only[All_years_park_type_date_only$REPORT_AC==0,]

Создание временных рядов для отчетов и мероприятий

 bwts_reports <- ts(reports_only, start=c(2018, 2, 5), end=c(2021, 7, 21), frequency=24)

bwts_AC <- ts(AC_only, start=c(2018, 2, 5), end=c(2021, 7, 21), frequency=24)

ts.plot(bwts_reports, bwts_AC, gpars = list(xlab = "Date", ylab = "Count", 
                                            col=c(1:5)))

Сравнение временных рядов с использованием критерия причинности Грейнджера

 library(lmtest)

grangertest(bwts_reports ~ bwts_AC)

grangertest(bwts_AC ~ bwts_reports)

Автокорреляция

 ccf_reports_ac <- ccf(bwts_AC, bwts_reports, ylab = "cross-correlation")
Метки: