' eViews program to perform leave-one-out validation
' The scalars AC, NB, NL, can be changed to select the desired model
' koyck1_no indicates the specific model (Model A - I)

'Indicate which model to run
scalar AC	= 1											' 1 = no autocorrelation, 2 = unrestricted autocorrelation, 3 = restricted autocorrelation
scalar NB	= 1											' 1 = one beta, 3 = three betas
scalar NL	= 1											' 1 = one lambda, 3 = three lambdas

' Select correct model to estimate
' Koyck model with one prom_cost
if AC = 1 and NB = 1 and NL = 1 then					' Koyck model with 1 prom_cost and no autocorrelation				
for !i= 5 to 236 
    smpl 5 !i !i+5 241
    koyck1_no.ml(showopts, m=1000, c=1e-5)
    show koyck1_no.output

    forecasterror(!i+1) = attendance(!i+1)  - ( cc(1) + cc(2) * prom_cost_sum(!i+1) + cc(5) * attendance(!i) + cc(15) * Season11(!i+1) + cc(16) * season12(!i+1) + cc(17) * opening11(!i+1) + cc(18) * opening12(!i+1) + cc(23) * day5(!i+1) + cc(24) * day6(!i+1) + cc(27) * day7(!i+1) + cc(28) * july(!i+1) + cc(29) * august(!i+1) + cc(30) * sep(!i+1) + cc(31) * perf_wp(!i+1) + cc(32) * opp_distance(!i+1) + cc(50) * temp(!i+1) +  cc(60) * nattv(!i+1) + cc(80) * Control_entertainment(!i+1) + cc(81) * control_giveaway(!i+1) + cc(82) * control_kids(!i+1) - cc(5) * ( cc(15) * Season11(!i) + cc(16) * season12(!i) + cc(17) * opening11(!i) + cc(18) * opening12(!i) + cc(23) * day5(!i) + cc(24) * day6(!i) + cc(27) * day7(!i) + cc(28) * july(!i) + cc(29) * august(!i) + cc(30) * sep(!i) + cc(31) * perf_wp(!i) + cc(32) * opp_distance(!i) + cc(50) * temp(!i) + cc(60) * nattv(!i) + cc(80) * Control_entertainment(!i) + cc(81) * control_giveaway(!i) + cc(82) * control_kids(!i) ) )

       next
endif

if AC = 2 and NB = 1 and NL = 1	then					' Koyck model with 1 prom_cost and unrestricted autocorrelation	
  for !i= 5 to 236 
    smpl 5 !i !i+5 241
    koyck1_unres.ml(showopts, m=1000, c=1e-5)
    show koyck1_unres.output

   for !j = 1 to 4 
      resma(!i + !j) = 0
   next

    forecasterror(!i+1) = attendance(!i+1)  - ( cc(1) + cc(2) * prom_cost_sum(!i+1) + cc(5) * attendance(!i) - ( cc(5) +cc(8) ) * resma(!i) + cc(15) * Season11(!i+1) + cc(16) * season12(!i+1) + cc(17) * opening11(!i+1) + cc(18) * opening12(!i+1) + cc(23) * day5(!i+1) + cc(24) * day6(!i+1) + cc(27) * day7(!i+1) + cc(28) * july(!i+1) + cc(29) * august(!i+1) + cc(30) * sep(!i+1) + cc(31) * perf_wp(!i+1) + cc(32) * opp_distance(!i+1) + cc(50) * temp(!i+1) +  cc(60) * nattv(!i+1) + cc(80) * Control_entertainment(!i+1) + cc(81) * control_giveaway(!i+1) + cc(82) * control_kids(!i+1) - cc(5) * ( cc(15) * Season11(!i) + cc(16) * season12(!i) + cc(17) * opening11(!i) + cc(18) * opening12(!i) + cc(23) * day5(!i) + cc(24) * day6(!i) + cc(27) * day7(!i) + cc(28) * july(!i) + cc(29) * august(!i) + cc(30) * sep(!i) + cc(31) * perf_wp(!i) + cc(32) * opp_distance(!i) + cc(50) * temp(!i) + cc(60) * nattv(!i) + cc(80) * Control_entertainment(!i) + cc(81) * control_giveaway(!i) + cc(82) * control_kids(!i) ) )

  next
