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

Changeset - 5e8f0e15e05e

Parent rev.

Child rev.

[Not reviewed]

0 2 0

Tom Bannink - 8 years ago 2017-07-13 16:28:06
tom.bannink@cwi.nl

Update exponent plot

2 files changed with 10 insertions and 1 deletions:

plots/triangle_exponent.pdf

bin+mod

triangle_exponent_plots.m

0 comments (0 inline, 0 general)

plots/triangle_exponent.pdf

➞

Show inline comments

binary diff not shown

triangle_exponent_plots.m

➞

Show inline comments

@@ @@ -17,12 +17,13 @@ Needs["ErrorBarPlots`"] @@
 (* When importing from exponent-only-data file or property-data file *)
 (* graphdata_exponent_mix32.m *)
 (* graphdata_exponent_highN.m *)
 (* graphdata_properties2.m *)
 (* graphdata_canonical_properties.m *)
 (* graphdata_exponent_accurate1.m *)
 gsraw=Import[NotebookDirectory[]<>"data/graphdata_exponent_hightau.m"];
 gsraw=SortBy[gsraw,#[[1,1]]&]; (* Sort by n *)
 averagesGrouped=GatherBy[gsraw,{#[[1,2]]&,#[[1,1]]&}];
 gsraw2=Import[NotebookDirectory[]<>"data/graphdata_canonical_properties2.m"];
 mediansFits=Map[Fit[#,{1,logn},logn]&,mediansLoglogdata];
 mediansFitsExtra=Map[LinearModelFit[#,logn,logn]&,mediansLoglogdata];
 canonicalloglogdata=Log[canonicalDatapoints[[All,nRange]]];
 canonicalFitsExtra=Map[LinearModelFit[#,logn,logn]&,canonicalloglogdata];
 canonicalNLMs=Map[NonlinearModelFit[#,factor*n^exponent,{factor,exponent},n]&,canonicalDatapoints];
 averagesFitsExtra[[1]]["ParameterConfidenceIntervalTable"]
 averagesFitsExtra[[1]]["BestFitParameters"]
 averagesFitsExtra[[1]]["ParameterErrors"]
 averagesFitsExtra[[1]]["ParameterConfidenceIntervals"]
@@ @@ -91,27 +93,34 @@ Export[NotebookDirectory[]<>"plots/avgtris_n.pdf",plot1] @@
 tauValues=averagesGrouped[[All,1,1,1,2]];
 averagesExponents=Transpose[{tauValues,averagesFits[[All,2,1]]}];
 mediansExponents=Transpose[{tauValues,mediansFits[[All,2,1]]}];
 Show[ListPlot[{averagesExponents,mediansExponents},Joined->True,PlotMarkers->Automatic,AxesLabel->{"tau","exponent"},PlotRange->{{2,3},{0,1.6}}],Plot[3/2(3-tau),{tau,2,3}]]
 canonicalNLMs[[1]]["ParameterConfidenceIntervals"]
 (* ::Subsection:: *)
 (*T(\[Tau]) including error bars*)
 (* For visual, shift the tau values slightly left or right to distinguish the two datasets *)
 tauValues=averagesGrouped[[All,1,1,1,2]];
 canonicalTauValues=averagesGrouped2[[All,1,1,1,2]];
 averagesExponentsErrorBars=Map[{{#[[1]],#[[2]]["BestFitParameters"][[2]]},ErrorBar[#[[2]]["ParameterConfidenceIntervals"][[2]]-#[[2]]["BestFitParameters"][[2]]]}&,
 Transpose[{tauValues-0.003,averagesFitsExtra}]];
 mediansExponentsErrorBars=Map[{{#[[1]],#[[2]]["BestFitParameters"][[2]]},ErrorBar[#[[2]]["ParameterConfidenceIntervals"][[2]]-#[[2]]["BestFitParameters"][[2]]]}&,
 Transpose[{tauValues+0.003,mediansFitsExtra}]];
 canonicalExponentsErrorBars=Map[{{#[[1]],#[[2]]["BestFitParameters"][[2]]},ErrorBar[#[[2]]["ParameterConfidenceIntervals"][[2]]-#[[2]]["BestFitParameters"][[2]]]}&,
 Transpose[{tauValues+0.000,canonicalFitsExtra}]];
 Transpose[{canonicalTauValues+0.000,canonicalFitsExtra}]];
 canonicalExponentsErrorBars2=Map[{{#[[1]],exponent/.#[[2]]["BestFitParameters"][[2]]},ErrorBar[#[[2]]["ParameterConfidenceIntervals"][[2]]-exponent/.#[[2]]["BestFitParameters"]]}&,
 Transpose[{canonicalTauValues+0.000,canonicalNLMs}]];
 plot2=Show[
 ErrorListPlot[{averagesExponentsErrorBars,mediansExponentsErrorBars,canonicalExponentsErrorBars},
 Joined->True,PlotMarkers->Automatic,
 PlotLegends->Placed[{"average","median","canonical"},{Left,Bottom}],
 Frame->True,FrameLabel->{"tau","triangle exponent"},
 PlotRange->{{2,3},{0,1.6}},
 ImageSize->300],
 Plot[3/2(3-tau),{tau,2,3},PlotStyle->{Black,Dashed}]]

0 comments (0 inline, 0 general)