lmomco (version 2.4.14)

tttlmomco: Total Time on Test Transform of Distributions

Description

This function computes the Total Time on Test Transform Quantile Function for a quantile function \(x(F)\) (par2qua, qlmomco). The TTT is defined by Nair et al. (2013, p. 171--172, 176) has several expressions $$T(u) = \mu - (1 - u) M(u)\mbox{,}$$ $$T(u) = x(u) - u R(u)\mbox{,}$$ $$T(u) = (1-u) x(u) + \mu L(u)\mbox{,}$$ where \(T(u)\) is the total time on test for nonexceedance probability \(u\), \(M(u)\) is the residual mean quantile function (rmlmomco), \(x(u)\) is a constant for \(x(F = u)\), \(R(u)\) is the reversed mean residual quantile function (rrmlmomco), \(L(u)\) is the Lorenz curve (lrzlmomco), and \(\mu\) as the following definitions $$\mu \equiv \lambda_1(u=0)\mbox{\ first L-moment of residual life for\ }u=0\mbox{,}$$ $$\mu \equiv \lambda_1(x(F))\mbox{\ first L-moment of the quantile function}\mbox{,}$$ $$\mu \equiv \mu(0)\mbox{\ conditional mean for\ }u=0\mbox{.}$$ The definitions imply that within numerical tolerances that \(\mu(0)\) (cmlmomco) should be equal to \(T(1)\), which means that the conditional mean that the 0th percentile in life has been reached equals that total time on test for the 100th percentile. The later can be interpreted as meaning that each of realization of the lifetime distribution for the respective sample size lived to its expected ordered lifetimes.

Usage

tttlmomco(f, para)

Value

Total time on test value for \(F\).

Arguments

f

Nonexceedance probability (\(0 \le F \le 1\)).

para

The parameters from lmom2par or vec2par.

Author

W.H. Asquith

References

Nair, N.U., Sankaran, P.G., and Balakrishnan, N., 2013, Quantile-based reliability analysis: Springer, New York.

See Also

qlmomco, rmlmomco, rrmlmomco, lrzlmomco

Examples

# It is easiest to think about residual life as starting at the origin, units in days.
A <- vec2par(c(0.0, 2649, 2.11), type="gov") # so set lower bounds = 0.0
tttlmomco(0.5, A)  # The median lifetime = 859 days

f <- c(0.25,0.75) # All three computations report: 306.2951 and 1217.1360 days.
Tu1 <- cmlmomco(f=0, A) - (1-f)* rmlmomco(f, A)
Tu2 <-    par2qua(f, A) -    f * rrmlmomco(f, A)
Tu3 <- (1-f)*par2qua(f, A) + cmlmomco(f=0, A)*lrzlmomco(f, A)

if(abs(cmlmomco(0,A) - tttlmomco(1,A)) < 1E-4) {
   print("These two quantities should be nearly identical.\n")
}