Вмешательство с дифференцированием

10

При проведении анализа вмешательства с данными временного ряда (также известного как Прерванный временной ряд), как обсуждалось здесь, например, одно из требований, которое у меня есть, - оценить общий выигрыш (или убыток) от вмешательства - то есть количество единиц, полученных или потерянных (переменная Y ).

Не совсем понимая, как оценивать функцию вмешательства с использованием функции фильтра в R, я делал это грубо, надеясь, что это достаточно обобщенно, чтобы работать в любой ситуации.

Скажем так, учитывая данные

 cds<- structure(c(2580L, 2263L, 3679L, 3461L, 3645L, 3716L, 3955L, 
    3362L, 2637L, 2524L, 2084L, 2031L, 2256L, 2401L, 3253L, 2881L, 
    2555L, 2585L, 3015L, 2608L, 3676L, 5763L, 4626L, 3848L, 4523L, 
    4186L, 4070L, 4000L, 3498L), .Dim = c(29L, 1L), .Dimnames = list(
        NULL, "CD"), .Tsp = c(2012, 2014.33333333333, 12), class = "ts")

мы решили, что наилучшая подходящая модель выглядит следующим образом с

гдеXt- импульс в октябре 2013 г.mt=ω0(1δB)XtXt

fit4 <- arimax(log(cds), order = c(1,1,0),include.mean=FALSE, 
               xtransf = data.frame(Oct13 = 1*(seq_along(cds)==22)),
               transfer = list(c(1,0))
               ,xreg=1*(seq_along(cds)==3))
fit4

#    ARIMA(1,1,0)                    

#    Coefficients:
#              ar1    xreg  Oct13-AR1  Oct13-MA0
#          -0.0184  0.2718     0.4295     0.4392
#    s.e.   0.2124  0.1072     0.3589     0.1485

#    sigma^2 estimated as 0.02176:  log likelihood=13.85
#    AIC=-19.71   AICc=-16.98   BIC=-13.05

У меня есть два вопроса:

1) Несмотря на то, что мы дифференцировали ошибки ARIMA, для оценки функции вмешательства, которая затем была технически подобрана с использованием дифференцированного ряда нужно ли что-то сделать, чтобы «изменить» оценку ω 0 или δ с используя X т к X т ?ИксTω0δИксTИксT

2) Правильно ли это: чтобы определить усиление вмешательства, я построил вмешательство из параметров. Получив m t, я сравниваю подгонянные значения из модели fit4 (exp (), чтобы перевернуть журнал) с exp (подогнанные значения минус m t ) и определяю, что в течение наблюдаемого периода вмешательство привело к 3342,37 дополнительных единиц.мTмTмT

Является ли этот процесс правильным для определения выгоды в целом от анализа вмешательства?

    int_vect1<-1*(seq_along(cds)==22)
    wo<- 0.4392
    delta<-0.4295


    mt<-rep(0,length(int_vect1))

    for (i in 1:length(int_vect1))
    {

      if (i>1)
      {
        mt[i]<-wo*int_vect1[i]+delta*mt[i-1]
      }

    }


    mt

sum(exp(fitted(fit4)) - (exp(fitted(fit4) - mt)))
B_Miner
источник
2
Хотите знать, если кто-нибудь может пролить свет на правильный способ оценки воздействия вмешательства - вообще, если процедура, которую я продемонстрировал, на самом деле является правильной?
B_Miner
мTзнак равноехр(α(T-я))еслияTеще0α

Ответы:

4

Предполагая, что это игрушечный пример:

Чтобы ответить на ваш первый вопрос:

1) Несмотря на то, что мы дифференцировали ошибки ARIMA, чтобы оценить функцию вмешательства, которая затем была технически подобрана с использованием разностного ряда, ▽ Xt есть что-то, что нам нужно сделать, чтобы «изменить» оценку ω0 или δ от использования ▽ Xt к Xt?

Когда вы различаете данные, вы должны различать переменные ответа / вмешательства. Когда вы возвращаете разницу (трансформируете) после моделирования, это автоматически учитывает разницу ** Я знаю, что это очень легко, когда вы используете SAS Proc ARIMA. Я не знаю, как это сделать R.

Второй вопрос:

2) Верно ли это: чтобы определить усиление вмешательства, я построил вмешательство mt из параметров. Получив mt, я сравниваю подогнанные значения из модели fit4 (exp (), чтобы перевернуть журнал) с exp (подогнанные значения минус mt) и определяю, что в течение наблюдаемого периода вмешательство привело к 3342,37 дополнительных единиц.

