A novel numerical method for solution of fractional partial differential equations involving the $\psi$-Caputo fractional derivative

1.   Introduction

Statistical process monitoring (SPM) tools are a set of instruments utilized for monitoring the quality of products in manufacturing and production processes. In these instruments, the charts are more vital tools that detect changes in process parameter(s). Shewhart ^[1] introduced the basic charts, which are referred to as the Shewhart charts. The Shewhart charts for attribute data, such as $c$, $p$ and $u$ charts, are used to monitor the nonconforming events. Despite being user-friendly, the Shewhart attribute charts tend to be insensitive in monitoring nonconforming events, particularly in high-quality processes. High-quality processes consistently produce products or services that meet or exceed customer expectations while minimizing waste and maximizing efficiency and have very low numbers of nonconforming events. In such processes, apart from monitoring nonconforming events, an alternate strategy is to utilize a chart to monitor the inter-arrival time of nonconforming events, which may follow a skewed distribution such as the gamma distribution. Time between events (TBE) charts are typically preferred over the Shewhart charts for high-quality processes where nonconforming events rarely occur ^[2,3].

In literature, many studies have been conducted that investigate TBE charts. For example, Lucas ^[4] proposed the cumulative sum (CUSUM) chart for exponentially distributed TBE data. Later, Vardeman and Ray ^[5] used integral equations to compute the average run length (ARL) for the exponential CUSUM chart. Gan ^[6] introduced the exponentially weighted moving average (EWMA) chart based on the exponential distribution. Chan et al. ^[7] designed the Shewhart-type chart, i.e., the cumulative quantity control (CQC) chart, or $t$-chart, to monitor the exponentially distributed TBE. Xie et al. ^[8] designed the Shewhart‐type chart, i.e., the ${t}_{r}$-chart for monitoring the $r$th $(r\ge 1)$ failures in the gamma distribution. Pehlivan and Testik ^[9] investigated the robustness of the exponential EWMA chart, where the TBE mean is assumed to be known; but, TBE observations are modeled by Weibull or lognormal distribution. Chakraborty et al. ^[10] proposed the gamma distribution-based lower-sided generally weighted moving average (GWMA) chart, called the GWMA TBE chart, that monitors downward changes in TBE data. Alevizakos et al. ^[2] improved the performance of the GWMA TBE chart and proposed the lower-sided double GWMA chart, referred to as the DGWMA TBE chart. Alevizakos and Koukouvinos ^[11] designed the double EWMA chart, referred to as the DEWMA TBE chart, that outperforms the EWMA TBE chart in monitoring downward changes in the process. Alevizakos et al. ^[12] developed one- and two-sided triple EWMA chart to monitor the TBE data, referred to as the TEWMA TBE chart.

Roberts ^[13] introduced the traditional EWMA chart, which monitors the mean level of the process. In general, the EWMA chart is more sensitive than the Shewhart-type chart in identifying small changes in the process as it uses both current and previous process observations. However, the EWMA chart allocates relatively more weight to recent process observations than previous observations ^[14]. As a result, Abbas ^[15] proposed the homogeneously weighted moving average (HWMA) chart for monitoring process mean, which outperforms the EWMA chart. Abid et al. ^[16] extended the concept of the HWMA chart and proposed the double HWMA (DHWMA) chart. Riaz et al. ^[17] introduced the triple HWMA (THWMA) chart for the mean level of the process that improves the detection ability of the HWMA and DHWMA charts. Other charts based on the HWMA statistic are proposed in the studies of Rasheed et al. ^[18], Rasheed et al. ^[19] and Rasheed et al. ^[20].

As stated previously, Chakraborty et al. ^[10], Alevizakos et al. ^[2] and Alevizakos and Koukouvinos ^[11] introduced the one-sided GWMA TBE, DGWMA TBE and DEWMA TBE charts, respectively, for monitoring TBE observations following the gamma distribution. Aslam et al. ^[3] proposed the HWMA chart to monitor the TBE, referred to as the HWMA TBE chart. These charts are lower-sided charts that effectively monitor only the process's downward shifts (i.e., process deterioration). In practice, however, both downward shifts and upward shifts (i.e., process improvement) in the process must be monitored because corrective actions must be taken in both cases, either to diagnose and remove the causes of the process deterioration or to detect and maintain the reasons of the process improvements. As a result, in high-quality processes, a two-sided chart that monitors both downward and upward shifts at the same time is required. In this regard, Alevizakos et al. ^[12], Alevizakos and Koukouvinos ^[21] and Alevizakos and Koukouvinos ^[22] suggested the gamma distribution-based two-sided TEWMA TBE, provisional mean TBE (PM TBE) and double PM TBE (DPM TBE) charts, which monitor both the process's downward and upward shifts. These charts may indicate that the process performance needs to be improved further. To enhance the performance of TBE charts, this study proposes the design structures of one-sided (upper and lower) and two-sided THWMA charts with time-varying control limits that efficiently monitor the gamma-distributed TBE observations. The proposed charts are referred toas THWMA TBE charts. The lower-sided THWMA TBE chart is designed to detect downward shifts in the process, while the upper-sided THWMA TBE chart is constructed to monitor upwards changes in the process. Similarly, the two-sided THWMA TBE chart is developed to detect both downwards and upwards shifts in the process. Extensive Monte Carlo simulations are used to compute the numerical finding associated with the THWMA TBE charts. In parallel to THWMA TBE charts, the one- and two-sided competing TBE charts, such as the DHWMA TBE, HWMA TBE, TEWMA TBE, DEWMA TBE and EWMA TBE charts, are also constructed for comparison purposes. The proposed charts are compared to competing charts. The comparison of the THWMA TBE charts to the competing TBE charts suggests that the competing TBE charts perform better than the competing TBE charts in detecting changes in the process. Lastly, two applications, i.e., boring machine failure data and traumatic brain injury patients' hospital stay time data applications are provided to assess the THWMA TBE charts practically.

The layout for the rest of the article is as follows: Section 2 describes the methodologies of the gamma-distributed THWMA TBE charts, along with several competing charts. Section 3 offers the design and implementation of the THWMA TBE charts. Performance and comparison analysis of the THWMA TBE charts are provided in Section 4. Section 5 presents the two real-life applications of the THWMA TBE chart. The last section addresses the overall summary, conclusions and recommendations of the study.

