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

@@ @@ -191,18 +191,27 @@ mixingTimesBars=Map[ @@
 Show[
 ErrorListPlot[mixingTimesBars,Joined->True,PlotMarkers->Automatic,AxesLabel->{"tau","~~mixing time divided by n, for n = 1000"},PlotRange->{{2,3},{0,30}}]
 ,Plot[(32-26(tau-2)),{tau,2,3}]]
 (* ::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]]];
 averagesErrorBars=Map[
 {{#[[1,1,1]],Mean[#[[All,2]]]},
 ErrorBar[StandardDeviation[#[[All,2]]]]
 (* ::Subsection:: *)
 (*Fitting the log-log-plot*)
 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]];
 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)*)
 plotRangeByTau[tau_]:=Piecewise[{{{0,30000},tau<2.3},{{0,4000},2.3<tau<2.7},{{0,800},tau>2.7}},Automatic]
 histograms=Map[Histogram[#[[All,2]],PlotRange->{plotRangeByTau[#[[1,1,2]]],Automatic}]&,averagesGrouped,{2}];

0 comments (0 inline, 0 general)