Чтобы определить выигрыш в вмешательстве, вам нужно взять показатель степени и затем вычесть -1, это даст пропорцию или дополнительный эффект. Чтобы продемонстрировать это в вашем случае, смотрите ниже. В первый месяц влияние составило 55% от первоначальных продаж и быстро затухало. В совокупности у вас есть 4580 единиц инкрементального эффекта (с 13 октября по февраль 2014 года. (Я упомянул Delurgio P: Принцип прогнозирования и его применения : 518. Существует превосходная объемная глава по анализу вмешательства).

Кто-нибудь, пожалуйста, исправьте, если эта методология верна?

В этом случае вмешательство импульса + ​​затухание явно недостаточны, я бы сделал сдвиг импульса + ​​постоянный уровень, как показано на диаграмме (e) ниже, из классической статьи Бокса и Тиао .

введите описание изображения здесь

введите описание изображения здесь

предсказатель
источник
Привет @ Forecaster. Как вы получили 3170 в качестве эффекта? Вот что я сделал, я посмотрел на соответствующие значения модели, которая была 8,64245833 (все еще в логарифмическом масштабе). Тогда exp (8,64245833) = 5667,244674. Затем я взял 8,64245833 - 0,4392 = 8,20325833. Так как exp (8.64245833) - exp (8.20325833) = 2014.411599, это эффект. exp (8.64245833) / exp (8.20325833) = 1.55, что показалось мне поддержкой этого.
B_Miner
Я полагаю, вы использовали фактические данные и смоделированный эффект, в отличие от моего подхода, который использовал модель для обоих. Я использовал идею о том, что говорит модель с эффектом и без него. Какой правильный?
B_Miner
Привет @B_miner, с логарифмической шкалой преобразования нам нужно посмотреть на скорость изменения. Подход, который я изложил, - это прямой подход к учебнику, который я цитировал. Тем не менее, ваш подход также разумен. Я сделаю скриншот страницы учебника в ближайшем будущем.
синоптик
Скорость изменения составила 0,55, что также является скоростью изменения в модельном подходе, который я выбрал. Интересно, какой подход более правильный? Я склоняюсь к своему, поскольку подход основан на модели (соответствуют действительным фактам). Если модель очень близка к реальной, два подхода будут выборкой. Я хотел бы видеть страницы. Я вижу, книга выходит из печати?
B_Miner
да, книга вышла из печати. Книжный пример - это постоянное изменение против импульсного вмешательства в вашем примере. Я думаю, что ваш подход прост и точен.
синоптик
0

После того, как @forecaster позволяя Autobox для идентификации 3 выбросов с использованием 29 значений (не в несоответствующих у опыта) полезная модель была найдена введите описание изображения здесьи здесь введите описание изображения здесь. График остаточного акф не предполагает недостаточно определенной модели введите описание изображения здесь. Сюжет Actual / Fit / Forecast введите описание изображения здесьс Fit / Forecast здесь введите описание изображения здесь. Предсказатель (правильно) ранее упомянул, как переменная импульса может превратиться в переменную уровня / шага, когда коэффициент знаменателя равен почти 1,0. При нахождении двух сдвигов уровней (самое последнее, начиная с 9/2013) и импульса в 10/2013, модель дает более четкую картину. С точки зрения воздействия импульса на 10/13 это просто значение коэффициента. НТН

IrishStat
источник
2
На какой из двух вопросов вы ответили?
B_Miner
Первый вопрос предполагал модель, которая предполагала логарифмическое преобразование, которое, я считаю, не оправдано. Пульс на 10/2013 = 1710, что является оценкой эффекта на 10/2013
IrishStat
@B_Miner, можно сказать, что сдвиг уровня на 9/2013 поднял вещи на 1480, таким образом, чистый подъем на 10/2013 будет 1710 + 1480 = 3190
IrishStat
Этот пост, кажется, является скорее расширенным комментарием к вопросу, чем ответом на любую часть вопроса. Возможно, это можно было бы дополнить для непосредственного решения вопросов?
whuber
Ложная предпосылка первого вопроса - это суть моего ответа: на мой взгляд, сомнение / неправильное ведение журнала и включение ненужных различий ... таким образом, мой «ответ» отчасти состоит в том, чтобы исправить предпосылку и предложить эффект воздействия на 10/2013 это просто сумма временных и постоянных изменений. Несмотря на то, что ФП принял другой ответ, я не получил.
IrishStat