endif

if AC = 3 and NB = 1 and NL = 1	then					' Koyck model with 1 prom_cost and restricted autocorrelation	
  for !i= 5 to 236 
    smpl 5 !i !i+5 241
    koyck1_res.ml(showopts, m=1000, c=1e-5)
    show koyck1_res.output

   for !j = 1 to 4 
      resma(!i + !j) = 0
   next

    forecasterror(!i+1) = attendance(!i+1)  - ( cc(1) + cc(2) * prom_cost_sum(!i+1) + cc(5) * attendance(!i) - cc(5) * resma(!i) + cc(15) * Season11(!i+1) + cc(16) * season12(!i+1) + cc(17) * opening11(!i+1) + cc(18) * opening12(!i+1) + cc(23) * day5(!i+1) + cc(24) * day6(!i+1) + cc(27) * day7(!i+1) + cc(28) * july(!i+1) + cc(29) * august(!i+1) + cc(30) * sep(!i+1) + cc(31) * perf_wp(!i+1) + cc(32) * opp_distance(!i+1) + cc(50) * temp(!i+1) +  cc(60) * nattv(!i+1) + cc(80) * Control_entertainment(!i+1) + cc(81) * control_giveaway(!i+1) + cc(82) * control_kids(!i+1) - cc(5) * ( cc(15) * Season11(!i) + cc(16) * season12(!i) + cc(17) * opening11(!i) + cc(18) * opening12(!i) + cc(23) * day5(!i) + cc(24) * day6(!i) + cc(27) * day7(!i) + cc(28) * july(!i) + cc(29) * august(!i) + cc(30) * sep(!i) + cc(31) * perf_wp(!i) + cc(32) * opp_distance(!i) + cc(50) * temp(!i) + cc(60) * nattv(!i) + cc(80) * Control_entertainment(!i) + cc(81) * control_giveaway(!i) + cc(82) * control_kids(!i) ) )

  next
endif

' Koyck model with three prom_costs, same carryover
if AC = 1 and NB = 3 and NL = 1 then					' Koyck model with 3 prom_costs and no autocorrelation
  for !i= 5 to 236
    smpl 5 !i !i+5 241
    koyck3_no_m3is1.ml(showopts, m=1000, c=1e-5)
    show koyck3_no_m3is1.output

    forecasterror(!i+1) = attendance(!i+1)  - ( cc(1) + cc(2) * prom_cost_entertainment(!i+1) + cc(3) * prom_cost_giveaway(!i+1) + cc(4) * prom_cost_kids(!i+1) + cc(5) * attendance(!i) + ( cc(12) - cc(5) * cc(4) ) * prom_cost_kids(!i) + cc(15) * Season11(!i+1) + cc(16) * season12(!i+1) + cc(17) * opening11(!i+1) + cc(18) * opening12(!i+1) + cc(23) * day5(!i+1) + cc(24) * day6(!i+1) + cc(27) * day7(!i+1) + cc(28) * july(!i+1) + cc(29) * august(!i+1) + cc(30) * sep(!i+1) + cc(31) * perf_wp(!i+1) + cc(32) * opp_distance(!i+1) + cc(50) * temp(!i+1) +  cc(60) * nattv(!i+1) + cc(80) * Control_entertainment(!i+1) + cc(81) * control_giveaway(!i+1) + cc(82) * control_kids(!i+1) - cc(5) * ( cc(15) * Season11(!i) + cc(16) * season12(!i) + cc(17) * opening11(!i) + cc(18) * opening12(!i) + cc(23) * day5(!i) + cc(24) * day6(!i) + cc(27) * day7(!i) + cc(28) * july(!i) + cc(29) * august(!i) + cc(30) * sep(!i) + cc(31) * perf_wp(!i) + cc(32) * opp_distance(!i) + cc(50) * temp(!i) + cc(60) * nattv(!i) + cc(80) * Control_entertainment(!i) + cc(81) * control_giveaway(!i) + cc(82) * control_kids(!i) ) )

  next