2.   Gamma-distributed TBE charts

In this section, an introduction to gamma-distributed TBE observations is provided in Subsection 2.1. The construction of the one- and two-sided TBE, DEWMA TBE, TEWMA TBE, HWMA TBE and DHWMA TBE charts is outlined in Subsections 2.2 to 2.6, respectively. Likewise, Subsection 2.7 explains the methodology and formulation of the proposed one- and two-sided THWMA TBE charts.

2.1.   Gamma-distributed TBE observations

Let ${T}_{1}^{*}, {T}_{2}^{*}, \dots, {T}_{k}^{*}$ represent the interarrival times between two successive occurrences in a homogeneous Poisson process having a rate of ${\theta }^{-1}$($\theta > 0$). In this case, each ${T}_{i}^{*}$ is known as an independently and identically distributed (iid) exponential random variable with scale parameter θ, i.e., ${T}_{i}^{*}\stackrel{iid}{\sim }\mathrm{E}\mathrm{x}\mathrm{p}\mathrm{o}\mathrm{n}\mathrm{e}\mathrm{n}\mathrm{t}\mathrm{i}\mathrm{a}\mathrm{l}\left(\theta \right)$ for $i = \mathrm{1, 2}, \dots k$. If $T$ denotes the sum of interarrival times until the first $k$ failures occur, i.e., $T = \sum _{i = 1}^{k}{T}_{i}^{*}$, then $T$ follows the gamma distribution, i.e., $T~\mathrm{G}\mathrm{a}\mathrm{m}\mathrm{m}\mathrm{a}(k, \theta)$. The probability density function (pdf) of $T$, is given as follows:

where $\mathrm{\Gamma }( \cdot )$ is called the gamma function, $k$ (assumed known) is a shape parameter and $\theta$ is a scale parameter of the gamma distribution. The mean and variance of the gamma distribution are $k\theta$ and $k{\theta }^{2}$, respectively. The gamma distribution defined by Eq (1) is also referred to as the Erlang distribution. It is important to note that if $k = 1,$ then $\mathrm{g}\mathrm{a}\mathrm{m}\mathrm{m}\mathrm{a}(k, \theta)$ reduces to $\mathrm{e}\mathrm{x}\mathrm{p}\mathrm{o}\mathrm{n}\mathrm{e}\mathrm{n}\mathrm{t}\mathrm{i}\mathrm{a}\mathrm{l}\left(\theta \right)$. In this case, the interest only lies in monitoring the time until one failure occurs.

Let ${T}_{t}\stackrel{iid}{\sim }\mathrm{G}\mathrm{a}\mathrm{m}\mathrm{m}\mathrm{a}(k, \theta)$ at time $t = 1, 2, \dots$ having mean and variance $k\theta$ and $k{\theta }^{2}$, respectively. The process is deemed to be in control (IC) when $\theta = {\theta }_{0}$; otherwise, if $\theta = {\theta }_{1}$ then the process is said to be out-of-control (OOC). In the OOC process, if ${\theta }_{1} < {\theta }_{0}$ then the process is considered to have deteriorated; however, if ${\theta }_{1} > {\theta }_{0}$ then the process may have improved. In both scenarios, corrective action should be taken to diagnose (if any ${\theta }_{1} < {\theta }_{0}$) or preserve (if any ${\theta }_{1} > {\theta }_{0}$) the reasons for change.

2.2.   EWMA TBE chart

The EWMA TBE chart is defined by the charting statistic, ${X}_{t}$ and can be defined as follows:

where $\varphi$ is a smoothing or sensitivity parameter, such that $0 < \varphi \le 1$. The starting value of ${X}_{t}$ is assumed to be $k{\theta }_{0}$, i.e., ${X}_{0} = k{\theta }_{0}$. The mean and variance of ${X}_{t}$ for the IC process, respectively, are given as follows:

To detect the two-sided changes in the process, the control limits $UC{L}_{t}$ and $LC{L}_{t}$ for the two-sided EWMA TBE chart are given as follows:

where $L$ denotes the width coefficient for the EWMA TBE chart. The two-sided EWMA TBE chart and statistic ${X}_{t}$ is plotted versus time $t,$ and the process is declared to be OOC if ${X}_{t}\ge UC{L}_{t}$ or ${X}_{t}\le LC{L}_{t}$; otherwise, the process is declared to be IC.

2.3.   DEWMA TBE chart

The DEWMA TBE chart introduced by Alevizakos and Koukouvinos ^[11], can be defined by the plotting statistic ${Y}_{t}$, defined as follows:

where the starting value ${Y}_{0} = k{\theta }_{0}$. Based on Eq (5), the mean and variance of ${Z}_{t}$ for the IC process, respectively, are defined as follows:

Based on $E\left({Y}_{t}\right)$ and $Var\left({Y}_{t}\right)$, the two-sided DEWMA TBE hart control limits $UC{L}_{t}$ and $LC{L}_{t}$ are given as follows:

To monitor the decreasing or increasing shifts in the process, the two-sided DEWMA TBE chart is constructed by plotting the statistic ${Y}_{t}$ against time $t$ and if ${Y}_{t}\ge UC{L}_{t}$ or ${Y}_{t}\le LC{L}_{t}$, then the process is deemed to be OOC; otherwise, the process is considered to be IC.

2.4.   TEWMA TBE chart

Alevizakos et al. ^[12] proposed the TEWMA TBE chart, which has the charting statistic ${Z}_{t}$, defined as follows:

The initial value of ${Z}_{t}$ is equal to $k{\theta }_{0}$, i.e., ${Z}_{0} = k{\theta }_{0}$. The IC mean of the statistic ${Z}_{t}$ is $E\left({Z}_{t}\right) = k{\theta }_{0}$, while its variance is given as follows:

where $\tau = {\left(1-\varphi \right)}^{2}$. To detect the decreasing or increasing changes in the process, via TEWMA TBE chart, the control limits, i.e., $UC{L}_{t}$ and $LC{L}_{t}$ are defined as follows:

In this case, the TEWMA TBE statistic ${Z}_{t}$ is plotted versus time $t$, and if ${Z}_{t}\ge UC{L}_{t}$ or ${Z}_{t}\le LC{L}_{t}$ then the process is regarded as OOC; otherwise, the process is considered to be IC.

2.5.   HWMA TBE chart

The HWMA chart for TBE introduced by Aslam et al. ^[3]is referred to as the HWMA TBE chart. Let ${U}_{t}$ be the charting statistic, then it can be defined as follows:

where ${\stackrel{-}{T}}_{t-1} = \frac{{\sum }_{j = 1}^{t-1}{T}_{j}}{t-1}$ is the mean of previous $t-1$ observations. The starting values ${U}_{0}$ and ${\stackrel{-}{T}}_{0}$ are set to $k{\theta }_{0}$, i.e., ${U}_{0} = {\stackrel{-}{T}}_{0} = k{\theta }_{0}$. The IC mean and variance of ${U}_{t}$ (cf. Appendices A, and B) are provided as follows:

The two-sided HWMA TBE chart control limits, i.e., $UC{L}_{t}$ and $LC{L}_{t}$ can be given as follows:

For monitoring the downward or upward shifts in the process, the two-sided HWMA TBE chart is formulated and the statistic ${U}_{t}$ is plotted versus time $t$. Whenever ${U}_{t}\ge UC{L}_{t}$ or ${U}_{t}\le LC{L}_{t}$, the process is considered to be OOC; otherwise, the process is deemed to be IC.

2.6.   DHWMA TBE chart

The gamma distribution-based DHWMA charts for monitoring TBE may is referred to as the DHWMA TBE charts. The DHWMA TBE statistic is symbolized by ${V}_{t}$, and is defined as follows:

The starting value of the statistic ${V}_{t}$ is equal to the IC $TBE$ mean, i.e., ${V}_{0} = k{\theta }_{0}$. The IC process mean and variance of the statistic ${V}_{t}$, (cf. Appendices A, B and C), are given as follows:

The control limits $UC{L}_{t}$ and $LC{L}_{t}$ of the two-sided DHWMA TBE chart to identify the downward or upward shifts, are given as follows:

The statistic ${V}_{t}$ is plotted against time $t$, and the process is said to be OOC it ${V}_{t}\ge UC{L}_{t}$ or ${V}_{t}\le LC{L}_{t}$; otherwise, the process is said to be IC.

2.7.   THWMA TBE chart

Riaz et al. ^[17] designed the THEMA chart for monitoring the mean level of the process. The methodology of the THWMA chart can be extended to design the develop one- and two-sided THWMA charts to monitor the gamma-distributed TBE observations. The proposed charts are called the THWMA TBE charts. If W_t is the charting statistic for the THWMA TBE chart, then it can be expressed by linear equations, defined as follows:

with starting value ${W}_{0} = k{\theta }_{0}$. The IC mean and variance of the statistic ${W}_{t}$, (cf. Appendices A, B, C and D) are given as follows:

To identify the decrease in interarrival time, the limit $LC{L}_{t}$ of the lower-sided THWMA TBE chart is defined as follows:

The statistic ${W}_{t}$ is plotted against time $t$, and if ${W}_{t}\le LC{L}_{t}$, the THWMA TBE chart states the process is OOC; otherwise, the process is IC. Similarly, to monitor the increase in interarrival time, the limit $UC{L}_{t}$ of the upper-sided THWMA TBE chart is given as follows:

The charting statistic ${W}_{t}$ is plotted against time $t$, and if ${W}_{t}\ge UC{L}_{t}$ the process under the THWMA TBE chart is OOC; otherwise, the process is IC. Likewise, to monitor the decrease or increase in interarrival time, the two-sided THWMA TBE chart control limits, i.e., $UC{L}_{t}$ and $LC{L}_{t}$, are presented as follows:

In this case, the charting statistic ${W}_{t}$ is plotted against $t$, and the process under the two-sided THWMA TBE chart is considered to be OOC when ${W}_{t}\ge UC{L}_{t}$ or ${W}_{t}\le LC{L}_{t}$; otherwise, the process is deemed to be IC.

3.   Design and implementation of THWMA TBE chart

The design structure of the THWMA TBE chart is based on the smoothing parameters $\varphi$ and chart width $L$. The combinations of $\varphi$ and $L$ are chosen so that the $AR{L}_{0}$ is close to the desired value. The $AR{L}_{0}$ is set to 370 in this study. Furthermore, THWMA TBE charts are used with time-varying control limits.

3.1.   Performance evaluation measures

The performance of the chart is often evaluated using run-length metrics such as ARL, SDRL, percentile points, etc. The ARL is referred to as the average number of sample points depicted on the chart until a point reveals an OOC signal ^[23]. ARLs are classified into two types: IC ARL $\left(AR{L}_{0}\right)$ and OOC ARL $\left(AR{L}_{1}\right)$. When a process is operating in an IC state, then the $AR{L}_{0}$ is anticipated to be large enough in order to avert false alarms, whereas when a process is functioning in an OOC state, then $AR{L}_{1}$ should remain small to identify the process change quickly ^[24]. To compare the efficiency of two or more charts, it is advisable to set the common $AR{L}_{0}$ for these charts. For a particular shift, the chart with the smaller $AR{L}_{1}$ value is deemed to be sensitive and may detect a specific shift faster than the other charts.

3.2.   Overall performance measures

The $ARL$ measures assess the performance of the chart for a single specified shift in the process parameters. Extra quadratic loss (EQL), relative average run length (RARL) and performance comparison index (PCI) are additional performance measures that evaluate a chart's overall performance over a certain range of shifts.

The EQL is the weighted ARL over the range of shifts, i.e., ${\delta }_{min}$ to ${\delta }_{max}$, using the square of the shift as a weight ^[25]. Mathematically, EQL can be defined as follows:

where ARL $\left(\delta \right)$ is the ARL of a particular chart at a shift $\delta$. A chart with a lower EQL value is regarded as more efficient.

The RARL is referred to as the ratio of the ARL of a certain chart, i.e., ARL $\left(\delta \right)$ to the ARL benchmark chart, i.e., ${ARL}_{bmk}\left(\delta \right)$ ^[26]. The RARL may be defined mathematically, as follows:

A benchmark chart might be considered as a chart with a minimum EQL or as one of the current standard charts.

According to Ou et al. ^[27], PCI is defined as the ratio of a chart's EQL to the EQL of the benchmark chart. It is mathematically given as follows:

The PCI value is equal to one for the benchmark chart and it is more than one for the other charts.

3.3.   Implementation of THWMA TBE chart

To implement the THWMA TBE chart, a simulation study is carried out. A computational methodology of Monte Carlo simulations is used as the design structure of the THWMA TBE chart and is more complicated than that of other charts. Furthermore, when time-varying control limits are considered, computational techniques such as the Markov chain and integral equations have been used limitedly.

To execute the simulation study, various combinations of $k = 1, 2, 3$ and $\varphi = 0.05, \mathrm{ }\mathrm{ }\mathrm{ }\mathrm{ }0.1, \mathrm{ }0.2$ are used to compute the chart coefficients values, i.e., $L$, so that $AR{L}_{0}$ is close to 370. For the sake of convenience, if the process is IC, then it is assumed that $\theta = {\theta }_{0} = 1$; however, if the process is OOC, then $\theta = {\theta }_{1} = \delta {\theta }_{0}$, where $\delta =$ 0.975, 0.95, 0.925, 0.9, 0.85, 0.8, 0.7, 0.5, 0.25 for the downward changes and $\delta =$ 1.025, 1.05, 1.08, 1.11, 1.18, 1.25, 1.43, 2, 4 for upward shifts. A simulation algorithm is developed in $\mathrm{R}$ software and can be described in the steps given as follows:

(ⅰ) Generate ${T}_{t}$ from $\mathrm{g}\mathrm{a}\mathrm{m}\mathrm{m}\mathrm{a}(k, \theta)$, i.e., ${T}_{t}~\mathrm{ }\mathrm{g}\mathrm{a}\mathrm{m}\mathrm{m}\mathrm{a}\left(k, {\theta }_{0}\right), t = 1, 2, \dots$ for each $k = 1, 2, 3$.

(ⅱ) Compute the statistic ${U}_{t}$ from Eq (11) using ${T}_{t}$.

(ⅲ) Using the statistic ${U}_{t}$ as input, compute the statistic ${V}_{t}$ in Eq (14).

(ⅳ) Using the statistic ${V}_{t}$ as input, compute the 𝑇𝐻𝑊𝑀𝐴 TBE statistic ${W}_{t}$ in Eq (17)

(ⅴ) Choose the design parameters $\left(\varphi, L\right)$ for the desired ${ARL}_{0}$.

(ⅵ) Calculate ${UCL}_{t}$ and ${LCL}_{t}$ in Eq (21) based on $L$ and $\varphi$.

(ⅶ) Plot the 𝑇𝐻𝑊𝑀𝐴 TBE statistic ${W}_{t}$ along ${UCL}_{t}$ and ${LCL}_{t}$.

(ⅷ) When ${W}_{t}\ge {UCL}_{t}$ or ${W}_{t}\le {LCL}_{t}$, the sequence order called run length ($\mathrm{R}\mathrm{L}$) is recorded.

(ⅸ) Replicate steps (ⅰ)-(ⅷ) 10⁶ times and record $\mathrm{R}\mathrm{L}\mathrm{s}$.

(ⅹ) Calculate the average and standard deviation of 10⁶ $\mathrm{R}\mathrm{L}$, i.e., ${ARL}_{0}$ and ${SDRL}_{0}$. Check whether the calculated ${ARL}_{0}$ is the desired ${ARL}_{0}$; otherwise, adjust $L$ and repeat steps (i)-(ⅸ) until the desired ${ARL}_{0}$ is achieved.

(ⅹⅰ) Compute ${ARL}_{1}$ and ${SDRL}_{1}$ (OOC SDRL) values by considering ${T}_{t}~\mathrm{g}\mathrm{a}\mathrm{m}\mathrm{m}\mathrm{a}\left(k, {\theta }_{1}\right), t = 1, \mathrm{2, 3}, \dots$, where ${\theta }_{1} = \delta {\theta }_{0}$ with $\delta \ne 1$ and replicate steps (ⅱ)-(ⅸ).

It should be noted that the ARL and SDRL values may be determined using their respective expressions, which are given as follows:

where $N$ is the number of simulation runs, which is specified to be in this study as $N = {10}^{6}$.

The above simulation algorithm is developed to compute the ARL and SDRL values for the two-sided THWMA TBE chart. However, to calculate the ARL and SDRL values for the one-sided chart, the control limit is specified in Eqs (19) and (20) are, respectively, required for lower- and upper-sided THWMA TBE charts. Moreover, a similar algorithm can also be used to determine the ARL and SDRL findings for the competing charts. The simulation results, i.e., the ARL and SDRL (specified in the parenthesis) values for the THWMA TBE charts, along with the competing charts, are shown in Tables 1–6. For each shift, the smallest ${ARL}_{1}$ values are shown in bold print.

4.   Performance and comparative study

This section covers the performance and comparative study for the THWMA TBE charts. The THWMA TBE charts are compared to competing charts, including the DHWMA TBE, HWMA TBE, TEWMA TBE, DEWMA TBE and EWMA TBE charts. Details of the performance and comparative analysis are as follows:

