**Deviance and Pearson Residuals-Based Control Charts with Different Link Functions for Monitoring Logistic Regression Profiles: An Application to COVID-19 Data**

**Maryam Cheema 1, Muhammad Amin 1, Tahir Mahmood 2,3, Muhammad Faisal 4, Kamel Brahim <sup>5</sup> and Ahmed Elhassanein 5,\***


**Abstract:** In statistical process control, the control charts are an effective tool to monitor the process. When the process is examined based on an exponential family distributed response variable along with a single explanatory variable, the generalized linear model (GLM) provides better estimates and GLM-based charts are preferred. This study is designed to propose GLM-based control charts using different link functions (i.e., logit, probit, c-log-log, and cauchit) with the binary response variable. The Pearson residuals (PR)- and deviance residuals (DR)-based control charts for logistic regression are proposed under different link functions. For evaluation purposes, a simulation study is designed to evaluate the performance of the proposed control charts. The results are compared based on the average run length (ARL). Moreover, the proposed charts are implemented on a real application for COVID-19 death monitoring. The Monte Carlo simulation study and real applications show that the performance of the model-based control charts with the c-log-log link function gives a better performance as compared to model-based control charts with other link functions.

**Keywords:** ARL; control charts; COVID-19 data; deviance residuals; link functions; logistic profiling; Pearson residuals

**MSC:** 62P30; 62J02; 62J12; 62J20

#### **1. Introduction**

Quality is one of the main requirements of any organization's product and service goodwill. The control charts are the main tools of statistical process control (SPC) for monitoring the quality of products to improve the process capability [1]. In the last few decades, different control charts have been introduced and implemented in various industries to monitor the online process. These control charts include the Shewhart control chart [2], the cumulative sum (CUSUM) control chart [3] and the exponential weighted moving average (EWMA) control chart [4]. The Shewhart-type control chart is the memoryless control chart and is applied to detect the large shift in the mean and standard deviation of the process. However, EWMA and CUSUM are memory-type structures used to monitor small changes in the process parameters.

A single quality characteristic of a process is monitored by control charts, which may be quantitative or qualitative. Sometimes, the quality characteristic depends upon the explanatory variable(s). When the quality characteristic follows a normal distribution and

**Citation:** Cheema, M.; Amin, M.; Mahmood, T.; Faisal, M.; Brahim, K.; Elhassanein, A. Deviance and Pearson Residuals-Based Control Charts with Different Link Functions for Monitoring Logistic Regression Profiles: An Application to COVID-19 Data. *Mathematics* **2023**, *11*, 1113. https://doi.org/10.3390/ math11051113

Academic Editors: Wen Zhang, Xiaofeng Xu, Jun Wu and Kaijian He

Received: 16 January 2023 Revised: 11 February 2023 Accepted: 20 February 2023 Published: 23 February 2023

**Copyright:** © 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

is linearly associated with the explanatory variable, it is termed a simple linear profile. In linear profiling, one may want to monitor the intercept, slope, and error variance [5].

In profiling, if the normality assumption is violated, we move towards the generalized linear model (GLM)based profiling. The term GLM refers to a large class of models popularized by Nelder and Wedderburn [6]. In the literature, several monitoring studies were designed based on the GLM approach, and such charts are termed model-based control charts. Skinner et al. [7] proposed a model-based control chart using deviance residuals while the response variable follows a Poisson distribution, and they used the square root link function. Jearkpaporn et al. [8] proposed a model-based scheme based on the deviance residual for the system in which the response variable follows a Gamma distribution, and they used the log link function. In addition, Skinner et al. [9] studied the effectiveness of GLM-based control charts on the semiconductor process data based on the deviance residuals. Koosha and Amiri [10] considered the effect of autocorrelation presence between the observations in different levels of the independent variable in a logistic regression profile on the monitoring procedure (*T*<sup>2</sup> control chart) and proposed two remedies to account for the autocorrelation within logistic profiles. Shu et al. [11] reviewed the literature on regression control charts based on their importance and practical applications. Asgari et al. [12] proposed the GLM-based control chart to monitor a twostage procedure under Poisson distribution. Amiri et al. [13] investigated the profiles with binomial and Poisson responses in phase I and monitored them using three methods, namely Hoteling *T*<sup>2</sup> statistic, the F method, and the likelihood ratio test (LRT). Amiri et al. [14] concentrated on Phase II monitoring and proposed procedures for monitoring multivariate linear and GLM regression profiles. Qi et al. [15] developed a control chart to monitor GLM profiles using the weighted likelihood ratio tests. Moheghi et al. [16] studied the robust estimation and monitoring of parameters in GLM profiles in the presence of outliers. Recently, Kinat et al. [17] proposed the Pearson and deviance residuals-based control charts for the inverse gaussian response. Moreover, recent studies on GLM-based control charts and their applications can be found in the literature [18–23].

