hide long namesshow long names
hide short namesshow short names
Integer type:  int32  int64  nag_int  show int32  show int32  show int64  show int64  show nag_int  show nag_int

PDF version (NAG web site, 64-bit version, 64-bit version)
Chapter Contents
Chapter Introduction
NAG Toolbox

NAG Toolbox: nag_tsa_multi_inputmod_estim (g13be)

Purpose

nag_tsa_multi_inputmod_estim (g13be) fits a multi-input model relating one output series to the input series with a choice of three different estimation criteria: nonlinear least squares, exact likelihood and marginal likelihood. When no input series are present, nag_tsa_multi_inputmod_estim (g13be) fits a univariate ARIMA model.

Syntax

[para, xxy, zsp, itc, sd, cm, s, d, ndf, res, sttf, nsttf, ifail] = g13be(mr, mt, para, xxy, 'nser', nser, 'npara', npara, 'kfc', kfc, 'nxxy', nxxy, 'kef', kef, 'nit', nit, 'zsp', zsp, 'iwa', iwa, 'imwa', imwa, 'kpriv', kpriv)
[para, xxy, zsp, itc, sd, cm, s, d, ndf, res, sttf, nsttf, ifail] = nag_tsa_multi_inputmod_estim(mr, mt, para, xxy, 'nser', nser, 'npara', npara, 'kfc', kfc, 'nxxy', nxxy, 'kef', kef, 'nit', nit, 'zsp', zsp, 'iwa', iwa, 'imwa', imwa, 'kpriv', kpriv)
Note: the interface to this routine has changed since earlier releases of the toolbox:
Mark 22: nxxy has been made optional
Mark 23: isttf, kzef, kzsp have been dropped, iwa, imwa, zsp, kfc, nit, kef, kpriv now optional
.

Description

The Multi-input Model

The output series ytyt, for t = 1,2,,nt=1,2,,n, is assumed to be the sum of (unobserved) components zi,tzi,t which are due respectively to the inputs xi,txi,t, for i = 1,2,,mi=1,2,,m.
Thus yt = z1,t + + zm,t + ntyt=z1,t++zm,t+nt where ntnt is the error, or output noise component.
A typical component ztzt may be either
(a) a simple regression component, zt = ωxtzt=ωxt (here xtxt is called a simple input), or
(b) a transfer function model component which allows for the effect of lagged values of the variable, related to xtxt by
zt = δ1 zt1 + δ2 zt2 + + δp ztp + ω0 xtb ω1 xtb1 ωq xtbq .
zt = δ1 zt-1 + δ2 zt-2 + + δp zt-p + ω0 xt-b - ω1 x t-b-1 - - ωq x t-b-q .
The noise ntnt is assumed to follow a (possibly seasonal) ARIMA model, i.e., may be represented in terms of an uncorrelated series, atat, by the hierarchy of equations
(i) d sD nt = c + wt d sD nt = c+wt
(ii) wt = Φ1 wts + Φ2 wt2 × s + + ΦP wtP × s + et Θ1 ets Θ2 et2 × s ΘQ etQ × s wt = Φ1 wt-s + Φ2 w t-2×s ++ ΦP w t-P×s + et - Θ1 et-s - Θ2 et-2×s -- ΘQ et-Q×s
(iii) et = φ1 et1 + φ2 et2 + + φp etp + at θ1 at1 θ2 at2 θq atq et = ϕ1 et-1 + ϕ2 e t-2 ++ ϕp et-p +at - θ1 at-1 - θ2 at-2 -- θq at-q
as outlined in Section [Description] in (g13ae).
Note:  the orders p,qp,q appearing in each of the transfer function models and the ARIMA model are not necessarily the same; dsDntdsDnt is the result of applying non-seasonal differencing of order dd and seasonal differencing of seasonality ss and order DD to the series ntnt: the differenced series is then of length N = nds × DN=n-d-s×D; the constant term parameter cc may optionally be held fixed at its initial value (usually, but not necessarily zero) rather than being estimated.
For the purpose of defining an estimation criterion it is assumed that the series atat is a sequence of independent Normal variates having mean 00 and variance σa2σa2. An allowance has to be made for the effects of unobserved data prior to the observation period. For the noise component an allowance is always made using a form of backforecasting.
For each transfer function input, you have to decide what values are to be assumed for the pre-period terms z0,z1,,z1pz0,z-1,,z1-p and x0,x1,,x1bqx0,x-1,,x1-b-q which are in theory necessary to re-create the component series z1,z2,,znz1,z2,,zn, during the estimation procedure.
The first choice is to assume that all these values are zero. In this case, in order to avoid undesirable transient distortion of the early values z1,z2,z1,z2,, you are advised first to correct the input series xtxt by subtracting from all the terms a suitable constant to make the early values x1,x2,x1,x2,, close to zero. The series mean xx- is one possibility, but for a series with strong trend the constant might be simply x1x1.
The second choice is to treat the unknown pre-period terms as nuisance parameters and estimate them along with the other parameters. This choice should be used with caution. For example, if p = 1p=1 and b = q = 0b=q=0, it is equivalent to fitting to the data a decaying geometric curve of the form AδtAδt, for t = 1,2,t=1,2,, along with the other inputs, this being the form of the transient. If the output ytyt contains a strong trend of this form, which is not otherwise represented in the model, it will have a tendency to influence the estimate of δδ away from the value appropriate to the transfer function model.
In most applications the first choice should be adequate, with the option possibly being used as a refinement at the end of the modelling process. The number of nuisance parameters is then max (p,b + q)max(p,b+q), with a corresponding loss of degrees of freedom in the residuals. If you align the input xtxt with the output by using in its place the shifted series xtbxt-b, then setting b = 0b=0 in the transfer function model, there is some improvement in efficiency. On some occasions when the model contains two or more inputs, each with estimation of pre-period nuisance parameters, these parameters may be co-linear and lead to failure of the function. The option must then be ‘switched off’ for one or more inputs.

