数量のシミュレーション
数量分布から待機時間のシミュレーションを行う.
In[1]:=
![Click for copyable input](assets.ja/simulate-quantities/In_6.png)
wait\[ScriptCapitalD] =
QuantityDistribution[ExponentialDistribution[.4], "Minutes"]
Out[1]=
![](assets.ja/simulate-quantities/O_5.png)
単一の観察の結果は,Quantityオブジェクトで表される.
In[2]:=
![Click for copyable input](assets.ja/simulate-quantities/In_7.png)
RandomVariate[wait\[ScriptCapitalD]]
Out[2]=
![](assets.ja/simulate-quantities/O_6.png)
複数の観察のサンプル結果は,QuantityArrayオブジェクトで表される.
In[3]:=
![Click for copyable input](assets.ja/simulate-quantities/In_8.png)
sample = RandomVariate[wait\[ScriptCapitalD], 100]
Out[3]=
![](assets.ja/simulate-quantities/O_7.png)
Histogramを使ってサンプルを可視化する.
In[4]:=
![Click for copyable input](assets.ja/simulate-quantities/In_9.png)
Histogram[sample, AxesLabel -> Automatic]
Out[4]=
![](assets.ja/simulate-quantities/O_8.png)
サンプルの記述統計の一部を計算する.
In[5]:=
![Click for copyable input](assets.ja/simulate-quantities/In_10.png)
stats = {Mean, Median, Min, Max};
TableForm[Map[#[sample] &, stats], TableHeadings -> {stats, None}]
Out[5]//TableForm=
![](assets.ja/simulate-quantities/O_9.png)