From Tables 1–6, it can be concluded that for a given value of $\varphi$, the performance behavior of the THWMA TBE charts is enhanced as the value of $k$ rises. For example, if $\varphi =$ 0.05 and $\delta =$ 0.95, then the lower-sided THWMA TBE chart provides ${ARL}_{1} =$ 22.7 when $k =$ 1, while it yields ${ARL}_{1} =$ 15.63 when $k =$ 2 and it has ${ARL}_{1} =$ 12.93 when $k =$ 3. A similar type of behavior can be observed with the upper- and two-sided THWMA TBE charts. In addition, for a specified $k$, the performance behavior of the two-sided THWMA TBE chart worsens when the values of $\varphi$ increase. For instance, when $k =$ 1 and $\delta =$ 0.9, then the two-sided THWMA TBE chart achieves ${ARL}_{1} =$ 40.01 with $\varphi =$ 0.05, whereas it holds ${ARL}_{1} =$ 48.88 with $\varphi =$ 0.1, and it receives ${ARL}_{1} =$ 123.12 with $\varphi =$ 0.2. Similarly, for $k =$ 1, 2. The performance of the upper-sided THWMA TBE chart deteriorates with increasing values of $\varphi$. For instance, at $k =$ 1 and $\delta =$ 1.025, for $\varphi =$ 0.05 the upper-sided THWMA TBE chart has ${ARL}_{1} =$ 54.62, for $\varphi =$ 0.1, it has ${ARL}_{1} =$ 59.52 and for $\varphi =$ 0.2, it has ${ARL}_{1} =$ 185.03. However, this is not observed for the lower-sided THWMA TBE chart with any $k$. For instance, when $k =$ 2 and $\delta =$ 0.925, the lower-sided THWMA TBE chart has ${ARL}_{1} =$ 8.24 if $\varphi =$ 0.05, ${ARL}_{1} =$ 8.02 if $\varphi =$ 0.1 and ${ARL}_{1} =$ 12.58 if $\varphi =$ 0.2. Furthermore, for a given value of $k$ and $\varphi$, the lower-sided THWMA TBE chart has better detection ability than the upper-sided THWMA TBE chart. For example, if $k =$ 3 and $\varphi =$ 0.1, the lower-sided THWMA TBE chart has ${ARL}_{1} =$ 12.61 when $\delta =$ 0.95 (1.05 time less than ${\theta }_{0}$), whereas the upper-sided THWMA TBE chart has ${ARL}_{1} =$ 15.15 when $\delta =$ 1.05 (1.05 time greater than ${\theta }_{0}$). Also, for $k = 1$, the two-sided THWMA TBE chart is just $ARL$-biased for very small downward shifts, i.e., $0.975\le \delta < 1$ when $\varphi = 0.2$. Moreover, as far as the SDRL values are concerned, it can be observed that when the value of $\varphi$ rises, the IC SDRL (${SDRL}_{0}$) values reduce. For example, when $k = 3$, then the two-sided THWMA TBE chart has ${SDRL}_{0} = 3368.86$ when $\varphi = 0.05$, though the two-sided THWMA TBE chart has ${SDRL}_{0} = 2444.97$ when $\varphi = 0.1$ and the two-sided THWMA TBE chart has ${SDRL}_{0} = 677.96$ when $\varphi = 0.2$. Similar inferences may be drawn for the one-sided THWMA TBE charts. It should be noted that the results regarding the competing charts, i.e., the DHWMA TBE, HWMA TBE, TEWMA TBE, DEWMA TBE, and EWMA TBE charts, those observed along with the THWMA TBE chart can be discussed in the same manner.

In comparing the THWMA TBE charts against the competing TBE charts, it is observed that the THWMA TBE charts gain excellent performance against the corresponding EWMA TBE charts for small and large shifts. For example, if $k = 1$, $\varphi = 0.1$ and $\delta =$ 0.975, 0.95, …, 0.25, then the lower-sided THWMA TBE chart has ${ARL}_{1} =$ 48.61, 21.54, …, 1.11, while the lower-sided EWMA TBE has ${ARL}_{1} =$ 296.22,237.58, …, 7.22. The THWMA TBE charts reveal superior performance over the DHWMA TBE charts for every shift value. For instance, when $k = 2$, $\varphi = 0.05$ and $\delta =$ 1.025, 1.05, …, 4, the upper-sided THWMA TBE chart has ${ARL}_{1} =$ 37.41, 18.15, …, 1.22, whereas the upper-sided DHWMA TBE has ${ARL}_{1} =$ 225.87,148.84, … 1.48. The THWMA TBE charts show better detection ability relative to the TEWMA TBE charts for each shift value. For example, when $k = 3$, $\varphi = 0.05$ and $\delta =$ 0.25, 0.5, …, 4, the two-sided THWMA TBE chart has ${ARL}_{1} =$ 4.06, 6.90, …, 2.26, but the upper-sided DHWMA TBE has ${ARL}_{1} =$ 8.51, 18.46, …, 2.77. The THWMA TBE charts achieve better detection ability over the TEWMA TBE charts for each shift value. For example, when $k = 3$, $\varphi = 0.05$ and $\delta =$ 0.25, 0.5, …, 4, the two-sided THWMA TBE chart has ${ARL}_{1} =$ 4.06, 6.90, …, 2.26, where the upper-sided DHWMA TBE has ${ARL}_{1} =$ 8.51, 18.46, …, 2.77. The THWMA TBE charts have better performance compared with the respective 𝐻𝑊𝑀𝐴 TBE charts for every single shift. For instance, if $k = 1$, $\varphi = 0.2$ and $\delta =$ 0.975, 0.95, …, 0.25, then the lower-sided THWMA TBE chart has ${ARL}_{1} =$ 58.14, 27.88, …, 1.17; however, the upper-sided HWMA TBE has ${ARL}_{1} =$ 105.69, 58.09, …, 2.43. The THWMA TBE charts have better performance compared with the DHWMA TBE charts for every shift. For example, when $k = 2$, $\varphi = 0.1$ and $\delta =$ 0.25, 0.5, …, 4, the ${ARL}_{1}$ values of the two-sided THWMA TBE chart are 4.57, 8.66, …, 1.58, on the other hand, the ${ARL}_{1}$ values the upper-sided DHWMA TBE are 6.85, 14.96, …, 1.72. Similar types of conclusions can also be drawn from the $ARL$ plots of the one- and two-sided THWMA TBE versus one- and two-sided competing charts, given in the Figures 1 to 3.