The Estimation Criterion

This is a measure of how well a proposed set of parameters in the transfer function and noise ARIMA models matches the data. The estimation function searches for parameter values which minimize this criterion. For a proposed set of parameter values it is derived by calculating
(i) the components z1,t,z2,t,,zm,tz1,t,z2,t,,zm,t as the responses to the input series x1,t,x2,t,xm,tx1,t,x2,t,xm,t using the equations (a) or (b) above,
(ii) the discrepancy between the output and the sum of these components, as the noise
nt = yt (z 1 , t + z 2 , t + + z m , t ) ,
n t = y t - ( z 1 , t + z 2 , t + + z m , t ) ,
(iii) the residual series atat from ntnt by reversing the recursive equations (i), (ii) and (iii) above.
This last step again requires treatment of the effect of unknown pre-period values of ntnt and other terms in the equations regenerating atat. This is identical to the treatment given in Section [Description] in (g13ae), and leads to a criterion which is a sum of squares function SS, of the residuals atat. It may be shown that the finite algorithm presented there is equivalent to taking the infinite set of past values n0,n1,n2,n0,n-1,n-2,, as (linear) nuisance parameters. The pre-period nuisance parameters for the input series are included in the reduction of dfdf, as is the constant if it is estimated.
The covariance matrix of the vector of model parameter estimates is given by
erv × H1
erv×H-1
where HH is the linearized least squares matrix taken from the final iteration of the algorithm of Marquardt. From this expression are derived the vector of standard deviations, and the correlation matrix of parameter estimates. These are approximations which are only valid asymptotically, and must be treated with great caution when the parameter estimates are close to their constraint boundaries.
The residual series atat is available upon completion of the iterations over the range t = 1 + d + s × D,,nt=1+d+s×D,,n corresponding to the differenced noise series wtwt.
Because of the algorithm used for backforecasting, these are only true residuals for t1 + q + s × Qps × Pds × Dt1+q+s×Q-p-s×P-d-s×D, provided this is positive. Estimation of pre-period terms for the inputs will also tend to reduce the magnitude of the early residuals, sometimes severely.
The model component series z1,t,,zm,tz1,t,,zm,t and ntnt may optionally be returned in place of the supplied series values, in order to assess the effects of the various inputs on the output.

Forecasting Information

