smpl @all

	'	Run this program after the data has been patched using gaps.prg

'	This program imports the load and temp data, and creates the summer models

'	*******************************************
'	* user enters the utility region below *
' 	*******************************************

%utility = "ae"

'	****************************************
'	* program fills in other information *
' 	****************************************

if %utility = "ae" then
	%weather = "acy"
	%city = "AtlanticCity,NJ" 
endif

if %utility = "bc" then 
	%weather = "bwi"
	%city = "Baltimore,MD"
endif

if %utility = "dpl" then 
	%weather = "ilg"
	%city = "Wilmington,DE"
endif

if %utility = "jc" then 
	%weather = "blm"
	%city = "Belmar,NJ"
endif

if %utility = "me" then 
	%weather = "thv"
	%city = "Reading,PA_RDG_York,PA"
endif

if %utility = "pe" then 
	%weather = "phl"
	%city = "Philadelphia,PA"
endif

if %utility = "pep" then 
	%weather = "dca"
	%city = "Washington,DC"
endif

if %utility = "pl" then 
	%weather = "abe"
	%city = "Allentown,PA"
endif

if %utility = "pn" then 
	%weather = "eri_ipt_jst"
	'this utility is a special case, see below
endif

if %utility = "ps" then 
	%weather = "ewr"
	'this utility is a special case, see below
endif

%data_path = "C:\EPA\Data\Working_Data\PJM_Utilities\"

	%data_file = %data_path + %utility +"_"+ %city +"_"+ %weather +"\"+ %utility +"_"+ %weather +".xls"
	%weather_sheet = %weather + "_Temp_WB"
	%output = "c:\epa\hourly_models\PJM_Utilities\"+ %utility +"_"+ %city +"_"+ %weather +"\" +%utility + "_summer_models"

'     *******************************
'     *  Special Case Utilities   *
'     *******************************

if %utility = "pn" then 

%data_file = %data_path + "PN_Erie,PA_ERI_Williamsport,PA_IPT_Johnstown,PA_JST\"+ %utility +"_"+ %weather +".xls"
	%weather_sheet = %weather + "_ave_Temp_WB"
	%output = "c:\epa\hourly_models\PJM_Utilities\"+ "PN_Erie,PA_ERI_Williamsport,PA_IPT_Johnstown,PA_JST\" +%utility + "_summer_models"

endif

if %utility = "ps" then 

%data_file = %data_path + "PS_Newark,NJ_EWR_Trenton,NJ_TTN\"+ %utility +"_"+ %weather +".xls"
	%weather_sheet = %weather + "_Temp_WB"
	%output = "c:\epa\hourly_models\PJM_Utilities\"+ "PS_Newark,NJ_EWR_Trenton,NJ_TTN\" +%utility + "_summer_models"

endif

'	c:\epa\hourly_models\PJM_Utilities\AE_AtlanticCity,NJ_ACY\ae_°Ftemp_hypoth1_cub

'     *******************************
'     *  Create the pool object  *
'     *******************************

pool pool_%utility _01 _02 _03 _04 _05 _06 _07 _08 _09 _10 _11 _12 _13 _14 _15 _16 _17 _18 _19 _20 _21 _22 _23 _24


'     *****************************
'     *  Import Calendar data  *
'     *****************************

read(d2) C:\EPA\Data\Calendar.xls 23

rename january jan
rename february feb
rename march mar
rename april apr
rename June jun
rename July jul
rename August aug
rename September sep
rename October oct
rename November nov
rename December dec

read(bh2) C:\EPA\Data\Calendar.xls Holiday

series dsum = jun + jul + aug + sep
series dwin = dec + jan + feb + mar 


'     *************************
'     *  Import Load data  *
'     *************************

read(c732, s=%utility) %data_file load_01 load_02 load_03 load_04 load_05 load_06 load_07 load_08 load_09 load_10 load_11 load_12 load_13 load_14 load_15 load_16 load_17 load_18 load_19 load_20 load_21 load_22 load_23 load_24

'     *************************
'     *  Import Temp data  *
'     *************************

read(d2, s=%weather_sheet) %data_file temp_01 temp_02 temp_03 temp_04 temp_05 temp_06 temp_07 temp_08 temp_09 temp_10 temp_11 temp_12 temp_13 temp_14 temp_15 temp_16 temp_17 temp_18 temp_19 temp_20 temp_21 temp_22 temp_23 temp_24

'*****************************
'* End of Setup Program *
'*****************************

'=======================================*** SUBROUTINE ***=================================================
subroutine matcount
'
'	Subroutine to initialize matrix for storing hourly lag indicators
'
for !a=1 to 27
'
	if !a=1 then
	mat_count(!a,1)=2
	mat_count(!a,2)=2
	endif
'
	if !a=2 then
	mat_count(!a,1)=2
	mat_count(!a,2)=3
	endif
'
	if !a=3 then
	mat_count(!a,1)=2
	mat_count(!a,2)=4
	endif
'
	if !a>3 and !a<13 then
	mat_count(!a,1)=0
	mat_count(!a,2)=!a-3
	endif
'
	if !a>13 then
	mat_count(!a,1)=1
	mat_count(!a,2)=!a-13
	endif
