# A Bounding Technique for Probabilistic PERT

Maksim Goman

► To cite this version:

Maksim Goman. A Bounding Technique for Probabilistic PERT. 13th International Conference on Research and Practical Issues of Enterprise Information Systems (CONFENIS), Dec 2019, Prague, Czech Republic. pp.112-122, 10.1007/978-3-030-37632-1_10. hal-03408414

HAL Id: hal-03408414

https://inria.hal.science/hal-03408414

Submitted on 29 Oct 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# A Bounding Technique for Probabilistic PERT

Maksim Goman

Johannes Kepler University, Altenberger Street 69, 4040 Linz, Austria

Maksim.Goman@jku.at

http://www.jku.at/ie

**Abstract.** Bounding time distributions has been an effective way of improvements of the original PERT method. Analytical enhancement of PERT is often an effective time bounding approach. However, one thing that is missing today is a combination of time distributions which parameters can be effectively obtained empirically and the effective bounding technique for them. We aim at addressing this gap and suggest Cornish-Fisher expansion (CFE) to compute time bounds in formal models like the classical PERT.

We argue that CFE allows us to evaluate analytically approximate time bounds easily without resort to simulations. This bounding approach is useful in case of complex distribution functions of task durations, because analytical derivation of project completion time distribution is tedious. Our example shows CFE usage for uniform time distributions and comparison with time bounds of classical PERT.

**Keywords:** Stochastic PERT · Bounding technique · Time bounds · Project management · Path duration.

## 1 Introduction

Project Evaluation and Review Technique (PERT) [1] (classical PERT) is known today as a method of project time evaluation. The principal idea is to perform a probabilistic analysis of the project completion time. Unfortunately, there has been a lot of criticism of the classical PERT due to its model assumptions, e.g. [2]. One of the most unfortunate assumptions is use of central limit theorem (CLT) to approximate path durations with Gaussian distribution, regardless of distribution function (DF) of individual activity times. CLT appeals to be an easy solution to the complex problem of project time evaluation in a stochastic activity network (SAN). The alternative is aggregation of random time variables along the paths. CLT introduces ambiguity about initial time distributions of tasks in classical PERT, i.e. they had initially beta DF, but going back from the resulting normal distribution of the whole project time we can assume that they could be marginally normal distributions. In this way, there is little use of constructing initial beta distribution of tasks. There are a number of equivalent conditions for data in order that CLT can be applied in the theory [3]. Goman [4] has analyzed the applicability of CLT for classical PERT problems and confirmed the need to verify the CLT conditions for the given data in order to have consistent and reasonable results of classical PERT analysis.
There have been many attempts at using times with Gaussian DF in PERT analysis [5]. Choice of the normal DF was due to simplified modelling and calculations, especially for multivariate case. Computational difficulties called attention to normal distributions many decades ago because manipulation of arbitrary distribution functions (DFs) was hard due to slow computers and lack of readily available software tools for that. Computation time for a project SAN is not an issue today for any realistic size of PERT SAN. Aggregation of known DFs can be done fast using modelling. Moreover, the choice of normal distribution needs more explanations for its negative range and infinite tails of time values.

However, this is not only about assumptions. Basic principles of management are even more important, e.g. what kind of the outcome should the analysis produce to a decision maker (DM) in project management? DM can not expect exact time estimation of possible project time due to unreliable initial time estimations, unknown real distribution types and other uncertainties in the project model. These distributions are only artefacts of the model and follow from some assumptions. In fact, the only thing possible about prospective time analysis is reduction of uncertainty, i.e. more or less accurate estimation of time. The resulting time distribution is obtained from complex aggregation of time distributions of all project activities and can not be verified until certain time in the future. Therefore, any mathematically precise derivations of the resulting distribution are not very vital. Instead, a good guess of the form of the distribution seems more appropriate.

In reality, DM does not need a good approximation of time DF, but a good estimation of time bounds (intervals) for activity start and end times. A good example is critical path method (CPM) that is considered useful today and referred to in textbooks. Although deterministic, it shows the information that a project manager needs and a schedule can be built on it. The same is true for different possible project events like milestones. Because the form of the resulting distribution can be diverse (skewed, non-continuous, discrete), bounds that reflect probability density function (PDF) concentration can be of more value for the DM. A bounding technique that returns quantiles of DF as bounds with certain reliability level can help the DM to manage project in the circumstances of uncertainty. In particular, it can be simple and useful in the pre-project estimations and in the very beginning of a project when lack of empirical data (related to the current project) is an issue. Additionally, a good bounding method should be simpler than derivation of a precise DF or application of a modelling technique.

Analytical bounding techniques have long history [6–11]. There are bounds for PDF and for expected value of the project time distribution. Classical PERT also produces a lower bound of expected value of project completion time. Moreover, classical PERT implies beta distributions to express uncertain time estimations. There is no practically proven method today that can easily operationally obtain time estimations with normal, exponential or beta time distribution. On the opposite, simpler distributions like uniform and triangular are possible to obtain from experts and verify the estimations (e.g. [12]). Analytical bounding is still appealing for its completeness. One of the goals of deterministic CPM
method is determination of time bounds of earliest and latest start and end
times for each activity in SAN. As DF assumptions are predefined and SAN is
given, it should allow computational complexity of O(n) to compute (or recom-
pute) n activity times like those of CPM without modelling for the stochastic
PERT problem.

Consideration of uniform or triangular distributions instead of beta distri-
bution are new in this context since they almost have not been considered for
PERT problems in the last 30 years. These DFs do not have many useful math-
ematical properties that simplify aggregation in PERT analysis. However, it is
easier to obtain their parameters in practice [12, 13]. Johnson [13] has shown
that simpler and intuitively obvious triangular distribution can be very close to
beta DF and proposed a procedure for its parameter estimation.

In order to set up a basis for classical PERT improvement, we would like to
address one piece of the aggregation task in this paper. This is aggregation of
a number of serial activity durations for known identical independent distribu-
tions (iid.) of activity times. We suggest an analytical bounds for random time
aggregation in order to solve this problem *in general* for different possible activ-
ity time distributions. This technique can be applied to *consecutive tasks* or *to
paths* for their comparison (paths need to be enumerated first). We believe, that
it is possible to extend this technique of time analysis to full project graph, but
this needs consideration of conditional DF behavior during collapsing of parallel
activities or paths. We leave this task for the future.

Our bounding approach presumes application of CFE [14] using known mo-
ment generating function (MGF) of activity DFs. This promises an approxi-
mation of time quantiles of the real convoluted DF of the last activity in the
sequence. In order to verify the quality of our bounds, we consider the case of
uniformly distributed activity times. Fortunately, there is a closed form PDF
expression for the sum of random variables for this case. Thus, we can resolve
any time quantile with PDF. However, we are searching for a better bounding
technique *in general*. Therefore, this particular case is convenient for verifica-
tion of the quality of our bounds because this DF type expresses the largest
uncertainty. Thus, we verify the quality of our CFE bounds (quantiles of the
sum) as an absolute error in comparison with known analytical solution for the
quantiles based on known PDF of the sum of uniform distributions (UD). The
CFE technique should work with any DF types, including discrete, mixed or non
identical DFs. In case of convolution of other DFs, verification of the bounding
error in general can be done with simulation techniques.

