Load data, assume Alpha is 0.33, forecast 20 points past the last value.
Uses MtxExpr, StatTimeSerAnalysis, Math387; procedure Example; var Data,YHat: Vector; T: Integer; NumPoints: Integer; Residuals: TVec; begin NumPoints := 20; Data.LoadFromFile('aerosol_particles.vec'); // last point period = Data.Length-1 + NumPoints T := Data.Length-1+NumPoints; SingleExpForecast(Data,YHat,0.33,T,0); // YHat now stores estimates for YHat[1,...Length-1] // so, if we need residuals, we have to subtract // these values from y[1,...,Length-1) Residuals.Size(YHat); Residuals.Sub(Y,YHat,1,0,0,YHat.Length); end;
#include "MtxExpr.hpp" #include "Math387.hpp" #include "StatTimeSerAnalysis.hpp" void __fastcall Example(); { sVector Data,YHat, Residuals; int NumPoints = 20; Data.LoadFromFile("aerosol_particles.vec"); // last point period = Data.Length-1 + NumPoints int T = Data.Length-1+NumPoints; SingleExpForecast(Data,YHat,0.33,T,0); // YHat now stores estimates for YHat[1,...Length-1] // so, if we need residuals, we have to subtract // these values from y[1,...,Length-1) Residuals.Size(YHat); Residuals.Sub(Y,YHat,1,0,0,YHat.Length); }
Copyright (c) 1999-2025 by Dew Research. All rights reserved.
|