endif

if AC = 2 and NB = 3 and NL = 1 then					' Koyck model with 3 prom_costs and unrestricted autocorrelation
  for !i= 5 to 236 
    smpl 5 !i !i+5 241
    koyck3_unres_m3is1.ml(showopts, m=1000, c=1e-5)
    show koyck3_unres_m3is1.output

   for !j = 1 to 4 
      resma(!i + !j) = 0
   next

    forecasterror(!i+1) = attendance(!i+1)  - ( cc(1) + cc(2) * prom_cost_entertainment(!i+1) + cc(3) * prom_cost_giveaway(!i+1) + cc(4) * prom_cost_kids(!i+1) + cc(5) * attendance(!i) - (cc(5) + cc(8)) * resma(!i) + ( cc(12) - cc(5) * cc(4) ) * prom_cost_kids(!i) + cc(15) * Season11(!i+1) + cc(16) * season12(!i+1) + cc(17) * opening11(!i+1) + cc(18) * opening12(!i+1) + cc(23) * day5(!i+1) + cc(24) * day6(!i+1) + cc(27) * day7(!i+1) + cc(28) * july(!i+1) + cc(29) * august(!i+1) + cc(30) * sep(!i+1) + cc(31) * perf_wp(!i+1) + cc(32) * opp_distance(!i+1) + cc(50) * temp(!i+1) +  cc(60) * nattv(!i+1) + cc(80) * Control_entertainment(!i+1) + cc(81) * control_giveaway(!i+1) + cc(82) * control_kids(!i+1) - cc(5) * ( cc(15) * Season11(!i) + cc(16) * season12(!i) + cc(17) * opening11(!i) + cc(18) * opening12(!i) + cc(23) * day5(!i) + cc(24) * day6(!i) + cc(27) * day7(!i) + cc(28) * july(!i) + cc(29) * august(!i) + cc(30) * sep(!i) + cc(31) * perf_wp(!i) + cc(32) * opp_distance(!i) + cc(50) * temp(!i) + cc(60) * nattv(!i) + cc(80) * Control_entertainment(!i) + cc(81) * control_giveaway(!i) + cc(82) * control_kids(!i) ) )


  next
endif

if AC = 3 and NB = 3 and NL = 1 then					' Koyck model with 3 prom_costs and restricted autocorrelation
  for !i= 5 to 236 
    smpl 5 !i !i+5 241
    koyck3_res_m3is1.ml(showopts, m=1000, c=1e-5)
    show koyck3_res_m3is1.output

   for !j = 1 to 4 
      resma(!i + !j) = 0
   next

    forecasterror(!i+1) = attendance(!i+1)  - ( cc(1) + cc(2) * prom_cost_entertainment(!i+1) + cc(3) * prom_cost_giveaway(!i+1) + cc(4) * prom_cost_kids(!i+1) + cc(5) * attendance(!i) - ( cc(5) ) * resma(!i) + ( cc(12) - cc(5) * cc(4) ) * prom_cost_kids(!i) + cc(15) * Season11(!i+1) + cc(16) * season12(!i+1) + cc(17) * opening11(!i+1) + cc(18) * opening12(!i+1) + cc(23) * day5(!i+1) + cc(24) * day6(!i+1) + cc(27) * day7(!i+1) + cc(28) * july(!i+1) + cc(29) * august(!i+1) + cc(30) * sep(!i+1) + cc(31) * perf_wp(!i+1) + cc(32) * opp_distance(!i+1) + cc(50) * temp(!i+1) +  cc(60) * nattv(!i+1) + cc(80) * Control_entertainment(!i+1) + cc(81) * control_giveaway(!i+1) + cc(82) * control_kids(!i+1) - cc(5) * ( cc(15) * Season11(!i) + cc(16) * season12(!i) + cc(17) * opening11(!i) + cc(18) * opening12(!i) + cc(23) * day5(!i) + cc(24) * day6(!i) + cc(27) * day7(!i) + cc(28) * july(!i) + cc(29) * august(!i) + cc(30) * sep(!i) + cc(31) * perf_wp(!i) + cc(32) * opp_distance(!i) + cc(50) * temp(!i) + cc(60) * nattv(!i) + cc(80) * Control_entertainment(!i) + cc(81) * control_giveaway(!i) + cc(82) * control_kids(!i) ) )

  next 