Using normal distribution and CLT is still popular and original PERT em-
ploys it. In fact, the kernel of the original PERT is nothing else as aggregation
of Gaussian random time variables along one critical path (CP). Therefore, we
consider it useful to compare CFE bounds and quantiles for a normal DF as
in the classical PERT analysis and for another perspective distribution, namely
uniform distribution. Nevertheless, we should repeat again that serial sub-paths
may not be long enough in order to approximate their duration with CLT.
The paper is organized as follows. In section 2 we define the problem and give necessary theoretical information required for CFE understanding and bounds derivation. In section 3 we provide expressions for determination of necessary DF parameters for CFE application, especially for a sum of iid. uniformly distributed time variables. Bounds are computed and compared to the classical PERT in an example problem in section 4. Conclusion summarizes our findings.

## 2 Theoretical Background

### Problem Definition

Let $\mathbf{A} = (A_1, A_2, \dots, A_n)$ be a vector of random time durations of $n$ project activities $t_1, t_2, \dots, t_n$ (activities on arcs) with defined precedence relations $t_i < t_j, i \in \overline{1, n-1}, j \in \overline{i+1, n}$ from the start node (event) $e_0$ to the finish node $e_n$. The activities are performed sequentially. These sequence of activities $\mathbf{A}$ is a subset of all activities of the project and they constitute one possible path or a part of it. We consider the case of UD activity durations $A_i \sim U(a_i, b_i)$ in this paper. The duration of their execution time is the completion time of the last activity $t_n$ before the end event $e_n$. Naturally, this time is a new random variable $Z = \sum_{i=1}^{n} A_i$.

We are interested in concentration of PDF of $Z$. We introduce the reliability threshold $\alpha \in [0.5, 1]$ for the bounds derivation. The lower bound (LB) and upper bound (UB), of project time, such that the real time can be below (above) the bound with respective probability $P_{LB} = (1-\alpha)/2 = \alpha'$ and $P_{UB} = 1-\alpha'$. Thus, the probability that the PDF of the value $Z$ gets into the interval [LB, UB] is $\alpha$. These bounds are related to one of the most important project management objectives, viz. the LB (UB) is bound of the earliest (latest) start or end times of any activity. The LB and UB of the start time of the specific activity $k$ is a random time variable $\sum_{i=1}^{k-1} A_i$ at the moment just before the activity $k$ begins. The LB and UB of the end time of the specific activity $k$ is a random time variable $\sum_{i=1}^{k} A_i$ at the moment when the activity $k$ has ended.

### Moment-generating Function

For $\forall w \in \mathbb{R}$, MGF of a random variable $X$ is defined as [15]:

$$M_X(w) = E(e^{wx}) = \int_{-\infty}^{\infty} e^{wx} f_X(x) dx, \quad (1)$$

Let $Y = bX + a$ be a new random variable for a random variable $X$ and $a, b \in \mathbb{R}$, then the MGF for $Y$ is related to the MGF of $X$ by $M_Y(w) = e^{aw}M_X(bw)$. For our case of a $n$-dimensional random vector $\mathbf{A}$ with independent components $A_i$, $n$-dimensional row vector $a, \mathbf{b} \in \mathbb{R}$, the random variable $Y = a + \mathbf{bX}$ has the MGF: $M_Y(w) = e^{aw}M_{A_1}(b_1w)M_{A_2}(b_2w)\dots M_{A_n}(b_nw)$.

MGF produces $k$-th moments $\mu_k' = E(X^k)$ of the convolution of $n$ random variables as the value of the $k$-th derivative at its parameter point $w = 0$.
## Cornish-Fisher Expansion

The CFE is a tool for random variables quantile approximation using only its first few cumulants. According to Stuart and Ord [14], the cumulants of the order r of a random variable X are values $κ_r$, such that $\forall t$

$$exp\left(\sum_{r=1}^{\infty} \frac{\kappa_r t^r}{r!}\right) = \sum_{r=1}^{\infty} \frac{E(X^r)t^r}{r!}. \quad (2)$$

Cumulants of X are connected to moments of X. Cumulant-generating function is $K(t) = \ln(M_X(w))$. Like MGF, a specific cumulant $κ_r$ is obtained from r-order differential of $K(t)$ at point zero: $κ_r = K^{(r)}(0)$.

Use of cumulants gives an advantage in the context of our problem. Cumulants of a sum of n independent random variables is the sum of their respective cumulants, e.g. for independent X and Y, $K_{X+Y}(t) = K_X(t) + K_Y(t)$ and this is true for all orders r of cumulants of the sum $κ_r(X + Y) = κ_r(X) + κ_r(Y)$.

Cumulants $κ_r$ of a random variable X can be expressed in terms of its mean value $μ$ and its central moments $μ_r = E[(X - μ)^r]$. Alternatively, the same cumulants $κ_r$ can be expressed in terms of only raw (noncentral) moments $μ'_r = E[X^r]$. Both the cases are summarized in Table 1 [14].

The CFE tries to approximate the quantile q of a target DF taking into consideration higher moments (skewness and kurtosis) of that DF to adjust for its non-normality. Thus, for a normally distributed random variable X with $μ = 0$ and $σ = 0$, Cornish and Fisher [16] derived an expansion that enables approximation of q-quantile $\Phi_X^{-1}(q)$ using the five cumulants of X and quantile function of Gaussian distribution $\Phi_Z^{-1}(q)$, $Z \sim N(0, 1)$. There are several versions of the expansion that use different number of cumulants. The formula that uses five cumulants is as follows (all terms are required) [15]:

$$x^* = \Phi_X^{-1}(q) = \Phi_Z^{-1}(q) + \frac{(\Phi_Z^{-1}(q))^2 - 1}{6}\kappa_3 + \frac{(\Phi_Z^{-1}(q))^3 - 3\Phi_Z^{-1}(q)}{24}\kappa_4 - \frac{2(\Phi_Z^{-1}(q))^3 - 5\Phi_Z^{-1}(q)}{36}\kappa_3^2 + \frac{(\Phi_Z^{-1}(q))^4 - 6(\Phi_Z^{-1}(q))^2 + 3}{120}\kappa_5 - \frac{(\Phi_Z^{-1}(q))^4 - 5(\Phi_Z^{-1}(q))^2 + 2}{24}\kappa_3\kappa_4 + \frac{12(\Phi_Z^{-1}(q))^4 - 53(\Phi_Z^{-1}(q))^2 + 17}{324}\kappa_3^3. \quad (3)$$

