模拟数量分布
根据数量分布模拟等待时间.
In[1]:=
![Click for copyable input](assets.zh/simulate-quantities/In_6.png)
wait\[ScriptCapitalD] =
QuantityDistribution[ExponentialDistribution[.4], "Minutes"]
Out[1]=
![](assets.zh/simulate-quantities/O_5.png)
一个观测即是一个 Quantity.
In[2]:=
![Click for copyable input](assets.zh/simulate-quantities/In_7.png)
RandomVariate[wait\[ScriptCapitalD]]
Out[2]=
![](assets.zh/simulate-quantities/O_6.png)
多条观测的样本组成一个 QuantityArray.
In[3]:=
![Click for copyable input](assets.zh/simulate-quantities/In_8.png)
sample = RandomVariate[wait\[ScriptCapitalD], 100]
Out[3]=
![](assets.zh/simulate-quantities/O_7.png)
用 Histogram 可视化样本.
In[4]:=
![Click for copyable input](assets.zh/simulate-quantities/In_9.png)
Histogram[sample, AxesLabel -> Automatic]
Out[4]=
![](assets.zh/simulate-quantities/O_8.png)
计算样本的几个选定的描述性统计量.
In[5]:=
![Click for copyable input](assets.zh/simulate-quantities/In_10.png)
stats = {Mean, Median, Min, Max};
TableForm[Map[#[sample] &, stats], TableHeadings -> {stats, None}]
Out[5]//TableForm=
![](assets.zh/simulate-quantities/O_9.png)