ISmKwiatkowskiPhillipsSchmidtShinTest.KPSSStat

Syntax

Description

The KPSSStat property returns the value of the Kwiatkowski-Phillips-Schmidt-Shin statistics.

Comments

To get the value of adjusted residual variance, use the ISmKwiatkowskiPhillipsSchmidtShinTest.F0 property.

Example

To execute the example, add a link to the Stat system assembly.

Sub UserProc;
Var
    KPSS: SmKwiatkowskiPhillipsSchmidtShinTest;
    KPSSStat: IUnitRootTestStatistic;
    can: Array[43] Of Double;
    i, res: Integer;
Begin
    KPSS := New SmKwiatkowskiPhillipsSchmidtShinTest.Create;
    // Set values for variables
    can[0] := 6209; can[1] := 6385; can[2] := 6752; can[3] := Double.Nan; can[4] := 6495;
    can[5] := Double.Nan; can[6] := 7349; can[7] := 7213; can[8] := 7061; can[9] := 7180;
    can[10] := 7132; can[11] := 7137; can[12] := 7473; can[13] := 7722; can[14] := 8088;
    can[15] := 8516; can[16] := 8941; can[17] := 9064; can[18] := 9380; can[19] := 9746;
    can[20] := 9907; can[21] := 10333; can[22] := 10863; can[23] := 11693; can[24] := 12242;
    can[25] := 12227; can[26] := 12910; can[27] := 13049; can[28] := 13384; can[29] := 14036;
    can[30] := 14242; can[31] := 14704; can[32] := 13802; can[33] := 14197; can[34] := 15010;
    can[35] := 15589; can[36] := 15932; can[37] := 16631; can[38] := 17394; can[39] := 17758;
    can[40] := 17308; can[41] := 16444; can[42] := 16413;
    //Select tested series
    KPSS.Serie.Value := can;
    // Type of tested series
    KPSS.TestedSeries := ADFTestedSeriesType.Level;
    // Method of missing data treatment
    KPSS.MissingData.Method := MissingDataMethod.LinTrend;
    // Model type
    KPSS.Equation := EquationType.Constant;
    // Method of calculation of adjusted residuals variance
    KPSS.F0SpectrumEstimation := F0SpectrumEstimationType.BartlettKernel;
    // Autoregression order
    KPSS.AutoRegressionOrder := 9;
    // Sample period
    KPSS.ModelPeriod.FirstPoint := 1;
    KPSS.ModelPeriod.LastPoint := 43;
    res := KPSS.Execute;
    For i := 0 To KPSS.WarningsCount - 1 Do
        Debug.WriteLine(KPSS.Warnings[i]);
    End For;
    KPSSStat := KPSS.KPSSStat;
    If res = 0 Then
        Debug.WriteLine("===The Kwiatkowski-Phillips-Schmidt-Shin statistic===");
Debug.WriteLine("Statistics value: " + KPSSStat.Statistic.ToString);
        Debug.WriteLine("Probability value: " + KPSSStat.Probability.ToString);
        Debug.WriteLine("Critical values");
        Debug.Indent;
        For i := 0 To KPSSStat.CriticalValues.Length - 1 Do
            Debug.Write(i.ToString + " ");
            Debug.WriteLine(KPSSStat.CriticalValues[i]);
        End For;
        Debug.Unindent;
        Debug.WriteLine("Adjusted residuals variance: " + KPSS.F0.ToString);
        Debug.WriteLine("Residuals variance: " + KPSS.S0.ToString);
        Debug.WriteLine("===Auxiliary regression===");
        Debug.WriteLine("Auxiliary regression coefficients:");
        Debug.Indent;
        For i := 0 To KPSS.ModelCoefficients.Estimate.Length - 1 Do
            Debug.WriteLine(KPSS.ModelCoefficients.Estimate[i]);
        End For;
        Debug.Unindent;
        Debug.WriteLine("Descriptive characteristics of auxiliary regression:");
        Debug.Indent;
        Debug.WriteLine("Adjusted determination coefficient: " + KPSS.SummaryStatistics.AdjR2.ToString);
        Debug.Unindent;
        Else
            Debug.WriteLine(KPSS.Errors);
    End If;
End Sub UserProc;

After executing the example the console window shows calculation results for the Kwiatkowski-Phillips-Schmidt-Shin generalized test.