One of the most important members of the GLM family is the logistic regression profile. When the quality characteristics follow the Bernoulli distribution, we use the logistic regression profile. Different link functions can be used for logistic regression profiling. These link functions include the logit, probit, log-log, complementary log-log (c-log-log), and cauchit link function. Koosha et al. [10] studied the effect of applying different link functions on the performance of the *T*<sup>2</sup> control chart in monitoring the parameters of logistic regression profiles, and they found that the logit link function is best for logistic regression. Yu et al. [24] analyzed the performance of the LRProb chart under the assumption that only a small number of predictable abnormal patterns are available. Hakimi et al. [19] proposed some robust approaches to estimate the logistic regression profile parameters to decrease the effects of outliers on the performance of the *T*<sup>2</sup> control chart. Khosravi et al. [25] proposed three self-starting control charts to monitor a logistic regression profile that models the relationship between a binomial response variable and explanatory variables. Alevizakos [26] proposed two indices, cp and Spmk, for logistic regression profiles using different link functions: the logit, probit, and the c-log-log. The value of each index is approximately the same regardless of the used link functions. Jahani et al. [27] developed two control charts based on Wald and Rao score test (RST) to monitor nominal logistic regression profiles in Phase II.

The available literature showed that most researchers focused on deviance residualsbased control charts with probit and logit link functions. In this study, we evaluate the performance of various link functions in logistic regression profiling based on Pearson and deviance residuals, and find out which one of the link functions shows better performance. The outline of the research work is described as: Section 2 involves the methodology. Section 3 presents the structure of the proposed control charts based on Pearson and deviance residuals. Section 4 consists of numerical evaluations, which include the simulation

study of the proposed control charts. Section 5 describes a real-life application, and finally, Section 6 consists of the conclusion and future recommendations.

#### **2. The GLM-Based Control Charts**

The GLM-based control chart is used to enhance the ability of linear profile when the variable of interest follows an exponential family distribution. In this study, GLM-based control charts are designed based on deviance residuals (DR) and Pearson residuals (PR) of the logistic regression. Figure 1 shows all the steps of our proposed approach.

**Figure 1.** Follow chart of the process.

Suppose that our variable of interest (*y*) follows the Bernoulli distribution with probability mass function given by:

$$p(y\_i|\pi\_i) = \pi\_i^{y\_i} (1 - \pi\_i)^{1 - y\_i} \text{ ; } \mathbf{y} = \mathbf{0} \text{ } \mathbf{1} \tag{1}$$

where *π<sup>i</sup>* is the probability of success. The logistic regression model is a subset of the binomial regression model. When the response variable (*y*) of a regression model belongs to the Bernoulli distribution, the logistic regression model is the most commonly used statistical model; i.e., *yi* ∼ *Be*(*πi*), where the probability *π<sup>i</sup>* is a function of *Xi* and is defined by:

$$\pi\_i = \frac{\exp\left(\boldsymbol{x}\_i^\prime \boldsymbol{\beta}\right)}{1 + \exp\left(\boldsymbol{x}\_i^\prime \boldsymbol{\beta}\right)},\tag{2}$$

