Re: [Gretl-users] Function for Automatic Modelling of ARIMA Models

On Mon, 23 Apr 2012, Henrique Andrade wrote: True, you can't save models by name inside a function. This is part of the "encapsulation" idea for gretl functions: their ability to produce "side effects" is strictly circumscribed. Basically, functions are allowed to print stuff and/or return a value, and nothing else. But if you want to return multiple objects (in a user-friendly form) you can use the bundle data-type for that purpose. (See also http://www.wfu.edu/~cottrell/tmp/extending.pdf ) I'm not seeing a problem with that here. This works OK: <hansl> function void testarima (series y) loop p=1..2 loop q=1..2 arima p 1 q ; y --nc fcast Y_hat_$p_1_$q print Y_hat_$p_1_$q endloop endloop end function open fedstl.bin data exbzus dataset addobs 12 testarima(exbzus) </hansl> Of course the series generated by "fcast" are local to the function -- they won't be available outside the function unless you return then in some way. Here are some comments on your function: 1. Using "set halt_on_error off" is a bad idea (perhaps we should get rid of it): if you do that then any errors will cascade and error messages will likely be hard to understand. Do use, instead, the "catch" modifier, and check the $error variable afterwards. This is important for functions that invoke commands such as arima, where failure of convergence is a live possibility. 2. There's no need to "set echo off" or "set messages off" inside a function; that's the default behavior. But you can set those things on for debugging. 3. It's good gretl programming style to use the string representation of loop indices only where a string is actually needed; otherwise just use the numerical value of the index. So: fcast Y_hat_$P_1_$Q # OK, strings needed arima P 1 Q # strings not needed, don't use "arima $P 1 $Q" 4. Consider limiting the sum of the AR and MA orders to something sane. Your function as written allows up to arima(12, 1, 12), which would surely be grossly over-parameterized. 5. There seem to be some problems with initialization of your arima models that include dummies for observations where the plain arima residuals are greater than two standard deviations. We use nonlinear least squares for initialization of arima with exogenous regressors: that works quite well in some cases but apparently does not work well here. You might consider using "set initvals" to specify your own initialization. For example, in the arima(1,1,1) case: <hansl> diff y ols d_y d_y(-1) x1 x2 matrix m = {$coeff[1], 0.001, $coeff[2], $coeff[3]} set initvals m arima 1 1 1 ; y x1 x2 --nc </hansl> Allin Cottrell