endif

' Koyck model with three prom_costs, different carryovers
if AC = 1 and NB = 3 and NL = 3 then					' Koyck model with 3 prom_costs and no autocorrelation
  for !i= 151 to 236 
    smpl 5 !i !i+5 241
    koyck33_no_m3is1.ml(showopts, m=1000, c=1e-5)
    show koyck33_no_m3is1.output

    forecasterror(!i+1) = attendance(!i+1)  - ( cc(1) + cc(2) * prom_cost_entertainment(!i+1) - cc(2) * (exp(cc(6)) +exp(cc(7))) * prom_cost_entertainment(!i) + cc(2) * exp(cc(6)) * exp(cc(7)) * prom_cost_entertainment(!i-1) + cc(3) * prom_cost_giveaway(!i+1) - cc(3) * (cc(5) +exp(cc(7))) * prom_cost_giveaway(!i) + cc(3) * cc(5) * exp(cc(7)) * prom_cost_giveaway(!i-1) + cc(4) * prom_cost_kids(!i+1) - cc(4) * (cc(5) +exp(cc(6))) * prom_cost_kids(!i) + cc(4) * cc(5) * exp(cc(6)) * prom_cost_kids(!i-1) +	( cc(5) + exp(cc(6)) + exp(cc(7)) ) * attendance(!i) - ( cc(5) * exp(cc(6)) + cc(5) * exp(cc(7)) + exp(cc(6)) * exp(cc(7)) ) * attendance(!i-1) + ( cc(5) * exp(cc(6)) * exp(cc(7)) ) * attendance(!i-2)  + ( cc(12) - exp(cc(7)) * cc(4) ) * prom_cost_kids(!i) - ( cc(5) + exp(cc(6)) ) * ( cc(12) - exp(cc(7)) * cc(4) ) * prom_cost_kids(!i-1) + cc(5) * exp(cc(6)) * ( cc(12) - exp(cc(7)) * cc(4) ) * prom_cost_kids(!i-2) + cc(15) * Season11(!i+1) + cc(16) * season12(!i+1) + cc(17) * opening11(!i+1) + cc(18) * opening12(!i+1) + cc(23) * day5(!i+1) + cc(24) * day6(!i+1) + cc(27) * day7(!i+1) + cc(28) * july(!i+1) + cc(29) * august(!i+1) + cc(30) * sep(!i+1) + cc(31) * perf_wp(!i+1) + cc(32) * opp_distance(!i+1) + cc(50) * temp(!i+1) +  cc(60) * nattv(!i+1) + cc(80) * Control_entertainment(!i+1) + cc(81) * control_giveaway(!i+1) + cc(82) * control_kids(!i+1) - (cc(5) + exp(cc(6)) + exp(cc(7))) * ( cc(15) * Season11(!i) + cc(16) * season12(!i) + cc(17) * opening11(!i) + cc(18) * opening12(!i) + cc(23) * day5(!i) + cc(24) * day6(!i) + cc(27) * day7(!i) + cc(28) * july(!i) + cc(29) * august(!i) + cc(30) * sep(!i) + cc(31) * perf_wp(!i) + cc(32) * opp_distance(!i) + cc(50) * temp(!i) + cc(60) * nattv(!i) + cc(80) * Control_entertainment(!i) + cc(81) * control_giveaway(!i) + cc(82) * control_kids(!i) ) + ( cc(5) * exp(cc(6)) + cc(5) * exp(cc(7)) + exp(cc(6)) * exp(cc(7)) ) * ( cc(15) * Season11(!i-1) + cc(16) * season12(!i-1) + cc(17) * opening11(!i-1) + cc(18) * opening12(!i-1) + cc(23) * day5(!i-1) + cc(24) * day6(!i-1) + cc(27) * day7(!i-1) + cc(28) * july(!i-1) + cc(29) * august(!i-1) + cc(30) * sep(!i-1) + cc(31) * perf_wp(!i-1) + cc(32) * opp_distance(!i-1) + cc(50) * temp(!i-1) + cc(60) * nattv(!i-1) + cc(80) * Control_entertainment(!i-1) + cc(81) * control_giveaway(!i-1) + cc(82) * control_kids(!i-1) ) - ( cc(5) * exp(cc(6)) * exp(cc(7)) ) * ( cc(15) * Season11(!i-2) + cc(16) * season12(!i-2) + cc(17) * opening11(!i-2) + cc(18) * opening12(!i-2) + cc(23) * day5(!i-2) + cc(24) * day6(!i-2) + cc(27) * day7(!i-2) + cc(28) * july(!i-2) + cc(29) * august(!i-2) + cc(30) * sep(!i-2) + cc(31) * perf_wp(!i-2)  + cc(32) * opp_distance(!i-2) + cc(50) * temp(!i-2) +  cc(60) * nattv(!i-2) + cc(80) * Control_entertainment(!i-2) + cc(81) * control_giveaway(!i-2) + cc(82) * control_kids(!i-2) ) )
  next 