where *xi* is the ith row of *X*, which is an *n* × (*p* + 1) data matrix with *p* explanatory variables, and *β* is a (*p* + 1) × 1 vector of regression coefficients. The mean and variance of Bernoulli distribution are, respectively. given by *E*(*yi*) = *π<sup>i</sup>* and *var*(*yi*) = *πi*(1 − *πi*).

To convert Equation (1) into an exponential format, we rewrite it as:

$$f(y \Big| \pi\_i) = \exp\left[\log\left(\pi\_i^{y\_i} (1 - \pi\_i)^{1 - y\_i}\right)\right] \tag{3}$$

$$\mathbf{x} = \exp(y\_i \log \pi\_i + (1 - y\_i) \log(1 - \pi\_i))\tag{4}$$

$$=\exp(y\_i \log \pi\_i + \log(1 - \pi\_i) - y\_i \log(1 - \pi\_i))\tag{5}$$

$$=\exp(y\_i \log\left(\frac{\pi\_i}{1-\pi\_i}\right) - y\_i \log(1-\pi\_i))\tag{6}$$

The logit link function is suitable for linking *π<sup>i</sup>* and *Xi* in the logistic regression profiles denoted by *g*(*πi*), which is defined by:

$$\log(\pi\_i) = \log\left(\frac{\pi\_i}{1-\pi\_i}\right) = \beta\_0 + \beta\_1 \mathbf{x}\_{i1} + \beta\_2 \mathbf{x}\_{i2} + \dots + \beta\_p \mathbf{x}\_{ip\prime} \tag{7}$$

We also assume some other link functions such as the probit, c-log-log, and the cauchit to fit the logistic regression model. These link functions are given below:

$$\text{prob}\qquad \mathfrak{g}(\pi\_i) = \mathfrak{g}^{-1}(\pi\_i)$$

$$\begin{aligned} \text{c-log-log} & \quad \mathfrak{g}(\pi\_i) = \ln(-\ln(1-\pi\_i)) \\ \text{c-auchit} & \quad \mathfrak{g}(\pi\_i) = \tan(\pi(\pi\_i - 0.5)) \end{aligned}$$

The maximum likelihood estimator (MLE), based on the iterative reweighted least square technique (IRLS), is the most often used approach for estimating the parameter *β*, where the following log-likelihood of Equation (5) should be maximized.

$$l = \sum\_{i=1}^{n} \left[ y\_i \log(\pi\_i) + (1 - y\_i) \log(1 - \pi\_i) \right]$$

$$l = \sum\_{i=1}^{n} \left[ y\_i \log(\pi\_i) + \log(1 - \pi\_i) - y\_i \log(1 - \pi\_i) \right]$$

$$l = \sum\_{i=1}^{n} \left[ \log(1 - \pi\_i) \right] + \sum\_{i=1}^{n} y\_i \left[ \log \pi\_i - \log(1 - \pi\_i) \right]$$

$$l = \sum\_{i=1}^{n} \left[ \log(1 - \pi\_i) \right] + \sum\_{i=1}^{n} y\_i \left[ \log \left( \frac{\pi\_i}{1 - \pi\_i} \right) \right]$$

Substituting Equations (1) and (7) in the above expression, we have:

$$l = \sum\_{i=1}^{n} \left[ \log \left( 1 - \frac{\exp \left( \mathbf{x}\_i^\prime \boldsymbol{\beta} \right)}{1 + \exp \left( \mathbf{x}\_i^\prime \boldsymbol{\beta} \right)} \right) \right] + \sum\_{i=1}^{n} y\_i \mathbf{x}\_i^\prime \boldsymbol{\beta} \mathbf{x}\_i$$

After some simplifications, the result is as follows:

$$l = -\sum\_{i=1}^{n} \left[ \log(1 + \exp(\mathbf{x}\_i \boldsymbol{\beta})) \right] + \sum\_{i=1}^{n} y\_i \mathbf{x}\_i' \boldsymbol{\beta} \tag{8}$$

