AENC/switchchain Changeset - d0cb8a9f1b0e · Centrum Wiskunde & Informatica (CWI)

Changeset - d0cb8a9f1b0e

Parent rev.

Child rev.

[Not reviewed]

0 1 0

Tom Bannink - 8 years ago 2017-05-08 14:21:48
tombannink@gmail.com

Add error bars to triangle exponent plot

1 file changed with 24 insertions and 2 deletions:

showgraphs.m

0 comments (0 inline, 0 general)

showgraphs.m

➞

Show inline comments

 (* ::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]]];
 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)*)

0 comments (0 inline, 0 general)