For the purpose of constructing forecasts of the output series at future time points t = n + 1,n + 2,t=n+1,n+2, using nag_tsa_multi_inputmod_forecast_state (g13bh), it is not necessary to use the whole set of observations ytyt and x1,t,x2,t,,xm,tx1,t,x2,t,,xm,t, for t = 1,2,,mt=1,2,,m. It is sufficient to retain a limited set of quantities constituting the ‘state set’ as follows: for each series which appears with lagged subscripts in equations (a), (b), (i), (ii) and (iii) above, include the values at times n + 1kn+1-k for k = 1k=1 up to the maximum lag associated with that series in the equations. Note that (i) implicitly includes past values of ntnt and intermediate differences of ntnt such as d1sDd-1sD.
If later observations of the series become available, it is possible to update the state set (without re-estimating the model) using nag_tsa_multi_inputmod_update (g13bg). If time series data is supplied with a previously estimated model, it is possible to construct the state set (and forecasts) using nag_tsa_multi_inputmod_forecast (g13bj).

References

Box G E P and Jenkins G M (1976) Time Series Analysis: Forecasting and Control (Revised Edition) Holden–Day
Marquardt D W (1963) An algorithm for least squares estimation of nonlinear parameters J. Soc. Indust. Appl. Math. 11 431

Parameters

Compulsory Input Parameters

1:     mr(77) – int64int32nag_int array
The orders vector (p,d,q,P,D,Q,s)(p,d,q,P,D,Q,s) of the ARIMA model for the output noise component.
pp, qq, PP and QQ refer respectively to the number of autoregressive (φ)(ϕ), moving average (θ)(θ), seasonal autoregressive (Φ)(Φ) and seasonal moving average (Θ)(Θ) parameters.
dd, DD and ss refer respectively to the order of non-seasonal differencing, the order of seasonal differencing and the seasonal period.
Constraints:
  • pp, dd, qq, PP, DD, QQ, s0s0;
  • p + q + P + Q > 0p+q+P+Q>0;
  • s1s1;
  • if s = 0s=0, P + D + Q = 0P+D+Q=0;
  • if s > 1s>1, P + D + Q > 0P+D+Q>0;
  • d + s × (P + D)nd+s×(P+D)n;
  • p + dq + s × (P + DQ)np+d-q+s×(P+D-Q)n.
2:     mt(44,nser) – int64int32nag_int array
The transfer function model orders bb, pp and qq of each of the input series. The order parameters for input series ii are held in column ii. Row 11 holds the value bibi, row 2 holds the value qiqi and row 3 holds the value pipi. For a simple input, bi = qi = pi = 0bi=qi=pi=0.
Row 4 holds the value riri, where ri = 1ri=1 for a simple input, ri = 2ri=2 for a transfer function input for which no allowance is to be made for pre-observation period effects, and ri = 3ri=3 for a transfer function input for which pre-observation period effects will be treated by estimation of appropriate nuisance parameters.
When ri = 1ri=1, any nonzero contents of rows 1, 2, and 3 of column ii are ignored.
Constraint: mt(4,i) = 1, ​2​ or ​3mt4i=1, ​2​ or ​3, for i = 1,2,,nser1i=1,2,,nser-1.
3:     para(npara) – double array
npara, the dimension of the array, must satisfy the constraint npara = p + q + P + Q + nser + (pi + qi)npara=p+q+P+Q+nser+(pi+qi), the summation being over all the piqipiqi supplied in mt. cc must be included, whether fixed or estimated.
Initial values of the multi-input model parameters. These are in order, firstly the ARIMA model parameters: pp values of φϕ parameters, qq values of θθ parameters, PP values of ΦΦ parameters and QQ values of ΘΘ parameters. These are followed by initial values of the transfer function model parameters ω0,ω1,,ωq1ω0,ω1,,ωq1, δ1,δ2,,δp1δ1,δ2,,δp1 for the first of any input series and similarly for each subsequent input series. The final component of para is the initial value of the constant cc, whether it is fixed or is to be estimated.
4:     xxy(ldxxy,nser) – double array
ldxxy, the first dimension of the array, must satisfy the constraint ldxxynxxyldxxynxxy.
The columns of xxy must contain the nxxy original, undifferenced values of each of the input series and the output series xtxt in that order.

Optional Input Parameters

1:     nser – int64int32nag_int scalar
Default: The second dimension of the arrays mt, xxy. (An error is raised if these dimensions are not equal.)
The total number of input and output series. There may be any number of input series (including none), but always one output series.
Constraints:
  • nser1nser1;
  • if there are no parameters in the model (that is, p = q = P = Q = 0p=q=P=Q=0 and kfc = 0kfc=0), nser > 1nser>1.