Differentiating Equation (3) with respect to *β* and equating to zero yield:

$$\frac{\partial l}{\partial \beta} = \sum\_{i=1}^{n} (y\_i - \pi\_i) x\_i' = 0 \tag{9}$$

Equation (4) is solved using the IRLS algorithm and obtained:

$$
\widehat{\boldsymbol{\beta}}\_{MLE} = \left(\boldsymbol{X}^{\prime}\widehat{\boldsymbol{W}}\boldsymbol{X}\right)^{-1}\boldsymbol{X}^{\prime}\widehat{\boldsymbol{W}}\boldsymbol{z},\tag{10}
$$

where

$$\hat{\mathcal{W}} = \text{diag}\left(\mathfrak{H}\_{\hat{\imath}}(1-\mathfrak{H}\_{\hat{\imath}})\right), \text{ and } \hat{z}\_{\hat{\imath}} = \log(\hat{\pi}\_{\hat{\imath}}) + \frac{y\_{\hat{\imath}}-\hat{\pi}\_{\hat{\imath}}}{\hat{\pi}\_{\hat{\imath}}(1-\hat{\pi}\_{\hat{\imath}})},\tag{11}$$

where *<sup>π</sup><sup>i</sup>* <sup>=</sup> exp(*x i β*ˆ *MLE*) 1+exp(*x i β*ˆ *MLE*) . The general form of the DR for the logistic regression model has the following form:

$$r\_i^D = \operatorname{sign}(y\_i - \hat{\pi}\_i) \sqrt{|2\{y\_i \log \hat{\pi}\_i + (1 - y\_i)\log(1 - \hat{\pi}\_i)\}|} / \sqrt{}$$

However, the PR for the logistic regression is defined as:

$$r\_i^P = \frac{y\_i - \hat{\pi}\_i}{\sqrt{\hat{\pi}\_i}(1 - \hat{\pi}\_i)}$$

#### **3. Structure of the Control Charts**

The Shewhart control chart includes a baseline as well as an upper control limit (UCL) and lower control limit (LCL), which are represented as dashed lines that are symmetric around the baseline. Measurements are plotted against a timeline on the chart. Measurements that exceed the limits are considered out of control (OOC). The control chart's central line (CL) is the acceptable value, which is an average of the historical check standard values. If our interest is to plot or monitor the ˆ *θ* statistic, then the UCL, CL, and LCL for ˆ *θ* are, respectively, given by:

$$\begin{aligned} \text{UCL} &= \text{E} \left( \hat{\theta} \right) + \text{kSD} \left( \hat{\theta} \right) \\ \text{CL} &= \text{E} \left( \hat{\theta} \right) \\ \text{LCL} &= \text{E} \left( \hat{\theta} \right) - \text{kSD} \left( \hat{\theta} \right) \end{aligned}$$

where E is the expectation of ˆ *θ*, *k* is the charting constant, and SD is the standard deviation of ˆ *θ*.

$$\text{1.3.1. Logistic Regression Model-based Control Chern Band on PR.}$$

In the model-based control chart, control limits of the PR are obtained as:

$$\begin{aligned} \text{UCL} &= \mathbb{E}\left(r\_i^P\right) + k\_1 \sqrt{\text{Var}\left(r\_i^P\right)} \\ &\text{CL} = \mathbb{E}\left(r\_i^P\right) \\ \text{LCL} &= \mathbb{E}\left(r\_i^P\right) - k\_1 \sqrt{\text{Var}\left(r\_i^P\right)} \end{aligned}$$

where *k*<sup>1</sup> is a constant that defines the size of control limits and is selected based on the fixed in-control (IC) average run length (ARL0). When any PR crosses the control limits, then the PR-logistic chart indicates the process is OOC otherwise, in IC condition.

#### *3.2. Logistic Regression Model-Based Control Chart Based on DR*

