Search Results for “quantile parameterized distribution”

Source: Quantile-parameterized distribution

A quantile-parameterized distribution (QPD) is a probability distributions that is directly parameterized by data. They were created to meet the need for easy-to-use continuous probability distributions flexible enough to represent a wide range of uncertainties, such as those commonly encountered in business, economics, engineering, and science. Because QPDs are directly parameterized by data, they have the practical advantage of avoiding the intermediate step of parameter estimation, a time-consuming process that typically requires non-linear iterative methods to estimate probability-distribution parameters from data. Some QPDs have virtually unlimited shape flexibility and closed-form moments as well.

History

quantile

Definition

where

L

0

=

lim

y
→

0

+

F

−
1

(
y
)

L

1

=

lim

y
→

1

−

F

−
1

(
y
)

{\displaystyle {\begin{array}{rcl}L_{0}&=&\lim _{y\rightarrow 0^{+}}F^{-1}(y)\\L_{1}&=&\lim _{y\rightarrow 1^{-}}F^{-1}(y)\end{array}}}

and the functions

g

i

(
y
)

{\displaystyle g_{i}(y)}

are continuously differentiable and linearly independent basis functions. Here, essentially,

L

0

{\displaystyle L_{0}}

and

L

1

{\displaystyle L_{1}}

are the lower and upper bounds (if they exist) of a random variable with quantile function

F

−
1

(
y
)

{\displaystyle F^{-1}(y)}

. These distributions are called quantile-parameterized because for a given set of quantile pairs

{
(

x

i

,

y

i

)
∣
i
=
1
,
…
,
n
}

{\displaystyle \{(x_{i},y_{i})\mid i=1,\ldots ,n\}}

, where

x

i

=

F

−
1

(

y

i

)

{\displaystyle x_{i}=F^{-1}(y_{i})}

, and a set of

n

{\displaystyle n}

basis functions

g

i

(
y
)

{\displaystyle g_{i}(y)}

, the coefficients

a

i

{\displaystyle a_{i}}

can be determined by solving a set of linear equations. If one desires to use more quantile pairs than basis functions, then the coefficients

a

i

{\displaystyle a_{i}}

can be chosen to minimize the sum of squared errors between the stated quantiles

x

i

{\displaystyle x_{i}}

and

F

−
1

(

y

i

)

{\displaystyle F^{-1}(y_{i})}

. Keelin and Powley illustrate this concept for a specific choice of basis functions that is a generalization of quantile function of the normal distribution,

x
=
μ
+
σ

Φ

−
1

(
y
)

{\displaystyle x=\mu +\sigma \Phi ^{-1}(y)}

, for which the mean

μ

{\displaystyle \mu }

and standard deviation

σ

{\displaystyle \sigma }

are linear functions of cumulative probability

y

{\displaystyle y}

:

μ
(
y
)
=

a

1

+

a

4

y

{\displaystyle \mu (y)=a_{1}+a_{4}y}

σ
(
y
)
=

a

2

+

a

3

y

{\displaystyle \sigma (y)=a_{2}+a_{3}y}

The result is a four-parameter distribution that can be fit to a set of four quantile/probability pairs exactly, or to any number of such pairs by linear least squares. Keelin and Powley call this the Simple Q-Normal distribution. Some skewed and symmetric Simple Q-Normal PDFs are shown in the figures below.

Properties

QPD’s that meet Keelin and Powley’s definition have the following properties.

= Probability density function

=
Differentiating

x
=

F

−
1

(
y
)
=

∑

i
=
1

n

a

i

g

i

(
y
)

{\displaystyle x=F^{-1}(y)=\sum _{i=1}^{n}a_{i}g_{i}(y)}

with respect to

y

{\displaystyle y}

yields

d
x

/

d
y

{\displaystyle dx/dy}

. The reciprocal of this quantity,

d
y

/

d
x

{\displaystyle dy/dx}

, is the probability density function (PDF)