The overall comparisons of the THWMA TBE charts against the rival TBE charts indicate that the THWMA TBE charts show better overall performance as compared to the EWMA TBE charts in term of minimum $EQL$, $RARL$ and $PCI$ values. For example, when $k = 1$ and $\varphi = 0.05,$ then the upper-sided THWMA TBE chart has $EQL = 17.35$, $RARL = 1.00$ and $PCI = 1.00$, whereas the upper-sided EWMA TBE chart has $EQL = 47.59$, $RARL = 2.74$ and $PCI =$ 3.35. The THWMA TBE charts demonstrate superior overall performance relative to the DEWMA TBE charts. For instance, when $k = 2$ and $\varphi = 0.1,$ the lower-sided THWMA TBE chart has $EQL = 9.08$, $RARL = 1.00$ and $PCI = 1.00$; whereas, the upper-sided DHWMA TBE chart has $EQL = 34.68$, $RARL = 3.82$ and $PCI =$ 7.53. The THWMA TBE charts reveal an improved overall performance against the TEWMA TBE charts. For instance, if $k = 3$ and $\varphi = 0.05,$ then the two-sided THWMA TBE chart has $EQL = 16.35$, $RARL = 1.00$ and $PCI = 1.00$, where the two-sided TEWMA TBE chart has $EQL = 30.33$, $RARL = 1.86$ and $PCI =$ 1.88. The THWMA TBE charts indicate dominant overall performance against the HWMA TBE charts. For instance, when $k = 2$ and $\varphi = 0.2,$ the lower-sided THWMA TBE chart has $EQL = 10.35$, $RARL = 1.00$ and $PCI = 1.00$, but the lower-sided HWMA TBE chart has $EQL = 18.19$, $RARL = 1.76$ and $PCI =$ 2.76. The THWMA TBE charts have enhanced overall performance as compared to the DHWMA TBE charts. For example, for $k = 1$ and $\varphi = 0.1,$ the lower-sided THWMA TBE chart has $EQL = 10.35$, $RARL = 1.00$ and $PCI = 1.00$; in contrast, the lower-sided DHWMA TBE chart has $EQL = 18.19$, $RARL = 1.76$ and $PCI =$ 2.76.

5.   Applications of THWMA TBE chart

This section provides two real-life examples to demonstrate how THWMA TBE charts can be used practically. Subsection 5.1 offers Example 5.1, which analyzes real-life data of a vertical boring machine's failure time. Additionally, the analysis of the real-life data representing the hospital stay time of the traumatic brain injury patients is given in Subsection 5.2.

5.1.   Example 5.1

Boring machines are extensively utilized in boring, drilling, milling, thread cutting, face-turning processes, etc. These machines may be horizontal or vertical, depending on the axis. The boring machines use a spinning tool known as a machine spindle, which generates the basic motion. By enlarging the existing holes in the workpiece, these machines generate a smooth and precise hole diameter. To accomplish this, boring machines employed a cutter, boring rod, drill or milling head to bore, thread, ream and mill surfaces, among other operations. A single steel cutting tip, a diamond, cemented carbide or a tiny grinding wheel can be used as the boring instrument. Special boring machines with many spindles are often utilized in large production plants ^[28].

In this example, the dataset is analyzed that reports the vertical boring machine failure times. This dataset was also utilized by Hossain et al. ^[29] and Hossain et al. ^[30] to construct various chart schemes. The TBE charts are not only implemented to monitor the processes of components or systems ^[8] in addition to being utilized to monitor the boring machine failure process. All possible actions must be taken to maintain the capability of the boring process to reduce the expense and hardship of boring machine failure.

Table 7 shows the data set, which consists of 32 vertical boring machine failure time (in hours) observations. Before implementing the THWMA TBE charts, the Kolmogorov-Smirnov ($\mathrm{K}\mathrm{S}$) test is used to test the goodness-of-fit for the gamma distribution. The $\text{KS}$ test provides a $p$-value of 0.3053, which confirms that the boring machine failure time follows the gamma distribution at the parametric values $k =$ 3 and $\theta$ = 959.3. To implement the THWMA TBE charts, the methodology of Alevizakos and Koukouvinos ^[11] is followed. The IC value ${\theta }_{0}$ of the process is considered to be 1200 for the downward shift, and for the upward shift, it is assumed to be 800. As a result, the shifts $\delta$ for downward and upward cases are approximately equal to 0.8 and 1.2, respectively. Because the earlier developments assume ${\theta }_{0} = 1$, so each observation in Table 7 is divided by ${\theta }_{0}$ to rescale the failure time. Setting ${ARL}_{0} = 370$, the lower-sided THWMA TBE, DHWMA TBE, HWMA TBE, TEWMA TBE, DEWMA TBE and EWMA TBE charts are constructed using the design parameters for each chart, respectively chosen as; $\left(\varphi, L\right) = \left(0.2, 0.453\right)$, $\left(\varphi, L\right) = \left(0.2, 1.027\right)$, $\left(\varphi, L\right) = \left(0.2, 1.498\right)$, $\left(\varphi, L\right) = \left(0.2, 1.885\right)$, $\left(\varphi, L\right) = \left(0.2, 2.027\right)$ and $\left(\varphi, L\right) = \left(0.2, 2.117\right)$ (see Table 5). Similarly, the upper-sided THWMA TBE, DHWMA TBE, HWMA TBE, TEWMA TBE, DEWMA TBE and EWMA TBE charts are formulated, based on the design parameters, respectively, given as; $\left(\varphi, L\right) = \left(0.2, 0.794\right)$, $\left(\varphi, L\right) = \left(0.2, 2.211\right)$, $\left(\varphi, L\right) = \left(0.2, 3.819\right)$, $\left(\varphi, L\right) = \left(0.2, 2.168\right)$, $\left(\varphi, L\right) = \left(0.2, 2.424\right)$ and $\left(\varphi, L\right) = \left(0.2, 3.128\right)$ (see Table 5). Likewise, the two-sided THWMA TBE, DHWMA TBE, HWMA TBE, TEWMA TBE, DEWMA TBE and EWMA TBE charts are developed, which used the design parameters provided as; $\left(\varphi, L\right) = \left(0.05, 1.226\right)$, $\left(\varphi, L\right) = \left(0.05, 1.342\right)$, $\left(\varphi, L\right) = \left(0.05, 2.261\right)$, $\left(\varphi, L\right) = \left(0.05, 1.742\right)$, $\left(\varphi, L\right) = \left(0.05, 1.945\right)$ and $\left(\varphi, L\right) = \left(0.05, 2.524\right)$, respectively (see Table 6). The lower, upper and two-sided THWMA TBE, DHWMA TBE, HWMA TBE, TEWMA TBE, DEWMA TBE and EWMA TBE charts are displayed in Figures 4–6, respectively.