In the model-based control chart, control limits of DR are obtained by the following expressions:

$$\text{UCL} = \text{E}\left(r\_i^D\right) + k\_2 \sqrt{\text{Var}\left(r\_i^D\right)}$$

$$\text{CL} = \text{E}\left(r\_i^D\right)$$

$$\text{LCL} = \text{E}\left(r\_i^D\right) - k\_2 \sqrt{\text{Var}\left(r\_i^D\right)}$$

where *k*<sup>2</sup> is a constant that defines the size of control limits and is obtained against the fixed ARL0. When any DR crosses the control limits, the DR-logistic chart indicates that the process is OOC otherwise, in IC condition.

#### **4. Monte Carlo Simulation Study**

In this section, we evaluate the performance of the proposed control charts under different link functions with the help of a simulation study.

#### *4.1. Performance Measures*

The performance of the proposed control charts is evaluated using average run length (ARL) (cf. [17,28]). The run length (RL) is described as the number of points till a signal is indicated, and ARL is the average of the RL. The ARL is divided into two categories: IC (ARL0) and OOC (ARL1). When the control chart is in a stable state, we use ARL0 as the chart's performance measure, but when the process is in an unstable state, we use ARL1. A chart is considered to be the best among those under discussion if, on the fixed ARL0, it has the smallest ARL1. Hence, in this study, we focused on the ARL to evaluate the performance of the control charts under different link functions.

#### *4.2. Data Generation*

For the simulation purpose, firstly, we generate the response variable from the Bernoulli distribution with parameter *<sup>π</sup><sup>i</sup>* i.e., *yi* <sup>∼</sup> *Be*(*πi*), where *<sup>π</sup>* <sup>=</sup> *<sup>e</sup>*(*β*0+*β*1*x*) <sup>1</sup>+*e*(*β*0+*β*1*x*) , *<sup>x</sup>* <sup>∼</sup> *uni f orm*(0, 100). The true values of *β*<sup>0</sup> and *β*<sup>1</sup> are set to be 0.05 and 0.06, respectively. These variables are generated for a fixed sample size *n* = 1000.

#### *4.3. Algorithm for Charting Constants*

The following algorithm is used to find out charting constants.


If specified ARL0 does not achieve, then adjust the previous arbitrary value and repeat steps i–viii until specified ARL0 is obtained. By using the above-stated algorithm, the control charting constants are reported in Table 1.


**Table 1.** Constants for the logistic regression profiling with ARL0 = 200 under various link functions.

#### *4.4. Results and Discussion*

This section discusses the simulated results of the logistic model-based control charts with different link functions. Furthermore, we evaluate the performance of control charts by inserting the additive shifts in parameters of the model, such as *β*0, *β*1, and *μx*.

#### Monitoring *β*0,*β*<sup>1</sup> and *μ<sup>x</sup>* in Terms of ARL

In this section, we evaluate the simulation results of the control charts based on PR and DR in terms of ARL. For the out-of-control scenario, we considered the shifts in the *μ<sup>i</sup>* by changing the *β*<sup>0</sup> into *β*<sup>0</sup> + *η*, *β*<sup>1</sup> into *β*<sup>1</sup> + *θ*, and by changing the *μ<sup>x</sup>* into *μ<sup>x</sup>* + *δ*. The OOC results in terms of ARL are reported in Tables 2–4.

**Table 2.** Performance of *ARL*<sup>0</sup> = 200 with shift in *β*<sup>0</sup> under considered link functions.


**Table 3.** Performance of *ARL*<sup>0</sup> = 200 with shift in *β*<sup>1</sup> under considered link functions.


