automatic — automatic general to specific regression
[results]=automatic(namey,varargin)
namey = a time series, a real (nx1) vector or a string equal to the name of a time series or a (nx1) real vector between quotes
varargin = arguments which can be:
* a time series
* a real (nxp) vector
* a string equal to the name of a time series or a (nx1) real vector between quotes
* the string 'prt=opt1,opt2, ,optn' where opti is one of the following available options:
- nothing (nothing printed !)
- initial (results of initial model printed)
- st0_mod (results of stage 0 model printed)
- st1_mod (results of stage 1 models printed)
- st1_union (results of stage 1 union model printed)
- st2_mod (results of stage 2 models printed)
- st2_union (results of stage 2 union model printed)
- final (results of final model printed)
- test_inter (results of specification tests for the intermediate results printed)
- test_final (results of specification tests for the final model printed)
- test (specification tests printed for every regression result printed)
- st1_path (results of stage 1 paths printed)
- st2_path (results of stage 2 paths printed)
- all (all results printed: can be very long !)
* the string 'test=opt1,opt2, ,optn' where opti is one of the following available options:
- chowtest(p)
- predfailin(p)
- doornhans
- arlm(p)
- arch(p)
- hetero_sq
* the string 'strategy=liberal' or 'strategy=conservative' if the user wants to use a predifined strategy (default = liberal)
* the string 'alpha=p' for the simplification significance level (default = set by the chosen strategy)
* the string 'f0_sig=p' for the top_down pre-test significance level (default = set by the chosen strategy)
* the string 'f0_bup=p' for the bottom-up pre-test significance level (default = set by the chosen strategy)
* the string 'eta=p' (default = 0.01)
* the string 'gam=p' for the F-tests significance level (default = set by the chosen strategy)
* the string 'groups_pval=m1;m2; ;mp' where the mi's are thresholds: coefficient whose p-values are greater than a given threshold are gathered to form a model included in stage 1 models (default = set by the chosen strategy)
* the string 'crit=x' where x= aic, bic, or hq (default bic)
* the string 'test=x1, ,xp' where xi is the name of a test function
* the string 'newname(x1,x2)' where x1 is the name of a test function and x2 is the corresponding name the user wants the program to display
* the string 'comp=x1; ;xn' where xi is the name of a variable that must be in the regression whatever significance it has
* the string 'estim=nwest' if the user wants to use the Newey-West correction on standard error Bartlett window (for the Newey-West estimation)
* the string 'win=xx' where xx is the length of the Bartlett window (for the Newey-West estimation)
a results tlist with the following fields:
- results('meth') = 'automatic';
- results('estim') = estimation method (ols or nwest)
- results('strategy') = estimation strategy (liberal, conservative or expert)
- results('alpha') = simplification significance level (default = set by the chosen strategy)
- results('f0_tdo') = top_down pre-test significance level (default = set by the chosen strategy)
- results('f0_bup') = bottom-up pre-test significance level (default = set by the chosen strategy)
- results('gam') = F-tests significance level (default = set by the chosen strategy)
- results('eta') = specification tests significance level (default = 0.01)
- results('comp') = matrix of values for the compulsory variables
- results('f_test') = the function used to perform the specification tests
- results('m_test') = the names of the specification tests
- results('initial model') = the estimation results of the unrestricted model
- results('ending reason') = the reason why the final model has been chosen
- results('final model') = the estimation results of the final model
- results('top-down model') = the top-down model
- results('bottom-up model') = the bottom-up model
- results('stage 0 model') = the estimation results of the stage 0 restricted model (all variables whose individual tstat level is lower than 1 and whose joint significance level is lower than f0_sig are withdrawn)
- results('selected stage 0 model') = the way the stage 0 model has been selected
- results('checking process model') = the model resulting from the group checking process
- results('stage i models') = the estimation results of the stage i (i=1,2) models and the corresponding paths
- results('stage i union model') = the estimation results of the model built from the union of stage i (i=1,2) models
load(GROCERDIR+'/data/bdhenderic.dat') ; bounds('1964q3','1989q2') [r1]=automatic('delts(lm1-lp)','lagts(lm1-lp-ly)','delts(lp)','rnet','delts(lagts(lm1-lp))', 'delts(lagts(2,lm1-lp))','delts(lagts(3,lm1-lp))','delts(lagts(4,lm1-lp))',&133; 'delts(ly)','delts(lagts(1,ly))', 'delts(lagts(2,ly))','delts(lagts(3,ly))','delts(lagts(4,ly))',&133; 'delts(delts(lp))','delts(delts(lagts(1,lp)))','delts(delts(lagts(2,lp)))','delts(delts(lagts(3,lp)))','delts(delts(lagts(4,lp)))',&133; 'delts(rnet)','delts(lagts(1,rnet))','delts(lagts(2,rnet))','delts(lagts(3,rnet))','delts(lagts(4,rnet))','const',&133; 'prt=initial,final','test=predfailin(0.5),predfailin(0.9),doornhans,arlm(5),hetero_sq','groups_pval=[]') // Starts from a rather extensive set of exogenous variables: from 'lagts(lm1-lp-ly)' to 'cte' in examples 1 and 2 and from 'lagts(lm1)' in example 3. // Here, the initial model is printed along with the final one; tests are given by the user and models of coefficients with p-value greater than some threshold are not examined. [r2]=automatic('delts(lm1-lp)','lagts(lm1-lp-ly)','delts(lp)','rnet','delts(lagts(lm1-lp))', 'delts(lagts(2,lm1-lp))','delts(lagts(3,lm1-lp))','delts(lagts(4,lm1-lp))',&133; 'delts(ly)','delts(lagts(1,ly))', 'delts(lagts(2,ly))','delts(lagts(3,ly))','delts(lagts(4,ly))',&133; 'delts(delts(lp))','delts(delts(lagts(1,lp)))','delts(delts(lagts(2,lp)))','delts(delts(lagts(3,lp)))','delts(delts(lagts(4,lp)))',&133; 'delts(rnet)','delts(lagts(1,rnet))','delts(lagts(2,rnet))','delts(lagts(3,rnet))','delts(lagts(4,rnet))','const',&133; 'prt=final','test=predfailin(0.5),predfailin(0.9),doornhans,arlm(5),hetero_sq','strategy=conservative') // Same set of variables, but only the final model is printed, tests are given by the user and the strategy is turned from default (liberal) to conservative. [r3]=automatic('delts(lm1-lp)','lagts(lm1-lp-ly)','delts(lp)','rnet','delts(lagts(lm1-lp))', 'delts(lagts(2,lm1-lp))','delts(lagts(3,lm1-lp))','delts(lagts(4,lm1-lp))',&133; 'delts(ly)','delts(lagts(1,ly))', 'delts(lagts(2,ly))','delts(lagts(3,ly))','delts(lagts(4,ly))',&133; 'delts(delts(lp))','delts(delts(lagts(1,lp)))','delts(delts(lagts(2,lp)))','delts(delts(lagts(3,lp)))','delts(delts(lagts(4,lp)))',&133; 'delts(rnet)','delts(lagts(1,rnet))','delts(lagts(2,rnet))','delts(lagts(3,rnet))','delts(lagts(4,rnet))','const',&133; 'prt=initial,final,st2_mod','test=predfailin(0.51),predfailin(0.862),doornhans,arlm(5),hetero_sq','alpha=0.01') // Same set of variables, but the initial model and all the stage 2 model (among which the final model is selected by the BIC information criterion) are printed along with the final one, tests are again given by the user and the significance level of individual coefficients is set to 0.01.