para.set_property_rights(scenario='OA')
para.solve_para()
home = '/home/nealbob'
folder = '/Dropbox/Model/results/chapter5/'
scenarios = ['CS', 'SWA', 'OA', 'NS', 'CS-SL', 'SWA-SL', 'CS-SWA']
results = {scen: 0 for scen in scenarios}
policies = {scen: 0 for scen in scenarios}
for i in range(1):
#try:
para.central_case(N = 100)
para.t_cost = 100000000000
para.aproximate_shares(nonoise=True)
#if i > 0:
# para.randomize(N = 100)
# para.aproximate_shares()
for scen in scenarios:
para.set_property_rights(scenario=scen)
res = {'paras' : [], 'stats' : [], 'VE': [], 'PE' : []}
pol = []
mod = model.Model(para)
VE, PE, stats, policy = mod.chapter5()
res['stats'].append(stats)
res['paras'].append(para.para_list)
res['VE'].append(VE)
out = '/Dropbox/Thesis/IMG/chapter6/'
para = Para()
scenarios = ['CS'] #['RS-HL-O', 'RS-HL', 'RS-O', 'RS', 'CS', 'CS-O', 'CS-HL', 'CS-HL-O', 'CS-U']
results = {scen: 0 for scen in scenarios}
Lambda = {scen: 0 for scen in scenarios}
LambdaK = {scen: 0 for scen in scenarios}
#==========================================
# Central case (with trade)
#==========================================
for scen in scenarios:
para.set_property_rights(scenario=scen)
para.aproximate_shares(nonoise=True)
mod = Model(para)
results[scen], Lambda[scen], LambdaK[scen] = mod.chapter6(sens=True)
del mod
#chapter6.tables(results, scenarios, Lambda, LambdaK, label='central')
#with open(home + folder + 'central_result.pkl', 'wb') as f:
# pickle.dump(results, f)
# f.close()
"""