In Table 2, we evaluate the performance of the proposed control charts for monitoring *β*<sup>0</sup> in terms of *ARL*<sup>1</sup> under various link functions. The outcomes revealed that raising shifts minimizes the *ARL*<sup>1</sup> values for control charts based on DR and PR by using considered link functions. The estimated *ARL*<sup>1</sup> values at maximum shift (*η* = 0.28) for the logit, probit, c-log-log, and the cauchit link functions are 2.58, 2.42, 2.41, and 2.49, respectively. On comparing the performance of link functions in the PR-based control charts, we found that the c-log-log link function outperforms other link functions. The estimated *ARL*<sup>1</sup> of the DR-based control chart at the maximum shift (*η* = 0.28) for the logit, probit, c-log-log, and the cauchit link functions are 2.48, 2.44, 2.41, and 2.44, respectively. It is clear from the results that the c-log-log link function performs better than the other link functions. Finally, it is also observed that the DR-based control chart with the c-log-log link function outperformed the PR-based control chart with the c-log-log link function.


**Table 4.** Performance of *ARL*<sup>0</sup> = 200 with shift in *μ<sup>x</sup>* under considered link functions.

In Table 3, we evaluate the performance of the proposed control charts for monitoring *β*<sup>1</sup> in terms of *ARL*<sup>1</sup> under various link functions. The outcomes revealed that an increase in shift minimizes the *ARL*<sup>1</sup> values for control charts. The estimated *ARL*<sup>1</sup> values at maximum shift (*θ* = 0.28) for the logit, probit, c-log-log, and the cauchit link functions are 10.47, 10.13, 9.99, and 10.81, respectively. On comparing the performance of link functions in PR-based control charts, we found that the c-log-log link function outperforms other link functions. The estimated *ARL*<sup>1</sup> of the DR-based control chart at the maximum shift (*θ* = 0.28) for the logit, probit, c-log-log, and the cauchit link functions are 10.67, 8.03, 11.10, and 12.70, respectively, which reveals that the probit link function performs better than the other link functions. On comparing the performance of the DR- and the PR-based control charts in monitoring *β*<sup>1</sup> under different link functions, we observed that the DR-based control charts with probit link function perform better than the PR-based control charts with other link functions.

In Table 4, we evaluate the performance of the proposed control charts for the shifts in *μ<sup>x</sup>* in terms of *ARL* <sup>1</sup>*s* under various link functions. The outcomes revealed that raising shifts minimizes the *ARL* <sup>1</sup>*s* values for control charts based on DR and PR by using considered link functions. Table 4 presents the estimated *ARL* <sup>1</sup> values of the PR-based control charts under different link functions. The estimated *ARL* <sup>1</sup> values report at maximum shift (*δ* = 0.28) for the logit, probit, c-log-log, and the cauchit link functions are 2.95, 2.89, 2.95, and 2.95, respectively. On comparing the performance of link functions in PR-based control charts, we found that the probit link function outperforms other link functions. Further, the estimated *ARL* <sup>1</sup> of the DR-based control chart at the maximum shift (*δ* = 0.28) for the logit, probit, c-log-log, and the cauchit link functions are 2.84, 2.94, 2.94, and 2.93, respectively. It is clear from the results that the logit link function performs better than the other link functions. Finally, on comparing the performance of DR- and PR-based control charts in monitoring *μ<sup>x</sup>* under different link functions, we observe that the DR-based control charts with the logit link function perform better than PR-based control charts as compared to the other link functions.

#### **5. Application: COVID-19 Deaths Profile Monitoring**

In this modern era, COVID-19 affects different people in different ways. Here, we monitor the mortality status of COVID-19 patients who were admitted to Benazir Bhutto Shaheed Hospital, Rawalpindi, Pakistan. This data had been already collected by Akhtar et al. [29], where they collected this data from three hospitals, and we are considering one of these hospital data sets. In this application, our response variable (*Y*) is binary (*y* = 1, if the COVID patient is discharged deceased; otherwise, *y* = 0). Therefore, we consider this application for the evaluation of our proposed control charts. The data about the demographic and vital signs of all adult COVID-19 patients who were discharged from Benazir Bhutto Shaheed Hospital, Rawalpindi, during the first wave of COVID-19 (February to August 2020) were retrospectively collected. The National Early Warning Score (NEWS) was calculated by following the work of the Royal College of Physicians (RCP) (2012, page 14), and considered an independent variable in this study. The NEWS is based on a simple aggregate scoring system, which is computed based on vital signs (see Table 5 [30]).