<table><caption>**Table 1.** Computation of cumulants $κ_r$ from moments of random variable X.</caption><thead><tr><th>Using central moments</th><th>Using raw moments</th></tr></thead><tbody><tr><td>κ<sub>1</sub> = μ</td><td>κ<sub>1</sub> = μ'<sub>1</sub></td></tr><tr><td>κ<sub>2</sub> = μ<sub>2</sub></td><td>κ<sub>2</sub> = μ'<sub>2</sub> - μ'<sub>1</sub><sup>2</sup></td></tr><tr><td>κ<sub>3</sub> = μ<sub>3</sub></td><td>κ<sub>3</sub> = μ'<sub>3</sub> - 3μ'<sub>2</sub>μ'<sub>1</sub> + 2μ'<sub>1</sub><sup>3</sup></td></tr><tr><td>κ<sub>4</sub> = μ<sub>4</sub> - 3μ'<sub>2</sub><sup>2</sup></td><td>κ<sub>4</sub> = μ'<sub>4</sub> - 4μ'<sub>3</sub>μ'<sub>1</sub> - 3μ'<sub>2</sub><sup>2</sup> + 12μ'<sub>2</sub>μ'<sub>1</sub><sup>2</sup> - 6μ'<sub>1</sub><sup>4</sup></td></tr><tr><td>κ<sub>5</sub> = μ<sub>5</sub> - 10μ<sub>3</sub>μ<sub>2</sub></td><td>κ<sub>5</sub> = μ'<sub>5</sub> - 5μ'<sub>4</sub>μ'<sub>1</sub> - 10μ'<sub>3</sub>μ'<sub>2</sub> + 20μ'<sub>3</sub>μ'<sub>1</sub><sup>2</sup> + 30μ'<sub>2</sub><sup>2</sup>μ'<sub>1</sub> - 60μ'<sub>2</sub>μ'<sub>1</sub><sup>3</sup> + 24μ'<sub>1</sub><sup>5</sup></td></tr></tbody></table>
### 3 CFE Application to Bounds of *n* Activity Distributions

Using concepts of CFE, MGF and basics of probability theory, we derive the simple analytical approximate bounds for random activity time durations. We can use CFE to approximate quantiles of $Z = \sum_{i=1}^{n} A_i$ with means $\mu_i$ and standard deviations $\sigma_i$. First, we normalize $A_i$ as required by CFE: $A'_i = (A_i - \mu_i)/\sigma_i$. Now $A'_i \sim N(0, 1)$ and central moments of $A'_i$ are obtained from the central moments of $A_i$ with the expression $\mu'_r = \mu_r/\sigma^r$ [15]. Because $A_i$ and consequently $A'_i$ are independent random variables, we compute individual cumulants of $A'_i$ and add the respective cumulants as it was explained in section 2 to get cumulants of $Z' = \sum_{i=1}^{n} A'_i$. Then, we resolve the required quantile $q$ approximation $z'^*$ of $Z'$ applying CFE expression (3). The approximated value of the original $z'^*$-quantile of $Z$ is obtained through de-normalization: $z' = z'^*\sigma + \mu$, where $\sigma = \sqrt{\sum_{i=1}^{n} \sigma_i^2}$ and $\mu = \sum_{i=1}^{n} \mu_i$.

In order to generalize the approach to moment derivation for other distributions in the prospective work, we will consider the ways of raw and central moment derivation in more details below. The central moments of $A_i$ are required for CFE application. However, we can express central moments $\mu_r = E[(X - \mu)^r]$ through non-central moments $\mu'_k = E(X^k)$. Denoting powers $k$ of the mean value $\mu^k$, the required formulas for the central moments are shown in Table 2.

In their turn, raw moments can be obtained from known MGF taking $\lim_{t \to 0}$. However, this is hard for UD. Alternatively, raw moments can be determined for a uniform distribution via straightforward integration of $\mu'_n = \int x^n f(x) dx$ with PDF in the terms of Heaviside step function $f(x) = \frac{H(x-a) - H(x-b)}{b-a}$ [17]:

$$ \mu'_{n} = \int_{-\infty}^{\infty} \frac{(H(x-a) - H(x-b))}{(b-a)x^n} dx = \frac{b^{n+1} - a^{n+1}}{(n+1)(b-a)}, \quad (4) $$

that enables to get the expressions of raw moments easier (see Table 2).

The respective central moments can be determined in the same way for UD via integration of $\mu_n = \int (x - \mu)^n f(x) dx$ (see $\mu = \mu'_1$ in Table 2) with PDF in the terms of Heaviside step function $f(x) = \frac{H(x-a) - H(x-b)}{b-a}$ [17]. This approach turns out $\mu_n = \frac{(a-b)^n - (b-a)^n}{2^{n+1}(n+1)}$ and gives the required central moments for the UD case:

$$ \mu_1 = \mu_3 = \mu_5 = 0; \quad \mu_2 = \frac{1}{12}(b-a)^2; \quad \mu_4 = \frac{1}{80}(b-a)^4. \quad (5) $$

<table><caption>**Table 2.** Computation of raw and central moments of random variable $X$.</caption><thead><tr><th>Raw moments</th><th>Central moments</th></tr></thead><tbody><tr><td>$\mu'_1 = \frac{1}{2}(a+b)$</td><td>$\mu_1 = 0$</td></tr><tr><td>$\mu'_2 = \frac{1}{3}(a^2 + ab + b^2)$</td><td>$\mu_2 = E[(X - \mu)^2] = \mu'_2 - \mu^2$</td></tr><tr><td>$\mu'_3 = \frac{1}{4}(a^3 + a^2b + ab^2 + b^3)$</td><td>$\mu_3 = E[(X - \mu)^3] = \mu'_3 + 2\mu^3 - 3\mu\mu'_2$</td></tr><tr><td>$\mu'_4 = \frac{1}{5}(a^4 + a^3b + a^2b^2 + ab^3 + b^4)$</td><td>$\mu_4 = E[(X - \mu)^4] = \mu'_4 - 4\mu\mu'_3 + 6\mu^2\mu'_2 - 3\mu^4$</td></tr><tr><td>$\mu'_5 = \frac{1}{6}(a^5 + a^4b + a^3b^2 + a^2b^3 + ab^4 + b^5)$</td><td>$\mu_5 = E[(X - \mu)^5] = \mu'_5 + 10\mu^2\mu'_3 - 10\mu^3\mu'_2 - 5\mu\mu'_4 + 4\mu^5$</td></tr></tbody></table>
## 4 Illustrative Example

We consider a sequential set of independent tasks with identical UD DFs and obtain approximations of the end time bounds with quantiles for probabilities from $q = 0.1$ to $q = 0.9$ with step 0.1 by means of CFE with five cumulants (3). Results of the classical PERT analysis are also computed. We evaluate the maximum error of CFE and PERT quantile approximation with available expression for PDF $f_Z(z)$ of the sum of $n$ UD. According to the original expression formula for $f_Z(z)$, the lower bound of DFs should strictly equal to zero $B_i \sim U(0, c_i)$ [18]:

$$f_n(z) = \frac{1}{(n-1)! \prod_{k=1}^{n} c_k} \left( z^{n-1} + \sum_{k=1}^{n} (-1)^k \sum_{l=1}^{k} \left[ (z - \sum_{j=1}^{k} c_{jl})_+ \right]^{n-1} \right), \quad (6)$$

where $J_k = \{(j_1, \dots, j_k); 1 \le j_1 < j_2 < \dots < j_k \le n\}; 0 \le z \le \sum_{i=1}^{n} c_i, n \in \mathbb{N}, z \in \mathbb{R}$ and $x_+ = \max(0, x) \forall x \in \mathbb{R}$.

We obtain $B_i \sim U(0, c_i)$ after introduction of new variables $c_i = b_i - a_i$: $B_i = A_i - a_i; i = 1, 2, \dots, n$. This is a shift of DFs by constants $a_i$. It is possible to process the random parts $B_i$ with known PDFs and then apply correction by adding $\sum_{i=1}^{n} a_i$.

By definition, a value $z^*$ such that cumulative distribution function (CDF) $F(z^*) = f(Z \le z^*) = q$ is called a quantile of order $q$, $q \in (0, 1)$. Quantiles for probability $q$ can be evaluated by solving the equality with known PDF $f(z)$ or the respective CDF $F(z)$: $q = F(z^*) = \int_{l}^{u} f(z)dz$, where in our case $l = 0$ is the lower bound and $u = z^*$ is the sought upper bound for $B_i$ ($q$ is specified).

We consider two sub problems. One is aggregation of two, five and ten identical DFs. The other is addition of random time variable $B \sim U(1, 3)$ $i = 2, 3, 4$ and several unequal time DFs $\sum_{j=1}^{i-1} A_j \sim U(1, c_j)$, $c_j = 2, 4, 6$. We consider that very long task sequences or paths without branches are very unlikely in real projects. Resulting quantiles as bounds for UD are shown in Table 3. We compute CFE approximation $z_{CFE}^*$ of quantiles $z^*$ for given probabilities $q = 0.1, \dots, 0.9$. For comparison, results of classical PERT (PERT$z^*$) are given as well. Absolute errors of $Z^*$ approximation can be seen by comparing the quantiles obtained with known PDF (6) in “Real $z^*$” rows and the respective values $z_{CFE}^*$ or PERT $z^*$.

We also verify the quality of approximations by substitution of the approximated quantile $z^*$ into the integral of the known PDF (6) and observing how different the probability $F_Z(z_{CFE}^*)$ from the initial $q$ value is. Absolute errors (abs.error $q$) of quantiles of orders $q$ between $F_Z(z_{CFE}^*)$ and $q$ are also given in the Table 3. The same verification is performed for classical PERT approximations $F_Z(\text{PERT } z^*)$ and absolute error is given.

It is obvious that classical PERT that uses CLT approximation is significantly worse than CFE approximation unless many identical iid. DFs are considered (and this makes the pattern closer to CLT applicability assumptions). PERT's probability estimation for $q$ quantiles is always worth than CFE for not very equal DFs (and mostly for equal ones) for the most important quantiles ($q =$
0.1, 0.2 for LB and $q = 0.8, 0.9$ for UB) and this difference is impressive, i.e. from 1.5 to 10 times and even more. The difference in time value $z^*$ approximation from the real one may not seem so dramatic, however, due to the scale of time units, the absolute difference can be substantial for the DM. Additionally, CFE underestimates the LB and overestimates the UB, i.e. produces real bounds, whereas classical PERT does vice versa, i.e. underestimates the quantile values from above and below that is undesirable.

