OrderDistribution[{dist,n},k]
分布 dist からの n 回の観測の k
次統計分布を表す.
OrderDistribution[{dist,n},{k1,k2,…}]
分布 dist からの n 回の観測の結合![]()
次統計分布を表す.
OrderDistribution[{dist1,…,distn},…]
独立分布 distiの順序統計分布を表す.
OrderDistribution[mdist,…]
多変量分布 mdist の順序統計分布を表す.
OrderDistribution
OrderDistribution[{dist,n},k]
分布 dist からの n 回の観測の k
次統計分布を表す.
OrderDistribution[{dist,n},{k1,k2,…}]
分布 dist からの n 回の観測の結合![]()
次統計分布を表す.
OrderDistribution[{dist1,…,distn},…]
独立分布 distiの順序統計分布を表す.
OrderDistribution[mdist,…]
多変量分布 mdist の順序統計分布を表す.
詳細
- OrderDistribution[{dist,n},k]は親分布 dist から取り出した n 個のサンプルをソートしたリスト中の k 番目の要素の分布を表す.
- 引数分布指定はProductDistributionにおけるものと等しい.これは,多変量分布を表し,すべての変数について連続的または離散的である必要がある.
- OrderDistribution[dists,1]はMin[x1,…,xn]の分布を表す.ただし,{x1,…,xn}は分布ProductDistribution[dists]に従う.
- OrderDistribution[dists,k]はRankedMin[{x1,…,xn},k]の分布を表す.ただし,{x1,…,xn}は分布ProductDistribution[dists]に従う.
- OrderDistributionはMean,CDF,RandomVariate等の関数で使うことができる.
例題
すべて開く すべて閉じる例 (3)
OrderDistribution[{UniformDistribution[], n}, k]OrderDistribution[{ExponentialDistribution[λ], n}, 1]1つの順序統計の順序分布は他の任意の一変量分布と同じように振る舞う:
𝒟 = OrderDistribution[{DiscreteUniformDistribution[{1, 10}], n}, n]PDF[𝒟, x]Block[{n = 3}, DiscretePlot[PDF[𝒟, x], {x, 0, 10}]]2つ以上の順序統計の結合順序分布は他の任意の多変量分布と同じように振る舞う:
𝒟 = OrderDistribution[{UniformDistribution[], n}, {i, j}]PDF[𝒟, {x, y}]Block[{n = 5, i = 2, j = 4}, Plot3D[%, {x, 0, 1}, {y, 0, 1}]]スコープ (35)
基本的な用法 (6)
𝒟min = OrderDistribution[{NormalDistribution[], 10}, 1];
𝒟max = OrderDistribution[{NormalDistribution[], 10}, 10];Plot[{Legended[PDF[𝒟min, x], "𝒟min"], Legended[PDF[𝒟max, x], "𝒟max"]}, {x, -4, 4}, Filling -> Axis]𝒟[k_] = OrderDistribution[{NormalDistribution[], 10}, k];pdf[k_] = PDF[𝒟[k], x]最小変数,7番目に小さい変数,最大変数の確率密度関数をプロットする:
Plot[{Legended[pdf[1], "min"], Legended[pdf[7], "7th"], Legended[pdf[10], "max"]}, {x, -3, 3}, Filling -> Axis, PlotStyle -> {Automatic, Automatic, Automatic}]quantile[k_] = Quantile[𝒟[k], q]Plot[{Legended[quantile[1], "min"], Legended[quantile[7], "k==7"], Legended[quantile[10], "max"]}, {q, 0, 1}, Filling -> Axis, PlotRange -> {-3, 3}, PlotStyle -> {Automatic, Automatic, Automatic}]𝒟 = OrderDistribution[{PoissonDistribution[μ], n}, n]PDF[𝒟, x]Block[{n = 10}, DiscretePlot[Evaluate[Table[Legended[PDF[𝒟, x], Row[{"μ", " = ", μ}]], {μ, {3, 5, 7}}]], {x, 1, 18}, PlotMarkers -> Automatic, Joined -> True, Filling -> None]]𝒟 = OrderDistribution[{ExponentialDistribution[λ], 5}, 3];PDF[𝒟, x]Plot[Evaluate[Table[Legended[PDF[𝒟, x], Row[{"λ", " = ", λ}]], {λ, {0.3, 0.7, 1}}]], {x, 0, 5}, Filling -> Axis]{Mean[𝒟], Variance[𝒟]}𝒟 = OrderDistribution[{BinomialDistribution[40, 1 / 6], 4}, {1, 4}];DiscretePlot3D[PDF[𝒟, {x, y}], {x, 0, 8}, {y, 5, 15}, ExtentSize -> 1 / 2]data = RandomVariate[OrderDistribution[{NormalDistribution[], 10}, 4], 10 ^ 2];EstimatedDistribution[data, OrderDistribution[{NormalDistribution[], n}, 4]]パラメトリック分布 (5)
ExponentialDistributionの k 次統計分布を定義する:
𝒟 = OrderDistribution[{ExponentialDistribution[λ], n}, k];PDF[𝒟, x]CDF[𝒟, x]Block[{n = 10, k = 3, λ = 2}, Table[Plot[fun[𝒟, x], {x, 0, 0.7}, Filling -> Axis, PlotLabel -> fun], {fun, {PDF, CDF, SurvivalFunction, HazardFunction}}]]𝒟 = OrderDistribution[{WeibullDistribution[3, 2], n}, n];PDF[𝒟, x]Plot[Table[Legended[PDF[𝒟, x], Row[{"n", " = ", n}]], {n, {5, 10, 20}}]//Evaluate, {x, 1, 4.5}, Filling -> Axis]𝒟 = OrderDistribution[{BenfordDistribution[b], n}, 1];sample = Block[{n = 4, b = 10}, RandomVariate[𝒟, 10 ^ 4]];Block[{n = 4, b = 10}, Show[Histogram[sample, {0.5, 5.5, 1}, "PDF"], DiscretePlot[PDF[𝒟, x], {x, 0, 10}, PlotStyle -> PointSize[Medium]]]]BinomialDistributionの k
次統計を求める:
𝒟 = OrderDistribution[{BinomialDistribution[100, 1 / 3], 10}, k];PDF[𝒟, x]最大統計,最小統計,4 次統計の確率密度関数をプロットする:
DiscretePlot[Table[Legended[PDF[𝒟, x], Row[{"k", " = ", k}]], {k, {1, 4, 10}}]//Evaluate, {x, 20, 50}, PlotMarkers -> Automatic]Block[{k = 4}, {Mean[𝒟], Variance[𝒟]}]//NBlock[{k = 4}, {Skewness[𝒟], Kurtosis[𝒟]}]//NSkellamDistributionからのサンプルの最大値の分布を求める:
𝒟 = OrderDistribution[{SkellamDistribution[20, 10], 20}, 20];DiscretePlot[PDF[𝒟, x], {x, 10, 30}]points = Table[N[CDF[𝒟, x]], {x, 10, 30, 2}];DiscretePlot[Quantile[𝒟, x], {x, points}, ExtentSize -> Right, ExtentMarkers -> {"Empty", "Filled"}]ノンパラメトリック分布 (3)
HistogramDistributionでサンプルの中の最大値の分布を求める:
hd = HistogramDistribution[RandomVariate[BirnbaumSaundersDistribution[1 / 2, 3], 10 ^ 3]];
𝒟 = OrderDistribution[{hd, 10}, 10];Plot[PDF[𝒟, x], {x, 0.2, 1.5}, PlotRange -> All, Exclusions -> None, Filling -> Axis]SmoothKernelDistributionで k 次統計の分布を求める:
skd = SmoothKernelDistribution[RandomVariate[NormalDistribution[], 10 ^ 3]];
𝒟 = OrderDistribution[{skd, 5}, k];Plot[Table[Legended[PDF[𝒟, x], Row[{"k", " = ", k}]], {k, {1, 3, 5}}]//Evaluate, {x, -3.5, 3.5}, PlotRange -> All, Exclusions -> None, Filling -> Axis]EmpiricalDistributionでサンプル中の最小値の分布を求める:
ed = EmpiricalDistribution[RandomVariate[PoissonDistribution[9], 20]];
𝒟 = OrderDistribution[{ed, 3}, 1];DiscretePlot[CDF[𝒟, x], {x, 0, 12}, ExtentSize -> Right]派生分布 (7)
TruncatedDistributionから2番目に大きい順序統計を求める:
𝒫 = TruncatedDistribution[{-1, 2}, NormalDistribution[]];
𝒟 = OrderDistribution[{𝒫, 10}, 8];CDF[𝒟, x]//FullSimplifyPlot[{Legended[Evaluate[PDF[𝒟, x]], "𝒟"], Legended[PDF[𝒫, x], "𝒫"]}, {x, -1.2, 2}, Filling -> Axis]ParameterMixtureDistributionからのサンプルの最大値の分布を求める:
ℳ = ParameterMixtureDistribution[RayleighDistribution[σ], σHalfNormalDistribution[θ]];
𝒟 = OrderDistribution[{ℳ, 10}, 10];PDF[𝒟, x]//FullSimplify重み分布の母数のさまざまな値についての確率密度関数をプロットする:
Plot[Evaluate[Table[Legended[PDF[𝒟, x], Row[{"θ", " = ", θ}]], {θ, {1 / 2, 1, 1.5}}]], {x, 0, 10}, Filling -> Axis]MixtureDistributionからの9要素のサンプルの中央要素の分布を求める:
ℳ = MixtureDistribution[{4, 3}, {NormalDistribution[], NormalDistribution[2, 1 / 2]}];
𝒟 = OrderDistribution[{ℳ, 9}, 5];Plot[{Legended[Evaluate[PDF[𝒟, x]], "𝒟"], Legended[PDF[ℳ, x], "ℳ"]}, {x, -2, 3.5}, Filling -> Axis]{Mean[ℳ//N], Mean[𝒟//N]}TransformedDistributionからの最小値の分布を求める:
𝒫 = TransformedDistribution[Log[u], uExponentialDistribution[1 / 10]];
𝒟 = OrderDistribution[{𝒫, 10}, 1];Plot[PDF[𝒟, x], {x, -3, 2}, Filling -> Axis]PDF[𝒟, x]CensoredDistributionからの最小値と最大値の結合分布を求める:
𝒞 = CensoredDistribution[{2, 10}, BinomialDistribution[20, 1 / 3]];
𝒟 = OrderDistribution[{𝒞, 4}, {1, 4}];DiscretePlot3D[Evaluate[PDF[𝒟, {x, y}]], {x, 0, 10}, {y, 5, 15}, ExtentSize -> Full]MarginalDistributionからの最大値の分布を求める:
ℳ = MarginalDistribution[ProbabilityDistribution[6( x - y / 2) ^ 2, {x, 0, 1}, {y, 0, 1}], 2];
𝒟 = OrderDistribution[{ℳ, 3}, 3];Plot[Evaluate[PDF[𝒟, x]], {x, 0, 1}, Filling -> Axis, PlotRange -> All]PDF[𝒟, x]QuantityDistributionからの順序統計分布を評価するとQuantityDistributionになる:
OrderDistribution[{QuantityDistribution[ExponentialDistribution[1 / 5], "Seconds"], 10}, 3]Plot[PDF[%, Quantity[x, "Seconds"]]//Evaluate, {x, 0, 10}, AxesLabel -> {"s"}, Filling -> 0]自動簡約 (14)
連続分布 (13)
BetaDistributionはUniformDistributionの順序分布である:
OrderDistribution[{UniformDistribution[{0, 1}], n}, k]DagumDistributionはMaxの下では閉じている:
OrderDistribution[{DagumDistribution[p, a, b], n}, n]ExponentialDistributionはMinの下では閉じている:
OrderDistribution[{ExponentialDistribution[λ], n}, 1]ExtremeValueDistributionはMaxの下では閉じている:
OrderDistribution[{ExtremeValueDistribution[α, β], n}, n]FrechetDistributionはMaxの下では閉じている:
OrderDistribution[{FrechetDistribution[α, β], n}, n]OrderDistribution[{FrechetDistribution[α, β, μ], n}, n]GompertzMakehamDistributionはMinの下では閉じている:
OrderDistribution[{GompertzMakehamDistribution[λ, ξ], n}, 1]OrderDistribution[{GompertzMakehamDistribution[λ, ξ, θ, α], n}, 1]GumbelDistributionはMinの下では閉じている:
OrderDistribution[{GumbelDistribution[α, β], n}, 1]MaxStableDistributionはMaxの下では閉じている:
OrderDistribution[{MaxStableDistribution[μ, σ, ξ], n}, n]OrderDistribution[{MaxStableDistribution[μ, σ, 0], n}, n]MinStableDistributionはMinの下では閉じている:
OrderDistribution[{MinStableDistribution[μ, σ, ξ], n}, 1]OrderDistribution[{MinStableDistribution[μ, σ, 0], n}, 1]ParetoDistributionはMinの下では閉じている:
OrderDistribution[{ParetoDistribution[k, α], n}, 1]PowerDistributionはMaxの下では閉じている:
n = 4;
OrderDistribution[Table[PowerDistribution[p, a[k]], {k, n}], n]SinghMaddalaDistributionはMinの下では閉じている:
OrderDistribution[{SinghMaddalaDistribution[q, a, b], n}, 1]WeibullDistributionはMinの下では閉じている:
OrderDistribution[{WeibullDistribution[α, β], n}, 1]OrderDistribution[{WeibullDistribution[α, β, μ], n}, 1]離散分布 (1)
GeometricDistributionはMinの下では閉じている:
OrderDistribution[{GeometricDistribution[p], n}, 1]一般化と拡張 (2)
Probability[k == 1, kOrderDistribution[MultinomialDistribution[4, {(1/5), (1/5), (1/5), (2/5)}], 1]]Mean[OrderDistribution[{ExponentialDistribution[Subscript[λ, 1]], ExponentialDistribution[Subscript[λ, 2]], ExponentialDistribution[Subscript[λ, 3]]}, 3]]//FullSimplifyアプリケーション (8)
六面のサイコロが4個投げられたとして,最小値の期待値を求める:
Expectation[x, xOrderDistribution[{DiscreteUniformDistribution[{1, 6}], 4}, 1]]N[%]Expectation[x, xOrderDistribution[{DiscreteUniformDistribution[{1, 6}], 4}, 4]]N[%]3つの最大値の合計の期待値を求める.恒等式
とExpectationから得られる線形性を使う:
Expectation[x1 + x2 + x3 + x4, {x1, x2, x3, x4}ProductDistribution[{DiscreteUniformDistribution[{1, 6}], 4}]] - Expectation[x, xOrderDistribution[{DiscreteUniformDistribution[{1, 6}], 4}, 1]]N[%]25人のあまり経験のないマネージャーの中で,最も有能なヘッジファンドのマネージャーが,10年のうち9年間マーケットで他を凌ぐ確率を求める.それぞれのパフォーマンスは互いに独立で,年から年への依存関係もないものとする:
managerSuccessDist = BinomialDistribution[10, 1 / 2];NProbability[goodYears ≥ 9, goodYearsOrderDistribution[{managerSuccessDist, 25}, 25]]でたらめに選ばれたマネージャーが同じ位よい結果を出す確率を求める:
NProbability[goodYears ≥ 9, goodYearsmanagerSuccessDist]連続分布
からのサイズ10のランダムサンプルが昇順にソートされている.新たな確率変量が生成された.11番目のサンプルがソートされたリストの4番目と5番目に小さい値の間に位置する確率を求める:
𝒟 = UniformDistribution[];n = 10;k = 4;Probability[x < z < y, {z𝒟, {x, y}OrderDistribution[{𝒟, n}, {k, k + 1}]}]Table[Probability[x < z < y, {z𝒟, {x, y}OrderDistribution[{𝒟, n}, {k, k + 1}]}], {k, 1, n - 1}]𝒟 = ExponentialDistribution[1];n = 6;Table[Probability[x < z < y, {z𝒟, {x, y}OrderDistribution[{𝒟, n}, {k, k + 1}]}], {k, 1, n - 1}]ExponentialDistributionのサンプルの中央値の期待値を計算する:
mOdd = Expectation[x, xOrderDistribution[{ExponentialDistribution[1], 2n + 1}, n + 1]]mEven = Expectation[(x + y) / 2, {x, y}OrderDistribution[{ExponentialDistribution[1], 2n}, {n, n + 1}]]//FullSimplifySeries[mEven, {n, Infinity, 4}]Series[mOdd, {n, Infinity, 4}]Median[ExponentialDistribution[1]]ExponentialDistributionに従うサイズ
のサンプルの分布範囲を求める:
cdf = Probability[max - min < z, {min, max}OrderDistribution[{ExponentialDistribution[λ], n}, {1, n}], Assumptions -> n > 1]dist[λ_, n_] = ProbabilityDistribution[{"CDF", cdf}, {z, 0, ∞}, Assumptions -> n > 1];pdf = PDF[dist[λ, n], z]Show[
Histogram[Max[#] - Min[#]& /@ RandomVariate[ExponentialDistribution[1], {10 ^ 4, 6}], 25, "PDF"],
Plot[PDF[dist[1, 6], z], {z, 0, 10}, PlotStyle -> Thick]]ErlangDistributionに従う確率変量が毎分新たに生成される.生成された値が先行する任意の値よりも大きい場合に記録されるとする.2番目の記録の値の分布を求める:
𝒟 = ErlangDistribution[k, β];最初の記録は最初の1分に起らなければならない.2番目の記録が
分目に起ったとすると,その確率密度は以下で与えられる:
condpdf[n_, x_] = PDF[OrderDistribution[{𝒟, n}, n], x]2番目の記録が
分目に起る確率は,最初と最後の要素を固定して置換した回数を総置換回数で割ったものに等しい:
pr[n_] := (n - 2)! / n!pdf[k_, β_, x_] = Sum[condpdf[n, x]pr[n], {n, 2, ∞}]Integrate[pdf[2, 1, x], {x, 0, ∞}]2番目の記録の確率密度関数をErlangの確率変量の確率密度関数と比較する:
Plot[{Legended[pdf[2, 1, x], "pdf[x]"], Legended[PDF[𝒟 /. {k -> 2, β -> 1}, x], "PDF[𝒟, x]"]}, {x, 0, 10}, Filling -> Axis]Mean[ProbabilityDistribution[pdf[2, 1, x], {x, 0, ∞}]]StandardDeviation[ProbabilityDistribution[pdf[2, 1, x], {x, 0, ∞}]]サイズ
の標準正規分布に従うサンプル中の最大要素の分布を求める.ただし,
それ自身が
でシフトされたGeometricDistributionに従う乱数であるとする:
𝒟 = ParameterMixtureDistribution[OrderDistribution[{NormalDistribution[], n + 1}, n + 1], nGeometricDistribution[1 / 4]];pdf = PDF[𝒟, x]data = RandomVariate[𝒟, 10 ^ 4];Show[Histogram[data, 25, "PDF"], Plot[pdf, {x, -4, 4}, PlotStyle -> Thick]]NExpectation[x, x𝒟]Mean[data]寿命の分布で定義された3つの同一部品からなるシステムがある.3つの部品のうち2つがダウンするとこのシステムは故障するという.このシステムの寿命分布を求める:
componentLifetime = WeibullDistribution[2, Quantity[3, "Years"]];
systemLifetime = OrderDistribution[{componentLifetime, 3}, 2];このシステムの密度関数といずれかの部品の密度分布を比較する:
Plot[{PDF[componentLifetime, Quantity[x, "Years"]], PDF[systemLifetime, Quantity[x, "Years"]]}, {x, 0, 8}, Filling -> Axis, PlotLegends -> {"Component Lifetime", "System Lifetime"}, AxesLabel -> {"yr"}]Mean[systemLifetime]//NMedian[systemLifetime]N[%]特性と関係 (3)
OrderDistributionはランダムサンプルのRankedMinの分布である:
𝒹 = GammaDistribution[2, 3];
data = RandomVariate[ProductDistribution[{𝒹, 10}], 10 ^ 3];
𝒟 = OrderDistribution[{𝒹, 10}, 2];データのヒストグラムを対応する順序分布の確率密度関数と比較する:
Show[Histogram[RankedMin[#, 2]& /@ data, Automatic, "PDF"], Plot[PDF[𝒟, x], {x, 0, 6}, PlotStyle -> Thick]]OrderDistributionはTransformedDistributionの特殊ケースである:
Table[PDF[OrderDistribution[{NormalDistribution[], 3}, k], w], {k, 3}]Table[PDF[TransformedDistribution[RankedMin[{x, y, z}, k], {x, y, z}ProductDistribution[{NormalDistribution[], 3}]], w], {k, 3}]FullSimplify[%% - %]{RankedMin[{x, y, z}, 1] - Min[x, y, z], RankedMin[{x, y, z}, 3] - Max[x, y, z]}//FullSimplifyExponentialDistributionは
の極限分布である.ただし
はUniformDistributionに従う:
mindist = OrderDistribution[{UniformDistribution[{0, 1}], n}, 1];CDF[mindist, x / n]Limit[CDF[mindist, x / n], n -> ∞]CDF[ExponentialDistribution[1], x]% - %%//FullSimplifyおもしろい例題 (1)
Table[Plot3D[PDF[OrderDistribution[{NormalDistribution[], 3}, k], {x, y}], {x, -3, 3}, {y, -3, 3}, Mesh -> None, PlotRange -> {0, 0.6}, PlotLabel -> k], {k, Subsets[Range[3], {2}]}]Table[Plot3D[PDF[OrderDistribution[{NormalDistribution[], 4}, k], {x, y}], {x, -3, 3}, {y, -3, 3}, Mesh -> None, PlotRange -> {0, 1}, PlotLabel -> k], {k, Subsets[Range[4], {2}]}]テキスト
Wolfram Research (2010), OrderDistribution, Wolfram言語関数, https://reference.wolfram.com/language/ref/OrderDistribution.html (2016年に更新).
CMS
Wolfram Language. 2010. "OrderDistribution." Wolfram Language & System Documentation Center. Wolfram Research. Last Modified 2016. https://reference.wolfram.com/language/ref/OrderDistribution.html.
APA
Wolfram Language. (2010). OrderDistribution. Wolfram Language & System Documentation Center. Retrieved from https://reference.wolfram.com/language/ref/OrderDistribution.html
BibTeX
@misc{reference.wolfram_2026_orderdistribution, author="Wolfram Research", title="{OrderDistribution}", year="2016", howpublished="\url{https://reference.wolfram.com/language/ref/OrderDistribution.html}", note=[Accessed: 04-November-2026]}
BibLaTeX
@online{reference.wolfram_2026_orderdistribution, organization={Wolfram Research}, title={OrderDistribution}, year={2016}, url={https://reference.wolfram.com/language/ref/OrderDistribution.html}, note=[Accessed: 04-November-2026]}