概要

Source code in Mathematica

g = Table[
  ListPlot[
    Table[{k, Binomial[k + r - 1, r - 1] (r/(10 + r))^r (1 - r/(10 + r))^k}, {k, 0, 25}],
    PlotRange -> 0.118, 
    Filling -> Axis,
    PlotStyle -> PointSize[0.02], 
    FillingStyle -> Directive[Hue[0.6], Thickness[0.01], Opacity[0.3]],
    Epilog -> Inset[Style["\!\(\*StyleBox[\"r\",\nFontSlant->\"Italic\"]\) = " <> ToString[r], 16], {23, 0.08}], 
    Prolog -> {Opacity[0.5], Green, Thin, Line[{{10 - Sqrt[(10 (10 + r))/r], .1}, {10 + Sqrt[(10 (10 + r))/r], .1}}], 
               Dashed, Thick, Orange, Line[{{10, 0}, {10, 1}}]},
    ImageSize -> 200,
    LabelStyle -> Directive[Larger],
    AxesLabel -> {"\!\(\*StyleBox[\"k\",\nFontSlant->\"Italic\"]\)"},
    AxesOrigin -> {-1, 0}
    ], 
  {r, {1, 2, 3, 4, 5, 10, 20, 40}}
  ]
Export["c:/users/stpasha/documents/negbinomial.gif", g, "DisplayDurations" -> {.5, .5, .5, .5, .5, .5, .5, .15}, "TransparentColor" -> White, ImageSize -> 250]

ライセンス

この著作物の著作権者である私は、この著作物における権利を放棄しパブリックドメインとします。これは全世界で適用されます。 一部の国では、これが法的に可能ではない場合があります。その場合は、次のように宣言します。 私は、あらゆる人に対して、法により必要とされている条件を除き、如何なる条件も課すことなく、あらゆる目的のためにこの著作物を使用する権利を与えます。