'
	if !a=13 then
	mat_count(!a,1)=1
	mat_count(!a,2)=0
	endif
'
	if !a>22 then
	mat_count(!a,1)=2
	mat_count(!a,2)=!a-23
	endif
'
next
endsub
'==============================================
subroutine makescalar
'
'	Subroutine to initialize scalars for hypothesis tests
'
!i =01
while !i < 55
	scalar lr01r!i = 0
	scalar lr01r!i_p=0
	scalar lr01r!i_r=0
	!i=!i+1
wend
'
endsub
'
'==============================================
'
subroutine lr_test
'
'	Subroutine LR_TEST
'
'	Likelhood Ratio Test for Restrictions
'
'	%i is the restricted equation identifier
'
'	scalar lr01r!i      = 
'	scalar lr01r!i_p = 
'	scalar lr01r!i_r  = 
	mat_logl(!i,!j) = eq_h!b!cr!i.@logl
	mat_hypoth(!i,!j*2-1) =  -2*(eq_h!b!cr!i.@logl - eqm_!b!c.@logl)

'	THE FOLLOWING IF STATEMENT GUARDS AGAINST NEGATIVE ARGUMENTS BEING PASSED TO THE @CCHISQ

	if mat_hypoth(!i,!j*2-1) < 0 then
		mat_hypoth(!i,!j*2) = 1 - @cchisq(0,@rows(eqm_!b!c.@coefs)-@rows(eq_h!b!cr!i.@coefs))
		else mat_hypoth(!i,!j*2) = 1 - @cchisq(mat_hypoth(!i,!j*2-1),@rows(eqm_!b!c.@coefs)-@rows(eq_h!b!cr!i.@coefs))
	endif

	mat_rest(!i,!j) = @rows(eqm_!b!c.@coefs)-@rows(eq_h!b!cr!i.@coefs)
	mat_aic(!i,!j) = @aic
	mat_schwarz(!i,!j) = @schwarz
'
endsub
'
'==============================================
'

'==================================================================================================
'
'	THIS PROGRAM IS DESIGNED TO CALCULATE COOLING DEGREE ELASTICITY ESTIMATES FOR HOURLY MODELS 
'
'
'==================================================================================================

smpl 1/1/1995 9/30/2003
output(t) %output
'	
'	Construct indicator matrix for modeling and variables
'	Create Storage Matrix for all Statistics
'
matrix(27,2) mat_count
'
'	Call subroutine to initialize scalars for hypothesis tests
'
'	call makescalar
'
'	Initialize matrix to store hypothesis test results
'
matrix(55,24) mat_logl
matrix(55,48) mat_hypoth
matrix(55,24) mat_rest
matrix(55,24) mat_aic
matrix(55,24) mat_schwarz
'

'	Why are the following two lines included?
'
'	EXTRACT MEAN AND MAXIMUM °F VALUES
'

'smpl 1/1/1998 3/31/1998

'
'
'	Construct indicator matrix for modeling and variables
'
call matcount

'=======================================*** Hours 01-24 Models ***====================================
'
'	Loop for Hours 1 to 24
'
for !z=4 to 27

	!b=mat_count(!z,1)
	!c=mat_count(!z,2)
	!m=mat_count(!z-1,1)
	!n=mat_count(!z-1,2)
	!o=mat_count(!z-2,1)
	!p=mat_count(!z-2,2)
	!q=mat_count(!z-3,1)
	!r=mat_count(!z-3,2)
'
'	Hourly Model Unrestricted Model
'
smpl 6/1/1995 9/30/1995  6/1/1996 9/30/1996  6/1/1997 9/30/1997  6/1/1998 9/30/1998  6/1/1999 9/30/1999   6/1/2000 9/30/2000  6/1/2001 9/30/2001  6/1/2002 9/30/2002  6/1/2003 9/30/2003  

'
'	Initialize Row (Eqn/Hypothesis) Hour Counter
'
!i=1
!j=!z-3
'
'	Adjust Lagged Hours for "beginning of day" effects (Hour1, Hour2, Hour3)
'
!day3=0
!day2=0
!day1=0

if !b=0 and !c=1 then
	!day3=-1
	!day2=-1
	!day1=-1
endif

if !b=0 and !c=2 then
	!day3=-1
	!day2=-1
endif

if !b=0 and !c=3 then
	!day3=-1
endif


'
'	Hourly Model Unrestricted Model
'
table(1,1) model_desc
setcolwidth(model_desc,1,80)
setcell(model_desc,1,1, "Model Number: 0")
setcell(model_desc,2,1,"Hourly Unrestricted Model")
print model_desc

equation eqm_!b!c.ls(n,p) log(load_!b!c) c tuesday wednesday thursday friday saturday sunday holiday jul aug sep  jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c  jun*temp_!b!c^2 jul*temp_!b!c^2 aug*temp_!b!c^2 sep*temp_!b!c^2 jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))


	mat_logl(!i,!j) = eqm_!b!c.@logl
	mat_aic(!i,!j) = @aic
	mat_schwarz(!i,!j) = @schwarz

