PascalDistribution[n,p]
母数が n と p のパスカル(Pascal)分布を表す.
PascalDistribution
PascalDistribution[n,p]
母数が n と p のパスカル(Pascal)分布を表す.
詳細
- パスカル分布は二項待ち時間分布とも呼ばれる.
- パスカル分布における整数値
の確率は,
では
であり,その他の場合はゼロである. - PascalDistribution[n,p]は非零の成功確率 p で n 回成功するまでの試行回数の分布を与える.
- PascalDistributionでは,n は正の整数でよく,p は1以下の正の実数でよい.
- PascalDistributionでは,n と p は無次元量でよい. »
- PascalDistributionは,Mean,CDF,RandomVariate等の関数とともに用いることができる.
予備知識
- PascalDistribution[n,p]は,整数値
について定義され,試行の成功数と成功確率をそれぞれ表す整数母数 n (
)と実数母数 p (
) によって決定された離散統計分布を表す.パスカル分布は離散的で単峰性の確率密度関数(PDF)を持つ.パスカル分布は,「負の二項分布」としてまとめられる数多くの分布の一つであるが,負の二項分布(NegativeBinomialDistribution)と混同してはならない. - パスカル分布は確率分布の中で最も早く研究されたものの一つであり,ブレーズ パスカル(Blaise Pascal)の1670年代の研究まで遡ることができる.パスカル分布は,負の二項分布族の一員として見ると,各取り出しが互いに独立で成功確率が p のときに,ある特定の色のビー玉を n 個壺から取り出すためのビー玉の取り出し回数を説明する壺モデルとして実現することができる.この分布は,誕生以来,NegativeBinomialDistribution[n,p]の整数 n の場合として実現されており,負の二項分布族に含まれているために,数多く応用されている(例:事故統計および無線通信).この分布のその他の応用例には,人口統計のモデル化,心理学データ,品質管理,待ち行列理論等がある.
- RandomVariateを使って,パスカル分布から,1つあるいは複数の機械精度あるいは任意精度(後者はWorkingPrecisionオプションを介す)の擬似乱数変量を得ることができる.Distributed[x,PascalDistribution[n,p]](より簡略な表記では xPascalDistribution[n,p])を使って,確率変数 x がパスカル分布に従って分布していると宣言することができる.このような宣言は,Probability,NProbability,Expectation,NExpectation等の関数で使うことができる.
- 確率密度関数および累積分布関数は,PDF[PascalDistribution[n,p],x]およびCDF[PascalDistribution[n,p],x]を使って得られることがある.平均,中央値,分散,原点の周りのモーメント,中心モーメントは,それぞれMean,Median,Variance,Moment,CentralMomentを使って計算することができる.これらの数量はDiscretePlotを使って可視化することができる.
- DistributionFitTestを使って,与えられたデータ集合がパスカル分布と一致するかどうかを検定することが,EstimatedDistributionを使って与えられたデータからパラメトリックパスカル分布を推定することが,FindDistributionParametersを使ってデータをパスカル分布にフィットすることができる.ProbabilityPlotを使って記号パスカル分布のCDFに対する与えられたデータのCDFのプロットを生成することが,QuantilePlotを使って記号パスカル分布の変位値に対する与えられたデータの変位値のプロットを生成することができる.
- TransformedDistributionを使って変換されたパスカル分布を表すことが,CensoredDistributionを使って上限値と下限値の間で切り取られた値の分布を表すことが,TruncatedDistributionを使って上限値と下限値の間で切断された値の分布を表すことができる.CopulaDistributionを使ってパスカル分布を含む高次元分布を構築することが,ProductDistributionを使ってパスカル分布を含む独立成分分布の結合分布を計算することができる.
- PascalDistributionは,他の数多くの統計分布と関連している.上述の通り,パルカル分布は負の二項分布族の一員であり,したがって,(非整数 n についてのポリヤ分布としても知られる)NegativeBinomialDistributionと定性的に関係している.GeometricDistributionはPascalDistributionの変換分布(TransformedDistribution)である.ここで,変量 u-1の累積分布関数は,uPascalDistribution[1,p]であるときは常に,GeometricDistribution[p]のそれと厳密に等しい.PascalDistributionはNegativeBinomialDistributionの変換でもある.PascalDistribution[n,p]の確率密度関数は,n→∞のときに,μ と σ についてNormalDistribution[μ,σ]の確率密度関数に収束するのに対し,平均(Mean)と標準偏差(StandardDeviation)はPascalDistribution[n,p]のそれにそれぞれ収束する.PascalDistributionはPoissonDistribution,PoissonConsulDistribution,BinomialDistribution,NegativeBinomialDistribution,MultinomialDistribution,NegativeMultinomialDistributionとも関係がある.
例題
すべて開く すべて閉じる例 (3)
DiscretePlot[Table[PDF[PascalDistribution[10, p], x], {p, {0.4, 0.5, 0.6}}]//Evaluate, {x, 10, 40}, PlotMarkers -> Automatic]DiscretePlot[Table[PDF[PascalDistribution[n, 0.4], x], {n, {5, 8, 10}}]//Evaluate, {x, 5, 40}, PlotMarkers -> Automatic]PDF[PascalDistribution[n, p], x]DiscretePlot[Table[CDF[PascalDistribution[10, p], x], {p, {0.4, 0.5, 0.6}}]//Evaluate, {x, 10, 40}, ExtentSize -> Right]DiscretePlot[Table[CDF[PascalDistribution[n, 0.4], x], {n, {5, 8, 10}}]//Evaluate, {x, 5, 45}, ExtentSize -> Right]CDF[PascalDistribution[n, p], x]Mean[PascalDistribution[n, p]]Variance[PascalDistribution[n, p]]スコープ (8)
data = RandomVariate[PascalDistribution[10, 0.3], 10 ^ 4];Show[
Histogram[data, {2}, "PDF"],
DiscretePlot[PDF[PascalDistribution[10, 0.3], x], {x, 10, 56}, PlotStyle -> PointSize[Medium]]]sample = RandomVariate[PascalDistribution[5, .3], 10 ^ 3];edist = EstimatedDistribution[sample, PascalDistribution[n, p]]サンプルの密度ヒストグラムを推定分布の確率密度関数と比較する:
Show[Histogram[sample, {1}, "PDF"], DiscretePlot[PDF[edist, x], {x, 5, 40}, PlotStyle -> PointSize[Medium]]]Plot[Table[Skewness[PascalDistribution[n, p]], {n, r = {5, 10, 15, 20}}]//Evaluate, {p, 0, 1}, PlotLegends -> (Row[{"n = ", #}]& /@ r)]Skewness[PascalDistribution[n, p]]Limit[Skewness[PascalDistribution[n, p]], n -> ∞]Limit[Skewness[PascalDistribution[n, p]], p -> 0]Limit[Skewness[PascalDistribution[n, p]], p -> 1, Direction -> 1, Assumptions -> n > 0]Plot[Table[Kurtosis[PascalDistribution[n, p]], {n, r = {5, 10, 15, 20}}]//Evaluate, {p, 0, 1}, PlotLegends -> (Row[{"n = ", #}]& /@ r)]Kurtosis[PascalDistribution[n, p]]大きい n の極限の値は標準的なNormalDistributionの尖度である:
Limit[Kurtosis[PascalDistribution[n, p]], n -> ∞]Kurtosis[NormalDistribution[]]Limit[Kurtosis[PascalDistribution[n, p]], p -> 0]Limit[Kurtosis[PascalDistribution[n, p]], p -> 1, Direction -> 1, Assumptions -> n > 0]FormulaGrid[list_, type_] := Grid[...]FormulaGrid[Table[Moment[PascalDistribution[n, p], k]//Factor, {k, 3}], M]FormulaGrid[Table[CentralMoment[PascalDistribution[n, p], k]//Factor, {k, 3}], CM]FormulaGrid[Table[FactorialMoment[PascalDistribution[n, p], k]//Factor, {k, 3}], FM]FactorialMoment[PascalDistribution[n, p], r]FormulaGrid[Table[Cumulant[PascalDistribution[n, p], k]//Factor, {k, 3}], C]DiscretePlot[Table[HazardFunction[PascalDistribution[10, p], x], {p, {0.3, 0.45, 0.6}}]//Evaluate, {x, 10, 30}, PlotMarkers -> Automatic]DiscretePlot[Table[HazardFunction[PascalDistribution[n, 0.4], x], {n, {5, 7, 10}}]//Evaluate, {x, 10, 30}, PlotMarkers -> Automatic]HazardFunction[PascalDistribution[n, p], x]𝒟 = PascalDistribution[2, .3];
points = Table[CDF[𝒟, x], {x, 0, 12}];DiscretePlot[Quantile[𝒟, p], {p, points}, ExtentSize -> Right, ExtentMarkers -> {"Empty", "Filled"}]無次元のQuantityを使ってPascalDistributionを定義する:
PascalDistribution[Quantity[3 + 1 / 2, "Dozens"], Quantity[379, "PerMil"]]アプリケーション (7)
PascalDistributionのCDFは右連続関数の例である:
DiscretePlot[CDF[PascalDistribution[2, .3], x], {x, 0, 10}, ExtentSize -> Right, ExtentMarkers -> {"Filled", "Empty"}]heads3 = PascalDistribution[3, 1 / 2];DiscretePlot[PDF[heads3, x], {x, 3, 10}]Probability[x ≤ 6, xheads3]N[%]Mean[heads3]表が3回出るまでのコイントスの回数のシミュレーションを行う:
ListPlot[RandomVariate[heads3, 20], Filling -> Axis]10回目のコイントスで7回目の表が出た.コインが公正な場合にそのような事象が起こる確率を求める:
fair𝒟 = PascalDistribution[7, 1 / 2];Probability[x == 10, xfair𝒟]N[%]コインが不正であるかもしれないと仮定して,p が最も取りそうな値を求める:
𝒟 = PascalDistribution[7, p];Plot[PDF[𝒟, 10], {p, 0, 1}, Filling -> Axis]FindMaximum[PDF[𝒟, 10], {p, 0.9}]バスケットボールの選手が4回成功するまでフリースローを続ける.この選手がフリースローで得点する確率は0.7である.このフリースローのシミュレーションを行う:
throws = RandomVariate[PascalDistribution[4, 0.7], 50]ListPlot[throws, Filling -> Axis]Mean[PascalDistribution[4, 0.7]]Probability[x == 4, xPascalDistribution[4, 0.7]]ファウルする確率が独立で1分間につき0.1だと仮定して,30分間のファウルのシミュレーションを行う:
RandomVariate[BernoulliDistribution[0.1], 30]ListPlot[%, Filling -> Axis]バスケットボール選手は6回ファウルすると退場になる.ファウルで退場になるまでの期待されるプレー時間を求める:
Mean[PascalDistribution[6, 0.1]]4つのデータパクレットを含むデータストリームがデータ順序に関する情報なしで繰り返し送信されている.全パクレットが正しい順序であるデータストリームが2回目に届くまでの送信回数の分布を求める:
𝒟 = PascalDistribution[2, 1 / 4!];PDF[𝒟, x]DiscretePlot[PDF[𝒟, x], {x, 0, 100}]送信回数が20回以下でパクレットが正しい順序で2回目に受信される確率を求める:
Probability[x ≤ 20, x𝒟]N[%]2回目に正しい順序でデータストリームが届くまでの平均試行回数を求める:
Mean[𝒟]パケットの順序が正しいデータストリームが2回目に届くまでの回数のシミュレーションを行う:
flows = RandomVariate[𝒟, 50]ListPlot[{flows, {{0, Mean[𝒟]}, {50, Mean[𝒟]}}}, Joined -> {False, True}, Filling -> Axis]2つのチャンネルを通してメッセージが送信.成功率は等しい.各チャンネルが少なくとも3つのメッセージを受信までに必要な送信回数が,5回を超えない確率を求める:
Probability[Max[n1, n2] ≤ 5, {n1PascalDistribution[3, p], n2PascalDistribution[3, p]}]Plot[%, {p, 0, 1}]特性と関係 (5)
PascalDistribution[n,p]は
のとき正規分布に収束する:
{μ, σ} = {Mean[PascalDistribution[n, p]], StandardDeviation[PascalDistribution[n, p]]}Block[{p = 1 / 3},
Table[Show[DiscretePlot[PDF[PascalDistribution[n, p], k], {k, Ceiling[μ - 2σ], Floor[μ + 2σ]}, PlotRange -> {{μ - 3σ, μ + 2σ}, All}], Plot[PDF[NormalDistribution[μ, σ], x], {x, μ - 2σ, μ + 2σ}], Ticks -> {{{μ - 2σ, HoldForm[μ - 2σ]}, {μ, HoldForm[μ]}, {μ + 2σ, HoldForm[μ + 2σ]}}, Automatic}, PlotLabel -> n], {n, {5, 15, 25, 50}}]]TransformedDistribution[u + v, {uPascalDistribution[Subscript[n, 1], p], vPascalDistribution[Subscript[n, 2], p]}]
GeometricDistributionはパスカル分布を変換したものである:
𝒟 = TransformedDistribution[u - 1, uPascalDistribution[1, p]];CDF[𝒟, x]CDF[GeometricDistribution[p], x]FullSimplify[% - %%]NegativeBinomialDistributionとパスカル分布の違いはシフトである:
TransformedDistribution[u + n, uNegativeBinomialDistribution[n, p]]テクニカルノート
-
▪
- 離散分布
関連するガイド
-
▪
- 離散一変量分布
テキスト
Wolfram Research (2010), PascalDistribution, Wolfram言語関数, https://reference.wolfram.com/language/ref/PascalDistribution.html (2016年に更新).
CMS
Wolfram Language. 2010. "PascalDistribution." Wolfram Language & System Documentation Center. Wolfram Research. Last Modified 2016. https://reference.wolfram.com/language/ref/PascalDistribution.html.
APA
Wolfram Language. (2010). PascalDistribution. Wolfram Language & System Documentation Center. Retrieved from https://reference.wolfram.com/language/ref/PascalDistribution.html
BibTeX
@misc{reference.wolfram_2026_pascaldistribution, author="Wolfram Research", title="{PascalDistribution}", year="2016", howpublished="\url{https://reference.wolfram.com/language/ref/PascalDistribution.html}", note=[Accessed: 19-July-2026]}
BibLaTeX
@online{reference.wolfram_2026_pascaldistribution, organization={Wolfram Research}, title={PascalDistribution}, year={2016}, url={https://reference.wolfram.com/language/ref/PascalDistribution.html}, note=[Accessed: 19-July-2026]}