diff --git a/triangle_exponent_plots.m b/triangle_exponent_plots.m index b88398c79d7059660338a64c3be02f5d1493305d..230af3640f86d2ac3426dabf6366eaff7a62f246 100644 --- a/triangle_exponent_plots.m +++ b/triangle_exponent_plots.m @@ -15,7 +15,10 @@ Needs["ErrorBarPlots`"] (*#triangles = const\[Cross]n^(T(\[Tau])) where T(\[Tau])=3/2 (3-\[Tau])*) -(* When importing from exponent-only-data file *) +(* When importing from exponent-only-data file or property-data file *) +(* graphdata_exponent_mix32.m *) +(* graphdata_exponent_highN.m *) +(* graphdata_properties2.m *) gsraw=Import[NotebookDirectory[]<>"data/graphdata_exponent_mix32.m"]; gsraw=SortBy[gsraw,#[[1,1]]&]; (* Sort by n *) averagesGrouped=GatherBy[gsraw,{#[[1,2]]&,#[[1,1]]&}]; @@ -40,7 +43,9 @@ ListLogLogPlot[averagesErrorBars[[All,All,1]],Joined->True,PlotMarkers->Automati (*Fitting the log-log-plot*) -loglogdata=Log[averagesErrorBars[[All,All,1]]]; +nRange=5;;-1; +nRange=All; +loglogdata=Log[averagesErrorBars[[All,nRange,1]]]; fits=Map[Fit[#,{1,logn},logn]&,loglogdata]; fitsExtra=Map[LinearModelFit[#,logn,logn]&,loglogdata];