endif


if AC = 2 and NB = 3 and NL = 3 then					' Koyck model with 3 prom_costs and unrestricted autocorrelation
  for !i= 5 to 236 
    smpl 5 !i !i+5 241
    koyck33_unres_m3is1.ml(showopts, m=1000, c=1e-9)
    show koyck33_unres_m3is1.output

   for !j = 1 to 4 
      resma(!i + !j) = 0
   next

    forecasterror(!i+1) = attendance(!i+1)  - ( cc(1) + cc(2) * prom_cost_entertainment(!i+1) - cc(2) * (cc(6) +cc(7)) * prom_cost_entertainment(!i) + cc(2) * cc(6) * cc(7) * prom_cost_entertainment(!i-1) + cc(3) * prom_cost_giveaway(!i+1) - cc(3) * (cc(5) +cc(7)) * prom_cost_giveaway(!i) + cc(3) * cc(5) * cc(7) * prom_cost_giveaway(!i-1) + cc(4) * prom_cost_kids(!i+1) - cc(4) * (cc(5) +cc(6)) * prom_cost_kids(!i) + cc(4) * cc(5) * cc(6) * prom_cost_kids(!i-1) +	( cc(5) + cc(6) + cc(7) ) * attendance(!i) - ( cc(5) * cc(6) + cc(5) * cc(7) + cc(6) * cc(7) ) * attendance(!i-1) + ( cc(5) * cc(6) * cc(7) ) * attendance(!i-2)  + ( cc(12) - cc(7) * cc(4) ) * prom_cost_kids(!i) - ( cc(5) + cc(6) ) * ( cc(12) - cc(7) * cc(4) ) * prom_cost_kids(!i-1) + cc(5) * cc(6) * ( cc(12) - cc(7) * cc(4) ) * prom_cost_kids(!i-2) - ( ( cc(5) + cc(6) + cc(7) + cc(8) ) * resma(!i) - ( cc(5) * cc(6) + cc(5) * cc(7) + cc(6) * cc(7) + cc(9) ) * resma(!i-1) + ( cc(5) * cc(6) * cc(7) + cc(10) ) * resma(!i-2) )  + cc(15) * Season11(!i+1) + cc(16) * season12(!i+1) + cc(17) * opening11(!i+1) + cc(18) * opening12(!i+1) + cc(23) * day5(!i+1) + cc(24) * day6(!i+1) + cc(27) * day7(!i+1) + cc(28) * july(!i+1) + cc(29) * august(!i+1) + cc(30) * sep(!i+1) + cc(31) * perf_wp(!i+1) + cc(32) * opp_distance(!i+1) + cc(50) * temp(!i+1) +  cc(60) * nattv(!i+1) + cc(80) * Control_entertainment(!i+1) + cc(81) * control_giveaway(!i+1) + cc(82) * control_kids(!i+1) - (cc(5) + cc(6) + cc(7)) * ( cc(15) * Season11(!i) + cc(16) * season12(!i) + cc(17) * opening11(!i) + cc(18) * opening12(!i) + cc(23) * day5(!i) + cc(24) * day6(!i) + cc(27) * day7(!i) + cc(28) * july(!i) + cc(29) * august(!i) + cc(30) * sep(!i) + cc(31) * perf_wp(!i) + cc(32) * opp_distance(!i) + cc(50) * temp(!i) + cc(60) * nattv(!i) + cc(80) * Control_entertainment(!i) + cc(81) * control_giveaway(!i) + cc(82) * control_kids(!i) ) + ( cc(5) * cc(6) + cc(5) * cc(7) + cc(6) * cc(7) ) * ( cc(15) * Season11(!i-1) + cc(16) * season12(!i-1) + cc(17) * opening11(!i-1) + cc(18) * opening12(!i-1) + cc(23) * day5(!i-1) + cc(24) * day6(!i-1) + cc(27) * day7(!i-1) + cc(28) * july(!i-1) + cc(29) * august(!i-1) + cc(30) * sep(!i-1) + cc(31) * perf_wp(!i-1) + cc(32) * opp_distance(!i-1) + cc(50) * temp(!i-1) + cc(60) * nattv(!i-1) + cc(80) * Control_entertainment(!i-1) + cc(81) * control_giveaway(!i-1) + cc(82) * control_kids(!i-1) ) - ( cc(5) * cc(6) * cc(7) ) * ( cc(15) * Season11(!i-2) + cc(16) * season12(!i-2) + cc(17) * opening11(!i-2) + cc(18) * opening12(!i-2) + cc(23) * day5(!i-2) + cc(24) * day6(!i-2) + cc(27) * day7(!i-2) + cc(28) * july(!i-2) + cc(29) * august(!i-2) + cc(30) * sep(!i-2) + cc(31) * perf_wp(!i-2)  + cc(32) * opp_distance(!i-2) + cc(50) * temp(!i-2) +  cc(60) * nattv(!i-2) + cc(80) * Control_entertainment(!i-2) + cc(81) * control_giveaway(!i-2) + cc(82) * control_kids(!i-2) ) )

   next 
