matrix(4,28) pv
matrix(4,28) coef 
matrix(4,28) tstat
matrix(28,15) RESULTS
vector(28) r2s
vector(28) aic
vector(28) rbar

!maxlags =2
!bestlag_a = 0

!bestlag_b = 0


 for !i=1 to 3
 !aic_1 = 99999999
!aic_2 = 99999999
!aic_3 = 99999999
!pval_a = 99999999
!pval_b = 99999999
!pval_C = 0.1
!pval_D = 0.1
for !lag_a = 0 to !maxlags
   for !lag_b = 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)

if (eq{!i}.@pval(2) <=!pval_C AND eq{!i}.@pval(3) <=!pval_D AND  eq{!i}.@aic <!aic_1  ) or ( eq{!i}.@pval(2) <=!pval_C AND eq{!i}.@pval(3) <=!pval_D  ) then   ', CHECKS IF THE PVALUES ARE LESS THAN 0.1 AND IF MANY PAIRS EXISTS SELECTS WITH BEST AI 
 		!bestlag_a =(!lag_a)  ', then store this lag as !bestlag.
		!bestlag_b =(!lag_b) 


			else if [eq{!i}.@pval(2) <!pval_C or eq{!i}.@pval(3)<!pval_D ] or [{eq{!i}.@pval(2) <!pval_C OR eq{!i}.@pval(3)<!pval_D} and { eq{!i}.@aic <!aic_2 }] then   ', SELECT EQUATION IF  XA OR XB HAS  PV LESS THAN 0.1 AND SELECT MINIMUM AIC EQUATION
				 !bestlag_a =(!lag_a)  ', then store this lag as !bestlag.
				!bestlag_b =(!lag_b)
							else if  eq{!i}.@aic < !aic_3 then 
							!bestlag_a =(!lag_a)  ', then store this lag as !bestlag.
							!bestlag_b =(!lag_b)

 ',

 













ENDIF
endif
endif

show eq{!i}.ls   Y{!i} c X{!i}A(-!bestlag_a) X{!i}B(-!bestlag_b ) Y{!i}(-1)
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