Shifted Gompertz distribution

Shifted Gompertz

Probability density function
Cumulative distribution function
Parameters	[math]\displaystyle{ b \geq 0 }[/math] scale (real) [math]\displaystyle{ \eta\geq 0 }[/math] shape (real)
Support	[math]\displaystyle{ x \in [0, \infty)\! }[/math]
PDF	[math]\displaystyle{ b e^{-bx} e^{-\eta e^{-bx}}\left[1 + \eta\left(1 - e^{-bx}\right)\right] }[/math]
CDF	[math]\displaystyle{ \left(1 - e^{-bx}\right)e^{-\eta e^{-bx}} }[/math]
Mean	[math]\displaystyle{ (-1/b)\{\mathrm{E}[\ln(X)] - \ln(\eta)\}\, }[/math] where [math]\displaystyle{ X = \eta e^{-bx}\, }[/math] and [math]\displaystyle{ \begin{align}\mathrm{E}[\ln(X)] =& [1 {+} 1 / \eta]\!\!\int_0^\eta \!\!\!\! e^{-X}[\ln(X)]dX\\ &- 1/\eta\!\! \int_0^\eta \!\!\!\! X e^{-X}[\ln(X)] dX \end{align} }[/math]
Mode	[math]\displaystyle{ 0 \text{ for }0 \lt \eta \leq 0.5 }[/math] [math]\displaystyle{ (-1/b)\ln(z^\star)\text{, for } \eta \gt 0.5 }[/math] [math]\displaystyle{ \text{ where }z^\star = [3 + \eta - (\eta^2 + 2\eta + 5)^{1/2}]/(2\eta) }[/math]
Variance	[math]\displaystyle{ (1/b^2)(\mathrm{E}\{[\ln(X)]^2\} - (\mathrm{E}[\ln(X)])^2)\, }[/math] where [math]\displaystyle{ X = \eta e^{-bx}\, }[/math] and [math]\displaystyle{ \begin{align}\mathrm{E}\{[\ln(X)]^2\} =& [1 {+} 1 / \eta]\!\!\int_0^\eta \!\!\!\! e^{-X}[\ln(X)]^2 dX\\ &- 1/\eta \!\!\int_0^\eta \!\!\!\! X e^{-X}[\ln(X)]^2 dX \end{align} }[/math]

The shifted Gompertz distribution is the distribution of the larger of two independent random variables one of which has an exponential distribution with parameter [math]\displaystyle{ b }[/math] and the other has a Gumbel distribution with parameters [math]\displaystyle{ \eta }[/math] and [math]\displaystyle{ b }[/math]. In its original formulation the distribution was expressed referring to the Gompertz distribution instead of the Gumbel distribution but, since the Gompertz distribution is a reverted Gumbel distribution, the labelling can be considered as accurate. It has been used as a model of adoption of innovations. It was proposed by Bemmaor (1994). Some of its statistical properties have been studied further by Jiménez and Jodrá (2009) and Jiménez Torres (2014).

It has been used to predict the growth and decline of social networks and on-line services and shown to be superior to the Bass model and Weibull distribution (Bauckhage and Kersting 2014).

Specification

Probability density function

The probability density function of the shifted Gompertz distribution is:

[math]\displaystyle{ f(x;b,\eta) = b e^{-bx} e^{-\eta e^{-bx}}\left[1 + \eta\left(1 - e^{-bx}\right)\right] \text{ for }x \geq 0. \, }[/math]

where [math]\displaystyle{ b \geq 0 }[/math] is a scale parameter and [math]\displaystyle{ \eta \geq 0 }[/math] is a shape parameter. In the context of diffusion of innovations, [math]\displaystyle{ b }[/math] can be interpreted as the overall appeal of the innovation and [math]\displaystyle{ \eta }[/math] is the propensity to adopt in the propensity-to-adopt paradigm. The larger [math]\displaystyle{ b }[/math] is, the stronger the appeal and the larger [math]\displaystyle{ \eta }[/math] is, the smaller the propensity to adopt.

