update results by starting from previous results

code/estimark/estimation.py

@@ -22,6 +22,8 @@
 )
 from time import time  # Timing utility
 
+from pathlib import Path
+
 import matplotlib.pyplot as plt
 import numpy as np  # Numeric Python
 
@@ -562,7 +564,18 @@ def estimate(
 
     targeted_moments = get_targeted_moments()
 
-    initial_guess = [parameters.DiscFacAdj_start, parameters.CRRA_start]
+    # start from previous estimation results if available
+    try:
+        csv_file_path = (
+            Path(__file__).resolve().parent
+            / ".."
+            / "tables"
+            / (estimation_agent + "_estimate_results.csv")
+        )
+
+        initial_guess = np.genfromtxt(csv_file_path, skip_header=1, delimiter=",")
+    except:
+        initial_guess = [parameters.DiscFacAdj_start, parameters.CRRA_start]
 
     # Estimate the model using Nelder-Mead
     if estimate_model:

code/tables/IndShock_estimate_results.csv

@@ -1,2 +1,2 @@
 DiscFacAdj,CRRA
-0.8783572029816282,3.5160942263466355
+0.8689011612165243,3.796066262763537

code/tables/Portfolio_estimate_results.csv

@@ -1,2 +1,2 @@
-DiscFacAdj,CRRA
-0.866131986228827,3.7557510682554494
+DiscFacAdj,CRRA
+0.8193095309029956,4.5574205165556725

code/tables/WarmGlowPortfolio_estimate_results.csv

@@ -1,2 +1,2 @@
-DiscFacAdj,CRRA
-0.8642073648151234,3.8055353724199676
+DiscFacAdj,CRRA
+0.8246562094500894,4.4642643726332185

code/tables/WarmGlow_estimate_results.csv

@@ -1,2 +1,2 @@
-DiscFacAdj,CRRA
-0.8764224065671837,3.505877643073127
+DiscFacAdj,CRRA
+0.8755702129442734,3.6158511660969577

code/tables/WealthPortfolio_estimate_results.csv

@@ -1,2 +1,2 @@
-DiscFacAdj,CRRA
-0.6013593197374101,5.031896157976007
+DiscFacAdj,CRRA
+0.6326006490726698,5.943630057453893