2:     npara – int64int32nag_int scalar
Default: The dimension of the array para.
The exact number of φ,θ,Φ,Θϕ,θ,Φ,Θ, ω,δω,δ and cc parameters.
Constraint: npara = p + q + P + Q + nser + (pi + qi)npara=p+q+P+Q+nser+(pi+qi), the summation being over all the piqipiqi supplied in mt. cc must be included, whether fixed or estimated.
3:     kfc – int64int32nag_int scalar
Must be set to 00 if the constant cc is to remain fixed at its initial value, and 11 if it is to be estimated.
Default: 11
Constraint: kfc = 0kfc=0 or 11.
4:     nxxy – int64int32nag_int scalar
Default: The first dimension of the array xxy.
The (common) length of the original, undifferenced input and output time series.
5:     kef – int64int32nag_int scalar
Indicates the likelihood option.
kef = 1kef=1
Gives least squares.
kef = 2kef=2
Gives exact likelihood.
kef = 3kef=3
Gives marginal likelihood.
Default: 22
Constraint: kef = 1kef=1, 22 or 33.
6:     nit – int64int32nag_int scalar
The maximum required number of iterations.
nit = 0nit=0
No change is made to any of the model parameters in array para except that the constant cc (if kfc = 1kfc=1) and any ωω relating to simple input series are estimated. (Apart from these, estimates are always derived for the nuisance parameters relating to any backforecasts and any pre-observation period effects for transfer function inputs.)
Default: 10001000
Constraint: nit0nit0.
7:     zsp(44) – double array
If kzsp = 1kzsp=1, then zsp must contain the four values used to control the strategy of the search procedure.
zsp(1)zsp1
Contains αα, the value used to constrain the magnitude of the search procedure steps.
zsp(2)zsp2
Contains ββ, the multiplier which regulates the value of αα.
zsp(3)zsp3
Contains δδ, the value of the stationarity and invertibility test tolerance factor.
zsp(4)zsp4
Contains γγ, the value of the convergence criterion.
If kzsp1kzsp1 before entry, default values of zsp are supplied by the function. These are 0.010.01, 10.010.0, 1000.01000.0 and max (100 × machine precision,0.0000001)max(100×machine precision,0.0000001), respectively.
Constraint: if kzsp = 1kzsp=1, zsp(1) > 0.0zsp1>0.0, zsp(2) > 1.0zsp2>1.0, zsp(3)1.0zsp31.0, 0zsp(4) < 1.00zsp4<1.0.
8:     iwa – int64int32nag_int scalar
The dimension of the array wa as declared in the (sub)program from which nag_tsa_multi_inputmod_estim (g13be) is called.
It is not practical to outline a method for deriving the exact minimum permissible value of iwa, but the following gives a reasonably good conservative approximation. (It should be noted that if iwa is too small (but not grossly so) then the exact minimum is printed if kpriv0kpriv0.)
Let q = q + (Q × s)q=q+(Q×s) and d = d + (D × s)d=d+(D×s) where the orders of the output noise model are pp, dd, qq, PP, DD, QQ, ss.
Let there be ll input series, where l = nser1l=nser-1.
Let
mxi = max ( bi + qi ,pi) , if ​ ri = 3 , for ​ i = 1,2,l
mxi = 0 , if ​ ri3 , for ​ i = 1,2,l
mxi = max( bi + qi ,pi) , if ​ ri=3 , for ​ i=1,2,l mxi = 0 , if ​ ri3 , for ​ i=1,2,l
where the transfer function model orders for input ii are given by bibi, qiqi, pipi, riri.
Let qx = max (q,mx1,mx2,,mxl)qx=max(q,mx1,mx2,,mxl).
Let ncd = npara + kfc + qx + i = 1lmxincd=npara+kfc+qx+i=1lmxi and nce = nxxy + d + 6 × qxnce=nxxy+d+6×qx.
Finally, let ncf = nserncf=nser, and then increment ncfncf by 11 every time any of the following conditions is satisfied. (The last six conditions should be applied separately to each input series, so that, for example, if we have two input series and if p1 > 0p1>0 and p2 > 0p2>0 then ncfncf is incremented by 22.)
The conditions are:
and ​q > 0​ and ​kef > 1. and ​kfc > 0​ and ​kef = 3. and ​ri = 1​ and ​kef > 3​ separately, for ​i = 1,2, … ,l.
p>0 q>0 P>0 Q>0 qx>0 kfc>0   p>0 q>0 P>0 Q>0 } and ​q>0​ and ​kef>1. p>0 q>0 P>0 Q>0 } and ​kfc>0​ and ​kef=3. mxi>0 pi>0 p>0 q>0 P>0 Q>0 } and ​ri=1​ and ​kef>3​ separately, for ​i=1,2,,l.
p > 0 q > 0 P > 0 Q > 0 qx > 0 kfc > 0  
p > 0 q > 0 P > 0 Q > 0 }
p > 0 q > 0 P > 0 Q > 0 }
mxi > 0 pi > 0 p > 0 q > 0 P > 0 Q > 0 }
Then iwa2 × (ncd)2 + (nce) × (ncf + 4)iwa2× (ncd) 2+(nce)×(ncf+4).
Default: 5 × (2 × ncd2 + nce × (ncf + 4))5×(2×ncd2+nce×(ncf+4)) 
9:     imwa – int64int32nag_int scalar
The dimension of the array mwa as declared in the (sub)program from which nag_tsa_multi_inputmod_estim (g13be) is called.
Default: 5 × (16 × nser + 7 × ncd + 3 × npara + 3 × kfc + 27)5×( 16×nser+7×ncd+3×npara+3×kfc+27 ) 
Constraint: imwa(16 × nser) + (7 × ncd) + (3 × npara) + (3 × kfc) + 27imwa(16×nser)+(7×ncd)+(3×npara)+(3×kfc)+27, where the derivation of ncdncd is shown under iwa.
If imwa is too small then if kpriv0kpriv0 it is printed.
10:   kpriv – int64int32nag_int scalar
Must not be set to 00, if it is required to monitor the course of the optimization or to print out the requisite minimum values of iwa or imwa in the event of an error of the type ifail = 6ifail=6 or 77. The course of the optimization is monitored by printing out at each iteration the iteration count (itc), the residual sum of squares (s), the objective function (d) and a description and value for each of the parameters in the para array. The descriptions are PHI for φϕ, THETA for θθ, SPHI for ΦΦ, STHETA for ΘΘ, OMEGA/SI for ωω in a simple input, OMEGA for ωω in a transfer function input, DELTA for δδ and CONSTANT for cc. In addition SERIES 1, SERIES 2, etc. indicate the input series relevant to the OMEGA and DELTA parameters.
kpriv must be set to 00 if the print-out of the above information is not required.
Default: 00