# 5 Conclusion

In this study, we attempted at improving only one problem of PERT analysis, namely obtaining time bounds on a sequential set of tasks. For the general case, we can determine upper and lower bounds using bounding techniques. This can be imprecise in case of CLT approximation in classical PERT or if bounding methods do not take into account specific properties of distributions. If time distributions are known, we can determine upper and lower bounds using CFE.

Using CFE, we derived bounds for iid. sum of uniform DFs of activity times and verified the absolute error with exact analytical solution of the quantile problem. The DF type is not used frequently, however it has an advantage of easily obtainable parameter estimation in practical application. We also compared our result with the results of classical PERT analysis. In fact, original PERT employs CLT approximation and thus uses quantiles of Gaussian distribution. Our experiments show that CFE approximation outperforms classical PERT in evaluation of time of sequential tasks with known DFs. This is absolutely important for bounds with high confidence (e.g. $\alpha \ge 0.8$). It seems that CFE application instead of CLT approximation can improve evaluation of CP in the classical PERT.

To the best of our knowledge, our bounding technique for UD based on CFE is pioneering in the area of PERT analysis. CPM like time bounds for the case of several consecutive activities with uniform DFs enable evaluation of duration of sub-critical paths in a SAN. Moreover, CFE is a universal tool that can be applied to the task with other DF types. The computational complexity is $O(n)$ for a known SAN structure without simulation modelling.