**Table 5.** The strategy to estimate national early warning score (NEWS).

We had data of a total of 916 COVID-19 patients, consisting of mortality status (*y*) and NEWS (*x*). Based on the available data, we run logistic regression by setting logit, probit, c-log-log, and cauchit link functions. The Pearson (PR) and deviance (DR) residuals were obtained for estimated logistic regression models under each link function, and their means and standard deviations (SD) are reported in Table 6. Further, we set *ARL*<sup>0</sup> = 200 and obtained the limits of PR and DR-based control charts, which are also given in Table 6. To make the OOC dataset, we have estimated a new response variable (*Y*1) by using *μ<sup>Y</sup>* − *uni f orm*(0, 0.2). Similarly, we estimated logistic regression models and obtained shifted PR's and DR's. The PR- and DR-based control charts under different link functions are presented in Figures 2–5. The points under the pink and white shaded areas belong to the IC and OOC situation, respectively. The blue color points indicate the IC PR's and DR's, while the red color shows the OOC points. The PR- and DR-based control charts under the logit function are presented in Figure 2, while proposed control charts based on the probit function are portrayed in Figure 3. Figures 4 and 5 consist of PR and DR-based control charts under c-log-log and cauchit link functions, respectively. The out-of-control points are counted for each chart and reported in Table 6. It is revealed that the PR- and DR-based control charts under the cauchit link function have captured a large number of OOC signals. The second largest OOC points were detected by the PR- and DR-based control charts for the c-log-log link function. Hence, PR- and DR-based control charts under cauchit link and c-log-log link functions outperform all their counterparts.

**Table 6.** For the COVID-19 data, the mean and standard deviation of logistic regression residuals, control limits and OOC signals under different link functions.


**Figure 2.** Control chart based on PR and DR with the logit link function.

**Figure 3.** Control chart based on PR and DR with the probit link function.

**Figure 4.** Control chart based on PR and DR with the c-log-log link function.

**Figure 5.** Control chart based on PR and Dr with the cauchit link function.

#### **6. Conclusions**

In this study, the control charts are constructed based on the LR residuals (Pearson residuals (PR) and deviance residuals (DR)) with different link functions (the logit, probit, c-log-log, and cauchit). We evaluate the performance of the control charts with the help of a simulation study, where ARL is considered as the evaluation criterion. The simulation study shows that in monitoring *β*0*,* the c-log-log link function showed good performance for both PR and DR control charts as compared to the other link functions. While monitoring *β*1, the DR-based control chart with the probit link function performs better as compared to the PR-based control charts and with other link functions due to its minimum ARL. We also conclude that in monitoring *μx*, the logit link function gives a better performance for the DR-based control chart as compared to the PR-based control charts with other link functions. For the COVID-19 death data, the PR- and DR-based control charts were implemented and showed that charts based on the cauchit link function indicated a large number of signals. The 2nd link function which gives the 2nd largest OOC is the c-log-log link function. Therefore, it is concluded from the results that mostly the DR-based control charts with the c-log-log link function give a better performance as compared to the other link functions.

**Author Contributions:** Methodology, M.C., M.A., T.M., M.F., K.B. and A.E.; Formal analysis, M.C., M.A., T.M., M.F., K.B. and A.E.; Writing—original draft, M.C.; Writing—review & editing, M.A., T.M., M.F., K.B. and A.E. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the Deanship of Scientific Research at University of Bisha, grant number UB-GRP-42-1444.

**Data Availability Statement:** Not Applicable.

**Acknowledgments:** The authors extend their appreciation to the Deanship of Scientific Research at University of Bisha for funding this research through the general research project under grant number (UB-GRP-42-1444).

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **References**


**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
