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 + "_winter_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 + "_winter_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 + "_winter_models"

endif

'	c:\epa\hourly_models\PJM_Utilities\AE_AtlanticCity,NJ_ACY\ae_cddtemp_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 = jan + feb + mar + dec
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 1/1/1995 3/31/1995  12/1/1995 3/31/1996  12/1/1996 3/31/1997  12/1/1997 3/31/1998  12/1/1998 3/31/1999  12/1/1999 3/31/2000  12/1/2000 3/31/2001  12/1/2001 3/31/2002  12/1/2002 3/31/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 feb mar dec jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c  jan*temp_!b!c^2 feb*temp_!b!c^2 mar*temp_!b!c^2 dec*temp_!b!c^2 jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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    feb mar dec  jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c  jan*temp_!b!c^2 feb*temp_!b!c^2 mar*temp_!b!c^2 dec*temp_!b!c^2 jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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  feb mar dec  jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c  jan*temp_!b!c^2 feb*temp_!b!c^2 mar*temp_!b!c^2 dec*temp_!b!c^2 jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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 feb mar dec  jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c  jan*temp_!b!c^2 feb*temp_!b!c^2 mar*temp_!b!c^2 dec*temp_!b!c^2 jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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  feb mar dec  jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c  jan*temp_!b!c^2 feb*temp_!b!c^2 mar*temp_!b!c^2 dec*temp_!b!c^2 jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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   feb mar dec  jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c  jan*temp_!b!c^2 feb*temp_!b!c^2 mar*temp_!b!c^2 dec*temp_!b!c^2 jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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
'	Feb=Mar=Dec=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,"Feb = Mar = Dec = 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  jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c  jan*temp_!b!c^2 feb*temp_!b!c^2 mar*temp_!b!c^2 dec*temp_!b!c^2 jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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
'	Feb=Mar=Dec=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 = Feb = Mar = Dec = 0")
print model_desc

!i=8
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c   jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c  jan*temp_!b!c^2 feb*temp_!b!c^2 mar*temp_!b!c^2 dec*temp_!b!c^2 jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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
'	Feb=Mar=Dec=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 = Feb = Mar = Dec = 0")
print model_desc