This is the first step in the development of an improved bounding technique for a project SAN. Based on the current results for serial activities, an extension of the bounding technique should be developed for a SAN with converging sub-paths. We need to examine conceptual meaning of converging activities into an event before the next common activity after them. Mathematically, this is consideration of a maximum operator for two or more random time variables and choosing the best promising method of DF evaluation for this operator. Another prospective task is evaluation of the bounding approach with triangular time distributions that are also perspective in operational parameter estimation and have been used for project time analysis (e.g. [12]).
<table><caption>**Table 3.** Quality of CFE q-quantile evaluation for the number of tasks i.</caption><thead><tr><th rowspan="2">i</th><th rowspan="2">Method</th><th colspan="9">z<sup>*</sup> for q for the sum of i of identical DFs A<sub>i</sub> ~ U(1,3)</th></tr><tr><th>0.1</th><th>0.2</th><th>0.3</th><th>0.4</th><th>0.5</th><th>0.6</th><th>0.7</th><th>0.8</th><th>0.9</th></tr></thead><tbody><tr><td rowspan="5">2</td><td>Real z<sup>*</sup></td><td>2.894</td><td>3.265</td><td>3.549</td><td>3.789</td><td>4.0</td><td>4.211</td><td>4.451</td><td>4.735</td><td>5.106</td></tr><tr><td>z<sup>*</sup><sub>CFE</sub></td><td>2.812</td><td>3.155</td><td>3.455</td><td>3.732</td><td>4.0</td><td>4.268</td><td>4.545</td><td>4.845</td><td>5.188</td></tr><tr><td>F<sub>Z</sub>(z<sup>*</sup><sub>CFE</sub>)</td><td>0.082</td><td>0.167</td><td>0.265</td><td>0.375</td><td>0.5</td><td>0.625</td><td>0.735</td><td>0.833</td><td>0.918</td></tr><tr><td>abs.error q</td><td>-0.018</td><td>-0.033</td><td>-0.035</td><td>-0.025</td><td>0.0</td><td>0.025</td><td>0.035</td><td>0.033</td><td>0.018</td></tr><tr><td>PERT z<sup>*</sup></td><td>3.715</td><td>3.813</td><td>3.883</td><td>3.944</td><td>4.0</td><td>4.056</td><td>4.117</td><td>4.187</td><td>4.285</td></tr><tr><td rowspan="5">5</td><td>Real z<sup>*</sup></td><td>8.321</td><td>8.887</td><td>9.3034</td><td>9.663</td><td>10.0</td><td>10.337</td><td>10.696</td><td>11.113</td><td>12.679</td></tr><tr><td>z<sup>*</sup><sub>CFE</sub></td><td>7.784</td><td>8.291</td><td>8.862</td><td>9.433</td><td>10.0</td><td>10.567</td><td>11.138</td><td>11.709</td><td>12.216</td></tr><tr><td>F<sub>Z</sub>(z<sup>*</sup><sub>CFE</sub>)</td><td>0.082</td><td>0.167</td><td>0.265</td><td>0.375</td><td>0.5</td><td>0.625</td><td>0.735</td><td>0.833</td><td>0.918</td></tr><tr><td>abs.error q</td><td>-0.018</td><td>-0.033</td><td>-0.035</td><td>-0.025</td><td>0.0</td><td>0.025</td><td>0.035</td><td>0.033</td><td>0.018</td></tr><tr><td>PERT z<sup>*</sup></td><td>9.288</td><td>9.532</td><td>9.709</td><td>9.859</td><td>10.0</td><td>10.141</td><td>10.291</td><td>10.468</td><td>10.712</td></tr><tr><td rowspan="5">10</td><td>Real z<sup>*</sup></td><td>17.644</td><td>18.445</td><td>19.029</td><td>19.530</td><td>20</td><td>20.469</td><td>20.971</td><td>21.555</td><td>22.356</td></tr><tr><td>z<sup>*</sup><sub>CFE</sub></td><td>16.072</td><td>16.703</td><td>17.738</td><td>18.859</td><td>20.0</td><td>21.142</td><td>22.262</td><td>23.297</td><td>23.928</td></tr><tr><td>F<sub>Z</sub>(z<sup>*</sup><sub>CFE</sub>)</td><td>0.015</td><td>0.035</td><td>0.109</td><td>0.269</td><td>0.5</td><td>0.731</td><td>0.891</td><td>0.965</td><td>0.985</td></tr><tr><td>abs.error q</td><td>-0.085</td><td>-0.165</td><td>-0.191</td><td>-0.131</td><td>0.0</td><td>0.131</td><td>0.191</td><td>0.165</td><td>0.085</td></tr><tr><td>PERT z<sup>*</sup></td><td>18.576</td><td>19.065</td><td>19.417</td><td>19.719</td><td>20.0</td><td>20.282</td><td>20.583</td><td>20.935</td><td>21.424</td></tr><tr><td rowspan="5">4</td><td>Real z<sup>*</sup></td><td>4.557</td><td>5.022</td><td>5.380</td><td>5.698</td><td>6.0</td><td>6.302</td><td>6.620</td><td>6.978</td><td>7.443</td></tr><tr><td>z<sup>*</sup><sub>CFE</sub></td><td>4.334</td><td>4.778</td><td>5.202</td><td>5.606</td><td>6.0</td><td>6.394</td><td>6.798</td><td>7.222</td><td>7.666</td></tr><tr><td>F<sub>Z</sub>(z<sup>*</sup><sub>CFE</sub>)</td><td>0.065</td><td>0.143</td><td>0.248</td><td>0.370</td><td>0.5</td><td>0.630</td><td>0.752</td><td>0.857</td><td>0.935</td></tr><tr><td>abs.error q</td><td>-0.035</td><td>-0.057</td><td>-0.052</td><td>-0.030</td><td>0.0</td><td>0.030</td><td>0.052</td><td>0.057</td><td>0.035</td></tr><tr><td>PERT z<sup>*</sup></td><td>5.502</td><td>5.673</td><td>5.796</td><td>5.901</td><td>6.0</td><td>6.099</td><td>6.204</td><td>6.328</td><td>6.498</td></tr><tr><td rowspan="5">4</td><td>Real z<sup>*</sup></td><td>7.107</td><td>7.868</td><td>8.462</td><td>8.993</td><td>9.5</td><td>10.007</td><td>10.538</td><td>11.132</td><td>11.893</td></tr><tr><td>z<sup>*</sup><sub>CFE</sub></td><td>6.562</td><td>7.287</td><td>8.039</td><td>8.775</td><td>9.5</td><td>10.225</td><td>10.961</td><td>11.713</td><td>12.438</td></tr><tr><td>F<sub>Z</sub>(z<sup>*</sup><sub>CFE</sub>)</td><td>0.051</td><td>0.120</td><td>0.227</td><td>0.358</td><td>0.5</td><td>0.642</td><td>0.773</td><td>0.880</td><td>0.949</td></tr><tr><td>abs.error q</td><td>-0.049</td><td>-0.080</td><td>-0.073</td><td>-0.042</td><td>0.0</td><td>0.042</td><td>0.073</td><td>0.080</td><td>0.049</td></tr><tr><td>PERT z<sup>*</sup></td><td>8.112</td><td>8.588</td><td>8.932</td><td>9.226</td><td>9.5</td><td>9.774</td><td>10.068</td><td>10.412</td><td>10.888</td></tr><tr><td rowspan="5">4</td><td>Real z<sup>*</sup></td><td>7.107</td><td>7.868</td><td>8.462</td><td>8.993</td><td>9.5</td><td>10.007</td><td>10.538</td><td>11.132</td><td>11.893</td></tr><tr><td>z<sup>*</sup><sub>CFE</sub></td><td>6.562</td><td>7.287</td><td>8.039</td><td>8.775</td><td>9.5</td><td>10.225</td><td>10.961</td><td>11.713</td><td>12.438</td></tr><tr><td>F<sub>Z</sub>(z<sup>*</sup><sub>CFE</sub>)</td><td>0.051</td><td>0.120</td><td>0.227</td><td>0.358</td><td>0.5</td><td>0.642</td><td>0.773</td><td>0.880</td><td>0.949</td></tr><tr><td>abs.error q</td><td>-0.049</td><td>-0.080</td><td>-0.073</td><td>-0.042</td><td>0.0</td><td>0.042</td><td>0.073</td><td>0.080</td><td>0.049</td></tr><tr><td>PERT z<sup>*</sup></td><td>8.112</td><td>8.588</td><td>8.932</td><td>9.226</td><td>9.5</td><td>9.774</td><td>10.068</td><td>10.412</td><td>10.888</td></tr><tr><td rowspan="5">4</td><td>Real z<sup>*</sup></td><td>7.107</td><td>7.868</td><td>8.462</td><td>8.993</td><td>9.5</td><td>10.007</td><td>10.538</td><td>11.132</td><td>11.893</td></tr><tr><td>z<sup>*</sup><sub>CFE</sub></td><td>6.562</td><td>7.287</td><td>8.039</td><td>8.775</td><td>9.5</td><td>10.225</td><td>10.961</td><td>11.713</td><td>12.438</td></tr><tr><td>F<sub>Z</sub>(z<sup>*</sup><sub>CFE</sub>)</td><td>0.051</td><td>0.120</td><td>0.227</td><td>0.358</td><td>0.5</td><td>0.642</td><td>0.773</td><td>0.880</td><td>0.949</td></tr><tr><td>abs.error q</td><td>-0.049</td><td>-0.080</td><td>-0.073</td><td>-0.042</td><td>0.0</td><td>0.042</td><td>0.073</td><td>0.080</td><td>0.049</td></tr><tr><td>PERT z<sup>*</sup></td><td>8.112</td><td>8.588</td><td>8.932</td><td>9.226</td><td>9.5</td><td>9.774</td><td>10.068</td><td>10.412</td><td>10.888</td></tr><tr><td rowspan="5">4</td><td>Real z<sup>*</sup></td><td>7.107</td><td>7.868</td><td>8.462</td><td>8.993</td><td>9.5</td><td>10.007</td><td>10.538</td><td>11.132</td><td>11.893</td></tr><tr><td>z<sup>*</sup><sub>CFE</sub></td><td>6.562</td><td>7.287</td><td>8.039</td><td>8.775</td><td>9.5</td><td>10.225</td><td>10.961</td><td>11.713</td><td>12.438</td></tr><tr><td>F<sub>Z</sub>(z<sup>*</sup><sub>CFE</sub>)</td><td>0.051</td><td>0.120</td><td>0.227</td><td>0.358</td><td>0.5</td><td>0.642</td><td>0.773</td><td>0.880</td><td>0.949</td></tr><tr><td>abs.error q</td><td>-0.049</td><td>-0.080</td><td>-0.073</td><td>-0.042</td><td>0.0</td><td>0.042</td><td>0.073</td><td>0.080</td><td>0.049</td></tr><tr><td>PERT z<sup>*</sup></td><td>8.112</td><td>8.588</td><td>8.932</td><td>9.226</td><td>9.5</td><td>9.774</td><td>10.068</td><td>10.412</td><td>10.888</td></tr><tr><td rowspan="5">4</td><td>Real z<sup>*</sup></td><td>7.107</td><td>7.868</td><td>8.462</td><td>8.993</td><td>9.5</td><td>10.007</td><td>10.538</td><td>11.132</td><td>11.893</td></tr><tr><td>z<sup>*</sup><sub>CFE</sub></td><td>6.562</td><td>7.287</td><td>8.039</td><td>8.775</td><td>9.5</td><td>10.225</td><td>10.961</td><td>11.713</td><td>12.438</td></tr><tr><td>F<sub>Z</sub>(z<sup>*</sup><sub>CFE</sub>)</td><td>0.051</td><td>0.120</td><td>0.227</td><td>0.358</td><td>0.5</td><td>0.642</td><td>0.773</td><td>0.880</td><td>0.949</td></tr><tr><td>abs.error q</td><td>-0.049</td><td>-0.080</td><td>-0.073</td><td>-0.042</td><td>0.0</td><td>0.042</td><td>0.073</td><td>0.080</td><td>0.049</td></tr><tr><td>PERT z<sup>*</sup></td><td>8.112</td><td>8.588</td><td>8.932</td><td>9.226</td><td>9.5</td><td>9.774</td><td>10.068</td><td>10.412</td><td>10.888</td></tr><tr><td rowspan="5">4</td><td>Real z<sup>*</sup></td><td>7.107</td><td>7.868</td><td>8.462</td><td>8.993</td><td>9.5</td><td>10.007</td><td>10.538</td><td>11.132</td><td>11.893</td></tr><tr><td>z<sup>*</sup><sub>CFE</sub></td><td>6.562</td><td>7.287</td><td>8.039</td><td>8.775</td><td>9.5</td><td>10.225</td><td>10.961</td><td>11.713</td><td>12.438</td></tr><tr><td>F<sub>Z</sub>(z<sup>*</sup><sub>CFE</sub>)</td><td>0.051</td><td>0.120</td><td>0.227</td><td>0.358</td><td>0.5</td><td>0.642</td><td>0.773</td><td>0.880</td><td>0.949</td></tr><tr><td>abs.error q</td><td>-0.049</td><td>-0.080</td><td>-0.073</td><td>-0.042</td><td>0.0</td><td>0.042</td><td>0.073</td><td>0.080</td><td>0.049</td></tr><tr><td>PERT z<sup>*</sup></td><td>8.112</td><td>8.588</td><td>8.932</td><td>9.226</td><td>9.5</td><td>9.774</td><td>10.068</td><td>10.412</td><td>10.888</td></tr><tr><td rowspan="5">4</td><td>Real z<sup>*</sup></td><td>7.107</td><td>7.868</td><td>8.462</td><td>8.993</td><td>9.5</td><td>10.007</td><td>10.538</td><td>11.132</td><td>11.893</td></tr><tr><td>z<sup>*</sup><sub>CFE</sub></td><td>6.562</td><td>7.287</td><td>8.039</td><td>8.775</td><td>9.5</td><td>10.225</td><td>10.961</td><td>11.713</td><td>12.438</td></tr><tr><td>F<sub>Z</sub>(z<sup>*</sup><sub>CFE</sub>)</td><td>0.051</td><td>0.120</td><td>0.227</td><td>0.358</td><td>0.5</td><td>0.642</td><td>0.773</td><td>0.880</td><td>0.949</td></tr><tr><td>abs.error q</td><td>-0.049</td><td>-0.080</td><td>-0.073</td><td>-0.042</td><td>0.0</td><td>0.042</td><td>0.073</td><td>0.080</td><td>0.049</td></tr><tr><td>PERT z<sup>*</sup></td><td>8.112</td><td>8.588</td><td>8.932</td><td>9.226</td><td>9.5</td><td>9.774</td><td>10.068</td><td>10.412</td><td>10.888</td></tr><tr><td rowspan="5">4</td><td>Real z<sup>*</sup></td><td>7.107</td><td>7.868</td><td>8.462</td><td>8.993</td><td>9.5</td><td>10.007</td><td>10.538</td><td>11.132</td><td>11.893</td></tr><tr><td>z<sup>*</sup><sub>CFE</sub></td><td>6.562</td><td>7.287</td><td>8.039</td><td>8.775</td><td>9.5</td><td>10.225</td><td>10.961</td><td>11.713</td><td>12.438</td></tr><tr><td>F<sub>Z</sub>(z<sup>*</sup><sub>CFE</sub>)</td><td>0.051</td><td>0.120</td><td>0.227</td><td>0.358</td><td>0.5</td><td>0.642</td><td>0.773</td><td>0.880</td><td>0.949</td></tr><tr><td>abs.error q</td><td>-0.049</td><td>-0.080</td><td>-0.073</td><td>-0.042</td><td>0.0</td><td>0.042</td><td>0.073</td><td>0.080</td><td>0.049</td></tr><tr><td>PERT z<sup>*</sup></td><td>8.112</td><td>8.588</td><td>8.932</td><td>9.226</td><td>9.5</td><td>9.774</td><td>10.068</td><td>10.412</td><td>10.888</td></tr><tr><td rowspan="5">4</td><td>Real z<sup>*</sup></td><td>7.107</td><td>7.868</td><td>8.462</td><td>8.993</td><td>9.5</td><td>10.007</td><td>10.538</td><td>11.132</td><td>11.893</td></tr><tr><td>z<sup>*</sup><sub>CFE</sub></td><td>6.562</td><td>7.287</td><td>8.039</td><td>8.775</td><td>9.5</td><td>10.225</td><td>10.961</td><td>11.713</td><td>12.438</td></tr><tr><td>F<sub>Z</sub>(z<sup>*</sup><sub>CFE</sub>)</td><td>0.051</td><td>0.120</td><td>0.227</td><td>0.358</td><td>0.5</td><td>0.642</td><td>0.773</td><td>0.880</td><td>0.949</td></tr><tr><td>abs.error q</td><td>-0.049</td><td>-0.080</td><td>-0.073</td><td>-0.042</td><td>0.0</td><td>0.042</td><td>0.073</td><td>0.080</td><td>0.049</td></tr><tr><td>PERT z<sup>*</sup></td><td>8.112</td><td>8.588</td><td>8.932</td><td>9.226</td><td>9.5</td><td>9.774</td><td>10.068</td><td>10.412</td><td>10.888</td></tr><tr><td rowspan="5">4</td><td>Real z<sup>*</sup></td><td>7.107</td><td>7.868</td><td>8.462</td><td>8.993</td><td>9.5</td><td>10.007</td><td>10.538</td><td>11.132</td><td>11.893</td></tr><tr><td>z<sup>*</sup><sub>CFE</sub></td><td>6.562</td><td>7.287</td><td>8.039</td><td>8.775</td><td>9.5</td><td>10.225</td><td>10.961</td><td>11.713</td><td>12.438</td></tr><tr><td>F<sub>Z</sub>(z<sup>*</sup><sub>CFE</sub>)</td><td>0.051</td><td>0.120</td><td>0.227</td><td>0.358</td><td>0.5</td><td>0.642</td><td>0.773</td><td>0.880</td><td>0.949</td></tr><tr><td>abs.error q</td><td>-0.049</td><td>-0.080</td><td>-0.073</td><td>-0.042</td><td>0.0</td><td>0.042</td><td>0.073</td><td>0.080</td><td>0.049</td></tr><tr><td>PERT z<sup>*</sup></td><td>8.112</td><td>8.588</td><td>8.932</td><td>9.226</td><td>9.5</td><td>9.774</td><td>10.068</td><td>10.412</td><td>10.888</td></tr><tr><td rowspan="5">4</td><td>Real z<sup>*</sup></td><td>7.107</td><td>7.868</td><td>8.462</td><td>8.993</td><td>9.5</td><td>10.007</td><td>10.538</td><td>11.132</td><td>11.893</td></tr><tr><td>z<sup>*</sup><sub>CFE</sub></td><td>6.562</td><td>7.287</td><td>8.039</td><td>8.775</td><td>9.5</td><td>10.225</td><td>10.961</td><td>11.713</td><td>12.438</td></tr><tr><td>F<sub>Z</sub>(z<sup>*</sup><sub>CFE</sub>)</td><td>0.051</td><td>0.120</td><td>0.227</td><td>0.358</td><td>0.5</td><td>0.642</td><td>0.773</td><td>0.880</td><td>0.949</td></tr><tr><td>abs.error q</td><td>-0.049</td><td>-0.080</td><td>-0.073</td><td>-0.042</td><td>0.0</td><td>0.042</td><td>0.073</td><td>0.080</td><td>0.049</td></tr><tr><td>PERT z<sup>*</sup></td><td>8.112</td><td>8.588</td><td>8.932</td><td>9.226</td><td>9.5</td><td>9.774</td><td>10.068</td><td>10.412</td><td>10.888</td></tr><tr><td rowspan="5">4</td><td>Real z<sup>*</sup></td><td>7.107</td><td>7.868</td><td>8.462</td><td>8.993</td><td>9.5</td><td>10.007</td><td>10.538</td><td>11.132</td><td>11.893</td></tr><tr><td>z<sup>*</sup><sub>CFE</sub></td><td>6.562</td><td>7.287</td><td>8.039</td><td>8.775</td><td>9.5</td><td>10.225</td><td>10.961</td><td>11.713</td><td>12.438</td></tr><tr><td>F<sub>Z</sub>(z<sup>*</sup><sub>CFE</sub>)</td><td>0.051</td><td>0.120</td><td>0.227</td><td>0.358</td><td>0.5</td><td>0.642</td><td>0.773</td><td>0.880</td><td>0.949</td></tr><tr><td>abs.error q</td><td>-0.049</td><td>-0.080</td><td>-0.073</td><td>-0.042</td><td>0.0</td><td>0.042</td><td>0.073</td><td>0.080</td><td>0.049</td></tr><tr><td>PERT z<sup>*</sup></td><td>8.112</td><td>8.588</td><td>8.932</td><td>9.226</td><td>9.5</td><td>9.774</td><td>10.068</td><td>10.412</td><td>10.888</td></tr><tr><td rowspan="5">4</td><td>Real z<sup>*</sup></td><td>7.107</td><td>7.868</td><td>8.462</td><td>8.993</td><td>9.5</td><td>10.007</td><td>10.538</td><td>11.132</td><td>11.893</td></tr><tr><td>z<sup>*</sup><sub>CFE</sub></td><td>6.562</td><td>7.287</td><td>8.039</td><td>8.775</td><td>9.5</td><td>10.225</td><td>10.961</td><td>11.713</td><td>12.438</td></tr><tr><td>F<sub>Z</sub>(z<sup>*</sup><sub>CFE</sub>)</td><td>0.051</td><td>0.120</td><td>0.227</td><td>0.358</td><td>0.5</td><td>0.642</td><td>0.773</td><td>0.880</td><td>0.949</td></tr><tr><td>abs.error q</td><td>-0.049</td><td>-0.080</td><td>-0.073</td><td>-0.042</td><td>0.0</td><td>0.042</td><td>0.073</td><td>0.080</td><td>0.049</td></tr><tr><td>PERT z<sup>*</sup></td><td>8.112</td><td>8.588</td><td>8.932</td><td>9.226</td><td>9.5</td><td>9.774</td><td>10.068</td><td>10.412</td><td>10.888</td></tr><tr><td rowspan="5">4</td><td>Real z<sup>*</sup></td><td>7.107</td><td>7.868</td><td>8.462</td><td>8.993</td><td>9.5</td><td>10.007</td><td>10.538</td><td>11.132</td><td>11.893</td></tr><tr><td>z<sup>*</sup><sub>CFE</sub></td><td>6.562</td><td>7.287</td><td>8.039</td><td>8.775</td><td>9.5</td><td>10.225</td><td>10.961</td><td>11.713</td><td>12.438</td></tr><tr><td>F<sub>Z</sub>(z<sup>*</sup><sub>CFE</sub>)</td><td>0.051</td><td>0.120</td><td>0.227</td><td>0.358</td><td>0.5</td><td>0.642</td><td>0.773</td><td>0.880</td><td>0.949</td></tr><tr><td>abs.error q</td><td>-0.049</td><td>-0.080</td><td>-0.073</td><td>-0.042</td><td>0.0</td><td>0.042</td><td>0.073</td><td>0.080</td><td>0.049</td></tr><tr><td>PERT z<sup>*</sup></td><td>8.112</td><td>8.588</td><td>8.932</td><td>9.226</td><td>9.5</td><td>9.774</td><td>10.068</td><td>10.412</td><td>10.888</td></tr><tr><td rowspan="5">4</td><td>Real z<sup>*</sup></td><td>7.107</td><td>7.868</td><td>8.462</td><td>8.993</td><td>9.5</td><td>10.007</td><td>10.538</td><td>11.132</td><td>11.893</td></tr><tr><td>z<sup>*</sup><sub>CFE</sub></td><td>6.562</td><td>7.287</td><td>8.039</td><td>8.775</td><td>9.5</td><td>10.225</td><td>10.961</td><td>11.713</td><td>12.438</td></tr><tr><td>F<sub>Z</sub>(z<sup>*</sup><sub>CFE</sub>)</td><td>0.051</td><td>0.120</td><td>0.227</td><td>0.358</td><td>0.5</td><td>0.642</td><td>0.773</td><td>0.880</td><td>0.949</td></tr><tr><td>abs.error q</td><td>-0.049</td><td>-0.080</td><td>-0.073</td><td>-0.042</td><td>0.0</td><td>0.042</td><td>0.073</td><td>0.080</td><td>0.049</td></tr><tr><td>PERT z<sup>*</sup></td><td>8.112</td><td>8.588</td><td>8.932</td><td>9.226</td><td>9.5</td><td>9.774</td><td>10.068</td><td>10.412</td><td>10.888</td></tr><tr><td rowspan="5">4</td><td>Real z<sup>*</sup></td><td>7.107</td><td>7.868</td><td>8.462</td><td>8.993</td><td>9.5</td><td>10.007</td><td>10.538</td><td>11.132</td><td>11.893</td></tr><tr><td>z<sup>*</sup><sub>CFE</sub></td><td>6.562</td><td>7.287</td><td>8.039</td><td>8.775</td><td>9.5</td><td>10.225</td><td>10.961</td><td>11.713</td><td>12.438</td></tr><tr><td>F<sub>Z</sub>(z<sup>*</sup><sub>CFE</sub>)</td><td>0.051</td><td>0.120</td><td>0.227</td><td>0.358</td><td>0.5</td><td>0.642</td><td>0.773</td><td>0.880</td><td>0.949</td></tr><tr><td>abs.error q</td><td>-0.049</td><td>-0.080</td><td>-0.073</td><td>-0.042</td><td>0.0</td><td>0.042</td><td>0.073</td><td>0.080</td><td>0.049</td></tr><tr><td>PERT z<sup>*</sup></td><td>8.112</td><td>8.588</td><td>8.932</td><td>9.226</td><td>9.5</td><td>9.774</td><td>10.068</td><td>10.412</td><td>10.888</td></tr><tr><td rowspan="5">4</td><td>Real z<sup>*</sup></td><td>7.107</td><td>7.868</td><td>8.462</td><td>8.993</td><td>9.5</td><td>10.007</td><td>10.538</td><td>11.132</td><td>11.893</td></tr><tr><td>z<sup>*</sup><sub>CFE</sub></td><td>6.562</td><td>7.287</td><td>8.039</td><td>8.775</td><td>9.5</td><td>10.225</td><td>10.961</td><td>11.713</td><td>12.438</td></tr><tr><td>F<sub>Z</sub>(z<sup>*</sup><sub>CFE</sub>)</td><td>0.051</td><td>0.120</td><td>0.227</td><td>0.358</td><td>0.5</td><td>0.642</td><td>0.773</td><td>0.880</td><td>0.949</td></tr><tr><td>abs.error q</td><td>-0.049</td><td>-0.080</td><td>-0.073</td><td>-0.042</td><td>0.0</td><td>0.042</td><td>0.073</td><td>0.080</td><td>0.049</td></tr><tr><td>PERT z<sup>*</sup></td><td>8.112</td><td>8.588</td><td>8.932</td><td>9.226</td><td>9.5</td><td>9.774</td><td>10.068</td><td>10.412</td><td>10.888</td></tr></tbody></table>
References

