Weakest-Link Scaling and Extreme Events in Finite-Sized Systems

Weakest-link scaling is used in the reliability analysis of complex systems. It is characterized by the extensivity of the hazard function instead of the entropy. The Weibull distribution is the archetypical example of weakest-link scaling, and it describes variables such as the fracture strength of brittle materials, maximal annual rainfall, wind speed and earthquake return times. We investigate two new distributions that exhibit weakest-link scaling, i.e., a Weibull generalization known as the κ-Weibull and a modified gamma probability function that we propose herein. We show that in contrast with the Weibull and the modified gamma, the hazard function of the κ-Weibull is non-extensive, which is a signature of inter-dependence between the links. We also investigate the impact of heterogeneous links, modeled by means of a stochastic Weibull scale parameter, on the observed probability distribution.


Introduction
In statistical mechanics, the velocities of non-interacting particles are viewed as random variables that follow the Maxwell-Boltzmann distribution.Quantities of interest are often macroscopic (system) properties, such as the energy or the entropy, that involve averages over many particles.Evaluating these properties involves calculating averages over suitable probability distributions.Randomness is also a factor in the analysis of macroscopic variables.Environmental time series, energy and mineral reserves and earthquakes are just a few examples of phenomena that exhibit stochastic fluctuations.In addition, the lifetimes (times to failure) of various natural and engineered materials and devices are also random variables.In contrast with statistical mechanics, there is a strong interest not only in average quantities, but also in the extreme values of such distributions [1].The classical extreme value theory studies the probability distributions of extreme events [2,3].The Fisher-Tippet-Gnedenko theorem shows that the distribution of suitably-scaled transformations of the maxima or minima of a collection of independent and identically distributed (i.i.d.) variables follows the generalized extreme value (GEV) distributions, which include the Gumbel, Weibull and Fréchet distributions.The type of GEV obtained depends on the tail behavior of the probability distribution for the i.i.d.random variables.Extensions that involve stochastic scalings have been recently proposed in [4].
We investigate systems that obey the weakest-link scaling (WLS) principle.In the WLS framework, the system response (e.g., material strength) is controlled by the weakest link (unit), which implies a marked deviation from statistical independence.The Weibull model [3] is the archetypical probability distribution for WLS systems and is widely used in reliability modeling [5].It finds application as a model of the fracture strength of brittle and quasi-brittle technological materials [6,7] and geologic media [8], waiting times between earthquakes [9][10][11][12][13][14][15], wind speed [16] and annual hydrological maxima [17].We study the impact of the randomness (heterogeneity) of the Weibull link parameters on the system survival function by introducing doubly-stochastic representations analogous to superstatistics [18].We also investigate weakest-link scaling with non-exponential link survival functions.In this context, we introduce a modified gamma distribution, and we provide theoretical arguments for the κ-Weibull distribution, which exhibits power-law tails in finite-sized systems.

Mathematical Preliminaries
We will use the capital letter X to denote a random variable (e.g., tensile strength, wind speed, earthquake return time, maximum annual rainfall or river discharge) and the lowercase x to denote values of the random variable, where x ∈ [0, ∞).We will assume that the random variable represents the property of interest for a particular spatial or temporal window (support), the size of which we will denote by .The probability density function (pdf) will be denoted by f (x).It is customary to use f X (x) to denote the pdf of X when more than one random variable is involved.In this case, however, there is no risk of confusion, and thus, we replace X with the support scale, which is relevant in our analysis.The integral of the pdf, i.e., the cumulative distribution function (cdf), will be denoted by: Another useful probability function, which is mostly used in reliability analysis, is the survival function or complementary cumulative distribution function, defined by: R (x) = 1 − F (x). (1) The survival function R (x) gives the probability that the system has not "failed" when the variable X has reached the level x.The survival function is equal to one for x = x min and tends to zero for x → ∞.
A characteristic property of the survival function is that its integral equals the expectation of X, i.e., The above is based on dR (x)/dx = −f (x), which follows from (1) and integration by parts.
• In the case of fracture strength, R (x) is the probability that the system has not failed if the external load takes values that do not exceed x.
• In the case of earthquake return times, the survival function is the probability that there are no successive earthquakes separated by a time interval less than or equal to x.
• In the case of annual rainfall maxima, R (x) gives the probability that the maximum rainfall over the support area for the duration of one year does not exceed the value x.
Finally, the hazard rate or hazard function h (x) is the conditional probability that the system will fail at time x * > x within the infinitesimal time window x < x * ≤ x + δx, given that the system has not failed during the interval [0, x].Let A denote the event that the system has not failed at x and B the event that the system fails at x < x * ≤ x + δx.Then, h (x) = Prob(B|A) = Prob(B ∩ A)/Prob(A), which leads to the following: The asymptotic dependence of the hazard rate is a matter of interest.For example, in the case of earthquakes, it is often desired to estimate the probability of the time until the next big earthquake.If the particular hazard rate tends to infinity as x → ∞, the probability that an earthquake will occur increases with the time elapsed since the last earthquake [19].In contrast, a decreasing or constant hazard rate indicates that the probability of a new occurrence is reduced or remains constant, respectively, as time elapses, and no event takes place.