'
'
'	All Day Effects are Zero
'	Tue=Wed=Thu=Fri=Sat=Sun=Holiday=0
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Day or Holiday effects")
setcell(model_desc,3,1,"Tue=Wed=Thu=Fri=Sat=Sun=Holiday=0")
print model_desc

!i=2
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c    jul aug sep  jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c  jun*temp_!b!c^2 jul*temp_!b!c^2 aug*temp_!b!c^2 sep*temp_!b!c^2 jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test

'
'
'	All Weekend and Holiday Effects are Zero
'	Sat=Sun=Holiday=0
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Weekend or Holiday effects")
setcell(model_desc,3,1,"Sat = Sun = Holiday = 0")
print model_desc

!i=3
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c  tuesday wednesday thursday friday  jul aug sep  jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c  jun*temp_!b!c^2 jul*temp_!b!c^2 aug*temp_!b!c^2 sep*temp_!b!c^2 jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'
'	All Weekend Effects are Zero
'	Sat=Sun=0
'
table(1,1) model_desc
setcolwidth(model_desc,1,80)
 setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Weekend effects")
setcell(model_desc,3,1,"Sat = Sun = 0")
print model_desc

!i=4
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c    tuesday wednesday thursday friday holiday jul aug sep  jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c  jun*temp_!b!c^2 jul*temp_!b!c^2 aug*temp_!b!c^2 sep*temp_!b!c^2 jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'
'	All Sunday and Holiday Effects are Zero
'	Sun=Holiday=0
'
table(1,1) model_desc
setcolwidth(model_desc,1,80)
 setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Sunday or Holiday Effects")
setcell(model_desc,3,1,"Sun = Holiday = 0")
print model_desc

!i=5
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c  tuesday wednesday thursday friday saturday  jul aug sep  jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c  jun*temp_!b!c^2 jul*temp_!b!c^2 aug*temp_!b!c^2 sep*temp_!b!c^2 jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'
'	All Weekday Effects are the same There are no weekend effects
'	Mon=Tue=Wed=Thu=Fri,   Sat=Sun=0, Holiday unrestricted
'
table(1,1) model_desc
setcolwidth(model_desc,1,80)
 setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Weekend effects; Weekday effects are all the same")
setcell(model_desc,3,1,"Sat = Sun = 0; Tue = Wed = Thu = Fri")
print model_desc

!i=6
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c (tuesday+wednesday+thursday+friday) holiday   jul aug sep  jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c  jun*temp_!b!c^2 jul*temp_!b!c^2 aug*temp_!b!c^2 sep*temp_!b!c^2 jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'
'	There are no monthly dummies
'	Jul=Aug=Sep=0
'
table(1,1) model_desc
setcolwidth(model_desc,1,80)
 setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Month effects")
setcell(model_desc,3,1,"Jul = Aug = Sep = 0")
print model_desc

!i=7
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c  tuesday wednesday thursday friday saturday sunday holiday jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c  jun*temp_!b!c^2 jul*temp_!b!c^2 aug*temp_!b!c^2 sep*temp_!b!c^2 jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'
'	All Day and Month Effects are Zero
'	Tue=Wed=Thu=Fri=Sat=Sun=Holiday=0
'	All Month Effects are Zero
'	Jul=Aug=Sep=0
'
table(1,1) model_desc
setcolwidth(model_desc,1,80)
 setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Day or Holiday or Month effects")
setcell(model_desc,3,1,"Tue = Wed = Thu = Fri = Sat = Sun = Holiday = Jul = Aug = Sep = 0")
print model_desc

!i=8
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c  jun*temp_!b!c^2 jul*temp_!b!c^2 aug*temp_!b!c^2 sep*temp_!b!c^2 jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'
'	All Weekend and Holiday Effects are Zero
'	Sat=Sun=Holiday=0
'	All Month Effects are Zero
'	Jul=Aug=Sep=0
'
table(1,1) model_desc
setcolwidth(model_desc,1,80)
 setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Weekend or Holiday or Month effects")
setcell(model_desc,3,1,"Sat = Sun = Holiday = Jul = Aug = Sep = 0")
print model_desc

!i=9
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c  tuesday wednesday thursday friday  jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c  jun*temp_!b!c^2 jul*temp_!b!c^2 aug*temp_!b!c^2 sep*temp_!b!c^2 jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'
'	All Weekend Effects are Zero
'	Sat=Sun=0
'	All Month Effects are Zero
'	Jul=Aug=Sep=0
'
table(1,1) model_desc
setcolwidth(model_desc,1,80)
 setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Weekend or Month effects")
setcell(model_desc,3,1,"Sat = Sun = Jul = Aug = Sep = 0")
print model_desc

!i=10
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c    tuesday wednesday thursday friday holiday   jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c  jun*temp_!b!c^2 jul*temp_!b!c^2 aug*temp_!b!c^2 sep*temp_!b!c^2 jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'
'	All Sunday and Holiday Effects are Zero
'	Sun=Holiday=0
'	All Month Effects are Zero
'	Jul=Aug=Sep=0
'
table(1,1) model_desc
setcolwidth(model_desc,1,80)
 setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Sunday or Holiday or Month Effects")
setcell(model_desc,3,1,"Sun = Holiday = Jul = Aug = Sep = 0")
print model_desc

