diff --git a/triangle_exponent_plots.m b/triangle_exponent_plots.m
index 23440533c6d45c265afb060731a1374913c63a04..6c3db7cc1c48c7859cf557ae0ec1e19386e076fc 100644
--- a/triangle_exponent_plots.m
+++ b/triangle_exponent_plots.m
@@ -20,11 +20,17 @@ Needs["ErrorBarPlots`"]
 (* graphdata_exponent_highN.m *)
 (* graphdata_properties2.m *)
 (* graphdata_canonical_properties.m *)
-gsraw=Import[NotebookDirectory[]<>"data/graphdata_exponent_mix32.m"];
+gsraw=Import[NotebookDirectory[]<>"data/graphdata_exponent_hightau.m"];
 gsraw=SortBy[gsraw,#[[1,1]]&]; (* Sort by n *)
 averagesGrouped=GatherBy[gsraw,{#[[1,2]]&,#[[1,1]]&}];
 
 
+gsraw2=Import[NotebookDirectory[]<>"data/graphdata_canonical_properties2.m"];
+gsraw2=SortBy[gsraw2,#[[1,1]]&]; (* Sort by n *)
+averagesGrouped2=GatherBy[gsraw2,{#[[1,2]]&,#[[1,1]]&}];
+canonicalDatapoints=Map[{#[[1,1,1]],Mean[#[[All,2]]]}&,averagesGrouped2,{2}];
+
+
 (* averagesGrouped[[ tau index, n index, run index , {ntau, avgtri} ]] *)
 nlabels=Map["n = "<>ToString[#]&,averagesGrouped[[1,All,1,1,1]]];
 taulabels=Map["tau = "<>ToString[#]&,averagesGrouped[[All,1,1,1,2]]];
@@ -61,6 +67,10 @@ mediansFits=Map[Fit[#,{1,logn},logn]&,mediansLoglogdata];
 mediansFitsExtra=Map[LinearModelFit[#,logn,logn]&,mediansLoglogdata];
 
 
+canonicalloglogdata=Log[canonicalDatapoints[[All,nRange]]];
+canonicalFitsExtra=Map[LinearModelFit[#,logn,logn]&,canonicalloglogdata];
+
+
 averagesFitsExtra[[1]]["ParameterConfidenceIntervalTable"]
 averagesFitsExtra[[1]]["BestFitParameters"]
 averagesFitsExtra[[1]]["ParameterErrors"]
@@ -91,10 +101,18 @@ Show[ListPlot[{averagesExponents,mediansExponents},Joined->True,PlotMarkers->Aut
 (* For visual, shift the tau values slightly left or right to distinguish the two datasets *)
 tauValues=averagesGrouped[[All,1,1,1,2]];
 averagesExponentsErrorBars=Map[{{#[[1]],#[[2]]["BestFitParameters"][[2]]},ErrorBar[#[[2]]["ParameterConfidenceIntervals"][[2]]-#[[2]]["BestFitParameters"][[2]]]}&,
-Transpose[{tauValues-0.001,averagesFitsExtra}]];
+Transpose[{tauValues-0.003,averagesFitsExtra}]];
 mediansExponentsErrorBars=Map[{{#[[1]],#[[2]]["BestFitParameters"][[2]]},ErrorBar[#[[2]]["ParameterConfidenceIntervals"][[2]]-#[[2]]["BestFitParameters"][[2]]]}&,
-Transpose[{tauValues+0.001,mediansFitsExtra}]];
-plot2=Show[ErrorListPlot[{averagesExponentsErrorBars,mediansExponentsErrorBars},Joined->True,PlotMarkers->Automatic,Frame->True,FrameLabel->{"tau","triangle exponent"},PlotRange->{{2,3},{0,1.6}},ImageSize->300],Plot[3/2(3-tau),{tau,2,3},PlotStyle->{Dashed}]]
+Transpose[{tauValues+0.003,mediansFitsExtra}]];
+canonicalExponentsErrorBars=Map[{{#[[1]],#[[2]]["BestFitParameters"][[2]]},ErrorBar[#[[2]]["ParameterConfidenceIntervals"][[2]]-#[[2]]["BestFitParameters"][[2]]]}&,
+Transpose[{tauValues+0.000,canonicalFitsExtra}]];
+plot2=Show[
+ErrorListPlot[{averagesExponentsErrorBars,mediansExponentsErrorBars,canonicalExponentsErrorBars},
+Joined->True,PlotMarkers->Automatic,
+Frame->True,FrameLabel->{"tau","triangle exponent"},
+PlotRange->{{2,3},{0,1.6}},
+ImageSize->300],
+Plot[3/2(3-tau),{tau,2,3},PlotStyle->{Black,Dashed}]]
 
 
 Export[NotebookDirectory[]<>"plots/triangle_exponent.pdf",plot2]