Dagum分布を使って収入のシミュレーションを行う
DagumDistributionを州立大学からの収入データにフィットさせ,そのような大学で可能な給料のシミュレーションを行う.
universityA = Select[ExampleData[{"Statistics", "UniversitySalaries"}], #[[4]] == "A"&];
salaries = Select[Table[If[universityA[[i, 2]] > 0, universityA[[i, 3]] / universityA[[i, 2]], 0], {i, Length[universityA]}], # > 0&];
edist = EstimatedDistribution[salaries, DagumDistribution[p, a, b]];Show[Histogram[salaries, {0, 200000, 10000}, "PDF", ChartStyle -> "BeachColors"], Plot[PDF[TruncatedDistribution[{0, 200000}, edist], x], {x, 0, 200000}, PlotStyle -> Thick], BaseStyle -> {FontFamily -> "Verdana"}, Epilog -> Inset[Framed[Style[Grid[{{"推定分布:"}, {Round[#, .1]& /@ edist}}], 11], Background -> StandardBlue, RoundingRadius -> 3], {Right, Top}, {Right, Top}], ImageSize -> 400]ListPlot[{RandomVariate[edist, 100], Tooltip[{{0, Median[edist]}, {100, Median[edist]}}, "Median"]}, Filling -> {1 -> Axis}, Joined -> {False, True}, AxesOrigin -> {0, 20000}, PlotStyle -> {PointSize[Medium], Thick}, BaseStyle -> {FontFamily -> "Verdana"}, PlotLabel -> "無作為に選ばれた100人の従業員の収入のシミュレーション", Epilog -> Inset[Style["中央値", Bold, 10], {75, Median[edist]}, {Center, Center}], ImageSize -> 400]