Files
@ 39d31ff390d2
Branch filter:
Location: AENC/resampling_chain/montecarlo/plots.m - annotation
39d31ff390d2
816 B
application/vnd.wolfram.mathematica.package
Add mathematica notebook to plot montecarlo data
39d31ff390d2 39d31ff390d2 39d31ff390d2 39d31ff390d2 39d31ff390d2 39d31ff390d2 39d31ff390d2 39d31ff390d2 39d31ff390d2 39d31ff390d2 39d31ff390d2 39d31ff390d2 39d31ff390d2 39d31ff390d2 39d31ff390d2 39d31ff390d2 39d31ff390d2 39d31ff390d2 39d31ff390d2 39d31ff390d2 39d31ff390d2 39d31ff390d2 39d31ff390d2 39d31ff390d2 39d31ff390d2 39d31ff390d2 39d31ff390d2 39d31ff390d2 39d31ff390d2 39d31ff390d2 | (* ::Package:: *)
(* ::Section:: *)
(*P(Z^n | start 011..11) on the chain {1,...,n}*)
rawdata=Import[NotebookDirectory[]<>"data/data_reach_chain_end.m"];
data=GatherBy[rawdata,{#[[2]]}&];
(* data[[n index, prob index, p / n / prob / runs within timelimit / total runs]] *)
nlabels=Map["n = "<>ToString[#[[1,2]]]&,data];
plotZoomed=ListPlot[data[[All,All,{1,3}]],
Joined->True,
PlotMarkers->Automatic,
PlotRange->{{0.6,0.7},{0,0.5}},
Frame->True,
ImageSize->250,
Epilog->{Dashed,Line[{{1-1/E,0},{1-1/E,1}}]}
];
plotBoth=ListPlot[data[[All,All,{1,3}]],
Joined->True,
FrameLabel->{"p (probability of BAD)","P(\!\(\*SuperscriptBox[\(Z\), \(n\)]\) | start 011..11)"},
PlotMarkers->Automatic,
PlotLegends->nlabels,
PlotRange->{{0,1},{0,1}},
Frame->True,
ImageSize->500,
Epilog->Inset[plotZoomed,Scaled[{0.32,0.6}]]]
|