1. Malcolm, D.G. et al.: Application of a Technique for Research and Development Program Evaluation. Operations Research. 7, (1959). https://doi.org/10.1287/opre.7.5.646.

2. Roy, D., Roy, R.: Distribution of the Activity Time in Network Analysis: A Critical Revisit with a Gamma Alternative. Communications in Statistics - Simulation and Computation. 42(6), 1288-1297 (2013). https://doi.org/10.1080/03610918.2012.664231.

3. Spanos, A.: Probability Theory and Statistical Inference: Econometric Modeling with Observational Data. Cambridge University Press (1999).

4. Goman, M.: Towards unambiguous IT risk definition. In: Proceedings of the Central European Cybersecurity Conference 2018, CECC 2018, pp. 15:115:6 ACM (2018). https://doi.org/10.1145/3277570.3277586.

5. Udoumoh, E.F., Ebong, D.W.: A Review of Activity Time Distributions in Risk Analysis. American Journal of Operations Research. 07(06), 356-371 (2017). https://doi.org/10.4236/ajor.2017.76027.

6. Devroye, L.P.: Inequalities for the Completion Times of Stochastic PERT Networks. Mathematics of Operations Research. 4(4), 441-447 (1979). https://doi.org/10.1287/moor.4.4.441.

7. Dodin, B.: Bounding the Project Completion Time Distribution in PERT Networks. Operations Research. 33(4), 862-881 (1985). https://doi.org/10.1287/opre.33.4.862.