Input Parameters Omitted from the MATLAB Interface

ldxxy kzsp ldcm wa mwa

Output Parameters

1:     para(npara) – double array
The latest values of the estimates of these parameters.
2:     xxy(ldxxy,nser) – double array
ldxxynxxyldxxynxxy.
the columns of xxy hold the corresponding values of the input component series ztzt in place of xtxt and the output noise component ntnt in place of ytyt, in that order.
3:     zsp(44) – double array
Contains the values, default or otherwise, used by the function.
4:     itc – int64int32nag_int scalar
The number of iterations carried out.
itc = -1itc=-1
Indicates that the only estimates obtained up to this point have been for the nuisance parameters relating to backforecasts, unless the marginal likelihood option is used, in which case estimates have also been obtained for simple input coefficients ωω and for the constant cc (if kfc = 1kfc=1). This value of itc usually indicates a failure in a consequent step of estimating transfer function input pre-observation period nuisance parameters.
itc = 0itc=0
Indicates that estimates have been obtained up to this point for the constant cc (if kfc = 1kfc=1), for simple input coefficients ωω and for the nuisance parameters relating to the backforecasts and to transfer function input pre-observation period effects.
5:     sd(npara) – double array
The npara values of the standard deviations corresponding to each of the parameters in para. When the constant is fixed its standard deviation is returned as zero. When the values of para are valid, the values of sd are usually also valid. However, if an exit value of ifail = 3ifail=3, 88 or 1010, then the contents of sd will be indeterminate.
6:     cm(ldcm,npara) – double array
ldcmnparaldcmnpara.
The first npara rows and columns of cm contain the correlation coefficients relating to each pair of parameters in para. All coefficients relating to the constant will be zero if the constant is fixed. The contents of cm will be indeterminate under the same conditions as sd.
7:     s – double scalar
The residual sum of squares, SS, at the latest set of valid parameter estimates.
8:     d – double scalar
The objective function, DD, at the latest set of valid parameter estimates.
9:     ndf – int64int32nag_int scalar
The number of degrees of freedom associated with SS.
10:   res(nxxy) – double array
The values of the residuals relating to the differenced values of the output series. The remainder of the first nxxy terms in the array will be zero.
11:   sttf(isttf) – double array
The nsttf values of the state set array.
12:   nsttf – int64int32nag_int scalar
The number of values in the state set array sttf.
13:   ifail – int64int32nag_int scalar
ifail = 0ifail=0 unless the function detects an error (see [Error Indicators and Warnings]).

