matrix(5,28) pv
matrix(5,28) coef 
matrix(5,28) tstat
matrix(28,30) RESULTS
vector(28) r2s
vector(28) aic
vector(28) rbar

!maxlags =2
!bestlag_a = 0
!bestlag_b = 0
!bestlag_C = 0

 for !i=1 to 28
!aic = 99999999
!pval_C = 0.1
!pval_D =0.1
for !lag_a = 0 to !maxlags
   for !lag_b = 0 to !maxlags
		 
    for !lag_C = 0 to !maxlags
        
equation eq{!i}.ls Y{!i} c X{!i}A(0 to -!lag_a) X{!i}B(0 to -!lag_b) Y{!i}(-1) X{!i}C(0 to -!lag_C)

!fake_aic = eq{!i}.@aic - 1000 * ((eq{!i}.@pval(2) < !pval_C) + (eq{!i}.@pval(3) < !pval_D))
if !fake_aic < !aic then
   !aic = !fake_aic
   !bestlag_a = !lag_a
   !bestlag_b = !lag_b
  !bestlag_C = !lag_C
endif
show eq{!i}.ls   Y{!i} c X{!i}A(-!bestlag_a) X{!i}B(-!bestlag_b ) Y{!i}(-1) X{!i}C(-!bestlag_C )
colplace(pv, eq!i.@pval, !i)
colplace(coef, eq!i.@coef, !i)
colplace(tstat, eq!i.@tstats, !i)


 r2s(!i) = eq!i.@r2 
 aic(!i) = eq!i.@aic
 rbar(!i) = eq!i.@rbar2


matrix PVALUE = @transpose(pv)   
matrix coefs = @transpose(coef)
matrix tstats = @transpose(tstat)
results.setwidth 25
matplace(results,PVALUE,1,1)
matplace(results,coefs,1,5)
matplace(RESULTS,tstats,1,9)

colplace(RESULTS,R2S,13)
colplace(RESULTS,AIC,14)
colplace(RESULTS,RBAR,15)

SHOW RESULTS
NEXT
NEXT
next
   next