The distribution can be reparametrized according to the external versus internal influence paradigm with [math]\displaystyle{ p = f(0;b,\eta) = be^{-\eta } }[/math] as the coefficient of external influence and [math]\displaystyle{ q = b - p }[/math] as the coefficient of internal influence. Hence:

[math]\displaystyle{ f(x;p,q) = (p + q) e^{-(p + q)x} e^{-\ln(1 + q/p) e^{-(p+q)x}}\left[1 + \ln(1 + q/p)\left(1 - e^{-(p + q)x}\right)\right] \text{ for }x \geq 0, p, q \geq 0. \, }[/math]

[math]\displaystyle{ = (p + q) e^{-(p + q)x} {(1 + q/p)^{-e^{-(p+q)x}}}\left[1 + \ln(1 + q/p)\left(1 - e^{-(p + q)x}\right)\right] \text{ for }x \geq 0, p, q \geq 0. \, }[/math]

When [math]\displaystyle{ q = 0 }[/math], the shifted Gompertz distribution reduces to an exponential distribution. When [math]\displaystyle{ p = 0 }[/math], the proportion of adopters is nil: the innovation is a complete failure. The shape parameter of the probability density function is equal to [math]\displaystyle{ q/p }[/math]. Similar to the Bass model, the hazard rate [math]\displaystyle{ z(x;p,q) }[/math] is equal to [math]\displaystyle{ p }[/math] when [math]\displaystyle{ x }[/math] is equal to [math]\displaystyle{ 0 }[/math]; it approaches [math]\displaystyle{ p + q }[/math] as [math]\displaystyle{ x }[/math] gets close to [math]\displaystyle{ \infty }[/math]. See Bemmaor and Zheng for further analysis.

Cumulative distribution function

The cumulative distribution function of the shifted Gompertz distribution is:

[math]\displaystyle{ F(x;b,\eta) = \left(1 - e^{-bx}\right)e^{-\eta e^{-bx}} \text{ for }x \geq 0. \, }[/math]

Equivalently,

[math]\displaystyle{ F(x;p, q) = \left(1 - e^{-(p + q)x}\right)e^{-\ln(1 + q/p)e^{-(p+q)x}} \text{ for }x \geq 0. \, }[/math]

[math]\displaystyle{ = \left(1 - e^{-(p + q)x}\right){(1 + q/p)^{-e^{-(p+q)x}}} \text{ for }x \geq 0. \, }[/math]

Properties

The shifted Gompertz distribution is right-skewed for all values of [math]\displaystyle{ \eta }[/math]. It is more flexible than the Gumbel distribution. The hazard rate is a concave function of [math]\displaystyle{ F(x;b,\eta) }[/math] which increases from [math]\displaystyle{ be^{-\eta} }[/math] to [math]\displaystyle{ b }[/math]: its curvature is all the steeper as [math]\displaystyle{ \eta }[/math] is large. In the context of the diffusion of innovations, the effect of word of mouth (i.e., the previous adopters) on the likelihood to adopt decreases as the proportion of adopters increases. (For comparison, in the Bass model, the effect remains the same over time). The parameter [math]\displaystyle{ q = b(1-e^{-\eta}) }[/math] captures the growth of the hazard rate when [math]\displaystyle{ x }[/math] varies from [math]\displaystyle{ 0 }[/math] to [math]\displaystyle{ \infty }[/math].

Shapes

The shifted Gompertz density function can take on different shapes depending on the values of the shape parameter [math]\displaystyle{ \eta }[/math]:

• [math]\displaystyle{ 0 \lt \eta \leq 0.5\, }[/math] the probability density function has its mode at 0.
• [math]\displaystyle{ \eta \gt 0.5\, }[/math] the probability density function has its mode at

[math]\displaystyle{ \text{mode}=-\frac{\ln(z^\star)}{b}\, \qquad 0 \lt z^\star \lt 1 }[/math]

where [math]\displaystyle{ z^\star\, }[/math] is the smallest root of

