diff --git a/triangle_exponent_plots.m b/triangle_exponent_plots.m index af61d1223c4fb6a47a3c62c48e505bd2fd7567fe..b88398c79d7059660338a64c3be02f5d1493305d 100644 --- a/triangle_exponent_plots.m +++ b/triangle_exponent_plots.m @@ -16,7 +16,7 @@ Needs["ErrorBarPlots`"] (* When importing from exponent-only-data file *) -gsraw=Import[NotebookDirectory[]<>"data/graphdata_exponent.m"]; +gsraw=Import[NotebookDirectory[]<>"data/graphdata_exponent_mix32.m"]; gsraw=SortBy[gsraw,#[[1,1]]&]; (* Sort by n *) averagesGrouped=GatherBy[gsraw,{#[[1,2]]&,#[[1,1]]&}]; @@ -51,7 +51,10 @@ 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}]] +plot1=Show[ListLogLogPlot[averagesErrorBars[[All,All,1]],Joined->False,PlotMarkers->Automatic,AxesLabel->{"n","\[LeftAngleBracket]triangles\[RightAngleBracket]"},PlotLegends->taulabels],Plot[fits,{logn,1,2000}]] + + +Export[NotebookDirectory[]<>"plots/avgtris_n.pdf",plot1] (* ::Subsection:: *) @@ -60,7 +63,7 @@ Show[ListLogLogPlot[averagesErrorBars[[All,All,1]],Joined->True,PlotMarkers->Aut tauValues=averagesGrouped[[All,1,1,1,2]]; exponents=Transpose[{tauValues,fits[[All,2,1]]}]; -Show[ListPlot[exponents,Joined->True,PlotMarkers->Automatic,AxesLabel->{"tau","T(\[Tau])"},PlotRange->{{2,3},{0,1.6}}],Plot[3/2(3-tau),{tau,2,3}]] +Show[ListPlot[exponents,Joined->True,PlotMarkers->Automatic,AxesLabel->{"tau","exponent"},PlotRange->{{2,3},{0,1.6}}],Plot[3/2(3-tau),{tau,2,3}]] (* ::Subsection:: *) @@ -70,4 +73,7 @@ Show[ListPlot[exponents,Joined->True,PlotMarkers->Automatic,AxesLabel->{"tau","T 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","T(\[Tau])"},PlotRange->{{2,3},{0,1.6}}],Plot[3/2(3-tau),{tau,2,3}]] +plot2=Show[ErrorListPlot[exponentsErrorBars,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}]] + + +Export[NotebookDirectory[]<>"plots/triangle_exponent.pdf",plot2]