8. Elmaghraby, S.E.: The Estimation of Some Network Parameters in the Pert Model of Activity Networks: Review and Critique. In: Advances in Project Scheduling. pp. 371-432 Elsevier (1989). https://doi.org/10.1016/B978-0-444-87358-3.50021-3.

9. Kamburowski, J.: An upper bound on the expected completion time of PERT networks. European Journal of Operational Research. 21(2), 206-212 (1985). https://doi.org/10.1016/0377-2217(85)90032-3.

10. Kleindorfer, G.B.: Bounding Distributions for a Stochastic Acyclic Network. Operations Research. 19(7), 1586-1601 (1971). https://doi.org/10.1287/opre.19.7.1586.

11. Shogan, A.W.: Bounding distributions for a stochastic pert network. Networks. 7(4), 359-381 (1977). https://doi.org/10.1002/net.3230070407.

12. Grey, S.: Practical risk assessment for project management. Wiley, Chichester: New York (1995).

13. Johnson, D.: The triangular distribution as a proxy for the beta distribution in risk analysis. Journal of the Royal Statistical Society: Series D (The Statistician). 46(3), 387-398 (1997). https://doi.org/10.1111/1467-9884.00091.

14. Stuart, A., Ord, J.K.: Kendall's Advanced Theory of Statistics, Volume 1: Distribution Theory, London: Arnold (1994).

15. Holton, G.A.: Value-at-Risk: Theory and Practice, 2nd Edition, e-book (2014). https://www.value-at-risk.net

16. Cornish, E.A., Fisher, R.A.: Moments and cumulants in the specification of distributions, Review of the International Statistical Institute, 5, 307-320 (1937).

17. Weisstein, E.W.: Uniform Distribution. From MathWorld-A Wolfram Web Resource. http://mathworld.wolfram.com/UniformDistribution.html. Last accessed 24 Aug 2019

18. Sadooghi-Alvandi, S.M. et al.: On the distribution of the sum of independent uniform random variables. Statistical Papers. 50(1), 171-175 (2007). https://doi.org/10.1007/s00362-007-0049-4.