diff --git a/showgraphs.m b/showgraphs.m
index 593de9225e4d6c401516efe9871c65edc9e4a200..79fbfec0387c8a75edf66b7e06e28162b73742a2 100644
--- a/showgraphs.m
+++ b/showgraphs.m
@@ -194,12 +194,21 @@ ErrorListPlot[mixingTimesBars,Joined->True,PlotMarkers->Automatic,AxesLabel->{"t
 
 
 (* ::Subsection:: *)
-(*Plot average #triangles vs n*)
+(*Triangle exponent: Plot average #triangles vs n*)
 
 
+(* When importing from exponent-only-data file *)
+gsraw=Import[NotebookDirectory[]<>"data/graphdata_partial.m"];
+gsraw=SortBy[gsraw,#[[1,1]]&]; (* Sort by n *)
+averagesGrouped=GatherBy[gsraw,{#[[1,2]]&,#[[1,1]]&}];
+
+
+(* When importing from general file *)
 averages=Map[{#[[1]],Mean[#[[2,1;;-1]]]}&,gsraw];
 (* averages=SortBy[averages,#[[1,1]]&]; (* Sort by n *) *)
 averagesGrouped=GatherBy[averages,{#[[1,2]]&,#[[1,1]]&}]; (* Split by n,tau *)
+
+
 (* 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]]];
@@ -221,9 +230,16 @@ ListLogLogPlot[averagesErrorBars[[All,All,1]],Joined->True,PlotMarkers->Automati
 
 loglogdata=Log[averagesErrorBars[[All,All,1]]];
 fits=Map[Fit[#,{1,logn},logn]&,loglogdata];
+fitsExtra=Map[LinearModelFit[#,logn,logn]&,loglogdata];
 
 
-Show[ListLogLogPlot[averagesErrorBars[[All,All,1]],PlotMarkers->Automatic,AxesLabel->{"n","\[LeftAngleBracket]triangles\[RightAngleBracket]"},PlotLegends->taulabels],Plot[fits,{logn,1,2000}]]
+fitsExtra[[1]]["ParameterConfidenceIntervalTable"]
+fitsExtra[[1]]["BestFitParameters"]
+fitsExtra[[1]]["ParameterErrors"]
+fitsExtra[[1]]["ParameterConfidenceIntervals"]
+
+
+Show[ListLogLogPlot[averagesErrorBars[[All,All,1]],Joined->True,PlotMarkers->Automatic,AxesLabel->{"n","\[LeftAngleBracket]triangles\[RightAngleBracket]"},PlotLegends->taulabels],Plot[fits,{logn,1,2000}]]
 
 
 tauValues=averagesGrouped[[All,1,1,1,2]];
@@ -231,6 +247,12 @@ exponents=Transpose[{tauValues,fits[[All,2,1]]}];
 Show[ListPlot[exponents,Joined->True,PlotMarkers->Automatic,AxesLabel->{"tau","triangle-law-exponent"},PlotRange->{{2,3},{0,1.6}}],Plot[3/2(3-tau),{tau,2,3}]]
 
 
+tauValues=averagesGrouped[[All,1,1,1,2]];
+exponentsErrorBars=Map[{{#[[1]],#[[2]]["BestFitParameters"][[2]]},ErrorBar[#[[2]]["ParameterConfidenceIntervals"][[2]]-#[[2]]["BestFitParameters"][[2]]]}&,
+Transpose[{tauValues,fitsExtra}]];
+Show[ErrorListPlot[exponentsErrorBars,Joined->True,PlotMarkers->Automatic,AxesLabel->{"tau","triangle-law-exponent"},PlotRange->{{2,3},{0,1.6}}],Plot[3/2(3-tau),{tau,2,3}]]
+
+
 (* ::Subsection:: *)
 (*Plot #triangles distribution for specific (n,tau)*)