数量データからの母数推定
北米の主要な河川の長さを単位をマイルとして得る.
In[1]:=
![Click for copyable input](assets.ja/parameter-estimation-from-quantity-data/In_34.png)
data = Quantity[ExampleData[{"Statistics", "RiverLengths"}], "Miles"];
In[2]:=
![Click for copyable input](assets.ja/parameter-estimation-from-quantity-data/In_35.png)
h = Histogram[data, 30, PDF, Frame -> True, FrameLabel -> Automatic,
PlotTheme -> "Detailed"]
Out[2]=
![](assets.ja/parameter-estimation-from-quantity-data/O_31.png)
EstimatedDistributionは数量データを取り,QuantityDistributionをデータの単位とともに返す.
In[3]:=
![Click for copyable input](assets.ja/parameter-estimation-from-quantity-data/In_36.png)
edist1 = EstimatedDistribution[data,
InverseGammaDistribution[a, b, c, d]]
Out[3]=
![](assets.ja/parameter-estimation-from-quantity-data/O_32.png)
異なる単位でモデルを推定する.
In[4]:=
![Click for copyable input](assets.ja/parameter-estimation-from-quantity-data/In_37.png)
edist2 = EstimatedDistribution[data,
QuantityDistribution[InverseGammaDistribution[a, b, c, d],
"Kilometers"]]
Out[4]=
![](assets.ja/parameter-estimation-from-quantity-data/O_33.png)
最初のモデルをキロメートルに変換する.
In[5]:=
![Click for copyable input](assets.ja/parameter-estimation-from-quantity-data/In_38.png)
UnitConvert[edist1, "Kilometers"]
Out[5]=
![](assets.ja/parameter-estimation-from-quantity-data/O_34.png)
データヒストグラムをフィットされたモデルの確率密度関数とともにプロットする.
In[6]:=
![Click for copyable input](assets.ja/parameter-estimation-from-quantity-data/In_39.png)
Show[h, Plot[PDF[edist1, Quantity[x, "Miles"]], {x, 0, 2000}]]
Out[6]=
![](assets.ja/parameter-estimation-from-quantity-data/O_35.png)