Error Indicators and Warnings

Note: nag_tsa_multi_inputmod_estim (g13be) may return useful information for one or more of the following detected errors or warnings.
Errors or warnings detected by the function:
  ifail = 1ifail=1
On entry,kfc < 0kfc<0,
orkfc > 1kfc>1,
orldxxy < nxxyldxxy<nxxy,
orldcm < nparaldcm<npara,
orkef < 1kef<1,
orkef > 3kef>3,
ornit < 0nit<0,
ornser < 1nser<1,
ornser = 1nser=1 and there are no parameters in the model (p = q = P = Q = 0p=q=P=Q=0 and kfc = 0kfc=0).
  ifail = 2ifail=2
On entry, there is inconsistency between npara and kfc on the one hand and the orders in arrays mr and mt on the other, or one of the riri, stored in mt(4,i)1mt4i1, 22 or 33.
  ifail = 3ifail=3
On entry or during execution, one or more sets of δδ parameters do not satisfy the stationarity or invertibility test conditions.
  ifail = 4ifail=4
On entry,when kzsp = 1kzsp=1, zsp(1)0.0zsp10.0,
orzsp(2)1.0zsp21.0,
orzsp(3) < 1.0zsp3<1.0,
orzsp(4) < 0.0zsp4<0.0,
orzsp(4)1.0zsp41.0.
  ifail = 5ifail=5
On entry, iwa is too small by a considerable margin. No information is supplied about the requisite minimum size.
  ifail = 6ifail=6
On entry, iwa is too small, but the requisite minimum size is returned in mwa(1)mwa1, which is printed if kpriv0kpriv0.
  ifail = 7ifail=7
On entry, imwa is too small, but the requisite minimum size is returned in mwa(1)mwa1, which is printed if kpriv0kpriv0.
  ifail = 8ifail=8
This indicates a failure in nag_linsys_real_posdef_solve_1rhs (f04as) which is used to solve the equations giving the latest estimates of the parameters.
  ifail = 9ifail=9
This indicates a failure in the inversion of the second derivative matrix. This is needed in the calculation of the correlation matrix and the standard deviations of the parameter estimates.
  ifail = 10ifail=10
On entry or during execution, one or more sets of the ARIMA (φϕ, θθ, ΦΦ or ΘΘ) parameters do not satisfy the stationarity or invertibility test conditions.
  ifail = 11ifail=11
On entry, isttf is too small. The state set information will not be produced and if kzef0kzef0 array xxy will remain unchanged. All other parameters will be produced correctly.
  ifail = 12ifail=12
The function has failed to converge after nit iterations. If steady decreases in the objective function, DD, were monitored up to the point where this exit occurred, then the exit probably occurred because nit was set too small, so the calculations should be restarted from the final point held in para.
  ifail = 13ifail=13
On entry, isttf is too small (see ifail = 11ifail=11) and nit iterations were carried out without the convergence conditions being satisfied (see ifail = 12ifail=12).

Accuracy

The computation used is believed to be stable.

Further Comments

The time taken by nag_tsa_multi_inputmod_estim (g13be) is approximately proportional to nxxy × itc × npara2nxxy×itc×npara2.

Example

function nag_tsa_multi_inputmod_estim_example
mr = [int64(1);0;0;0;0;1;4];
mt = [int64(1),0; ...
             0,0; ...
             1,0; ...
             3,0];