f
(
y
)
=

(

∑

i
=
1

n

a

i

d

g

i

(
y
)

d
y

)

−
1

{\displaystyle f(y)=\left(\sum _{i=1}^{n}a_{i}{{dg_{i}(y)} \over {dy}}\right)^{-1}}

where

0
<
y
<
1

{\displaystyle 0
. Note that this PDF is expressed as a function of cumulative probability

y

{\displaystyle y}

rather than

x

{\displaystyle x}

. To plot it, as shown in the figures, vary

y
∈
(
0
,
1
)

{\displaystyle y\in (0,1)}

parametrically. Plot

x
=

F

−
1

(
y
)

{\displaystyle x=F^{-1}(y)}

on the horizontal axis and

f
(
y
)

{\displaystyle f(y)}

on the vertical axis.

= Feasibility

=
A function of the form of

F

−
1

(
y
)

{\displaystyle F^{-1}(y)}

is a feasible probability distribution if and only if

f
(
y
)
>
0

{\displaystyle f(y)>0}

for all

y
∈
(
0
,
1
)

{\displaystyle y\in (0,1)}

. This implies a feasibility constraint on the set of coefficients

a

=
(

a

1

,
…
,

a

n

)
∈

R

n

{\displaystyle {\boldsymbol {a}}=(a_{1},\ldots ,a_{n})\in \mathbb {R} ^{n}}

:

∑

i
=
1

n

a

i

d

g

i

(
y
)

d
y

>
0

{\displaystyle \sum _{i=1}^{n}a_{i}{{dg_{i}(y)} \over {dy}}>0}

for all

y
∈
(
0
,
1
)

{\displaystyle y\in (0,1)}

In practical applications, feasibility must generally be checked rather than assumed.

= Convexity

=
A QPD’s set of feasible coefficients

S

a

=
{

a

∈

R

n

∣

∑

i
=
1

n

a

i

d

g

i

(
y
)

/

d
y
>
0

{\displaystyle S_{\boldsymbol {a}}=\{{\boldsymbol {a}}\in \mathbb {R} ^{n}\mid \sum _{i=1}^{n}a_{i}dg_{i}(y)/dy>0}

for all

y
∈
(
0
,
1
)
}

{\displaystyle y\in (0,1)\}}

is convex. Because convex optimization requires convex feasible sets, this property simplifies optimization applications involving QPDs.

= Fitting to data

=
The coefficients

a

{\displaystyle {\boldsymbol {a}}}

can be determined from data by linear least squares. Given

m

{\displaystyle m}

data points

(

x

i

,

y

i

)

{\displaystyle (x_{i},y_{i})}

that are intended to characterize the CDF of a QPD, and

m
×
n

{\displaystyle m\times n}

matrix

Y

{\displaystyle {\boldsymbol {Y}}}

whose elements consist of

g

j

(

y

i

)

{\displaystyle g_{j}(y_{i})}

, then, so long as

Y

T

Y

{\displaystyle {\boldsymbol {Y}}^{T}{\boldsymbol {Y}}}

is invertible, coefficients' column vector

a

{\displaystyle {\boldsymbol {a}}}

can be determined as

a

=
(

Y

T

Y

)

−
1

Y

T

x

{\displaystyle {\boldsymbol {a}}=({\boldsymbol {Y}}^{T}{\boldsymbol {Y}})^{-1}{\boldsymbol {Y}}^{T}{\boldsymbol {x}}}

, where

m
≥
n

{\displaystyle m\geq n}

and column vector

x

=
(

x

1

,
…
,

x

m

)

{\displaystyle {\boldsymbol {x}}=(x_{1},\ldots ,x_{m})}

. If

m
=
n

{\displaystyle m=n}

, this equation reduces to

a

=

Y

−
1

x

{\displaystyle {\boldsymbol {a}}={\boldsymbol {Y}}^{-1}{\boldsymbol {x}}}

, where the resulting CDF runs through all data points exactly. An alternate method, implemented as a linear program, determines the coefficients by minimizing the sum of absolute distances between the CDF and the data subject to feasibility constraints.

= Shape flexibility

=
A QPD with

n

{\displaystyle n}

terms, where

n
≥
2

{\displaystyle n\geq 2}

, has

n
−
2

{\displaystyle n-2}

shape parameters. Thus, QPDs can be far more flexible than the Pearson distributions, which have at most two shape parameters. For example, ten-term metalog distributions parameterized by 105 CDF points from 30 traditional source distributions (including normal, student-t, lognormal, gamma, beta, and extreme value) have been shown to approximate each such source distribution within a K–S distance of 0.001 or less.

= Transformations

=
QPD transformations are governed by a general property of quantile functions: for any quantile function

x
=
Q
(
y
)

{\displaystyle x=Q(y)}

and increasing function

t
(
x
)
,
x
=

t

−
1

(
Q
(
y
)
)

{\displaystyle t(x),x=t^{-1}(Q(y))}

is a quantile function. For example, the quantile function of the normal distribution,

x
=
μ
+
σ

Φ

−
1

(
y
)

{\displaystyle x=\mu +\sigma \Phi ^{-1}(y)}

, is a QPD by the Keelin and Powley definition. The natural logarithm,

t
(
x
)
=
ln
⁡
(
x
−

b

l

)

{\displaystyle t(x)=\ln(x-b_{l})}

, is an increasing function, so

x
=

b

l

+

e

μ
+
σ

Φ

−
1

(
y
)

{\displaystyle x=b_{l}+e^{\mu +\sigma \Phi ^{-1}(y)}}

is the quantile function of the lognormal distribution with lower bound

b

l

{\displaystyle b_{l}}

. Importantly, this transformation converts an unbounded QPD into a semi-bounded QPD. Similarly, applying this log transformation to the unbounded metalog distribution yields the semi-bounded (log) metalog distribution; likewise, applying the logit transformation,

t
(
x
)
=
ln
⁡
(
(
x
−

b

l

)

/

(

b

u

−
x
)
)

{\displaystyle t(x)=\ln((x-b_{l})/(b_{u}-x))}

, yields the bounded (logit) metalog distribution with lower and upper bounds

b

l

{\displaystyle b_{l}}

and

b

u

{\displaystyle b_{u}}

, respectively. Moreover, by considering

t
(
x
)

{\displaystyle t(x)}

to be

F

−
1

(
y
)

{\displaystyle F^{-1}(y)}

distributed, where

F

−
1

(
y
)

{\displaystyle F^{-1}(y)}

is any QPD that meets Keelin and Powley’s definition, the transformed variable maintains the above properties of feasibility, convexity, and fitting to data. Such transformed QPDs have greater shape flexibility than the underlying

F

−
1

(
y
)

{\displaystyle F^{-1}(y)}

, which has

n
−
2

{\displaystyle n-2}

shape parameters; the log transformation has

n
−
1

{\displaystyle n-1}

shape parameters, and the logit transformation has

n

{\displaystyle n}

shape parameters. Moreover, such transformed QPDs share the same set of feasible coefficients as the underlying untransformed QPD.

= Moments

=
The

k

t
h

{\displaystyle k^{th}}

moment of a QPD is:

E
[

x

k

]
=

∫

0

1

(

∑

i
=
1

n

a

i

g

i

(
y
)

)

k

d
y

{\displaystyle E[x^{k}]=\int _{0}^{1}\left(\sum _{i=1}^{n}a_{i}g_{i}(y)\right)^{k}dy}

Whether such moments exist in closed form depends on the choice of QPD basis functions

g

i

(
y
)

{\displaystyle g_{i}(y)}

. The unbounded metalog distribution and polynomial QPDs are examples of QPDs for which moments exist in closed form as functions of the coefficients

a

i

{\displaystyle a_{i}}

.

= Simulation

=
Since the quantile function

x
=

F

−
1

(
y
)

{\displaystyle x=F^{-1}(y)}

is expressed in closed form, Keelin and Powley QPDs facilitate Monte Carlo simulation. Substituting in uniformly distributed random samples of

y

{\displaystyle y}

produces random samples of

x

{\displaystyle x}

in closed form, thereby eliminating the need to invert a CDF expressed as

y
=
F
(
x
)

{\displaystyle y=F(x)}

.

Related distributions

The following probability distributions are QPDs according to Keelin and Powley’s definition:

The quantile function of the normal distribution,

x
=
μ
+
σ

Φ

−
1

(
y
)

{\displaystyle x=\mu +\sigma \Phi ^{-1}(y)}

.
The quantile function of the Gumbel distribution,

x
=
μ
−
β
ln
⁡
(
−
ln
⁡
(
y
)
)

{\displaystyle x=\mu -\beta \ln(-\ln(y))}

.
The quantile function of the Cauchy distribution,

x
=

x

0

+
γ
tan
⁡
[
π
(
y
−
0.5
)
]

{\displaystyle x=x_{0}+\gamma \tan[\pi (y-0.5)]}

.
The quantile function of the logistic distribution,

x
=
μ
+
s
ln
⁡
(
y

/

(
1
−
y
)
)

{\displaystyle x=\mu +s\ln(y/(1-y))}

.
The unbounded metalog distribution, which is a power series expansion of the

μ

{\displaystyle \mu }

and

s

{\displaystyle s}

parameters of the logistic quantile function.
The semi-bounded and bounded metalog distributions, which are the log and logit transforms, respectively, of the unbounded metalog distribution.
The SPT (symmetric-percentile triplet) unbounded, semi-bounded, and bounded metalog distributions, which are parameterized by three CDF points and optional upper and lower bounds.
The Simple Q-Normal distribution
The metadistributions, including the meta-normal
Quantile functions expressed as polynomial functions of cumulative probability

y

{\displaystyle y}

, including Chebyshev polynomial functions.
Like the SPT metalog distributions, the Johnson Quantile-Parameterized Distributions (JQPDs) are parameterized by three quantiles. JQPDs do not meet Keelin and Powley’s QPD definition, but rather have their own properties. JQPDs are feasible for all SPT parameter sets that are consistent with the rules of probability.

Applications

The original applications of QPDs were by decision analysts wishing to conveniently convert expert-assessed quantiles (e.g., 10th, 50th, and 90th quantiles) into smooth continuous probability distributions. QPDs have also been used to fit output data from simulations in order to represent those outputs (both CDFs and PDFs) as closed-form continuous distributions. Used in this way, they are typically more stable and smoother than histograms. Similarly, since QPDs can impose fewer shape constraints than traditional distributions, they have been used to fit a wide range of empirical data in order to represent those data sets as continuous distributions (e.g., reflecting bimodality that may exist in the data in a straightforward manner). Quantile parameterization enables a closed-form QPD representation of known distributions whose CDFs otherwise have no closed-form expression. Keelin et al. (2019) apply this to the sum of independent identically distributed lognormal distributions, where quantiles of the sum can be determined by a large number of simulations. Nine such quantiles are used to parameterize a semi-bounded metalog distribution that runs through each of these nine quantiles exactly. QPDs have also been applied to assess the risks of asteroid impact, cybersecurity, biases in projections of oil-field production when compared to observed production after the fact, and future Canadian population projections based on combining the probabilistic views of multiple experts. See metalog distributions and Keelin (2016) for additional applications of the metalog distribution.

External links

The Metalog Distributions, www.metalogs.org

History

Definition

Properties

= Probability density function

= Feasibility

= Convexity

= Fitting to data

= Shape flexibility

= Transformations

= Moments

= Simulation

Related distributions

Applications

External links

References

Kata Kunci Pencarian:

Recent Movies

Recent Movies

Categories

Recent Movies