根据数量数据进行参数估计
获取北美主要河流的长度,单位为英里.
In[1]:=
data = Quantity[ExampleData[{"Statistics", "RiverLengths"}], "Miles"];
In[2]:=
h = Histogram[data, 30, PDF, Frame -> True, FrameLabel -> Automatic,
PlotTheme -> "Detailed"]
Out[2]=
EstimatedDistribution 接受数量数据,返回带有数据的单位的 QuantityDistribution.
In[3]:=
edist1 = EstimatedDistribution[data,
InverseGammaDistribution[a, b, c, d]]
Out[3]=
采用不同的单位,估计模型.
In[4]:=
edist2 = EstimatedDistribution[data,
QuantityDistribution[InverseGammaDistribution[a, b, c, d],
"Kilometers"]]
Out[4]=
把第一个模型的单位转换成公里.
In[5]:=
UnitConvert[edist1, "Kilometers"]
Out[5]=
同时绘制数据的直方图和拟合模型的 PDF.
In[6]:=
Show[h, Plot[PDF[edist1, Quantity[x, "Miles"]], {x, 0, 2000}]]
Out[6]=