!i=9
equation eq_h!b!cr!i.ls(n,p) log(load_!b!c) c  tuesday wednesday thursday friday  jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c  jan*temp_!b!c^2 feb*temp_!b!c^2 mar*temp_!b!c^2 dec*temp_!b!c^2 jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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
'	Feb=Mar=Dec=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 = Feb = Mar = Dec = 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   jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c  jan*temp_!b!c^2 feb*temp_!b!c^2 mar*temp_!b!c^2 dec*temp_!b!c^2 jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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
'	Feb=Mar=Dec=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 = Feb = Mar = Dec = 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   jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c  jan*temp_!b!c^2 feb*temp_!b!c^2 mar*temp_!b!c^2 dec*temp_!b!c^2 jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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
'	Feb=Mar=Dec=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 = Feb = Mar = Dec = 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   jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c  jan*temp_!b!c^2 feb*temp_!b!c^2 mar*temp_!b!c^2 dec*temp_!b!c^2 jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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) feb mar dec  jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c  jan*temp_!b!c^2 feb*temp_!b!c^2 mar*temp_!b!c^2 dec*temp_!b!c^2 jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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
'	Feb=Mar=Dec=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 = Feb = Mar = Dec = 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)   jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c  jan*temp_!b!c^2 feb*temp_!b!c^2 mar*temp_!b!c^2 dec*temp_!b!c^2 jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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
'	Jan=Feb=Mar=Dec=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,"Jan °F = Feb °F = Mar °F = Dec °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 feb mar dec  jan*temp_!b!c^2 feb*temp_!b!c^2 mar*temp_!b!c^2 dec*temp_!b!c^2 jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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 Jan or Dec Level Effects
'	Jan=Dec=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 January or December")
setcell(model_desc,3,1,"Jan °F = Dec °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 feb mar dec  feb*temp_!b!c mar*temp_!b!c   jan*temp_!b!c^2 feb*temp_!b!c^2 mar*temp_!b!c^2 dec*temp_!b!c^2 jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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 January Level Effects
'	Jan=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 January")
setcell(model_desc,3,1,"Jan °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 feb mar dec   feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c  jan*temp_!b!c^2 feb*temp_!b!c^2 mar*temp_!b!c^2 dec*temp_!b!c^2 jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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 December Level Effects
'	Dec=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 December")
setcell(model_desc,3,1,"Dec °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 feb mar dec  jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c   jan*temp_!b!c^2 feb*temp_!b!c^2 mar*temp_!b!c^2 dec*temp_!b!c^2 jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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 December Level Effects Others are the same
'	Jan=Feb=Mar,  Dec=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 December")
setcell(model_desc,3,1,"Temp Level effects for other months are all equal")
setcell(model_desc,4,1,"Dec °F = 0; Jan °F = Feb °F = Mar °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 feb mar dec  (jan+feb+mar)*temp_!b!c  jan*temp_!b!c^2 feb*temp_!b!c^2 mar*temp_!b!c^2 dec*temp_!b!c^2 jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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 January Level Effects, Others are the Same
'	Jan=0,    Feb=Mar=Dec
'
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 January")
setcell(model_desc,3,1,"Temp Level effects for other months are all equal")
setcell(model_desc,4,1,"Jan °F = 0; Feb °F = Mar °F = Dec °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 feb mar dec (feb+mar+dec)*temp_!b!c  jan*temp_!b!c^2 feb*temp_!b!c^2 mar*temp_!b!c^2 dec*temp_!b!c^2 jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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 January and December Level Effects, February and March are the Same
'	Jan=0,  Feb=Mar, Dec=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 January or December")
setcell(model_desc,3,1,"Temp Level effects for February and March are equal")
setcell(model_desc,4,1,"Jan °F = Dec °F = 0; Feb °F = Mar °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 feb mar dec  (feb+mar)*temp_!b!c   jan*temp_!b!c^2 feb*temp_!b!c^2 mar*temp_!b!c^2 dec*temp_!b!c^2 jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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
'	January and December Effects and February and March are the Same
'	Jan=Dec and Feb=Mar
'
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 January and December are equal")
setcell(model_desc,3,1,"Temp Level effects for February and March are equal")
setcell(model_desc,4,1,"Jan °F = Dec °F; Feb °F = Mar °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 feb mar dec  (jan+dec)*temp_!b!c (feb+mar)*temp_!b!c  jan*temp_!b!c^2 feb*temp_!b!c^2 mar*temp_!b!c^2 dec*temp_!b!c^2 jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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
'	Jan=Feb=Mar=Dec
'
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,"Jan °F = Feb °F = Mar °F = Dec °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 feb mar dec  (jan+feb+mar+dec)*temp_!b!c   jan*temp_!b!c^2 feb*temp_!b!c^2 mar*temp_!b!c^2 dec*temp_!b!c^2 jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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
'	Jan=Feb=Mar=Dec=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,"Jan °F^2 = Feb °F^2 = Mar °F^2 = Dec °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 feb mar dec jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c  jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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 Jan or Dec SquaresEffects
'	Jan=Dec=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 January or December")
setcell(model_desc,3,1,"Jan °F^2 = Dec °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 feb mar dec  jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c  feb*temp_!b!c^2 mar*temp_!b!c^2 jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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 January SquaresEffects
'	Jan=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 January")
setcell(model_desc,3,1,"Jan °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 feb mar dec jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c  feb*temp_!b!c^2 mar*temp_!b!c^2 dec*temp_!b!c^2 jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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 December SquaresEffects
'	Dec=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 December")
setcell(model_desc,3,1,"Dec °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 feb mar dec jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c  jan*temp_!b!c^2 feb*temp_!b!c^2 mar*temp_!b!c^2  jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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 December SquaresEffects Others are the same
'	Jan=Feb=Mar,  Dec=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 December")
setcell(model_desc,3,1,"Temp Squared effects for other months are all equal")
setcell(model_desc,4,1,"Dec °F^2 = 0; Jan °F^2 = Feb °F^2 = Mar °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 feb mar  jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c  (jan+feb+mar)*temp_!b!c^2  jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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 January SquaresEffects, Others are the Same
'	Jan=0,    Feb=Mar=Dec
'
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 January")
setcell(model_desc,3,1,"Temp Squared effects for other months are all equal")
setcell(model_desc,4,1,"Jan °F^2 = 0; Feb °F^2 = Mar °F^2 = Dec °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 feb mar jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c  (feb+mar+dec)*temp_!b!c^2  jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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 January and December SquaresEffects, February and March are the Same
'	Jan=0,  Feb=Mar, Dec=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 January or December")
setcell(model_desc,3,1,"Temp Squared effects for February and March are equal")
setcell(model_desc,4,1,"Jan °F^2 = Dec °F^2 = 0; Feb °F^2 = Mar °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 feb mar dec jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c   (feb+mar)*temp_!b!c^2  jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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
'	January and December Effects and February and March are the Same
'	Jan=Dec, and Feb=Mar
'
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 January and December are equal")
setcell(model_desc,3,1,"Temp Squared effects for February and March are equal")
setcell(model_desc,4,1,"Jan °F^2 = Dec °F^2; Feb °F^2 = Mar °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 feb mar dec jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c   (jan+dec)*temp_!b!c^2 (feb+mar)*temp_!b!c^2 jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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
'	Jan=Feb=Mar=Dec
'
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,"Jan °F^2 = Feb °F^2 = Mar °F^2 = Dec °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 feb mar dec jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c  (jan+feb+mar+dec)*temp_!b!c^2 jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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
'	Jan=Feb=Mar=Dec=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,"Jan °F = Feb °F = Mar °F = Dec °F = ")
setcell(model_desc,4,1,"Jan °F^2 = Feb °F^2 = Mar °F^2 = Dec °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 feb mar dec   jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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 Jan or Dec Effects
'	Jan=Dec=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 January or December")
setcell(model_desc,3,1,"Jan °F = Dec °F = Jan °F^2 = Dec °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 feb mar dec  feb*temp_!b!c mar*temp_!b!c   feb*temp_!b!c^2 mar*temp_!b!c^2 jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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 January Effects
'	Jan=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 January")
setcell(model_desc,3,1,"Jan °F = Jan °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 feb mar dec feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c  feb*temp_!b!c^2 mar*temp_!b!c^2 dec*temp_!b!c^2 jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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 December Effects
'	Dec=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 December")
setcell(model_desc,3,1,"Dec °F = Dec °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 feb mar dec jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c  jan*temp_!b!c^2 feb*temp_!b!c^2 mar*temp_!b!c^2  jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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 December Effects Others are the same
'	Jan=Feb=Mar,  Dec=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 December")
setcell(model_desc,3,1,"Temp effects for other months are all equal")
setcell(model_desc,4,1,"Dec °F = 0; Jan °F = Feb °F = Mar °F")
setcell(model_desc,5,1,"Dec °F^2 = 0; Jan °F^2 = Feb °F^2 = Mar °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 feb mar  (jan+feb+mar)*temp_!b!c  (jan+feb+mar)*temp_!b!c^2  jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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 January Effects, Others are the Same
'	Jan=0,    Feb=Mar=Dec
'
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 January")
setcell(model_desc,3,1,"Temp effects for other months are all equal")
setcell(model_desc,4,1,"Jan °F = 0; Feb °F = Mar °F = Dec °F")
setcell(model_desc,5,1,"Jan °F^2 = 0; Feb °F^2 = Mar °F^2 = Dec °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 feb mar (feb+mar+dec)*temp_!b!c   (feb+mar+dec)*temp_!b!c^2  jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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 January and December Effects, February and March are the Same
'	Jan=0,  Feb=Mar, Dec=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 January or December")
setcell(model_desc,3,1,"Temp effects for February and March are equal")
setcell(model_desc,4,1,"Jan °F = Dec °F = 0; Feb °F = Mar °F")
setcell(model_desc,5,1,"Jan °F^2 = Dec °F^2 = 0; Feb °F^2 = Mar °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 feb mar dec (feb+mar)*temp_!b!c    (feb+mar)*temp_!b!c^2  jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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
'	January and December Effects and February and March Effects are the Same
'	Jan=Dec, and Feb=Mar
'
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 January and December are equal; Temp effects for February and March are equal")
setcell(model_desc,3,1,"Jan °F = Dec °F; Feb °F = Mar °F")
setcell(model_desc,4,1,"Jan °F^2 = Dec °F^2; Feb °F^2 = Mar °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 feb mar dec (jan+dec)*temp_!b!c (feb+mar)*temp_!b!c    (jan+dec)*temp_!b!c^2 (feb+mar)*temp_!b!c^2 jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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
'	Jan=Feb=Mar=Dec
'
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,"Jan °F = Feb °F = Mar °F = Dec °F")
setcell(model_desc,4,1,"Jan °F^2 = Feb °F^2 = Mar °F^2 = Dec °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 feb mar dec (jan+feb+mar+dec)*temp_!b!c   (jan+feb+mar+dec)*temp_!b!c^2 jan*temp_!b!c^3 feb*temp_!b!c^3 mar*temp_!b!c^3 dec*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,"Jan °F^3 = Feb °F^3 = Mar °F^3 = Dec °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 feb mar dec  jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c  jan*temp_!b!c^2 feb*temp_!b!c^2 mar*temp_!b!c^2 dec*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 Jan or Dec 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 Jan/Sept Cubic Effects")
setcell(model_desc,3,1,"Jan °F^3 = Dec °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 feb mar dec  jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c  jan*temp_!b!c^2 feb*temp_!b!c^2 mar*temp_!b!c^2 dec*temp_!b!c^2  feb*temp_!b!c^3 mar*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 Jan = Dec; Feb = Mar
'
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 Jan = Dec; Feb = Mar")
setcell(model_desc,3,1,"Jan °F^3 = Dec °F^3; Feb °F^3 = Mar °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 feb mar dec  jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c  jan*temp_!b!c^2 feb*temp_!b!c^2 mar*temp_!b!c^2 dec*temp_!b!c^2 (jan+dec)*temp_!b!c^3 (feb+mar)*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,"Jan °F^3 = Dec °F^3 = Feb °F^3 = Mar °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 feb mar dec  jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c  jan*temp_!b!c^2 feb*temp_!b!c^2 mar*temp_!b!c^2 dec*temp_!b!c^2 (jan+feb+mar+dec)*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,"Jan °F^2 = Feb °F^2 = Mar °F^2 = Dec °F^2 = ")
setcell(model_desc,4,1,"Jan °F^3 = Feb °F^3 = Mar °F^3 = Dec °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 feb mar dec jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c dec*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 Jan or Dec 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 Jan or Dec Squared/Cubic Effects")
setcell(model_desc,3,1,"Jan °F^2 = Dec °F^2 = Jan °F^3 = Dec °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 feb mar dec  jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c  feb*temp_!b!c^2 mar*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 Jan = Dec; Feb = Mar
'
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 Jan = Dec; Feb = Mar")
setcell(model_desc,3,1,"Jan °F^2 = Dec °F^2; Feb °F^2 = Mar °F^2")
setcell(model_desc,4,1,"Jan °F^3 = Dec °F^3; Feb °F^3 = Mar °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 feb mar dec  jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c  (jan+dec)*temp_!b!c^2 (feb+mar)*temp_!b!c^2 (jan+dec)*temp_!b!c^3 (feb+mar)*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,"Jan °F^2 = Feb °F^2 = Mar °F^2 = Dec °F^2")
setcell(model_desc,4,1,"Jan °F^3 = Feb °F^3 = Mar °F^3 = Dec °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 feb mar dec  jan*temp_!b!c feb*temp_!b!c mar*temp_!b!c dec*temp_!b!c  (jan+feb+mar+dec)*temp_!b!c^2 (jan+feb+mar+dec)*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,"Jan °F = Feb °F = Mar °F = Dec °F = ")
setcell(model_desc,4,1,"Jan °F^2 = Feb °F^2 = Mar °F^2 = Dec °F^2 = ")
setcell(model_desc,5,1,"Jan °F^3 = Feb °F^3 = Mar °F^3 = Dec °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 feb mar dec 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 Jan or Dec 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 Jan or Dec Level/Squared/Cubic Effects")
setcell(model_desc,3,1,"Jan °F = Dec °F = Jan °F^2 = Dec °F^2 = Jan °F^3 = Dec °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 feb mar dec  feb*temp_!b!c mar*temp_!b!c feb*temp_!b!c^2 mar*temp_!b!c^2 feb*temp_!b!c^3 mar*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 Jan = Dec; Feb = Mar
'
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 Jan = Dec; Feb = Mar")
setcell(model_desc,3,1,"Jan °F = Dec °F; Feb °F = Mar °F")
setcell(model_desc,4,1,"Jan °F^2 = Dec °F^2; Feb °F^2 = Mar °F^2")
setcell(model_desc,5,1,"Jan °F^3 = Dec °F^3; Feb °F^3 = Mar °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 feb mar dec  (jan+dec)*temp_!b!c (feb+mar)*temp_!b!c (jan+dec)*temp_!b!c^2 (feb+mar)*temp_!b!c^2 (jan+dec)*temp_!b!c^3 (feb+mar)*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,"Jan °F = Feb °F = Mar °F = Dec °F")
setcell(model_desc,4,1,"Jan °F^2 = Feb °F^2 = Mar °F^2 = Dec °F^2")
setcell(model_desc,5,1,"Jan °F^3 = Feb °F^3 = Mar °F^3 = Dec °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 feb mar dec  (jan+feb+mar+dec)*temp_!b!c (jan+feb+mar+dec)*temp_!b!c^2 (jan+feb+mar+dec)*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