ISmx12arima.OutliersARIMAao

Копировать пример

Sub UserProc;
Var
    cens: ISmX12ARIMA;
    y: Array[62] Of Double;
    AR, MA: Array[1] Of Integer;
    outlier: ICensus2PeriodBegin;
    i, res: Integer;    
Begin
    cens := New SmX12ARIMA.Create;
     // Значения объясняемого ряда

    y[0] :=  284; y[1] := 277;  y[2] := 338;  y[3] := 363; y[4] := Double.Nan; y[5] := 388;
    y[6] :=  407; y[7] := 427;  y[8] := 368;  y[9] := Double.Nan; y[10] := 353; y[11] := 375;
    y[12] := 248; y[13] := 263; y[14] := 322; y[15] := 396; y[16] := 412; y[17] := 451;
    y[18] := 457; y[19] := 457; y[20] := 463; y[21] := 443; y[22] := 411; y[23] := 398;
    y[24] := 335; y[25] := 321; y[26] := 393; y[27] := 418; y[28] := 431; y[29] := 499;
    y[30] := 583; y[31] := 578; y[32] := 497; y[33] := 519; y[34] := 528; y[35] := 508;
    y[36] := 422; y[37] := 439; y[38] := 494; y[39] := 539; y[40] := 559; y[41] := 603;
    y[42] := 627; y[43] := 612; y[44] := 548; y[45] := 480; y[46] := 385; y[47] := 428;
    y[48] := 307; y[49] := 273; y[50] := 334; y[51] := 390; y[52] := 415; y[53] := 453;
    y[54] := 450; y[55] := 479; y[56] := 420; y[57] := 387; y[58] := 381; y[59] := 414;
    y[60] := 335; y[61] := 326;
     /// Общие настройки ///
     //Период расчета

    cens.ModelPeriod.FirstPoint := 1;
    cens.ModelPeriod.LastPoint := 62;
     //Выбор объясняемой переменной
    cens.Serie.Value := y;
     //Кварталы или месяцы
    cens.SeasonalComponentCycleType := SeasonalityCycleType.Month;
    cens.StartPeriod.MonthOrQuarter := 2;
    cens.StartPeriod.Year := 2000;
     //Обработка пропусков
    cens.MissingData.Method := MissingDataMethod.SampleAverage;
     // Порядок авторегрессии
    AR[0] := 2;
    cens.OrderAR := AR;
     // Порядок скользящего среднего
    MA[0] := 1;
    cens.OrderMA := MA;
     // Порядок сезонной авторегрессии
    AR[0] := 2;
    cens.OrderARSeas := AR;
     // Порядок сезонного скользящего среднего

    MA[0] := 1;
    cens.OrderMASeas := MA;
     /// Настройки сезонности ///
     //Тип сезонности
    cens.SeasonalAdjustmentMode := SeasonalityTypeX12.Additive;
     //Сезонный фильтр
    cens.SeasonalFilter := SeasonalFilterType.Auto; 
     /// Поправки на рабочие дни и праздники ///
     // Есть ли поправки
    cens.AdjustmentOptions := AdjustmentOptionsX12Type.X11Step;
     // Использовать поправки, только если необходимо
    cens.AICtest := False;
     // Поправка на рабочие дни
    cens.TradingDayEffects := TradingDayEffectsX12Type.TdStock;
    cens.TdstockValue := 1;
     // Преобразование исходного ряда
    cens.DataTransformation := DataTransformationX12Type.BoxCoxPower;
    cens.BoxCoxPowerTransform := 1;
     /// Выбросы ///

    outlier := cens.OutliersARIMAao.Add;
    outlier.Year := 2000;
    outlier.MonthOrQuarter := 2;
    outlier := cens.OutliersARIMArpbegin.Add;
    outlier.Year := 2002;
    outlier.MonthOrQuarter := 1;
    outlier := cens.OutliersARIMArpend.Add;
    outlier.Year := 2002;
    outlier.MonthOrQuarter := 5;
     // Регрессоры
    cens.PredefinedVariableConst := False;
    cens.PredefinedVariableSeasonal := False;
     /// Диагностика /// 
     // Анализ стабильности сезонных компонент
    cens.StabilityAnalysisofSeasonals := StabilityAnalysisofSeasonalsX12Type.None;
     //Прочее
    cens.ResidualsDiagnostic := False;
    cens.OutliersDetection := True;
     /// Расчет модели ///
    res := cens.Execute;
    Debug.WriteLine(cens.Errors);

    For i := 0 To cens.WarningsCount - 1 Do
        Debug.WriteLine(cens.Warnings[i]);
    End For;
    If (res = 0) Then
        Debug.WriteLine("=== Сезонная составляющая ===");
        Debug.Indent;
        For i := 0 To cens.SeasonallyAdjusted.Length - 1 Do
            Debug.WriteLine(cens.SeasonallyAdjusted[i]);
        End For;
        Debug.Unindent;
        Debug.WriteLine("=== Сезонная корректировка ===");
        Debug.Indent;
        For i := 0 To cens.SeasonalFactors.Length - 1 Do
            Debug.WriteLine(cens.SeasonalFactors[i]);
        End For;
        Debug.Unindent;
        Debug.WriteLine("=== Сезонная компонента ===");
        Debug.Indent;
        For i := 0 To cens.CombinedSeasonalFactors.Length - 1 Do
            Debug.WriteLine(cens.CombinedSeasonalFactors[i]);
        End For;
        Debug.Unindent;
    End If;
End Sub UserProc;