ARIMA Procedure

Initial Values

The syntax for giving initial values to transfer function parameters in the INITVAL= option parallels the syntax of the INPUT= option. For each transfer function in the INPUT= option, the INITVAL= option should give an initialization specification followed by the input series name. The initialization specification for each transfer function has the form

upper C normal dollar-sign left-parenthesis upper V Subscript 1 comma 1 Baseline comma upper V Subscript 1 comma 2 Baseline comma ellipsis right-parenthesis left-parenthesis upper V Subscript 2 comma 1 Baseline comma ellipsis right-parenthesis ellipsis slash left-parenthesis upper V Subscript i comma 1 Baseline comma ellipsis right-parenthesis ellipsis

where C is the lag 0 term in the first numerator factor of the transfer function (or the overall scale factor if the ALTPARM option is specified) and upper V Subscript i comma j is the coefficient of the upper L Subscript i comma j element in the transfer function.

To illustrate, suppose you want to fit the model

upper Y Subscript t Baseline equals mu plus StartFraction left-parenthesis omega 0 minus omega 1 upper B minus omega 2 upper B squared right-parenthesis Over left-parenthesis 1 minus delta 1 upper B minus delta 2 upper B squared minus delta 3 upper B cubed right-parenthesis EndFraction upper X Subscript t minus 3 Baseline plus StartFraction 1 Over left-parenthesis 1 minus phi 1 upper B minus phi 2 upper B cubed right-parenthesis EndFraction a Subscript t

and start the estimation process with the initial values mu=10, omega 0=1, omega 1=0.5, omega 2=0.03, delta 1=0.8,
delta 2=–0.1, delta 3=0.002, phi 1=0.1, phi 2=0.01. (These are arbitrary values for illustration only.) You would use the following statements:

identify var=y crosscorr=x;
estimate p=(1,3) input=(3$(1,2)/(1,2,3)x)
         mu=10 ar=.1 .01
         initval=(1$(.5,.03)/(.8,-.1,.002)x);

Note that the lags specified for a particular factor are sorted, so initial values should be given in sorted order. For example, if the P= option had been entered as P=(3,1) instead of P=(1,3), the model would be the same and so would the AR= option. Sorting is done within all factors, including transfer function factors, so initial values should always be given in order of increasing lags.

Here is another illustration, showing initialization for a factored model with multiple inputs. The model is

StartLayout 1st Row 1st Column upper Y Subscript t Baseline equals mu 2nd Column plus 3rd Column StartFraction omega Subscript 1 comma 0 Baseline Over left-parenthesis 1 minus delta Subscript 1 comma 1 Baseline upper B right-parenthesis EndFraction upper W Subscript t plus left-parenthesis omega Subscript 2 comma 0 Baseline minus omega Subscript 2 comma 1 Baseline upper B right-parenthesis upper X Subscript t minus 3 2nd Row 1st Column Blank 2nd Column plus 3rd Column StartFraction 1 Over left-parenthesis 1 minus phi 1 upper B right-parenthesis left-parenthesis 1 minus phi 2 upper B Superscript 6 Baseline minus phi 3 upper B Superscript 12 Baseline right-parenthesis EndFraction a Subscript t EndLayout

and the initial values are mu=10, omega Subscript 1 comma 0=5, delta Subscript 1 comma 1=0.8, omega Subscript 2 comma 0=1, omega Subscript 2 comma 1=0.5, phi 1=0.1, phi 2=0.05, and phi 3=0.01. You would use the following statements:

identify var=y crosscorr=(w x);
estimate p=(1)(6,12) input=(/(1)w, 3$(1)x)
         mu=10 ar=.1 .05 .01
         initval=(5$/(.8)w 1$(.5)x);
Last updated: June 19, 2025