!i=11
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c  tuesday wednesday thursday friday saturday   jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c  jun*temp_!b!c^2 jul*temp_!b!c^2 aug*temp_!b!c^2 sep*temp_!b!c^2 jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'
'	All Weekday Effects are the same There are no weekend effects
'	Mon=Tue=Wed=Thu=Fri,   Sat=Sun=0, Holiday unrestricted
'	All Month Effects are Zero
'	Jul=Aug=Sep=0
'
table(1,1) model_desc
setcolwidth(model_desc,1,80)
 setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Weekend or Month effects; Weekday effects are all the same")
setcell(model_desc,3,1,"Sat = Sun = Jul = Aug = Sep = 0; Tue = Wed = Thu = Fri")
print model_desc

!i=12
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c (tuesday+wednesday+thursday+friday) holiday   jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c  jun*temp_!b!c^2 jul*temp_!b!c^2 aug*temp_!b!c^2 sep*temp_!b!c^2 jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'	All Weekday Effects are Zero,  Mon=Tue=Wed=Thu=Fri=0,   
'	Sat=Sun=Holiday unrestricted
'
table(1,1) model_desc
setcolwidth(model_desc,1,80)
 setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Weekday effects")
setcell(model_desc,3,1,"Tue = Wed = Thu = Fri = 0")
print model_desc

!i=13
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c (weekend+holiday) jul aug sep  jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c  jun*temp_!b!c^2 jul*temp_!b!c^2 aug*temp_!b!c^2 sep*temp_!b!c^2 jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'	All Weekday Effects are the Zero , 	Mon=Tue=Wed=Thu=Fri=0  
'	Sat=Sun=Holiday unrestricted
'	All Month Effects are Zero
'	Jul=Aug=Sep=0
'
table(1,1) model_desc
setcolwidth(model_desc,1,80)
 setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Weekday or Month effects; Weekend and Holiday effects are the same")
setcell(model_desc,3,1,"Tue = Wed = Thu = Fri = Jul = Aug = Sep = 0; Sat = Sun = Holiday")
print model_desc

!i=14
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c (weekend+holiday)   jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c  jun*temp_!b!c^2 jul*temp_!b!c^2 aug*temp_!b!c^2 sep*temp_!b!c^2 jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'=========================*** Elasticity Levels Models =========================================================
'
'	Elasticity Levels
'	No Level Effects
'	Jun=Jul=Aug=Sep=0
'
table(1,1) model_desc
setcolwidth(model_desc,1,80)
setcell(model_desc,1,1,"Restrictions on Elasticity Levels")
print model_desc

table(1,1) model_desc
setcolwidth(model_desc,1,80)
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Temp Level effects")
setcell(model_desc,3,1,"Jun °F = Jul °F = Aug °F = Sep °F = 0")
print model_desc

!i=15
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   tuesday wednesday thursday friday saturday sunday holiday jul aug sep  jun*temp_!b!c^2 jul*temp_!b!c^2 aug*temp_!b!c^2 sep*temp_!b!c^2 jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test

'
'	Elasticity Levels
'	No Jun or Sep Level Effects
'	Jun=Sep=0
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Temp Level effects for June or September")
setcell(model_desc,3,1,"Jun °F = Sep °F = 0")
print model_desc

!i=16
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   tuesday wednesday thursday friday saturday sunday holiday jul aug sep  jul*temp_!b!c aug*temp_!b!c   jun*temp_!b!c^2 jul*temp_!b!c^2 aug*temp_!b!c^2 sep*temp_!b!c^2 jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'	Elasticity Levels
'	No June Level Effects
'	Jun=0
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Temp Level effects for June")
setcell(model_desc,3,1,"Jun °F = 0")
print model_desc

!i=17
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   tuesday wednesday thursday friday saturday sunday holiday jul aug sep   jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c  jun*temp_!b!c^2 jul*temp_!b!c^2 aug*temp_!b!c^2 sep*temp_!b!c^2 jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'	Elasticity Levels
'	No September Level Effects
'	Sep=0
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Temp Level effects for September")
setcell(model_desc,3,1,"Sep °F = 0")
print model_desc

!i=18
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   tuesday wednesday thursday friday saturday sunday holiday jul aug sep  jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c   jun*temp_!b!c^2 jul*temp_!b!c^2 aug*temp_!b!c^2 sep*temp_!b!c^2 jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'	Elasticity Levels
'	No September Level Effects Others are the same
'	Jun=Jul=Aug,  Sep=0
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Temp Level effects for September")
setcell(model_desc,3,1,"Temp Level effects for other months are all equal")
setcell(model_desc,4,1,"Sep °F = 0; Jun °F = Jul °F = Aug °F")
print model_desc

!i=19
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   tuesday wednesday thursday friday saturday sunday holiday jul aug sep  (jun+jul+aug)*temp_!b!c  jun*temp_!b!c^2 jul*temp_!b!c^2 aug*temp_!b!c^2 sep*temp_!b!c^2 jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'	Elasticity Levels
'	No June Level Effects, Others are the Same
'	Jun=0,    Jul=Aug=Sep
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Temp Level effects for June")
setcell(model_desc,3,1,"Temp Level effects for other months are all equal")
setcell(model_desc,4,1,"Jun °F = 0; Jul °F = Aug °F = Sep °F")
print model_desc

