BinomialDistribution[n,p]
n 回の試行での成功確率が p である二項分布を表す.
BinomialDistribution
BinomialDistribution[n,p]
n 回の試行での成功確率が p である二項分布を表す.
詳細
- 二項分布における値
の確率は,0から n までの整数について
である. » - BinomialDistributionでは,n は任意の非負の整数でよい.
- BinomialDistributionでは,n と p は無次元量でよい. »
- BinomialDistributionは,Mean,CDF,RandomVariate等の関数で使うことができる. »
予備知識
- BinomialDistribution[n,p]は,整数値
のところで定義され,非負の実数 p (
)でパラメータ化された離散統計分布を表す.二項分布は,単峰性で峰が
の平均にある,離散確率密度関数(PDF)を持つ.母数 n および p は,PDFの,高さ,水平位置,および歪みを決定する. - 二項分布の研究は,18世紀初頭のヤコブ・ベルヌーイの業績まで遡ることができ,研究対象となる最古の分布の一つとなっている.二項分布は,(公正あるいは不正な)n 個のコイントスの行為をモデル化のために設計された.トスするコインは,独立で等しく,独立的に順次置換えたものでサンプルを取られる.伝統的に,p は実験が「成功する」確率であり,
は「失敗する」確率であると考えられる.コイントスの例では,
は公正なコインのトスに相当する. - 非常に単純に定義されているにもかかわらず,二項分布は,より複雑な多くの数学概念の基礎になっている.例えば,二項分布は,成功確率が p の有限の n 段階ベルヌーイ過程(つまり,それぞれがBernoulliDistribution[p]によって独立どう分布に従っている,有限な一連の確率変数からなる離散時間確率過程)から生じた成功の分布として考えることができる.同様に,二項分布は,二項過程(BinomialProcess)のスライス分布(SliceDistribution)である.この二項過程は,二項分布に従う独立同分布の確率変数の有限列からなる離散時間・離散状態の確率過程であり,確率変数間の時間は幾何分布 (GeometricDistribution)に従う.確率変数の有限列からなる)離散時間・離散状態の確率過程,さらに,例えば n 個の(公正または不正な)サイコロを振って特定の値が出る確率等,多くの現実世界の事象が二項過程としてモデル化できる.移民パターンやある種の待ち行列モデル等で二項分布間の驚くべき関係も発見されている.
- RandomVariateを使って,二項分布から,1つあるいは複数の機械精度あるいは任意精度(後者はWorkingPrecisionオプションを介す)の擬似乱数変量を得ることができる.Distributed[x,BinomialDistribution[n,p]](より簡略すると xBinomialDistribution[n,p])を使って,確率変数 x が,二項分布に従って分布していると宣言することができる.このような宣言は,Probability,NProbability,Expectation,NExpectation等の関数で使うことができる.
- 確率密度関数および累積分布関数は,PDF[BinomialDistribution[n,p],x]およびCDF[BinomialDistribution[n,p],x]を使って得られることがある.平均,中央値,分散,原点の周りのモーメント,中心モーメントは,それぞれMean,Median,Variance,Moment,CentralMomentを使って計算することができる.これらの数量はDiscretePlotを使って可視化することができる.
- DistributionFitTestを使って,与えられたデータ集合が二項分布と一致するかどうかを検定することが,EstimatedDistributionを使って与えられたデータから二項パラメトリック分布を推定することが,FindDistributionParametersを使ってデータを二項分布にフィットすることができる.ProbabilityPlotを使って記号的な二項分布のCDFに対する与えられたデータのCDFのプロットを生成することが,QuantilePlotを使って記号的な二項分布の変位値に対する与えられたデータの変位値のプロットを生成することができる.
- TransformedDistributionを使って変換された二項分布を表すことが,CensoredDistributionを使って上限値と下限値の間で切り取られた値の分布を表すことが,TruncatedDistributionを使って上限値と下限値の間で切断された値の分布を表すことができる.CopulaDistributionを使って二項分布を含む高次元分布を構築することが,ProductDistributionを使って二項分布を含む独立成分分布の結合分布を計算することができる.
- BinomialDistributionは,数多くの他の統計分布と関係がある.例えば,BinomialDistribution[1,p]は,値
および
において,BernoulliDistribution[p]と厳密に等しく,BernoulliDistribution[p]に従って分布している n 個の独立変数の合計はBinomialDistribution[n,p]に従って分布している.同様に,BinomialDistribution[t,p]はBinomialProcess[p][t]と同じPDFを持つ.この特性は,BinomialDistribution[t,p]が厳密にSliceDistribution[BinomialProcess[p],t]であることによる.BinomialDistribution[n,p]はまた,いくつかの分布についての極限分布である.特に,BinomialDistribution[n,p]は,固定した p の値について,n がInfinityに近付く(ただし,
,
)につれてNormalDistribution[μ,σ]に収束し,n がInfinityに,p が0に近付くにつれて離散化された PoissonDistribution[p]に収束する.これに対し,HypergeometricDistribution[n,p]は,p がInfinityに近付くにつれてBinomialDistribution[n,p]に限定する.BinomialDistributionはMultinomialDistributionの2変数の形であり,BetaBinomialDistributionの構成成分であり,NegativeBinomialDistributionと自然に関係している.
例題
すべて開く すべて閉じる例 (3)
DiscretePlot[Table[PDF[BinomialDistribution[40, p], k], {p, {0.1, 0.5, 0.7}}]//Evaluate, {k, 36}, PlotRange -> All, PlotMarkers -> Automatic]DiscretePlot[Table[PDF[BinomialDistribution[n, .6], k], {n, {5, 20, 40}}]//Evaluate, {k, 32}, PlotRange -> All, PlotMarkers -> Automatic]PDF[BinomialDistribution[n, p], k]DiscretePlot[Table[CDF[BinomialDistribution[40, p], k], {p, {0.1, 0.5, 0.7}}]//Evaluate, {k, 36}, ExtentSize -> Right]DiscretePlot[Table[CDF[BinomialDistribution[n, .6], k], {n, {5, 20, 40}}]//Evaluate, {k, 32}, ExtentSize -> Right]CDF[BinomialDistribution[n, p], k]Mean[BinomialDistribution[n, p]]Variance[BinomialDistribution[n, p]]スコープ (8)
data = RandomVariate[BinomialDistribution[10, 0.4], 10 ^ 4];Show[
Histogram[data, {1}, "PDF"],
DiscretePlot[PDF[BinomialDistribution[10, 0.4], x], {x, 0, 10}, PlotStyle -> PointSize[Medium]]]sample = RandomVariate[BinomialDistribution[20, .3], 10 ^ 3];edist = EstimatedDistribution[sample, BinomialDistribution[n, p]]サンプルの密度ヒストグラムを推定分布の確率密度関数と比較する:
Show[Histogram[sample, {1}, "PDF"], DiscretePlot[PDF[edist, x], {x, 0, 15}, PlotStyle -> PointSize[Medium]]]Plot[Table[Skewness[BinomialDistribution[n, p]], {n, r = {1, 2, 5, 20}}]//Evaluate, {p, 0, 1}, PlotLegends -> (Row[{"n = ", #}]& /@ r)]Skewness[BinomialDistribution[n, p]]Solve[Skewness[BinomialDistribution[n, p]] == 0, p]Limit[Skewness[BinomialDistribution[n, p]], n -> ∞]Plot[Table[Kurtosis[BinomialDistribution[n, p]], {n, r = {1, 2, 5, 20}}]//Evaluate, {p, 0, 1}, PlotLegends -> (Row[{"n = ", #}]& /@ r)]Kurtosis[BinomialDistribution[n, p]]極限値は標準NormalDistributionの尖度の値である:
Limit[Kurtosis[BinomialDistribution[n, p]], n -> ∞]FormulaGrid[list_, type_] := Grid[...]FormulaGrid[Table[Moment[BinomialDistribution[n, p], k]//FullSimplify, {k, 5}], M]FormulaGrid[Table[CentralMoment[BinomialDistribution[n, p], k]//FullSimplify, {k, 5}], CM]FormulaGrid[Table[FactorialMoment[BinomialDistribution[n, p], k]//Simplify, {k, 5}], FM]FactorialMoment[BinomialDistribution[n, p], r]FormulaGrid[Table[Cumulant[BinomialDistribution[n, p], k], {k, 5}], C]DiscretePlot[Table[HazardFunction[BinomialDistribution[20, p], x], {p, {0.2, .5, 0.8}}]//Evaluate, {x, 20}, PlotMarkers -> Automatic]HazardFunction[BinomialDistribution[n, p], x]𝒟 = BinomialDistribution[10, .6];
points = Table[CDF[𝒟, x], {x, 0, 10}]//N;DiscretePlot[Quantile[𝒟, p], {p, points}, ExtentSize -> Right, ExtentMarkers -> {"Empty", "Filled"}]無次元のQuantityを使ってBinomialDistributionを定義する:
BinomialDistribution[Quantity[1, "Gross"], Quantity[30, "Percent"]]アプリケーション (12)
BinomialDistributionのCDFは右連続関数の例である:
DiscretePlot[CDF[BinomialDistribution[5, .6], x], {x, 0, 5}, ExtentSize -> Right, ExtentMarkers -> {"Filled", "Empty"}]CDF[BinomialDistribution[5, .6], x]//PiecewiseExpandフリースローの成功確率が75%のバスケットボール選手がいる.10回のフリースローのシミュレーションを行う:
RandomVariate[BernoulliDistribution[0.75], 10] /. {0 -> "Miss", 1 -> "Hit"}この選手が試合のフリースローで3回に2回成功する確率を求める:
Probability[x == 2, xBinomialDistribution[3, 0.75]]この選手が5回のフリースローで最後の2回成功する確率を求める:
Probability[x == 0, xBinomialDistribution[3, 0.75]]Probability[x == 2, xBinomialDistribution[2, 0.75]]n 回のフリースローがある試合での期待される成功回数を求める:
Expectation[x, xBinomialDistribution[n, 0.75]]打率0.300の野球選手がいる.5打席のシミュレーションを行う:
RandomVariate[BernoulliDistribution[0.300], 5]この選手が3回打席に立った場合の期待されるヒット数を求める:
Expectation[x, xBinomialDistribution[3, 0.300]]ある薬が30%のケースで有効であると証明された,患者4人のうち少なくとも3人に有効である確率はどうなるだろうか:
Probability[x ≥ 3, xBinomialDistribution[4, 0.3]]Expectation[x, xBinomialDistribution[500, 0.3]]公正なコインを n 回はじいた場合に表が出る回数はBinomialDistributionでモデル化することができる:
heads[n_] := BinomialDistribution[n, 1 / 2]DiscretePlot[PDF[heads[100], k], {k, 20, 80}, PlotRange -> All]100回コインをはじいた場合に表が60回から80回出る確率を計算する:
Probability[60 ≤ x ≤ 80, xheads[100]]N[%]不正なコインで表が出る確率が0.6だとした場合は以下のようになる:
uheads[n_] := BinomialDistribution[n, 0.6];DiscretePlot[PDF[uheads[100], k], {k, 20, 80}, PlotRange -> All]Probability[60 ≤ x ≤ 80, xuheads[100]]defective[n_] := BinomialDistribution[n, 1 / 10]Probability[x ≤ 1, xdefective[5]]N[%]航空機のエンジンに不具合が生じる確率は p である.2/4以下の割合で不具合が生じる場合を計算する:
prob1 = Probability[x ≤ 2, xBinomialDistribution[4, p]]prob2 = Probability[x ≤ 1, xBinomialDistribution[2, p]]4つのエンジンを選んだ方が2つのエンジンよりも優れている場合を考える:
Reduce[prob1 > prob2 && 0 < p < 1, p]Plot[{prob1, prob2}, {p, 0, 1}, PlotLegends -> {"4 engines", "2 engines"}]3つのマイクロプロセッサを有しそのうちの1つが機能する限り使えるように設計された三重冗長性を使うシステムがある.1つのマイクロプロセッサが
秒後に使用可能な確率は
である.
秒たってもシステムが使用可能な確率を求める:
Probability[x ≥ 1, xBinomialDistribution[3, Exp[-λ t]]]各プロセッサの故障までの平均時間が
であるとして,システムが99%未満の確率で動作するようになる時期を求める:
Reduce[(% /. λ -> 10^-8) ≤ 0.99, t, Reals]//QuietLast[%] / (3600 24 365)チェスのチャンピオンであるガルリ・カスパロフがトーナメントで同時に100人のアマチュアを相手にしてチェスを行う.このような場合に彼が負ける確率は約1%であると見積もることができる.彼が0,2,5,10ゲームに負ける確率を求める:
Table[Probability[x == k, xBinomialDistribution[100, 0.01]], {k, {0, 2, 5, 10}}]Table[Probability[x == k, xPoissonDistribution[100 0.01]], {k, {0, 2, 5, 10}}]チャンピオンが5ゲームを行う場合について同じ計算をする.ただし,対戦相手はより強く,チャンピオンが負ける確率は10%である:
Table[Probability[x == k, xBinomialDistribution[5, 0.10]], {k, {0, 2, 5, 10}}]Table[Probability[x == k, xPoissonDistribution[5 0.1]], {k, {0, 2, 5, 10}}]n 個の記号の文字列からなるパケットがノイズの多い通信路を通して送られた.各記号が正しく送られない確率は
である.パケットが正しく送られない確率が
より小さくなる n を求める:
Probability[x ≥ 1, xBinomialDistribution[n, 10^-4]]//Simplify[#, n > 0]&Reduce[0 ≤ % ≤ 10^-3, n, Integers]Probability[x ≥ 1, xPoissonDistribution[n 10^-4]]Reduce[0 ≤ % ≤ 10^-3, n, Integers]各顧客があるサービスを確率 p で使っている場合に n 人のうち
人の顧客がサービスを必要とする確率を求める:
PDF[BinomialDistribution[n, p], k]同時に
(キャパシティ)を超えるサービスのリクエストがある確率を計算する:
Probability[x > c, xBinomialDistribution[n, p], Assumptions -> c∈Integers && 0 ≤ c ≤ n]
で
のときにサービスが受けられる確率をさまざまなキャパシティ
について計算する:
DiscretePlot[CDF[BinomialDistribution[100, 0.1], c], {c, 30}]99.9%の確率でサービスが受けられる最小のキャパシティ
を求める:
Minimize[{c, CDF[BinomialDistribution[100, 1 / 10], c] ≥ 999 / 1000 && 0 < c ≤ 100}, c, Integers]{CDF[BinomialDistribution[100, 0.1], 19], CDF[BinomialDistribution[100, 0.1], 20]}2人のプレーヤーがサイコロを振っている.両方の目の合計が10に満たない場合は2番目のプレーヤーが4セントを受け取り,それ以外場合は最初のプレーヤーが9セント受け取る.このゲームは公平と言えるだろうか:
p = With[{𝒟 = DiscreteUniformDistribution[{1, 6}]}, Probability[k1 + k2 >= 10, k1𝒟 && k2𝒟]]1ゲームあたりの平均スコアが等しくないのでこのゲームは公平ではない:
{pay1, pay2} = {9, 4};
{pay1 p, pay2(1 - p)}//Nn 回ゲームを行った段階で不利な立場にあるプレーヤーのスコアが高くなる確率を求める:
f = FunctionExpand[Probability[pay1 * k > pay2 * (n - k), kBinomialDistribution[n, p]], n∈Integers && n ≥ 1]DiscretePlot[f, {n, 50}, PlotRange -> All]Maximize[{f, 1 < n < 50}, n, Integers]特性と関係 (9)
TransformedDistribution[x + y, {xBinomialDistribution[Subscript[n, 1], p], yBinomialDistribution[Subscript[n, 2], p] }]BinomialDistribution[n,p]は
で正規分布に収束する:
{μ, σ} = {Mean[BinomialDistribution[n, p]], StandardDeviation[BinomialDistribution[n, p]]}Block[{p = 1 / 3},
Table[Show[DiscretePlot[PDF[BinomialDistribution[n, p], k], {k, Ceiling[μ - 2σ], Floor[μ + 2σ]}, PlotRange -> {{μ - 3σ, μ + 2σ}, All}, AxesOrigin -> {μ - 3σ, 0}], Plot[PDF[NormalDistribution[μ, σ], x], {x, μ - 2σ, μ + 2σ}], Ticks -> {{{μ - 2σ, HoldForm[μ - 2σ]}, {μ, HoldForm[μ]}, {μ + 2σ, HoldForm[μ + 2σ]}}, Automatic}, PlotRange -> {{μ - 3σ, μ + 2σ}, All}, PlotLabel -> n], {n, {5, 25, 50, 100}}]]
のBinomialDistributionはBernoulliDistributionに等しい:
Table[PDF[BinomialDistribution[1, p], k], {k, 0, 1}]Table[PDF[BernoulliDistribution[p], k], {k, 0, 1}]BernoulliDistributionを示す n 個の独立変数の和は二項分布を示す:
TransformedDistribution[u + v + w, {u, v, w}ProductDistribution[{BernoulliDistribution[p], 3}]]BinomialDistributionはHypergeometricDistributionの母集団の上限を無限にしたものである:
Assuming[n > 0 && 0 < p < 1, Limit[PDF[HypergeometricDistribution[n, p * ntot, ntot], k], ntot -> Infinity]]PDF[BinomialDistribution[n, p], k]//FunctionExpandSimplify[% - %%]n が大きく p が小さい場合,BinomialDistributionはPoissonDistributionに近付く:
Limit[PDF[BinomialDistribution[n, μ / n], k], n -> ∞]PDF[PoissonDistribution[μ], k]FullSimplify[% - %%, k∈Integers]MarginalDistribution[MultinomialDistribution[n, {p, 1 - p}], 1]PDF[MultinomialDistribution[n, {p, 1 - p}], {x, n - x}]PDF[BinomialDistribution[n, p], x]FullSimplify[% - %%]BetaBinomialDistributionはBinomialDistributionとBetaDistributionを混合したものである:
ParameterMixtureDistribution[BinomialDistribution[n, p], pBetaDistribution[α, β]]考えられる問題 (3)
BinomialDistributionは,p が0と1の間にない場合は定義されない:
Mean[BinomialDistribution[8, 1.1]]BinomialDistributionは,n が正の整数ではない場合は定義されない:
Mean[BinomialDistribution[3 / 2, 0.5]]Mean[BinomialDistribution[n, p]] /. {n -> -90, p -> 3}おもしろい例題 (1)
0から1までの区間でBernstein分布に従う
の多項式近似を構築する:
f[x_] := Sin[Pi x] / (1 + x ^ 2)
の近似は
の期待値として構築できる.
は,
の平均が
と等しくなるような母数が
と
の二項確率変数である:
g[n_, x_] := Expectation[f[k / n], kBinomialDistribution[n, x]]Show[Plot[f[x], {x, 0, 1}, PlotStyle -> Directive[Magenta, Dashed]], Plot[Table[g[n, x], {n, r = Reverse@{2, 5, 10, 20, 50}}]//Evaluate, {x, 0, 1}, PlotLegends -> (Row[{"n = ", #}]& /@ r)]]テクニカルノート
-
▪
- 離散分布
関連するガイド
-
▪
- 離散一変量分布 ▪
- 有界領域分布 ▪
- 壷モデル分布 ▪
- 信頼性解析の分布 ▪
- 通信システムの分布 ▪
- パラメトリック統計分布 ▪
- 統計学で使用される関数
テキスト
Wolfram Research (2007), BinomialDistribution, Wolfram言語関数, https://reference.wolfram.com/language/ref/BinomialDistribution.html (2016年に更新).
CMS
Wolfram Language. 2007. "BinomialDistribution." Wolfram Language & System Documentation Center. Wolfram Research. Last Modified 2016. https://reference.wolfram.com/language/ref/BinomialDistribution.html.
APA
Wolfram Language. (2007). BinomialDistribution. Wolfram Language & System Documentation Center. Retrieved from https://reference.wolfram.com/language/ref/BinomialDistribution.html
BibTeX
@misc{reference.wolfram_2026_binomialdistribution, author="Wolfram Research", title="{BinomialDistribution}", year="2016", howpublished="\url{https://reference.wolfram.com/language/ref/BinomialDistribution.html}", note=[Accessed: 13-August-2026]}
BibLaTeX
@online{reference.wolfram_2026_binomialdistribution, organization={Wolfram Research}, title={BinomialDistribution}, year={2016}, url={https://reference.wolfram.com/language/ref/BinomialDistribution.html}, note=[Accessed: 13-August-2026]}