Premium monitor added.

Author: Juan Sabuco

insurancesimulation.py

@@ -138,6 +138,8 @@ def __init__(self, override_no_riskmodels, replic_ID, simulation_parameters):
         self.history_total_reincash = []
         self.history_total_reincontracts = []
         self.history_total_reinoperational = []
+
+        self.history_market_premium = []
 
 
     def build_agents(self, agent_class, agent_class_string, parameters, agent_parameters):
@@ -241,6 +243,7 @@ def save_data(self):
         self.history_total_reincash.append(total_reincash_no)
         self.history_total_reincontracts.append(total_reincontracts_no)
         self.history_total_reinoperational.append(reinoperational_no)
+        self.history_market_premium.append(self.market_premium)
 
         individual_contracts_no = [len(insurancefirm.underwritten_contracts) for insurancefirm in self.insurancefirms]
         for i in range(len(individual_contracts_no)):
@@ -433,6 +436,8 @@ def replication_log_prepare(self):
         to_log.append(("data/two_reinoperational.dat", self.history_total_reinoperational, "a"))
         to_log.append(("data/two_reincontracts.dat", self.history_total_reincontracts, "a"))
         to_log.append(("data/two_reincash.dat", self.history_total_reincash, "a"))
+        to_log.append(("data/two_premium.dat", self.history_market_premium, "a"))
+
         return to_log
 
     def replication_log_prepare_oneriskmodel(self):
@@ -443,6 +448,8 @@ def replication_log_prepare_oneriskmodel(self):
         to_log.append(("data/one_reinoperational.dat", self.history_total_reinoperational, "a"))
         to_log.append(("data/one_reincontracts.dat", self.history_total_reincontracts, "a"))
         to_log.append(("data/one_reincash.dat", self.history_total_reincash, "a"))
+        to_log.append(("data/one_premium.dat", self.history_market_premium, "a"))
+
         return to_log
 
     def single_log_prepare(self):
@@ -453,6 +460,8 @@ def single_log_prepare(self):
         to_log.append(("data/reinoperational.dat", self.history_total_reinoperational, "w"))
         to_log.append(("data/reincontracts.dat", self.history_total_reincontracts, "w"))
         to_log.append(("data/reincash.dat", self.history_total_reincash, "w"))
+        to_log.append(("data/premium.dat", self.history_market_premium, "w"))
+
         return to_log
 
 #if __name__ == "__main__":

plotter.py

@@ -25,18 +25,26 @@
 reincash = [eval(k) for k in rfile]
 rfile.close()
 
+rfile = open("data/premium.dat","r")
+premium = [eval(k) for k in rfile]
+rfile.close()
+
 c_s = []
 
 o_s = []
 
 h_s = []
 
+p_s = []
+
 c_re = []
 
 o_re= []
 
 h_re = []
 
+p_e = []
+
 for i in range(len(contracts[0])):
     cs = np.mean([item[i] for item in contracts])
     os = np.median([item[i] for item in op])
@@ -53,24 +61,30 @@
     o_re.append(ore)
     h_re.append(hre)
 
+    p_s = np.median([item[i] for item in premium])
+    p_e.append(p_s)
+
 fig = plt.figure()
-ax0 = fig.add_subplot(611)
+ax0 = fig.add_subplot(711)
 ax0.plot(range(len(c_s)), c_s,"b")
 ax0.set_ylabel("Contracts")
-ax1 = fig.add_subplot(612)
+ax1 = fig.add_subplot(712)
 ax1.plot(range(len(o_s)), o_s,"b")
-ax1.set_ylabel("Active firms (out of initially 20)")
-ax0 = fig.add_subplot(613)
+ax1.set_ylabel("Active firms")
+ax0 = fig.add_subplot(713)
 ax0.plot(range(len(h_s)), h_s,"b")
 ax0.set_ylabel("Cash")
-ax0 = fig.add_subplot(614)
+ax0 = fig.add_subplot(714)
 ax0.plot(range(len(c_re)), c_re,"r")
 ax0.set_ylabel("Contracts")
-ax1 = fig.add_subplot(615)
+ax1 = fig.add_subplot(715)
 ax1.plot(range(len(o_re)), o_re,"r")
-ax1.set_ylabel("Active reinfirms (out of initially 1)")
-ax0 = fig.add_subplot(616)
+ax1.set_ylabel("Active reinfirms")
+ax0 = fig.add_subplot(716)
 ax0.plot(range(len(h_re)), h_re,"r")
 ax0.set_ylabel("Cash")
+ax0 = fig.add_subplot(717)
+ax0.plot(range(len(p_e)), p_e,"k")
+ax0.set_ylabel("Premium")
 ax1.set_xlabel("Time")
 plt.show()