[math]\displaystyle{ \eta^2z^2 - \eta(3 + \eta)z + \eta + 1 = 0\,, }[/math]

which is

[math]\displaystyle{ z^\star = [3 + \eta - (\eta^2 + 2\eta + 5)^{1/2}]/(2\eta). }[/math]

Related distributions

When [math]\displaystyle{ \eta }[/math] varies according to a gamma distribution with shape parameter [math]\displaystyle{ \alpha }[/math] and scale parameter [math]\displaystyle{ \beta }[/math] (mean = [math]\displaystyle{ \alpha\beta }[/math]), the distribution of [math]\displaystyle{ x }[/math] is Gamma/Shifted Gompertz (G/SG). When [math]\displaystyle{ \alpha }[/math] is equal to one, the G/SG reduces to the Bass model (Bemmaor 1994). The three-parameter G/SG has been applied by Dover, Goldenberg and Shapira (2009) and Van den Bulte and Stremersch (2004) among others in the context of the diffusion of innovations. The model is discussed in Chandrasekaran and Tellis (2007).Similar to the shifted Gompertz distribution, the G/SG can either be represented according to the propensity-to-adopt paradigm or according to the innovation-imitation paradigm. In the latter case, it includes three parameters: [math]\displaystyle{ p, q }[/math] and [math]\displaystyle{ \alpha }[/math] with [math]\displaystyle{ p = f(0;b,\beta, \alpha) = b/(1+\beta)^{\alpha } }[/math] and [math]\displaystyle{ q = b - p }[/math]. The parameter [math]\displaystyle{ \alpha }[/math] modifies the curvature of the hazard rate as expressed as a function of [math]\displaystyle{ F(x;p,q, \alpha) }[/math]: when [math]\displaystyle{ \alpha }[/math] is less than 0.5, it decreases to a minimum prior to increasing at an increasing rate as [math]\displaystyle{ F(x;p,q, \alpha \lt 1/2) }[/math] increases, it is convex when [math]\displaystyle{ \alpha }[/math] is less than one and larger or equal to 0.5, linear when [math]\displaystyle{ \alpha }[/math] is equal to one, and concave when [math]\displaystyle{ \alpha }[/math] is larger than one. Here are some special cases of the G/SG distribution in the case of homogeneity (across the population) with respect to the likelihood to adopt at a given time:

                        [math]\displaystyle{ F(x;p,q, \alpha = 0) }[/math] = Exponential[math]\displaystyle{ (p + q) }[/math]
                        [math]\displaystyle{ F(x;p,q, \alpha = 1/2) }[/math] = Left-skewed two-parameter distribution[math]\displaystyle{ (p,q) }[/math]
                        [math]\displaystyle{ F(x;p,q, \alpha = 1) }[/math]  = Bass model[math]\displaystyle{ (p,q) }[/math]
                        [math]\displaystyle{ F(x;p,q, \alpha = \infty) }[/math] = Shifted Gompertz[math]\displaystyle{ (p,q) }[/math]
 

with:

              [math]\displaystyle{  F(x;p, q,\alpha = 1/2) = \left(1 - e^{-(p + q)x}\right)/{(1 + (q/p)(2+q/p)e^{-(p+q)x})^{1/2}} \text{ for }x \geq 0,p, q \geq 0. \, }[/math]

One can compare the parameters [math]\displaystyle{ p }[/math] and [math]\displaystyle{ q }[/math] across the values of [math]\displaystyle{ \alpha }[/math] as they capture the same notions. In all the cases, the hazard rate is either constant or a monotonically increasing function of [math]\displaystyle{ F(x;p,q, \alpha) }[/math] (positive word of mouth). As the diffusion curve is all the more skewed as [math]\displaystyle{ \alpha }[/math] becomes large, we expect [math]\displaystyle{ q }[/math] to decrease as the level of right-skew increases.

See also

• Gumbel distribution
• Generalized extreme value distribution
• Mixture model
• Bass model
• Gompertz distribution