!i=20
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   tuesday wednesday thursday friday saturday sunday holiday jul aug sep (jul+aug+sep)*temp_!b!c  jun*temp_!b!c^2 jul*temp_!b!c^2 aug*temp_!b!c^2 sep*temp_!b!c^2 jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'	Elasticity Levels
'	No June and September Level Effects, July and August are the Same
'	Jun=0,  Jul=Aug, Sep=0
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Temp Level effects for June or September")
setcell(model_desc,3,1,"Temp Level effects for July and August are equal")
setcell(model_desc,4,1,"Jun °F = Sep °F = 0; Jul °F = Aug °F ")
print model_desc

!i=21
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   tuesday wednesday thursday friday saturday sunday holiday jul aug sep  (jul+aug)*temp_!b!c   jun*temp_!b!c^2 jul*temp_!b!c^2 aug*temp_!b!c^2 sep*temp_!b!c^2 jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'	Elasticity Levels
'	June and September Effects and July and August are the Same
'	Jun=Sep and Jul=Aug
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"Temp Level effects for June and September are equal")
setcell(model_desc,3,1,"Temp Level effects for July and August are equal")
setcell(model_desc,4,1,"Jun °F = Sep °F; Jul °F = Aug °F ")
print model_desc

!i=22
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   tuesday wednesday thursday friday saturday sunday holiday jul aug sep  (jun+sep)*temp_!b!c (jul+aug)*temp_!b!c  jun*temp_!b!c^2 jul*temp_!b!c^2 aug*temp_!b!c^2 sep*temp_!b!c^2 jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'	Elasticity Levels
'	Level Effects are the Same
'	Jun=Jul=Aug=Sep
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"Temp Level effects for all months are equal")
setcell(model_desc,3,1,"Jun °F = Jul °F = Aug °F = Sep °F")
print model_desc

!i=23
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   tuesday wednesday thursday friday saturday sunday holiday jul aug sep  (jun+jul+aug+sep)*temp_!b!c   jun*temp_!b!c^2 jul*temp_!b!c^2 aug*temp_!b!c^2 sep*temp_!b!c^2 jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test

'=========================*** Elasticity Squares Models =========================================================
'
'	Elasticity Squares
'	No SquaresEffects
'	Jun=Jul=Aug=Sep=0
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Restrictions on Elasticity Squares")
print model_desc

table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Temp Squared effects")
setcell(model_desc,3,1,"Jun °F^2 = Jul °F^2 = Aug °F^2 = Sep °F^2 = 0")
print model_desc

!i=24
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   tuesday wednesday thursday friday saturday sunday holiday jul aug sep jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c  jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'	Elasticity Squares
'	No Jun or Sep SquaresEffects
'	Jun=Sep=0
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Temp Squared effects for June or September")
setcell(model_desc,3,1,"Jun °F^2 = Sep °F^2 = 0")
print model_desc

!i=25
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   tuesday wednesday thursday friday saturday sunday holiday jul aug sep  jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c  jul*temp_!b!c^2 aug*temp_!b!c^2 jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'	Elasticity Squares
'	No June SquaresEffects
'	Jun=0
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Temp Squared effects for June")
setcell(model_desc,3,1,"Jun °F^2 = 0")
print model_desc

!i=26
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   tuesday wednesday thursday friday saturday sunday holiday jul aug sep jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c  jul*temp_!b!c^2 aug*temp_!b!c^2 sep*temp_!b!c^2 jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'	Elasticity Squares
'	No September SquaresEffects
'	Sep=0
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Temp Squared effects for September")
setcell(model_desc,3,1,"Sep °F^2 = 0")
print model_desc

!i=27
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   tuesday wednesday thursday friday saturday sunday holiday jul aug sep jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c  jun*temp_!b!c^2 jul*temp_!b!c^2 aug*temp_!b!c^2  jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'	Elasticity Squares
'	No September SquaresEffects Others are the same
'	Jun=Jul=Aug,  Sep=0
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Temp Squared effects for September")
setcell(model_desc,3,1,"Temp Squared effects for other months are all equal")
setcell(model_desc,4,1,"Sep °F^2 = 0; Jun °F^2 = Jul °F^2 = Aug °F^2")
print model_desc

!i=28
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   tuesday wednesday thursday friday saturday sunday holiday jul aug  jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c  (jun+jul+aug)*temp_!b!c^2  jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'	Elasticity Squares
'	No June SquaresEffects, Others are the Same
'	Jun=0,    Jul=Aug=Sep
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Temp Squared effects for June")
setcell(model_desc,3,1,"Temp Squared effects for other months are all equal")
setcell(model_desc,4,1,"Jun °F^2 = 0; Jul °F^2 = Aug °F^2 = Sep °F^2")
print model_desc

