gsraw2=Import[NotebookDirectory[]<>"data/graphdata_exponent_hightau.m"];

gsraw2=SortBy[gsraw2,#[[1,1]]&]; (* Sort by n *)

averagesGrouped=GatherBy[gsraw2,{#[[1,2]]&,#[[1,1]]&}];

averagesLoglogdata=Map[{Log[#[[1,1,1]]],Log[Mean[#[[All,2]]]]}&,averagesGrouped[[All,nRange]],{2}];

averagesFitsExtra=Map[LinearModelFit[#,logn,logn]&,averagesLoglogdata];

avgTauValues=averagesGrouped[[All,1,1,1,2]];

averagesExponentsErrorBars=Map[{{#[[1]],#[[2]]["BestFitParameters"][[2]]},ErrorBar[#[[2]]["ParameterConfidenceIntervals"][[2]]-#[[2]]["BestFitParameters"][[2]]]}&,

Transpose[{avgTauValues-0.000,averagesFitsExtra}]];

tauValues=gdata[[All,1,1,2]];

(* For visual, shift the tau values slightly left or right to distinguish the two datasets *)

uniformExponents=Map[{{#[[1]],#[[2]]["BestFitParameters"][[2]]},ErrorBar[#[[2]]["ParameterConfidenceIntervals"][[2]]-#[[2]]["BestFitParameters"][[2]]]}&, Transpose[{tauValues+0.000,uniformFits}]];

ECMExponents    =Map[{{#[[1]],#[[2]]["BestFitParameters"][[2]]},ErrorBar[#[[2]]["ParameterConfidenceIntervals"][[2]]-#[[2]]["BestFitParameters"][[2]]]}&, Transpose[{tauValues+0.000,ECMFits}]];

Needs["MaTeX`"]

plot2=Show[

ErrorListPlot[{ECMExponents,uniformExponents,averagesExponentsErrorBars},

Joined->True,PlotMarkers->Automatic,

PlotLegends->Placed[{"ECM","canonical","average"},{Left,Bottom}],

PlotRange->{{2,3},{0,1.6}},

ImageSize->300],

Plot[3/2(3-tau),{tau,2,3},PlotStyle->{Black,Dashed},PlotLegends->Placed[LineLegend[{MaTeX["\\frac{3}{2}(3-\\tau)"]},LegendMarkerSize->20],{Left,Bottom}]]]

Export[NotebookDirectory[]<>"plots/triangle_exponent.pdf",plot2]