Novos Recursos do
Mathematica
8: Estimativa de Parâmetros e Testes
◄
anterior
|
próximo
►
Principais Algoritmos
Compare dois modelos de velocidades de vento
Ajuste dados de velocidade de vento a modelos de Rayleigh e de valores extremos, compare visualmente os dados com as estimativas de densidade de kernel, e compare quantitativamente os cuantis empírico e ajustado.
In[1]:=
X
maxWinds = WeatherData["Boston", "MaxWindSpeed", {{1950, 1, 1}, {2009, 12, 31}, "Month"}][[ All, -1]];
In[2]:=
X
rdist = EstimatedDistribution[maxWinds, RayleighDistribution[\[Theta]]]; evdist = EstimatedDistribution[maxWinds, ExtremeValueDistribution[a, b]];
In[3]:=
X
GraphicsGrid[{Table[ Show[Histogram[maxWinds, Automatic, "PDF", PlotRange -> {{0, 140}, All}, ChartStyle -> "Aquamarine"], Plot[PDF[i, x], {x, 0, 140}, PlotRange -> {0, .045}, PlotStyle -> {Blend[{Blue, Black}], Thick}], SmoothHistogram[maxWinds, PlotStyle -> {Blend[{Black, Purple}], Thick}, FillingStyle -> Directive[Opacity[.5], ColorData[1, 2]]], PlotLabel -> Style[StringDrop[ToString[Head[i]], -12], FontFamily -> "Verdana"], Axes -> {True, False}], {i, {rdist, evdist}}], Table[QuantilePlot[maxWinds, i, PlotStyle -> {Blend[{Blue, Black}], Thick}], {i, {rdist, evdist}}]}]
Out[3]=