!i=29
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   tuesday wednesday thursday friday saturday sunday holiday jul aug jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c  (jul+aug+sep)*temp_!b!c^2  jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'	Elasticity Squares
'	No June and September SquaresEffects, July and August are the Same
'	Jun=0,  Jul=Aug, Sep=0
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Temp Squared effects for June or September")
setcell(model_desc,3,1,"Temp Squared effects for July and August are equal")
setcell(model_desc,4,1,"Jun °F^2 = Sep °F^2 = 0; Jul °F^2 = Aug °F^2 ")
print model_desc

!i=30
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   tuesday wednesday thursday friday saturday sunday holiday jul aug sep jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c   (jul+aug)*temp_!b!c^2  jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'	Elasticity Squares
'	June and September Effects and July and August are the Same
'	Jun=Sep, and Jul=Aug
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"Temp Squared effects for June and September are equal")
setcell(model_desc,3,1,"Temp Squared effects for July and August are equal")
setcell(model_desc,4,1,"Jun °F^2 = Sep °F^2; Jul °F^2 = Aug °F^2")
print model_desc

!i=31
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   tuesday wednesday thursday friday saturday sunday holiday jul aug sep jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c   (jun+sep)*temp_!b!c^2 (jul+aug)*temp_!b!c^2 jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'	Elasticity Squares
'	SquaresEffects are the Same
'	Jun=Jul=Aug=Sep
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"Temp Squared effects for all months are equal")
setcell(model_desc,3,1,"Jun °F^2 = Jul °F^2 = Aug °F^2 = Sep °F^2")
print model_desc

!i=32
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   tuesday wednesday thursday friday saturday sunday holiday jul aug sep jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c  (jun+jul+aug+sep)*temp_!b!c^2 jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test

'=========================*** Elasticity Joint Levels and Squares Models =========================================================
'
'	Elasticity Levels and Squares
'	No Level and Squares Effects
'	Jun=Jul=Aug=Sep=0
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Restrictions on Elasticity Levels and Squares")
print model_desc

table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Temp Level or Squared effects")
setcell(model_desc,3,1,"Jun °F = Jul °F = Aug °F = Sep °F = ")
setcell(model_desc,4,1,"Jun °F^2 = Jul °F^2 = Aug °F^2 = Sep °F^2 = 0")
print model_desc

!i=33
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   tuesday wednesday thursday friday saturday sunday holiday jul aug sep   jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test

'
'	Elasticity Levels and Squares
'	No Jun or Sep Effects
'	Jun=Sep=0
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Temp Level or Squared effects for June or September")
setcell(model_desc,3,1,"Jun °F = Sep °F = Jun °F^2 = Sep °F^2 = 0")
print model_desc

!i=34
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   tuesday wednesday thursday friday saturday sunday holiday jul aug sep  jul*temp_!b!c aug*temp_!b!c   jul*temp_!b!c^2 aug*temp_!b!c^2 jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'	Elasticity Levels and Squares
'	No June Effects
'	Jun=0
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Temp Level or Squared effects for June")
setcell(model_desc,3,1,"Jun °F = Jun °F^2 = 0")
print model_desc

!i=35
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   tuesday wednesday thursday friday saturday sunday holiday jul aug sep jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c  jul*temp_!b!c^2 aug*temp_!b!c^2 sep*temp_!b!c^2 jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'	Elasticity Levels and Squares
'	No September Effects
'	Sep=0
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Temp Level or Squared effects for September")
setcell(model_desc,3,1,"Sep °F = Sep °F^2 = 0")
print model_desc

!i=36
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   tuesday wednesday thursday friday saturday sunday holiday jul aug sep jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c  jun*temp_!b!c^2 jul*temp_!b!c^2 aug*temp_!b!c^2  jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'	Elasticity Levels and Squares
'	No September Effects Others are the same
'	Jun=Jul=Aug,  Sep=0
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Temp Level or Squared effects for September")
setcell(model_desc,3,1,"Temp effects for other months are all equal")
setcell(model_desc,4,1,"Sep °F = 0; Jun °F = Jul °F = Aug °F")
setcell(model_desc,5,1,"Sep °F^2 = 0; Jun °F^2 = Jul °F^2 = Aug °F^2")
print model_desc

!i=37
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   tuesday wednesday thursday friday saturday sunday holiday jul aug  (jun+jul+aug)*temp_!b!c  (jun+jul+aug)*temp_!b!c^2  jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'	Elasticity Levels and Squares
'	No June Effects, Others are the Same
'	Jun=0,    Jul=Aug=Sep
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"Temp effects for other months are all equal")
setcell(model_desc,3,1,"Temp effects for other months are all equal")
setcell(model_desc,4,1,"Jun °F = 0; Jul °F = Aug °F = Sep °F")
setcell(model_desc,5,1,"Jun °F^2 = 0; Jul °F^2 = Aug °F^2 = Sep °F^2")
print model_desc

!i=38
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   tuesday wednesday thursday friday saturday sunday holiday jul aug (jul+aug+sep)*temp_!b!c   (jul+aug+sep)*temp_!b!c^2  jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'	Elasticity Levels and Squares
'	No June and September Effects, July and August are the Same
'	Jun=0,  Jul=Aug, Sep=0
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"There are no Temp Level or Squared effects for June or September")
setcell(model_desc,3,1,"Temp effects for July and August are equal")
setcell(model_desc,4,1,"Jun °F = Sep °F = 0; Jul °F = Aug °F")
setcell(model_desc,5,1,"Jun °F^2 = Sep °F^2 = 0; Jul °F^2 = Aug °F^2")
print model_desc