para = [0;
     0;
     2;
     0.5;
     0];
xxy = [8.075, 105;
     7.819, 119;
     7.366, 119;
     8.113, 109;
     7.38, 117;
     7.134, 135;
     7.222, 126;
     7.768, 112;
     7.386, 116;
     6.965, 122;
     6.478, 115;
     8.105, 115;
     8.06, 122;
     7.684, 138;
     7.58, 135;
     7.093, 125;
     6.129, 115;
     6.026, 108;
     6.679, 100;
     7.414, 96;
     7.112, 107;
     7.762, 115;
     7.645, 123;
     8.639, 122;
     7.667, 128;
     8.08, 136;
     6.678, 140;
     6.739, 122;
     5.569, 102;
     5.049, 103;
     5.642, 89;
     6.808, 77;
     6.636, 89;
     8.241, 94;
     7.968, 104;
     8.044, 108;
     7.791, 119;
     7.024, 126;
     6.102, 119;
     6.053, 103];
kef = int64(3);
[paraOut, xxyOut, zspOut, itc, sd, cm, s, d, ndf, res, sttf, nsttf, ifail] = ...
     nag_tsa_multi_inputmod_estim(mr, mt, para, xxy, 'kef', kef)
 

paraOut =

    0.3809
   -0.2578
    8.9561
    0.6596
  -75.4355


xxyOut =

  180.5668  -75.5668
  191.4297  -72.4297
  196.3025  -77.3025
  195.4597  -86.4597
  201.5940  -84.5940
  199.0756  -64.0756
  195.2111  -69.2111
  193.4501  -81.4501
  197.1785  -81.1785
  196.2167  -74.2167
  191.8117  -76.8117
  184.5444  -69.5444
  194.3221  -72.3221
  200.3689  -62.3689
  200.9901  -65.9901
  200.4685  -75.4685
  195.7627  -80.7627
  184.0250  -76.0250
  175.3598  -75.3598
  175.4922  -79.4922
  182.1623  -75.1623
  183.8574  -68.8574
  190.7970  -67.7970
  194.3268  -72.3268
  205.5576  -77.5576
  204.2605  -68.2605
  207.1038  -67.1038
  196.4229  -74.4229
  189.9237  -87.9237
  175.1579  -72.1579
  160.7606  -71.7606
  156.5745  -79.5745
  164.2560  -75.2560
  167.7826  -73.7826
  184.4834  -80.4834
  193.0549  -85.0549
  199.3897  -80.3897
  201.3025  -75.3025
  195.6949  -76.6949
  183.7384  -80.7384


zspOut =

   1.0e+03 *

    0.0000
    0.0100
    1.0000
    0.0000


itc =

                   11


sd =

    0.1664
    0.1782
    0.9481
    0.0602
   33.5053


cm =

    1.0000   -0.1839   -0.1775   -0.0340    0.1394
   -0.1839    1.0000    0.0518    0.2547   -0.2860
   -0.1775    0.0518    1.0000   -0.3070   -0.2926
   -0.0340    0.2547   -0.3070    1.0000   -0.8185
    0.1394   -0.2860   -0.2926   -0.8185    1.0000


s =

   1.1980e+03


d =

   1.2866e+03


ndf =

                   34


res =

    0.3973
    3.0856
   -2.8185
   -9.9406
   -5.0615
   14.0532
    2.6237
   -5.8231
   -2.1471
   -0.2162
   -2.5168
    7.9165
    1.4228
   11.9364
    5.1168
   -5.6717
   -5.6814
   -1.6372
   -1.0188
   -2.6235
    3.2831
    6.8961
    5.3954
    0.8753
   -4.1526
    6.2056
    4.2078
   -2.3868
  -11.8034
    6.4349
    1.3417
   -4.9236
    4.7989
   -0.0743
   -6.0234
   -6.4273
   -2.5270
    2.0394
    0.2427
   -3.1663


sttf =

    6.0530
  183.7384
   -5.7855
   -0.1645
    0.1800
   -3.0977


nsttf =

                    6


ifail =

                    0


function g13be_example
mr = [int64(1);0;0;0;0;1;4];
mt = [int64(1),0; ...
             0,0; ...
             1,0; ...
             3,0];
