diff --git a/analyze_sayma_data.py b/analyze_sayma_data.py
index f00d753..95cb8e5 100644
--- a/analyze_sayma_data.py
+++ b/analyze_sayma_data.py
@@ -53,6 +53,7 @@ def report(*args, **kwargs):
     rpt.append("  Minimum:                 {:>10.4f} ps".format(prm["phase_ps_min"]))
     rpt.append("  Maximum:                 {:>10.4f} ps".format(prm["phase_ps_max"]))
     rpt.append("  Standard Deviation:      {:>10.4f} ps".format(prm["phase_ps_std"]))
+    rpt.append("  Max Absolute Deviation:  {:>10.4f} ps".format(prm["phase_ps_maxabsdev"]))
     rpt.append("  Mean Absolute Deviation: {:>10.4f} ps".format(prm["phase_ps_meanabsdev"]))
     rpt.append("")
     # TODO: Use jinja2 to produce a report
@@ -110,6 +111,8 @@ def main():
         rpt_params["phase_ps_min"] = data_phase_ps.min()
         rpt_params["phase_ps_max"] = data_phase_ps.max()
         rpt_params["phase_ps_std"] = data_phase_ps.std()
+        rpt_params["phase_ps_maxabsdev"] = \
+            np.absolute(data_phase_ps - data_phase_ps.max()).mean()
         rpt_params["phase_ps_meanabsdev"] = \
             np.absolute(data_phase_ps - data_phase_ps.mean()).mean()
         # Generate the report