Figure 4 reveals that the lower-sided THWMA TBE chart detects the first OOC point at failure number 7, while the lower-sided DHWMA TBE, HWMA TBE and TEWMA TBE charts diagnose the first OOC point at failure numbers 11, 23 and 20, respectively; however, the lower-sided DHWMA TBE and EWMA TBE charts fail to identify any OOC point. Similarly, Figure 5 indicates that the upper-sided THWMA TBE chart detects the first OOC point at failure number 3, while the upper-sided DHWMA TBE and TEWMA TBE charts identify the first OOC point at failure numbers 27 and 3, respectively, and the upper-sided HWMA TBE and EWMA TBE charts remain insensitive and do not detect any OOC signal. Likewise, Figure 6 suggests that with decreasing shift, the two-sided THWMA TBE chart provides the first OOC signal at failure number 16, while the two-sided DHWMA TBE chart triggers the first OOC point at failure number 19, but the HWMA TBE, TEWMA TBE, DEWMA TBE and EWMA TBE charts detect no OOC point. Equally, Figure 6 also shows that, with increasing shift, both the two-sided THWMA TBE and DHWMA TBE charts produce the first OOC point at failure number 4, whereas both the two-sided HWMA TBE and DEWMA TBE charts diagnose the OOC point at failure number 27 and the TEWMA TBE and EWMA TBE charts identify no OOC signal. The above discussion suggests that the THWMA TBE charts are more efficient in detecting changes in the process than the DHWMA TBE, HWMA TBE, TEWMA TBE, DEWMA TBE and EWMA TBE charts.

5.2.   Example 5.2

Trauma is a multifaceted psychiatric state that is the third leading source of mortality globally, after only cardiovascular disease and cancer. As a result of traumatic injury, over sixty million individuals are believed to be affected by trauma every year, with approximately 16,000 people dying each day. In trauma epidemiology statistics, traumatic brain injury (TBI) is one of the main roots of death and disability ^[31]. The vast majority of TBI patients are young and financially active adults who have a higher risk of being involved in road accidents ^[32]. As a result, TBI is viewed as a public health issue that results in large healthcare costs as well as many economic and social impacts. Therefore, a monitoring scheme is needed to assist in developing public health policies, reduce the TBI burden and enhance survival. In the current example, the one- and two-sided THWMA TBE with the competing charts, i.e., DHWMA TBE, HWMA TBE, TEWMA TBE, DEWMA TBE and EWMA TBE charts are applied to the data set of male patients suffering from TBI related to car accidents, hospitalized in the Emergency Department of Ribeirão Preto Medical School, University of São Paulo, Brazil.

The data set considered by Ramos et al. ^[33] representing the hospital stay time for male patients suffering from TBI related to car accidents is offered in Table 8. The KS goodness-of-fit test is applied, which provides a $p$-value of 0.918, confirming that the hospital stay times of the TBI patients follow the gamma distribution having parameters $k =$ 1 and $\theta$ = 32.5. In accordance with Alevizakos et al. ^[2], the IC process value ${\theta }_{0}$ is supposed to be 50 for the decrease in interarrival time; however, for the increase in interarrival time, it is considered to be 15. As a result, the data set may be regarded as the OOC process observations having a decreasing shift of size 0.65 and an increasing shift of size 2.17. All the previous developments are executed at ${\theta }_{0} = 1$, thus each observation of the data set is divided by the IC value ${\theta }_{0}$, i.e., ${\theta }_{0}$ = 50 for the downward shift and ${\theta }_{0}$ = 15.

The THWMA TBE, DHWMA TBE, HWMA TBE, TEWMA TBE, DEWMA TBE and EWMA TBE charts are constructed with a common value of ${ARL}_{0}$, i.e., ${ARL}_{0}\cong 370$. The design parameters; $\left(\varphi, L\right) = \left(0.2, 0.365\right)$, $\left(\varphi, L\right) = \left(\mathrm{0.2, 0.727}\right)$, $\left(\varphi, L\right) = \left(0.2, 0.9537\right)$, $\left(\varphi, L\right) = \left(0.2, 1.7765\right)$, $\left(\varphi, L\right) = \left(0.2, 1.8724\right)$ and $\left(\varphi, L\right) = \left(0.05, 1.813\right)$ are, respectively, chosen to formulate the lower-sided THWMA TBE, DHWMA TBE, HWMA TBE, TEWMA TBE, DEWMA TBE and EWMA TBE charts. The lower-sided charts are shown in Figure 7, which suggests that the lower-sided THWMA TBE chart diagnoses the first OOC signal after patient number 1, whereas the lower-sided 𝐷 DHWMA TBE, HWMA TBE, TEWMA TBE and DEWMA TBE charts, respectively, diagnose the first OOC point after patient numbers 2, 2, 10 and 13, but the EWMA TBE chart detects no OOC point. In addition, Figure 7 also indicates that the THWMA TBE chart overall detects 18 OOC signals, while the DHWMA TBE, HWMA TBE, TEWMA TBE and DEWMA TBE charts trigger 17, 17, 9 and 2 OOC signals, respectively.

Similarly, the upper-sided THWMA TBE, DHWMA TBE, HWMA TBE, TEWMA TBE, DEWMA TBE and EWMA TBE charts are developed using the design parameters, respectively, given as; $\left(\varphi, L\right) = \left(0.2, 0.878\right)$, $\left(\varphi, L\right) = \left(0.2, 2.608\right)$, $\left(\varphi, L\right) = \left(0.2, 4.601\right)$, $\left(\varphi, L\right) = \left(0.2, 2.258\right)$, $\left(\varphi, L\right) = \left(0.2, 2.548\right)$ and $\left(\varphi, L\right) = \left(0.2, 3.498\right)$. The upper-sided charts are depicted in Figure 8, which indicates that the upper-sided THWMA TBE chart identifies the first OOC point at patient number 7, while the upper-sided DHWMA TBE, TEWMA TBE, DEWMA TBE and EWMA TBE charts, respectively, trigger the first OOC point after patient numbers 15, 15, 16 and 15, but the HWMA TBE chart fails to detect any OOC signal. Furthermore, Figure 8 also demonstrates that the THWMA TBE chart diagnoses 14 OOC points in total; however, the DHWMA TBE, TEWMA TBE, DEWMA TBE and EWMA TBE charts detect 4, 4, 3 and 3 OOC points, respectively.