Weakest-Link Scaling Theory
The weakest-link scaling (WLS) theory has been developed to model the strength of disordered systems based on the strengths of "critical elements" or links [20].In WLS, the entire system fails if the link with the lowest mechanical strength breaks, hence the term "weakest-link".The foundations of weakest-link theory are found in the studies of Gumbel [2] and Weibull [3] on the statistics of extremes.The concept of links is straightforward in simple systems, e.g., one-dimensional chains.In higher dimensions, the links correspond to critical subsystems, the failure of which causes an avalanche, leading to the failure of the entire system.
Weakest-link scaling is used in statistical physics to investigate size effects in the strength of heterogeneous materials [6,[21][22][23].In reliability modeling, a weakest-link system is known as a series system [24].The same scaling ideas, however, can be used in other systems, as well.For example, in the case of maximum annual rainfall, the scale links could correspond to single days and the system scale to an entire year.Then, if "one link fails," i.e., if the rainfall on a specific day exceeds the specified threshold x, then the "entire system fails," i.e., the maximum annual daily rainfall also exceeds the threshold.
Let us consider a system with support N consisting of N links with support 0 .The WLS premise is then expressed as follows: where R i; 0 (x) refers to the survival function of the i-th link.If we assume that all of the links have a common survival function, then (4) can be expressed as: The above leads to an additive relation for the logarithms of the survival function, i.e., Now, if we take into account that: it follows from (6) that the WLS hazard function satisfies the additive property: For a system of N independent and identically distributed random variables, it holds that: Equation (5) implies that under WLS conditions, the link independence does not hold, since: which is markedly different from the link independence condition.
To compare general WLS distribution functions with the Weibull, we use the double logarithm of the reciprocal survival probability, If F N (x) follows the Weibull dependence, then Φ(x; N ) = m ln x xs .For a general WLS distribution, we obtain the following relation based on (5): The above scaling relation in combination with (10), replacing N in the latter with 0 , implies the following scaling property:

Systems Composed of Interacting Links
Above, we assumed an identical survival function for all of the links, a condition that is plausible in a mean-field framework.Herein, we consider systems in which the link properties vary.This is an interesting possibility, which can lead to new probability models for reliability engineering.In general, if R i; 0 (x) = Prob(Λ i ) is the probability that the i-th link survives, the system's survival probability under WLS conditions is given by: Using the law of joint probability it follows that: The conditional probability Prob(Λ 2 |Λ 1 ) may in general depend on the state of Λ 1 and the degree of interaction between the two links.This agrees with the observation of complex interaction patterns in geological fault systems, if the links are identified with single faults or fault groups [25].One possibility for expressing this dependence is by means of the following ansatz: where λ 1 > 0 is a scaling factor.Since R 0 (•) is a decreasing function, λ 1 > 1 implies that the conditional survival probability of Λ 2 is reduced compared to the probability R 0 (x).
More generally, for an N -link system, the survival function of the link Λ i can be expressed as: In light of the above, the system survival function is given by the following equation: which is determined by the coefficient vector λ = (λ 1 , λ 2 , . . ., λ N −1 ) T .The above system survival function is given by a product rule in agreement with the survival function ( 4) of the homogeneous-link system.In addition, all terms of the product involve the same function R 0 (•), albeit the argument of the latter is modified by the respective coefficients λ i , i = 1, . . ., N − 1.Hence, (13) maintains the product rule and the same functional form of the link survival function, which are central features of WLS, but it also captures interactions between the links via the coefficients λ.We illustrate the behavior of the system survival function for a fictitious two-link system in Figure 1, where the former is the survival function of the correlated system, whereas the latter is the standard WLS survival function.For the Weibull survival function, λ effectively renormalizes the scale parameter Hence, the ansatz (13) essentially accounts for the correlations by renormalizing a parameter of the "free distribution", which is an idea with deep roots in statistical physics.Figure 1 demonstrates that values λ > 1 (λ < 1) lead to lower (higher) survival probability for the same x than the uniform-link WLS function.
The proposed system survival function ( 13) is quite general, since it allows for different functional forms of the link survival function.In (13), it is possible to consider correlated random variables λ i as coefficients.In this case, the random variable X becomes doubly stochastic and is described by superstatistics [18,26].We expect that the variability of link parameters is a plausible scenario in heterogeneous physical or engineered systems.For example, a theoretical framework has been proposed that links the statistics of earthquake recurrence times to fracture strength in the Earth's crust [15,27].Since the crust is a heterogeneous material decorated with flaws (faults) of widely varying sizes and orientations, it is reasonable to expect that the fracture strength parameters exhibit local random fluctuations in addition to systematic changes with depth [25].

Hazard Function for Links with Variable Parameters
There is a distinct possibility that links have different statistical properties, even in the case of engineered materials.Thus, a link parameter θ can become a random variable with pdf g(θ).In light of the additivity property (4), the hazard rate of the system is given by the following sum: We can define the average hazard function over the links: If θ follows the pdf g(θ), and the support of θ is the finite interval [θ min , θ max ], we can define an effective link hazard rate h * 0 (x) by: If h 0 (x, θ) and g(θ) are continuous functions of θ, the requirements of the mean value theorem are fulfilled; thus, there exists a θ * ∈ [θ min , θ max ], such that: Equation ( 16) states that the effective link hazard rate maintains the original functional form and that the impact of random θ can be described by means of an "effective parameter" θ * , provided that the conditions of the mean value theorem are upheld.
If ergodic conditions are satisfied so that the events on the links sufficiently sample the space of probable states, one can assume that h 0 (x) = h * 0 (x).The observed survival function then is R 0 (x) = R 0 (x|θ * ), and: which implies that:

Observed Distributions for Non-ergodic Conditions
In the case of natural phenomena, however, the available sample is often finite, and the ergodic condition may not be justifiable.For example, an earthquake catalog for a particular area covers only a small portion of the area's geological history.This means that the full ensemble, involving both the distribution of the scale parameters and the distribution of x values for fixed parameter values, is not adequately sampled during the observation process.If g(θ) has a narrow mode (e.g., Weibull with m > 1), the estimated pdf ĝ(θ) will have a peak in the interval [θ min , θ max ] near the mode of g(θ).For disperse parameter distributions (e.g., Weibull with m ≤ 1), ĝ(θ) may exhibit more complicated behavior due to insufficient sampling.
In order to illustrate the effect of the sample size and the link parameter's variability on the observed link probability distribution, we conduct a numerical experiment.We simulate a random variable X that represents individual links, assuming that X follows the Weibull distribution f 0 (x; m, x s ) with a fixed shape parameter m and uniformly distributed scale parameter x s .The distribution of x s is g(x s ) = (x s;max − x s;min ) −1 if x s ∈ [x s;min , x s;max ] and g(x s ) = 0 otherwise.We generate a large set containing n total = 10 6 random numbers drawn from Weibull distributions with different values of x s .The latter correspond to different links drawn from the uniform distribution g(x s ).We define the mixing ratio as the number of sampled links over the number of samples per link, i.e., ρ = N n s , N = #number of links, n s = #samples per link The ratio ρ determines the mixing of random values from distributions with different scale parameters.Small values of ρ denote that a small number of links are sampled compared to the samples taken from each link, whereas high values of ρ denote the opposite.
In order to investigate different mixing scenarios, we consider two values of ρ equal to 0.1 and 10 −4 .The lower limit x s;min of the scale parameter distribution is set to one, whereas the width δ = x s;max − x s;min of the uniform distribution is set to δ = 1 and 100.From the initial set of n total random numbers, we randomly select n = 10 3 and 10 5 samples, for each of the four scenarios obtained by the different combinations of ρ and δ.The different values of n allow investigating the sample size effect.For each sample, we calculate the maximum likelihood estimates of the optimal-fit Weibull parameters, which are listed in Table 1.For m = 0.7, the sample-based estimates of the Weibull modulus are close to the true value, especially in the case of the narrower distribution of the scale parameter.For m = 2.5, the estimates of the Weibull modulus are close to the true value only for the narrow scale parameter distribution.We also show the estimated Φ(x; ) derived from the empirical (sample-based) cumulative distribution function in Figure 2.For reference purposes, if the data follow the Weibull distribution, Φ(x; ) is a straight line with slope m and intercept m ln xs .As shown in Figure 2(c), for m < 1 and large n, the Weibull distribution is obtained.For smaller n (Figure 2(a)), deviations from the Weibull emerge.For m > 1 and small n, the data cannot be described with a Weibull distribution (Figure 2(b)).For a larger sample size (Figure 2(d)), Φ(x; ) demonstrates Weibull tails with different slopes.Overall, for m = 0.7 (in general, for m < 1), the distribution of sampled values maintains, albeit approximately, the Weibull characteristic shape and the original modulus.For m > 1, the Weibull scaling is maintained only for a narrow scale parameter distribution.

Exponential and Non-Exponential Link Survival Functions
Whereas the Weibull distribution has proven quite successful in the analysis of extreme value data, in certain cases, deviations from the Weibull scaling have been noticed in the tails of empirical distributions [27,28].As shown above, such deviations could occur due to fluctuations of the link hazard rate parameters.Below, we investigate a different possibility, which is based on non-exponential link survival functions with deterministic parameters.Ideally, such functions should be obtained from microscopic theories that describe the dynamics of the relevant physical processes inside the links.Nevertheless, the definition of links and their respective probability functions is not always possible in terms of such detailed analysis [6,28].Hence, for practical purposes, it is useful to consider phenomenological definitions of the survival function.

Weibull Link Ansatz
In his original paper, Weibull assumed that R 0 (x) = e −φ(x) , where φ(x) ≥ 0 is a monotonically increasing function [3].In particular, he proposed that φ(x) = (x/x s ) m , where x s is the scale parameter and m is the shape parameter or Weibull modulus.More generally, the exponential dependence R n (x) = e −n φ(x) , implies the self-similarity relation: This self-similarity property implies that the survival function for any cluster comprising several links has the same functional form, except for a trivial size scaling, as the single-link survival function.The Weibull distribution can also be derived without invoking the exponential dependence of the link survival function.This approach is based on the additivity of the logarithmic survival functions, i.e., Since the system's failure is caused by the weakest link, we can assume that for N 1, this occurs at x, such that F 0 (x) 1. Thus, we can use the approximation: Next, one assumes that F 0 (x) ≈ (x/x s ) m , x < x c , where x c is a threshold above which the algebraic dependence of F 0 (x) ceases.This argument leads to the Weibull expression; however, it does not clarify the tail behavior of the distribution function for x x c .Based on the self-similarity relation (19) and the pioneering works of Daniels and Smith [29,30] that use a Gaussian strength distribution for the links, Curtin showed that the strength distribution of heterogeneous materials tends asymptotically to the Weibull form and that the system's strength parameters depend on both the system and the link sizes [21].

Non-Exponential Link Ansatz
An acceptable link survival function should satisfy the following properties: (i) it should take real values in the interval [0, 1]; (ii) it should be a monotonically decreasing function of x.Let us introduce the following simple functional form: which gives admissible survival functions if φ(x) satisfies (i) φ(0) = 0, (ii) lim x→∞ φ(x) → ∞ and (iii) φ(x) is a monotonically increasing function of x.Condition (iii) is necessary to ensure that dR 0 (x)/dx ≤ 0, because the latter is given by: The second term in the product on the right-hand side of the above is negative, and thus, dR 0 (x)/dx ≤ 0 only if dφ(x)/dx ≥ 0.
In light of (7), one can then show that the hazard function is given by: Hence, for x 1, i.e., φ(x) → ∞, it holds that h(x) ≈ d ln φ(x)/dx.The main appeal of ( 20) is its simplicity.On the other hand, the condition (19) requires solving the following set of nonlinear equations: , for all N = 2 2m , m = 1, 2, . . ., which is satisfied by: The above, however, does not lead to a self-similar expression for φ N/2 (x).
In the following sections, we investigate certain probability distributions for systems with a finite number of links that are based on non-exponential link survival functions.

The Gamma Distribution
First, we consider a slight modification of the exponential link survival function, which is based on the gamma distribution.
where x s and m are respectively the scale and shape parameters, γ(z, m) is the lower incomplete gamma function defined by means of the integral: and Γ(m) is the gamma function that is obtained from γ(z, m) at the asymptotic limit, Γ(m) = lim z→∞ γ(z, m).The pdf of the gamma distribution is given by: The hazard rate of the gamma distribution is respectively given by:

The Modified Gamma Distribution
We use the WLS scaling relation ( 5) and a gamma link probability function to define the modified gamma distribution, which applies to the entire system as follows: The quantile function (inverse survival function) of the modified gamma distribution is: where γ −1 (•) is the inverse of the lower incomplete gamma function.The quantile function of the modified gamma distribution is plotted in Figure 3.The median of the modified gamma distribution is Based on (5) and assuming uniform link parameters, Equation ( 9) leads to the following modified gamma pdf for the N -link system: Equation ( 29) can also be expressed as: where f γ (x) is given by ( 24) and g N (x) is a modulation function given by: The hazard rate of the modified gamma distribution is given according to (8) by: where h 0 (x) is given by ( 25).

Motivation and Properties
The gamma distribution is used in reliability engineering to model extreme events, such as the lifetime of devices, the intensity of rainfall and earthquake return times [1].The gamma distribution is also known in statistical physics as the stretched exponential distribution for m < 1, whereas it defaults to the exponential distribution for m = 1.
The modified gamma survival function is constructed assuming an WLS system that comprises N links that have an identical link failure probability that follows the gamma cdf.The modulation function g N (x) given by (30) takes values between N (at x = 0) and zero (at x → ∞).For m > 1 and increasing N , g N (x) shifts the mode of the pdf to lower values and reduces the width of the peak, as shown in Figure 4(a).Increasing N also suppresses the right tail of the pdf.The gamma distribution is obtained for N = 1, as it is easily verified from (29).For m < 1, as N increases, the g N (x) increases the values of f N (x) as x → 0 and reduces the right tail, as shown in Figure 4  The gamma hazard rate is plotted in Figure 5; it is an increasing function for x < x s , whereas at x 1, it is dominated by the exponential dependence, which leads to an asymptotically constant hazard function.In contrast, the Weibull distribution for m > 1 has a hazard rate that increases with x [13,19].In Figure 6 we plot the cdf and Φ(x; N ) of the modified gamma distribution for m = 0.7, m = 2.5 and for different values of N .We compare these functions to the respective functions for the Weibull distribution.For m > 1, the Φ(x; N ) curves are straight lines on the Weibull plot for ln(x/x s ) ≤ a c , where a c ≈ −1, but become concave as ln(x/x s ) increases further.This means that the right tail of the modified gamma distribution is heavier than the respective Weibull tail.In contrast, for m < 1, the Φ(x; N ) curve becomes convex upwards for ln(x/x s ) > a c , and the right tail of the modified gamma distribution is lighter than the respective Weibull tail.The Φ(x; N ) curves of the modified gamma distribution for different N with fixed m collapse onto a single curve by parallel shifting, in agreement with the scaling relation (12).The definition of the κ-exponential function for κ ≥ 0 is given by [31][32][33][34][35]: The κ-exponential tends asymptotically to the standard exponential at the limits z → 0 or κ → 0, i.e., The inverse function of the κ-exponential is the κ-logarithm, which is defined by:

The κ-Weibull Function
If κ = 1/N , where N is the number of links in a system, we can define a system survival function based on: The κ-Weibull admits explicit expressions for the cdf, the survival function and the pdf, which are given respectively by the following expressions: The κ-Weibull hazard function is obtained from (31) using ( 37) and ( 38) as follows: The κ-Weibull quantile function is defined by means of the survival function as follows:

Conclusions
We investigated the statistical implications of the principle of weakest-link scaling.This type of scaling is used to model the probability distributions of extreme events, including the failure strength of various heterogeneous materials, wind speed, earthquake recurrence times and hydrological maxima.We have focused on various statistical mechanisms leading to deviations from the Weibull distribution, which is commonly used in weakest-link scaling systems.Our motivation stems from empirical distributions that deviate from the Weibull in physical systems expected to follow weakest-link scaling, e.g., [7,15,27,28,36].Our results include a modified gamma distribution, the κ-Weibull distribution and weakest-link scaling distributions with non-homogeneous links that have random (possibly correlated) parameters.The distributions proposed above modify the tail behavior in comparison with the classical Weibull dependence.We thus believe that they are of interest to both statistical physicists and reliability or systems engineers.
The modified gamma distribution, which we introduce herein, is derived by applying weakest-link scaling to systems with a finite number of links.The modified gamma distribution has a heavier tail than the respective Weibull for m > 1 and a lighter tail for m < 1.The κ-Weibull distribution is also shown to follow from weakest-link scaling.A notable property of the κ-Weibull distribution is its ability to exhibit a power-law tail, whereas the modified gamma distribution is dominated by exponential decay at large values.The power-law dependence of the κ-Weibull is derived from the interaction between the system and the links, expressed in terms of a finite-valued scaling factor N in the link survival function.For N → ∞, the κ-Weibull tends to the Weibull model, and the power-law decay at x 1 is replaced by the exponential Weibull dependence.We have argued in [27] that the κ-Weibull distribution can be observed in the case of earthquake recurrence times if the monitored area comprises a finite number of faults.We have also linked the number of effective units in this system to the inverse of κ.This type of behavior can also be observed if the system under study does not include all of the units (e.g., faults) of an interdependent system.In this case, events that are outside the boundaries of the observed system are not taken into account, thus leading to longer recurrence times, which generate the power-law tail.The heavy tail of the κ-Weibull may also find applications in biology, e.g., in modeling the inter-spike interval distributions of cortical neurons, which is known to follow a power law [37].
We have also proposed a different approach for inserting correlations in weakest-link scaling systems by means of a variable scaling factor in the link survival function.If this factor is random (but possibly correlated between links), the observed random variable X is described by an appropriate superstatistics formulation.In the case of Weibull link scaling, we have shown that if the scale parameter is uniformly distributed and uncorrelated across links, the observable probability distribution is controlled by three factors: the width of the scale parameter distribution, the Weibull modulus and the mixing ratio, which determines the balance between the number of links and the values contributed by each link to the observed sample.In the limit of a very large sample, it is also shown that the system's survival function regains the Weibull form, albeit with a renormalized scale, provided that the distribution of the scale parameter has a compact support.In further research, we will investigate the impact of correlations between link distribution parameters on the observed system survival function, as well as applications of this formalism to environmental extreme events that exhibit weakest-link scaling.
Another extreme value distribution that fits within the weak-scaling formalism is the Gumbel distribution [38].The Gumbel distribution is of the general form F (x) = exp(− exp(−x)) and belongs to the family of generalized extreme value distributions.The Gumbel is characterized by a different mathematical expression and size scaling than the Weibull distribution.There is both theoretical and experimental evidence (e.g., in the fracture of brittle ceramics) that in certain cases, the Gumbel is a more suitable asymptotic form than the Weibull [22,38].One may therefore wonder if the same arguments that led to the κ-Weibull distribution also apply to the Gumbel distribution.At this point, however, we lack insight into how this construction should be implemented.Hence, we will leave this as an open subject for future research.

Figure 2 .
Figure 2. Sample estimates Φ(x; ) of data sets with sizes n = 10 3 (top) and n = 10 5 (bottom).The samples are randomly drawn from a total of n total = 10 6 simulated values generated from Weibull distributions with fixed m = 0.7 (left column) and m = 2.5 (right column).The x s are uniformly distributed in the interval [1, 1 + δ] with δ = 1, 100.The mixing ratio that determines the number of x s values (links) over the number of random values per link is ρ = 0.1, 10 −4 .

Figure 3 .
Figure 3. Quantile function for the modified gamma distribution (27) with m = 2.5 (a) and m = 0.7 (b) for different link sizes N (shown with different line types).

Figure 4 .
Figure 4. Modified gamma distribution pdf f N (x) with m = 2.5 (a) and m = 0.7 (b) for different link sizes N (shown with different line types).

Figure 6 .
Figure 6.Cumulative distribution function F N (x) (left) and associated Φ(x; N ) (right) for the modified gamma distribution with m = 2.5 (top) and m = 0.7 (bottom) and different values of N (shown by different line types), as well as the respective Weibull functions with the corresponding Weibull modulus.

7 Figure 7 .
Figure 7. Semi-logarithmic plots of κ-Weibull pdfs for x s = 1, different values of N (shown by different line types), and Weibull modulus equal to (a) m = 2.5 (b) m = 0.7.

Table 1 .
Estimates of Weibull distribution parameters for the simulated data of Figure2.