!i=39
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   tuesday wednesday thursday friday saturday sunday holiday jul aug sep (jul+aug)*temp_!b!c    (jul+aug)*temp_!b!c^2  jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'	Elasticity Levels and Squares
'	June and September Effects and July and August Effects are the Same
'	Jun=Sep, and Jul=Aug
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"Temp effects for June and September are equal; Temp effects for July and August are equal")
setcell(model_desc,3,1,"Jun °F = Sep °F; Jul °F = Aug °F")
setcell(model_desc,4,1,"Jun °F^2 = Sep °F^2; Jul °F^2 = Aug °F^2")
print model_desc

!i=40
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   tuesday wednesday thursday friday saturday sunday holiday jul aug sep (jun+sep)*temp_!b!c (jul+aug)*temp_!b!c    (jun+sep)*temp_!b!c^2 (jul+aug)*temp_!b!c^2 jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'	Elasticity Levels and Squares
'	Level Effects are the Same
'	Jun=Jul=Aug=Sep
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"Temp level and square effects for all months are equal")
setcell(model_desc,3,1,"Jun °F = Jul °F = Aug °F = Sep °F")
setcell(model_desc,4,1,"Jun °F^2 = Jul °F^2 = Aug °F^2 = Sep °F^2")
print model_desc

!i=41
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   tuesday wednesday thursday friday saturday sunday holiday jul aug sep (jun+jul+aug+sep)*temp_!b!c   (jun+jul+aug+sep)*temp_!b!c^2 jun*temp_!b!c^3 jul*temp_!b!c^3 aug*temp_!b!c^3 sep*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test


'=========================*** Elasticity Cubes Models =========================================================
'
'	Elasticity Cubes
'	No Cubic Effects
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Restrictions on Elasticity Cubes")
print model_desc

table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"No Cubic Effects")
setcell(model_desc,3,1,"Jun °F^3 = Jul °F^3 = Aug °F^3 = Sep °F^3 = 0")
print model_desc

!i=42

equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c tuesday wednesday thursday friday saturday sunday holiday jul aug sep  jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c  jun*temp_!b!c^2 jul*temp_!b!c^2 aug*temp_!b!c^2 sep*temp_!b!c^2  log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test

'
'	Elasticity Cubes
'	No Jun or Sep Cubic Effects
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"No Jun/Sept Cubic Effects")
setcell(model_desc,3,1,"Jun °F^3 = Sep °F^3 = 0")
print model_desc

!i=43

equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c tuesday wednesday thursday friday saturday sunday holiday jul aug sep  jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c  jun*temp_!b!c^2 jul*temp_!b!c^2 aug*temp_!b!c^2 sep*temp_!b!c^2  jul*temp_!b!c^3 aug*temp_!b!c^3  log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test

'
'	Elasticity Cubes
'	Cubic Effects for Jun = Sep; Jul = Aug
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"Cubic Effects for Jun = Sep; Jul = Aug")
setcell(model_desc,3,1,"Jun °F^3 = Sep °F^3; Jul °F^3 = Aug °F^3")
print model_desc

!i=44

equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c tuesday wednesday thursday friday saturday sunday holiday jul aug sep  jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c  jun*temp_!b!c^2 jul*temp_!b!c^2 aug*temp_!b!c^2 sep*temp_!b!c^2 (jun+sep)*temp_!b!c^3 (jul+aug)*temp_!b!c^3  log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test

'
'	Elasticity Cubes
'	Cubic Effects are all the same
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"Cubic Effects are the same across all months")
setcell(model_desc,3,1,"Jun °F^3 = Sep °F^3 = Jul °F^3 = Aug °F^3")
print model_desc

!i=45

equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c tuesday wednesday thursday friday saturday sunday holiday jul aug sep  jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c  jun*temp_!b!c^2 jul*temp_!b!c^2 aug*temp_!b!c^2 sep*temp_!b!c^2 (jun+jul+aug+sep)*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test



'=========================*** Elasticity Joint Squares and Cubes Models =========================================================
'
'	Elasticity Squares and Cubes
'	No Squared/Cubic Effects
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Restrictions on Elasticity Squares and Cubes")
print model_desc

table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"No Squared/Cubic Effect")
setcell(model_desc,3,1,"Jun °F^2 = Jul °F^2 = Aug °F^2 = Sep °F^2 = ")
setcell(model_desc,4,1,"Jun °F^3 = Jul °F^3 = Aug °F^3 = Sep °F^3 = 0")
print model_desc

!i=46
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   tuesday wednesday thursday friday saturday sunday holiday jul aug sep jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test
'
'	Elasticity Squares and Cubes
'	No Jun or Sep Squared/Cubic Effects
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"No Jun or Sep Squared/Cubic Effects")
setcell(model_desc,3,1,"Jun °F^2 = Sep °F^2 = Jun °F^3 = Sep °F^3 = 0")
print model_desc

!i=47

equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c tuesday wednesday thursday friday saturday sunday holiday jul aug sep  jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c  jul*temp_!b!c^2 aug*temp_!b!c^2  log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test

'
'	Elasticity Squares and Cubes
'	Squared/Cubic Effects for Jun = Sep; Jul = Aug
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"Squared/Cubic Effects for Jun = Sep; Jul = Aug")
setcell(model_desc,3,1,"Jun °F^2 = Sep °F^2; Jul °F^2 = Aug °F^2")
setcell(model_desc,4,1,"Jun °F^3 = Sep °F^3; Jul °F^3 = Aug °F^3")
print model_desc

!i=48

equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c tuesday wednesday thursday friday saturday sunday holiday jul aug sep  jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c  (jun+sep)*temp_!b!c^2 (jul+aug)*temp_!b!c^2 (jun+sep)*temp_!b!c^3 (jul+aug)*temp_!b!c^3  log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test

'
'	Elasticity Squares and Cubes
'	Squared/Cubic Effects are the same in all months
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"Squared/Cubic Effects are the same in all months")
setcell(model_desc,3,1,"Jun °F^2 = Jul °F^2 = Aug °F^2 = Sep °F^2")
setcell(model_desc,4,1,"Jun °F^3 = Jul °F^3 = Aug °F^3 = Sep °F^3")
print model_desc

!i=49

equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c tuesday wednesday thursday friday saturday sunday holiday jul aug sep  jun*temp_!b!c jul*temp_!b!c aug*temp_!b!c sep*temp_!b!c  (jun+jul+aug+sep)*temp_!b!c^2 (jun+jul+aug+sep)*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test





'=========================*** Elasticity Joint Levels, Squares and Cubes Models =========================================================
'
'	Elasticity Levels, Squares and Cubes
'	No Level/Squared/Cubic Effects
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Restrictions on Elasticity Levels, Squares and Cubes")
print model_desc

table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"No Level/Squared/Cubic Effects")
setcell(model_desc,3,1,"Jun °F = Jul °F = Aug °F = Sep °F = ")
setcell(model_desc,4,1,"Jun °F^2 = Jul °F^2 = Aug °F^2 = Sep °F^2 = ")
setcell(model_desc,5,1,"Jun °F^3 = Jul °F^3 = Aug °F^3 = Sep °F^3 = 0")
print model_desc

!i=50

equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   tuesday wednesday thursday friday saturday sunday holiday jul aug sep log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test

'
'	Elasticity Levels, Squares and Cubes
'	No Jun or Sep Level/Squared/Cubic Effects
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"No Jun or Sep Level/Squared/Cubic Effects")
setcell(model_desc,3,1,"Jun °F = Sep °F = Jun °F^2 = Sep °F^2 = Jun °F^3 = Sep °F^3 = 0")
print model_desc

!i=51

equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c tuesday wednesday thursday friday saturday sunday holiday jul aug sep  jul*temp_!b!c aug*temp_!b!c jul*temp_!b!c^2 aug*temp_!b!c^2 jul*temp_!b!c^3 aug*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test

'
'	Elasticity Levels, Squares and Cubes
'	Level/Squared/Cubic Effects for Jun = Sep; Jul = Aug
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"Level/Squared/Cubic Effects for Jun = Sep; Jul = Aug")
setcell(model_desc,3,1,"Jun °F = Sep °F; Jul °F = Aug °F")
setcell(model_desc,4,1,"Jun °F^2 = Sep °F^2; Jul °F^2 = Aug °F^2")
setcell(model_desc,5,1,"Jun °F^3 = Sep °F^3; Jul °F^3 = Aug °F^3")
print model_desc

!i=52

equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c tuesday wednesday thursday friday saturday sunday holiday jul aug sep  (jun+sep)*temp_!b!c (jul+aug)*temp_!b!c (jun+sep)*temp_!b!c^2 (jul+aug)*temp_!b!c^2 (jun+sep)*temp_!b!c^3 (jul+aug)*temp_!b!c^3  log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test


'	Elasticity Levels, Squares and Cubes
'	Level/Squared/Cubic Effects are the Same
'
table(1,1) model_desc
setcolwidth(model_desc,1,80) 
setcell(model_desc,1,1,"Model Number: " + @str(!i))
setcell(model_desc,2,1,"Level/Squared/Cubic Effects are the Same")
setcell(model_desc,3,1,"Jun °F = Jul °F = Aug °F = Sep °F")
setcell(model_desc,4,1,"Jun °F^2 = Jul °F^2 = Aug °F^2 = Sep °F^2")
setcell(model_desc,5,1,"Jun °F^3 = Jul °F^3 = Aug °F^3 = Sep °F^3")
print model_desc

!i=53
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   tuesday wednesday thursday friday saturday sunday holiday jul aug sep  (jun+jul+aug+sep)*temp_!b!c (jun+jul+aug+sep)*temp_!b!c^2 (jun+jul+aug+sep)*temp_!b!c^3 log(load_!m!n(!day1)) log(load_!o!p(!day2)) log(load_!q!r(!day3)) log(load_!b!c(-1))

call lr_test



'	NEXT Z:
next

delete model_desc
smpl @all

stop