para = [0;
     0;
     2;
     0.5;
     0];
xxy = [8.075, 105;
     7.819, 119;
     7.366, 119;
     8.113, 109;
     7.38, 117;
     7.134, 135;
     7.222, 126;
     7.768, 112;
     7.386, 116;
     6.965, 122;
     6.478, 115;
     8.105, 115;
     8.06, 122;
     7.684, 138;
     7.58, 135;
     7.093, 125;
     6.129, 115;
     6.026, 108;
     6.679, 100;
     7.414, 96;
     7.112, 107;
     7.762, 115;
     7.645, 123;
     8.639, 122;
     7.667, 128;
     8.08, 136;
     6.678, 140;
     6.739, 122;
     5.569, 102;
     5.049, 103;
     5.642, 89;
     6.808, 77;
     6.636, 89;
     8.241, 94;
     7.968, 104;
     8.044, 108;
     7.791, 119;
     7.024, 126;
     6.102, 119;
     6.053, 103];
kef = int64(3);
[paraOut, xxyOut, zspOut, itc, sd, cm, s, d, ndf, res, sttf, nsttf, ifail] = ...
     g13be(mr, mt, para, xxy, 'kef', kef)
 

paraOut =

    0.3809
   -0.2578
    8.9561
    0.6596
  -75.4355


xxyOut =

  180.5668  -75.5668
  191.4297  -72.4297
  196.3025  -77.3025
  195.4597  -86.4597
  201.5940  -84.5940
  199.0756  -64.0756
  195.2111  -69.2111
  193.4501  -81.4501
  197.1785  -81.1785
  196.2167  -74.2167
  191.8117  -76.8117
  184.5444  -69.5444
  194.3221  -72.3221
  200.3689  -62.3689
  200.9901  -65.9901
  200.4685  -75.4685
  195.7627  -80.7627
  184.0250  -76.0250
  175.3598  -75.3598
  175.4922  -79.4922
  182.1623  -75.1623
  183.8574  -68.8574
  190.7970  -67.7970
  194.3268  -72.3268
  205.5576  -77.5576
  204.2605  -68.2605
  207.1038  -67.1038
  196.4229  -74.4229
  189.9237  -87.9237
  175.1579  -72.1579
  160.7606  -71.7606
  156.5745  -79.5745
  164.2560  -75.2560
  167.7826  -73.7826
  184.4834  -80.4834
  193.0549  -85.0549
  199.3897  -80.3897
  201.3025  -75.3025
  195.6949  -76.6949
  183.7384  -80.7384


zspOut =

   1.0e+03 *

    0.0000
    0.0100
    1.0000
    0.0000


itc =

                   11


sd =

    0.1664
    0.1782
    0.9481
    0.0602
   33.5053


cm =

    1.0000   -0.1839   -0.1775   -0.0340    0.1394
   -0.1839    1.0000    0.0518    0.2547   -0.2860
   -0.1775    0.0518    1.0000   -0.3070   -0.2926
   -0.0340    0.2547   -0.3070    1.0000   -0.8185
    0.1394   -0.2860   -0.2926   -0.8185    1.0000


s =

   1.1980e+03


d =

   1.2866e+03


ndf =

                   34


res =

    0.3973
    3.0856
   -2.8185
   -9.9406
   -5.0615
   14.0532
    2.6237
   -5.8231
   -2.1471
   -0.2162
   -2.5168
    7.9165
    1.4228
   11.9364
    5.1168
   -5.6717
   -5.6814
   -1.6372
   -1.0188
   -2.6235
    3.2831
    6.8961
    5.3954
    0.8753
   -4.1526
    6.2056
    4.2078
   -2.3868
  -11.8034
    6.4349
    1.3417
   -4.9236
    4.7989
   -0.0743
   -6.0234
   -6.4273
   -2.5270
    2.0394
    0.2427
   -3.1663


sttf =

    6.0530
  183.7384
   -5.7855
   -0.1645
    0.1800
   -3.0977


nsttf =

                    6


ifail =

                    0



PDF version (NAG web site, 64-bit version, 64-bit version)
Chapter Contents
Chapter Introduction
NAG Toolbox

© The Numerical Algorithms Group Ltd, Oxford, UK. 2009–2013