endif

if AC = 3 and NB = 3 and NL = 3 then					' Koyck model with 3 prom_costs and restricted autocorrelation
  for !i= 5 to 236 
    smpl 5 !i !i+5 241
    koyck33_res_m3is1.ml(showopts, m=1000, c=1e-5)
    show koyck33_res_m3is1.output

   for !j = 1 to 4 
      resma(!i + !j) = 0
   next

     forecasterror(!i+1) = attendance(!i+1)  - ( cc(1) + cc(2) * prom_cost_entertainment(!i+1) - cc(2) * (cc(6) +cc(7)) * prom_cost_entertainment(!i) + cc(2) * cc(6) * cc(7) * prom_cost_entertainment(!i-1) + cc(3) * prom_cost_giveaway(!i+1) - cc(3) * (cc(5) +cc(7)) * prom_cost_giveaway(!i) + cc(3) * cc(5) * cc(7) * prom_cost_giveaway(!i-1) + cc(4) * prom_cost_kids(!i+1) - cc(4) * (cc(5) +cc(6)) * prom_cost_kids(!i) + cc(4) * cc(5) * cc(6) * prom_cost_kids(!i-1) +	( cc(5) + cc(6) + cc(7) ) * attendance(!i) - ( cc(5) * cc(6) + cc(5) * cc(7) + cc(6) * cc(7) ) * attendance(!i-1) + ( cc(5) * cc(6) * cc(7) ) * attendance(!i-2)  + ( cc(12) - cc(7) * cc(4) ) * prom_cost_kids(!i) - ( cc(5) + cc(6) ) * ( cc(12) - cc(7) * cc(4) ) * prom_cost_kids(!i-1) + cc(5) * cc(6) * ( cc(12) - cc(7) * cc(4) ) * prom_cost_kids(!i-2) - ( ( cc(5) + cc(6) + cc(7) ) * resma(!i) - ( cc(5) * cc(6) + cc(5) * cc(7) + cc(6) * cc(7) ) * resma(!i-1) + ( cc(5) * cc(6) * cc(7) ) * resma(!i-2) )  + cc(15) * Season11(!i+1) + cc(16) * season12(!i+1) + cc(17) * opening11(!i+1) + cc(18) * opening12(!i+1) + cc(23) * day5(!i+1) + cc(24) * day6(!i+1) + cc(27) * day7(!i+1) + cc(28) * july(!i+1) + cc(29) * august(!i+1) + cc(30) * sep(!i+1) + cc(31) * perf_wp(!i+1) + cc(32) * opp_distance(!i+1) + cc(50) * temp(!i+1) +  cc(60) * nattv(!i+1) + cc(80) * Control_entertainment(!i+1) + cc(81) * control_giveaway(!i+1) + cc(82) * control_kids(!i+1) - (cc(5) + cc(6) + cc(7)) * ( cc(15) * Season11(!i) + cc(16) * season12(!i) + cc(17) * opening11(!i) + cc(18) * opening12(!i) + cc(23) * day5(!i) + cc(24) * day6(!i) + cc(27) * day7(!i) + cc(28) * july(!i) + cc(29) * august(!i) + cc(30) * sep(!i) + cc(31) * perf_wp(!i) + cc(32) * opp_distance(!i) + cc(50) * temp(!i) + cc(60) * nattv(!i) + cc(80) * Control_entertainment(!i) + cc(81) * control_giveaway(!i) + cc(82) * control_kids(!i) ) + ( cc(5) * cc(6) + cc(5) * cc(7) + cc(6) * cc(7) ) * ( cc(15) * Season11(!i-1) + cc(16) * season12(!i-1) + cc(17) * opening11(!i-1) + cc(18) * opening12(!i-1) + cc(23) * day5(!i-1) + cc(24) * day6(!i-1) + cc(27) * day7(!i-1) + cc(28) * july(!i-1) + cc(29) * august(!i-1) + cc(30) * sep(!i-1) + cc(31) * perf_wp(!i-1) + cc(32) * opp_distance(!i-1) + cc(50) * temp(!i-1) + cc(60) * nattv(!i-1) + cc(80) * Control_entertainment(!i-1) + cc(81) * control_giveaway(!i-1) + cc(82) * control_kids(!i-1) ) - ( cc(5) * cc(6) * cc(7) ) * ( cc(15) * Season11(!i-2) + cc(16) * season12(!i-2) + cc(17) * opening11(!i-2) + cc(18) * opening12(!i-2) + cc(23) * day5(!i-2) + cc(24) * day6(!i-2) + cc(27) * day7(!i-2) + cc(28) * july(!i-2) + cc(29) * august(!i-2) + cc(30) * sep(!i-2) + cc(31) * perf_wp(!i-2)  + cc(32) * opp_distance(!i-2) + cc(50) * temp(!i-2) +  cc(60) * nattv(!i-2) + cc(80) * Control_entertainment(!i-2) + cc(81) * control_giveaway(!i-2) + cc(82) * control_kids(!i-2) ) )

  next 
endif