Likewise, the two-sided THWMA TBE, DHWMA TBE, HWMA TBE, TEWMA TBE, DEWMA TBE and EWMA TBE charts are constructed, which are based on the design parameters; $\left(\varphi, L\right) = \left(0.05, 1.196\right)$, $\left(\varphi, L\right) = \left(0.05, 1.342\right)$, $\left(\varphi, L\right) = \left(0.05, 2.998\right)$, $\left(\varphi, L\right) = \left(0.05, 1.724\right)$, $\left(\varphi, L\right) = \left(0.05, 1.9264\right)$ and $\left(\varphi, L\right) = \left(0.05, 2.536\right)$. The two-sided charts are portrayed in Figure 9, which demonstrates that both two-sided THWMA TBE and DHWMA TBE charts trigger the first OOC point after patient numbers 4 and 7 with downward and upward shift, respectively; however, the two-sided DEWMA TBE chart identifies the first OOC signal after patient number 13 with downward and after patient number 15 with upward shifts, respectively. On the other hand, the two-sided HWMA TBE and EWMA TBE charts diagnose the first OOC point after patient numbers 16 and 15, respectively, with an upward shift only, while the two-sided TEWMA TBE chart is unable to detect the OOC signal. Moreover, in overall detection, Figure 9 also shows that the two-sided THWMA TBE and DHWMA TBE charts detect 15 OOC signals with decreasing shift, whereas the two-sided DEWMA TBE chart trigger 6 OOC points. On the contrary, as the increasing shift is concerned, the two-sided THWMA TBE and DHWMA TBE charts detect 11 OOC points, whereas the two-sided HWMA TBE, DEWMA TBE and EWMA TBE charts identify 3, 4 and 3 OOC points, respectively. The above discussion suggests that, in terms of the ability to diagnose the changes in the process, the THWMA TBE appears to be more sensitive than other charts, i.e., the DHWMA TBE, HWMA TBE, TEWMA TBE, 𝐷DEWMA TBE, and EWMA TBE charts.

6.   Conclusions, summary and recommendations

This paper proposes triple homogeneously weighted moving average charts to monitor TBE observations following by the gamma distribution. The proposed charts are referred to as THWMA TBE charts. THWMA TBE charts are one- and two-sided charts. The one-sided THWMA TBE charts comprise the lower- and upper-sided THWMA TBE charts, which monitor the process's downward and upward shifts, respectively. Similarly, the two-sided THWMA TBE chart detects both the downward and upward changes in the process simultaneously. The Monte Carlo simulations methodology is used to compute the numerical findings for the performance measures, such as average run length (ARL), standard deviation run length (SDRL), extra quadratic loss (EQL), relative average run length ($RARL$) and performance comparison index (PCI). These numerical findings are obtained for the THWMA TBE charts, as well as for competing charts including the DHWMA TBE, HWMA TBE, TEWMA TBE, DEWMA TBE and EWMA TBE charts. A comparison study, based on aforementioned performance measures, is carried out and the THWMA TBE charts are compared to the competing charts. The comparison suggests that the THWMA TBE charts outperform the competing DHWMA TBE, HWMA TBE, TEWMA TBE, DEWMA TBE and EWMA TBE charts. In addition, two real-life data sets, i.e., one representing vertical boring machine failure time, and the other representing the hospital stay time for traumatic brain injury patients, are analyzed to implement the one- and two-sided THWMA TBE charts practically. Finally, it is pointed out that some of the issues may require further research. For instance, the performance of the THWMA TBE charts can be investigated in the case when the scale parameter of the gamma distribution is assumed to be unknown and thus estimated by the reference sample or historical data. Furthermore, the THWMA TBE methodology can be used in cases when TBE observations follow distributions such as Weibull or lognormal.

Use of AI tools declaration

The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

Acknowledgments

This research was supported by the Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2023R259), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

Conflict of interest

The authors reported that they do not have any conflicts of interest.

Availability of Data and Materials

Open source data sets are used and references are provided.

Supplementary

Appendix A: Mean, variance, and covariance of ${\mathit{T}}_{\mathit{t}}$

Suppose ${T}_{i}^{*}\stackrel{iid}{\sim }\mathrm{E}\mathrm{x}\mathrm{p}\mathrm{o}\mathrm{n}\mathrm{e}\mathrm{n}\mathrm{t}\mathrm{i}\mathrm{a}\mathrm{l}\left(\theta \right)$ for $i = \mathrm{1, 2}, \dots k$ and if $T = \sum _{i = 1}^{k}{T}_{i}^{*}$ then $T~\mathrm{G}\mathrm{a}\mathrm{m}\mathrm{m}\mathrm{a}(k, \theta)$. Let ${T}_{t}\stackrel{iid}{\sim }\mathrm{G}\mathrm{a}\mathrm{m}\mathrm{m}\mathrm{a}(k, \theta)$ at time $t = \mathrm{1, 2}, \dots$ then for the IC process, i.e., $\theta = {\theta }_{0}$,

For the IC process, the expected value of ${\stackrel{-}{T}}_{t}$ is given as follows:

Similarly, the variance of ${\stackrel{-}{T}}_{t}$ is given as follows:

Likewise, the covariance of ${T}_{t}$ and ${\stackrel{-}{T}}_{t-1}$, i.e., $Cov\left({T}_{t}, {\stackrel{-}{T}}_{t-1}\right)$ is given as follows:

Appendix B: Mean and variance of HWMA TBE chart

The charting statistic of the HWMA TBE chart is presented as follows:

For a stable process, the mean of ${U}_{t}$ is given as follows:

Similarly, the variance of ${U}_{t}$ is provided as follows:

or

Appendix C: Mean and variance of DHWMA TBE chart

The DHWMA TBE statistic, ${V}_{t}$ can be defined as follows:

For the IC process, the mean of ${V}_{t}$ is defined as follows:

In the same way, the variance of ${V}_{t}$ is defined as follows:

or

Appendix D: Mean and variance of THWMA TBE chart

The plotting statistic of the THWMA TBE chart is defined as follows:

The IC mean of ${W}_{t}$ is given as follows:

Similarly, the IC variance of ${W}_{t}$ is given as follows:

or