In[]:=
data=ParallelTable[SeedRandom[2524+i];DeathLossFunction[#,200]&/@NestList[First[MinimalBy[Append[Table[CARuleMutation[#],1],#],DeathLossFunction[#,200]&]]&,{0,3,1},300],{i,500}];
In[]:=
data5000=ParallelTable[SeedRandom[2524+i];DeathLossFunction[#,200]&/@NestList[First[MinimalBy[Append[Table[CARuleMutation[#],1],#],DeathLossFunction[#,200]&]]&,{0,3,1},300],{i,5000}];
In[]:=
ListLinePlot[data,PlotRange->All]
Out[]=
In[]:=
ListLinePlot[Mean/@Transpose[data]]
Out[]=
In[]:=
ListLinePlot[Table[Quantile[#,frac]&/@Transpose[data],{frac,0,1,.05}]]
Out[]=
In[]:=
Histogram[Flatten[Reverse[Transpose[With[{u=PDF[SmoothKernelDistribution[#],x]},Table[u,{x,-50,0,.5}]]&/@Transpose[data]]]]]
Out[]=
In[]:=
dxd=Flatten[Reverse[Transpose[With[{u=PDF[SmoothKernelDistribution[#],x]},Table[u,{x,-50,0,.5}]]&/@Transpose[data]]]];
In[]:=
Histogram[dxd^.1]
Out[]=
In[]:=
Show[ArrayPlot[Reverse[Transpose[With[{d=SmoothKernelDistribution[#]},Table[PDF[d,x]/PDF[d,Median[d]],{x,-100,0,1}]]&/@Transpose[data5000]]],AspectRatio->1/GoldenRatio,ColorFunction->(GrayLevel[1-#^.2]&),DataRange->{{0,300},{-100,0}}],ListLinePlot[Mean/@Transpose[data5000],PlotStyle->Red],ListLinePlot[Min/@Transpose[data5000],PlotStyle->Red],ListLinePlot[Transpose[Quartiles/@Transpose[data5000]],PlotStyle->Red]]
Out[]=
In[]:=
ListPlot3D[Reverse[Transpose[With[{d=SmoothKernelDistribution[#]},Table[PDF[d,x]/PDF[d,Median[d]],{x,-100,0,1}]]&/@Transpose[data5000]]]]
Out[]=
