Integral equation solutions for the average run length for monitoring shifts in the mean of a generalized seasonal ARFIMAX(P, D, Q, r)s process running on a CUSUM control chart

The CUSUM control chart is suitable for detecting small to moderate parameter shifts for processes involving autocorrelated data. The average run length (ARL) can be used to assess the ability of a CUSUM control chart to detect changes in a long-memory seasonal autoregressive fractionally integrated moving average with exogenous variable (SARFIMAX) process with underlying exponential white noise. Herein, new ARLs via an analytical integral equation (IE) solution as an analytical IE and a numerical IE method to test a CUSUM control chart’s ability to detect a wide range of shifts in the mean of a SARFIMAX(P, D, Q, r)s process with underlying exponential white noise are presented. The analytical IE formulas were derived by using the Fredholm integral equation of the second type while the numerical IE method for the approximate ARL is based on quadrature rules. After applying Banach’s fixed-point theorem to guarantee its existence and uniqueness, the precision of the proposed analytical IE ARL was the same as the numerical IE method. The sensitivity and accuracy of the ARLs based on both methods were assessed on a CUSUM control chart running a SARFIMAX(P, D, Q, r)s process with underlying exponential white noise. The results of an extensive numerical study comprising the examination of a wide variety of out-of-control situations and computational schemes reveal that none of the methods outperformed the IE. Specifically, the computational scheme is easier and can be completed in one step. Hence, it is recommended for use in this situation. An illustrative example based on real data is also provided, the results of which were found to be in accordance with the research results.

and will appear in the published article if the submission is accepted.Please make sure it is accurate.

Unfunded studies
Enter: The author(s) received no specific funding for this work.If the data are held or will be held in a public repository, include URLs, accession numbers or DOIs.If this information will only be available after acceptance, indicate this by ticking the box below.For example: All XXX files are available from the XXX database (accession number(s) XXX, XXX.).

•
If the data are all contained within the manuscript and/or Supporting Information files, enter the following: All relevant data are within the manuscript and its Supporting Information files.
• If neither of these applies but you are able to provide details of access elsewhere, with or without limitations, please do so.For example: Data cannot be shared publicly because of [XXX].Data are available from the XXX Institutional Data Access / Ethics Committee (contact via XXX) for

Introduction
The discipline of statistical process control (SPC) provides tools such as control charts for monitoring processes and detecting changes in a given in-control model.The Shewhart, 2 cumulative sum (CUSUM), and exponentially weighted moving average (EWMA) control charts are the three most commonly used ones.The CUSUM and EWMA control charts can effectively detect small shifts in the parameters of interest (usually location parameters) whereas the Shewhart control chart such as the well-known X -chart is used to detect large shifts in process parameters.Details of the charts can be found in Montgomery [1].The CUSUM control chart was used in the present study due to its superior capability of detecting small parameter shifts and to take advantage of the benefits of using an upper-sided control chart.
Small levels of autocorrelation between successive observations can have a significant effect on the statistical properties of control charts.Numerous researchers have considered the effect of autocorrelation on control chart performance.Johnson and Bagshaw [2] and Bagshaw and Johnson [3] concluded that even small levels of autocorrelation generated by using autoregressive AR (1) or moving average MA (1) time series models made it more difficult to detect out-of-control signals on a conventional CUSUM control chart.Nevertheless, a CUSUM control chart has been used for monitoring shifts in the mean of a stationary ARMA (p, q) model process, and its performance has been evaluated when applied to processes with autocorrelated data [4].
Observations from real-world stochastic processes frequently exhibit a time series component.In particular, economic observations can have AR and MA components in a time series model.The error or so-called white noise, which is defined as the difference between the actual and approximated values, should be kept to a minimum to maintain the highest possible accuracy rate.Although the white noise (also referred to as Gaussian white noise) generated by autocorrelated observations usually follows a normal distribution, this is not always the case.Indeed, non-Gaussian white noise processes are useful for modeling a wide variety of phenomena, including wind speed, and the oxygen concentration and flow rate of water.Numerous studies have been conducted on time series models with non-Gaussian white noise, and when it is exponentially distributed (referred to as exponential white noise) is particularly interesting.Jacob and Lewis [5] considered an ARMA(1,1) process with exponential white noise, while Mohamed and Hocine [6] conducted a Bayesian analysis of an AR(1) process with exponential white noise.Furthermore, Pereira and Turkman [7] conducted a Bayesian analysis of threshold AR models with exponential white noise, while Suparman [8] recently proposed parameter estimations for an AR model of unknown order with exponential white noise.
Time series models are capable of capturing process correlation.An important class of time series model is the stationary process in which it is assumed that the process remains stable around a constant mean.For monitoring autocorrelated processes, the type of model can provide a framework for establishing statistical control.The following basic time series models were used in this study.The conventional Box-Jenkins AR integrated MA (ARIMA) model can be generalized as the AR fractionally integrated MA (ARFIMA) model that allows noninteger (fractional) values for the differencing parameter.This is advantageous when modeling time series with inherent long-range dependence (i.e., long-memory) [9,10 ].Accordingly, ARFIMA models are appropriate for time series long-memory processes, and exogenous variables can be incorporated into the time series model (i.e., ARFIMAX) to improve the latter's performance.Ebens [11 ] first introduced the ARFIMAX model with k exogenous variables and used it to estimate the realized volatilities in a Dow Jones Industrial Average portfolio, while Degiannakis [12 ] estimated and forecasted intraday realized volatility using ARFIMAX and ARFIMAX-TARCH (threshold AR conditional heteroskedasticity) models.
Seasonality (seasonal periods over a single year) is a phenomenon that affects time series.Numerous scholars have discussed ARFIMA models with seasonal components, and it has been used to explain a large number of real-world cyclical phenomena in long-memory processes.For instance, it occurs in series of revenues, inflation rates, monetary aggregates, and gross national product.In the current study, we concentrate on seasonal ARFIMAX (SARFIMAX) models.In terms of economic forecasting and other fields, cases where an exogenous variable is included in the forecasting model are usually more accurate than ones without it.In addition, time series model processes have been monitored by using control charts.For instance, Ramjee [13] discovered that the performances of Shewhart and EWMA control charts with correlated observations in an ARFIMA model were not good, and so proposed the hyperbolic weighted MA (HWMA) control chart instead.Subsequently, Ramjee et al. [14] designed an HWMA forecast-based control chart particularly designed for a non-stationary ARFIMA model with autocorrelated data.When residual control charts with ARFIMA and ARIMA models were used to monitor the air quality Taiwan [15], the one with the ARFIMA model was more appropriate than the one with the ARIMA model.Rabyk and Schmid [16] recently introduced the EWMA control chart for detecting changes in the mean of a long-memory process with the control chart's design based on an ARFIMA( , , ) p d q process.The most widespread measure of evaluating the performance of the control chart is the average run length (ARL).The power of a control chart is defined as the probability of its out-ofcontrol signal detection, while the average number of samples required to signal an out-of-control situation in the process (or expected value of the run length) is represented by the ARL.The incontrol and out-of-control ARLs are denoted as ARL0 and ARL1, respectively: ARL0 should be as large as possible whereas ARL1 should be as small as possible when the process changes undesirably.To improve the out-of-control detection performance on a CUSUM control chart, two alternative schemes are proposed that help to minimize ARL1 when ARL0 is fixed.
Assessing the ARL of a control chart using randomness and numerical methods is possible using a variety of techniques, including Monte Carlo, Markov Chain, and integral equations (IEs).The Monte Carlo simulation technique can be used to determine the expected run length until the target ARL value approaches that of an IE under the assumption of a sufficiently long simulation time.For the Markov Chain approach, the range between the upper and lower control limits is divided into sub-intervals, while the probabilities of the runs on the control chart are presented as a transitional probability matrix [17] over the sub-interval so that the ARL can be evaluated by inverting the matrix.In the meantime, the IE method comprises two alternatives: analytical IE and numerical IE.The analytical IE method generates exact ARLs or explicit formulas.The numerical IE method is transformed into a system of linear equations that can be solved for the ARL by using certain linear algebra-based methods.Crowder [18] provided a numerical method by using an IE of the second kind and a computer program for ARL evaluation of a normally distributed process on a two-sided EWMA control chart.Moreover, Hawkins [19] solved IEs for determining the ARL of a CUSUM control chart.Later, ARLs by using analytical IEs and the Markov Chain approach were used on EWMA and CUSUM control charts for observations from an AR(1) process with additional random error [20 ].To evaluate the ARL for a SARMA(1, 1)s process with exponential white noise on a CUSUM control chart, Phanyaem [21 ] developed explicit formulas for IEs initially based on the SARMA model.In addition, analytical IEs for the solution for the ARL for a long-memory SARFIMA model on CUSUM control chart were presented by Peerajit et al. [22 ].Recently, complex data from seasonal and non-seasonal MA processes with exogenous variables were used to evaluate the ARL on a CUSUM control chart by Sunthornwat and Areepong [23 ].Moreover, it is important to apply Banach's fixed-point theorem to prove the existence and uniqueness of an ARL.In this paper, our focus was on establishing ARLs for monitoring changes in a long-memory SARFIMAX process with underlying exponential white noise running on a CUSUM control chart.
The remainder of the article is organized as follows.Section 2 contains a brief review of the ARFIMAX models taking long memory and seasonality into account and descriptions of the CUSUM control chart, as well as the characteristics of the ARL and Banach's fixed-point theorem.Theoretical proof for the existence and uniqueness of analytical IE IEs, along with its comparable accuracy to the numerical ARL are covered in Section 3. The next section comprises derivations of the proposed analytical and numerical IE methods for the ARL for monitoring changes in a long-memory SARFIMAX process with underlying exponential white noise on a CUSUM control chart.The results of an evaluation of the performances of the IE methods are provided in Section 5.In Section 6, application of the methods to the monthly stock price data for PTT Public Company Ltd with the EUR/THB exchange rate as the exogenous variable is reported.Conclusions and a discussion of possible future work are provided in Section 7. Finally, the appendix contains technical details of the computation of the ARL by using the proposed methods.

Preliminaries
The relevant fundamentals of the long-memory SARFIMAX process with underlying exponential white noise are described in this section.In addition, we determine whether detecting changes in the mean of a generalized SARFIMAX process running on a CUSUM control chart is suitable for this study.For simplicity, the second subsection is focused on the upper-sided CUSUM control chart only.Next, performance evaluation of a control chart and the characteristics of the ARL are covered.In the last subsection, definitions related to Banach's fixed-point theorem are covered.

Long-memory SARFIMAX processes with underlying exponential white noise
The ARFIMA model, a generalized form of the usual Box-Jenkins ARIMA model but allowing non-integer (fractional) values for the differencing parameter, is useful for modeling time series with inherent long-range dependence [9,10].Furthermore, including an exogenous variable in the ARFIMA model produces the ARFIMAX model with better performance [12 ].In practice, time series are often found in meteorology, economics, hydrology, and astronomy.Seasonality is an autocorrelation structure where the data vary according to specific periods in a year, a phenomenon that is included in our analysis.
A SARFIMAX model of order ( , , , ) s P D Q r , with AR order ( ), P fractional integration order ( ), D MA order ( ), Q and exogenous variables order ( ), r and period () s for seasonal data is an extension of the ARFIMAX( , , , ) p d q r model applied to a seasonal time series.
where initial value is normally the process mean, seasonal AR, seasonal MA and exogenous variables coefficient are, respectively, 2 , ,..., 1,

Process monitoring on a CUSUM and EWMA Control Charts
One of the simplest and efficient SPC techniques for detecting small and moderate changes in the mean of a process is the CUSUM and EWMA control charts.The statistic of the upper-sided CUSUM control chart [24] can be expressed by the following recursive equation: where CUSUM control chart parameter t Y is a sequence of long-memory SARFIMAX( , , , ) s P D Q r process with underlying exponential white noise.0   and 0 0 C  as the reference and starting value, respectively.The starting value of the CUSUM statistic is  (i.e., 0 , ), and parameter  is used to balance the series () t C that is commonly chosen to be larger than but close to the value of the in-control mean.Remark: an alarm is triggered that indicates that the process might be out-of-control.In view of an exact performance evaluation, we chose rational or even real number design parameter . For the above CUSUM sequence ( ), t C the corresponding stopping time with predetermined threshold b can be deployed: where b  is the stopping time, b is the upper control limit (UCL) of CUSUM chart.
The statistic of the upper-sided EWMA control chart [2 5] can be expressed by the following recursive equation: .
where  is an exponential smoothing parameter with 0 1.

 
The starting value of the EWMA statistic is 00 ,  ZY the target mean value is ,  t Y is a processes with mean  and variance 2 . The upper control limit (UCL) and lower control limit (LCL) to detect the sequence is given by, where L is the width of the control limits, where 0  is the target mean,  is the process standard deviation.

Characteristics of the ARL
Throughout the whole paper, the error term () t  is a sequence of continuous i.i.d.
random variables taken from an exponential distribution with distribution function.  ( , ) , Fx which is used in the next section concerning monitoring of the mean.We consider the following simplified change-point models: where 0  and 1  are known parameters 10 ;.   As mentioned previously, when considering the change point in Equation ( 6), the ARL defined with (.) E  can be characterized as follows: Fx for a given change-point time ( ).
   indicates no change in the statistical process and is referred to as the in-control ARL 0 (ARL ).Herein, 1   marks the first time point that a change takes place from 0  to  in the statistical process, which is called out-of-control ARL 1 (ARL ).

Definitions and theories for evaluating the ARL
Definitions, theories, and concepts necessary to prove the existence and uniqueness of the IEs in the proof based on the functional analysis are covered in this section.

Definition 2. Let () ,
Md be a metric space and let :  MM be a mapping.A point uM  is called a fixed point of  if .() uu  Remark: Fixed point theorems guarantee the existence of a fixed point under appropriate conditions on the map  and the set .M Definition 3. A normed vector space M is a Banach space if the metric space () , Mdis complete, where   where the norm .
is given by the Euclidean distance.The () C M is space of continuous functions where M is a compact interval and the norm . is given by Md be a metric space and let :  MM be a map, and  is called contraction if there exists a fixed constant [0,1) u v M 

Theorem 1 (Banach's fixed-point theorem or the contraction theorem) [26]
Let () , Md be a complete metric space, then each contraction map :  MM has a unique fixed point.

The existence and uniqueness of the analytical and numerical IE methods for the ARL on a CUSUM control chart
The analytical IE method for the ARL of a CUSUM control chart running a long-memory SARFIMAX process corresponds to a Fredholm integral equation of the second type.In this section, Banach's fixed-point theorem is applied to prove the existence and uniqueness of the analytical IE solution and demonstrate that it has the same accuracy as the numerical IE method.
Throughout, this paper, let ()  be the ARL conditioned on the initial value for a CUSUM control chart running a long-memory SARFIMAX process with underlying exponential white noise [17].The initial value for monitoring CUSUM statistic 0 and c be the probability measure and induced expectation corresponding to .
 Moreover, ARL is defined as a function of ( ) , () C b satisfies the solution for the analytical IE as follows: where indicator function is The IE for the ARL derived from the Fredholm integral equation of the second kind in Equation ( 3) can be written as follows: where (.) f is the probability density function (pdf) of an exponential distribution [17].
Hence, deriving the analytical and numerical IEs method via Equation ( 8) becomes b The right-hand side of Equation ( 9) becomes a continuous function, and so theoretically, it is clear that the ARL function () . is a continuous solution on [0, ].
b for the IE in Equation ( 9).After that, definitions and Banach's fixed-point theorem are applied to guarantee the existence and uniqueness of the solution for functional ( ).  Theorem 2. The existence and uniqueness of the solution for the analytical IE are derived for a long-memory SARFIMAX( , , , ) s P D Q r process with underlying exponential white noise.

nn nn
In this situation we have   where ()  is a fixed point of . Hence, the IE in Equation ( 9) can be rewritten as

Proof. (Uniqueness).
Let  be a contraction mapping on a complete metric space.By considering Banach's fixed-point theorem, it can be shown that where metric is the space of all continuous functions ()  on interval ; [0, ] b endowed with supremum norm Thus, as proved and guaranteed by applying Banach's fixed-point theorem, the analytical IE for a long-memory SARFIMAX( , , , ) s P D Q r process with underlying exponential white noise on the CUSUM control chart exists and is unique for solutions of functional ( ). 

Analytical and approximate ARLs for a long-memory SARFIMAX process with underlying exponential white noise on a CUSUM control chart
In this section, we apply the proposed analytical and numerical IEs corresponding to a Fredholm integral equation of the second type representing analytical and approximate ARLs for monitoring changes in the process mean of a long-memory SARFIMAX process with underlying exponential white noise on a CUSUM control chart.

The proposed analytical IE for a SARFIMAX( , , , ) s P D Q r process with underlying white noise on a CUSUM control chart
Theorem 3 Exact formula ()  for monitoring changes in the process mean of a long-memory SARFIMAX process with underlying exponential white noise on a CUSUM control charts can be expressed as   Proof: Let ()  be the ARL of the IE corresponding to the long-memory SARFIMAX( , , , ) can be substituted as follows: By replacing 0   in Equation ( 10), we obtain Subsequently, by substituting (0) into Equation (10), it follows that can be derived as


After rearranging, we obtain c  Finally, by substituting constant c from Equation ( 13) into Equation ( 12), we obtain ...
The proof is complete.
As shown in the previous equation, the in-control ARL 0 (ARL ) is assigned exponential parameter  0 ( ), as follows: ...
    , and the shift in the mean is determined by  and 0 1   as follows: These circumstances show that the computational scheme can be completed in one stage.

The numerical ARL for a SARFIMAX( , , , ) s P D Q r process with underlying white noise on a CUSUM control chart
An approximate ARL using the numerical IE method via Fredholm's integral equation of the second kind was previously formulated for a CUSUM control chart by Peerajit et al. [22].For this method, an approximate ARL was developed by using the Gauss-Legendre quadrature rules technique.

Definition 5
The quadrature rules usually applied to integral  0 () b f g dg can be approximated by the sum of the areas of a rectangle is as follows:   for monitoring changes in the process mean of a long-memory SARFIMAX process with underlying exponential white noise () t Y on a CUSUM control chart can be expressed as   as the approximate ARL from the IE solution in Equation ( 9) using the Gauss-Legendre quadrature rule technique and replacing  with , With the system of m linear equations with m unknowns exists, then the unique solution is ( ), ( ),..., ( ) , is a column vector of () Therefore, the approximation for the integral is obtained in the summation where  j is replaced by  in ( ), j   which completes the proof.

Performance evaluation and comparison
A popular benchmark to measure the performance of a control chart is the ARL for detecting shifts in the process mean as a measure of its sensitivity and practicability.Thus, the performances of the analytical and numerical IE methods for determining the ARL for monitoring shifts in the process mean for a long-memory SARFIMAX( , , , ) s P D Q r process on a CUSUM control chart were compared.We determined the number of division points 800 m  for the approximated ARL by using the numerical IE method.The in-control situation comprises exponential white noise with mean 0 1,   while the out-of-control process is signaled by a shift in the process mean from 0  to 1  , where 10 (1 ).

 
means that the process is incontrol whereas 0,   means that the process is out-of-control, denoted as 0 ARL and 1 ARL , respectively.Accordingly, 0 ARL was fixed at 370 or 500.For the CUSUM statistic, the value of  was determined as 3.0, 3.5, or 4.0.Thus, the CUSUM control limit () b was calculated by using Equation (15) to achieve the specified ARL, the results for which are reported in

 
For a fixed 0 ARL , the method achieving the smallest 1 ARL has delivered the best performance [1].The results for the sensitivity, accuracy (see Definition 6), and computational times of the two methods are provided in Tables 2-6 and Fig 1.

Definition 6
The percentage accuracy ( % Accuracy ) indicates the relative performances of the analytical and numerical IE methods defined as follows: ( ) ( ) % Accuracy 100 100%, () where ()  and ()  are the ARL values obtained using the analytical and numerical IE methods, respectively.
The following algorithm was developed by using the Mathematica program to calculate the values of the ARLs for a long-memory SARFIMAX( , , , ) s P D Q r process on a CUSUM control chart: Step 1: Determine the values of coefficients .
 for the in-control exponential white noise   Step 3: Determine the value of known parameter    and the initial value of CUSUM statistic .
Step 4: Calculate the upper control limit () b by using Equation ( 15 The above algorithm was executed to obtain ARLs for out-of-control situations as measurements of the performance of the proposed analytical and numerical IE methods.Table 1.Calculated b for corresponding reference parameter  values for a long- memory SARFIMAX( ,, )  , The results for computed parameter b and reference parameter  for ARL 0 = 370 or 500 using the proposed ARL schemes are reported in Table 1.b is inversely proportional to  for long- memory process models [1] to [3].For instance, for prespecified ARL 0 = 370 on long-memory SARFIMAX(2, 0.15, 1, 1)12 with coefficient values 1 0.01,     15)).The performance assessment of the analytical and numerical IE methods for determining the ARL for detecting shifts in the process means requires computing their ARL1 values, as described earlier.Tables 2-5 provide the numerical results obtain from both methods in which we report their ARL1 values, percentage accuracy, and computational times.We summarize our findings as follows: 1.The results of the analytical and numerical IE methods for detecting changes in the process mean tended to decrease when the shift size was increased in the order of long-memory process models [ 4-5).
3. The reference value ()  was inversely proportional to the upper control limit and directly proportional to ARL1 for both methods.
4. The percentage accuracy results were ~99% in all cases, meaning that the proposed analytical IE method was very accurate.5.The values in parentheses in each column are the computational times for calculating ARL1: 31-34 minutes for ARL0 = 370 and 32-38 minutes for ARL0 = 500 using the numerical IE but less than 1 second for both using the analytical IE.  .In addition, the lowest reference parameter ( ) 3.0   was better for detecting small to moderate shifts in the process mean than 3.5 and 4.Moreover, when D was increased from 0.15 to 0.30, detecting changes was more sensitive.The results also show that the efficacy of both IEs for detecting changes was best for models [3] and [-3] with reference parameter ( ) 3.0   and ARL0 = 500.
In summary, the results show that the analytical IE method is a good alternative to the numerical IE method for evaluating shifts in the process mean of a long-memory SARFIMAX( , , , ) s P D Q r process with exponential white noise on a CUSUM control chart, especially because of the marked reduction in computational time.

Illustrative example
In this section, the results of applying the proposed analytical and numerical IE methods for determining the ARL for monitoring changes in the mean of a long-memory , SARFIMAX( , ) , processes with exponential white noise on a CUSUM control chart involving real data are provided.The real dataset comprises monthly stock price data for the PTT Public Company Ltd. in Thailand with the EUR/THB exchange rate as the exogenous variable (https://th.investing.com).In this example, we used monthly data from January 2003 to August 2018 comprising 188 observations.A seasonal pattern exists when a series is influenced by a seasonal factor, which is a 12-month difference term in the model (S = 12).Accordingly,  The results in Table 6 reveal that the probability of variable X was significant (0.05; p-value = 0.0355), which indicates that exchange rate does impact the stock price data for PTT.It can also be seen that the dataset is suitable for a long-memory SARFIMAX(1, 0.496151, 1 1)12 model presenting statistically significant parameters with coefficients  After that, the distribution of the white noise was confirmed as asymptotic exponential by using the Kolmogorov-Smirnov test.The results show that the white noise significantly fitted an exponential distribution (P-value = 0.091 > 0.05) with a mean of 3.42526 when the process was in-control, as reported in Table 7.

Table 7. Testing the suitability of an exponential distribution for the exponential white noise.
Based on the results in Tables 6 and 7 , the SARFIMAX( Next, we calculated ARL0 using the analytical and numerical IEs and chose the value of parameter ()  = 6.5 after calculating the control limits for ARL0 = 370 and 500 using Equation ( 16) 17.37735 and 19.2241, respectively).The rationale behind the use of 370 and 500 of the in-control ARLs is that they give good results in practice.The performances of the two methods for determining the ARL for detecting shifts in the process mean of a long-memory SARFIMAX(1, 0.496151, 1 1)12 on a CUSUM control chart are presented in Table 8 and Fig


It can be seen that the ARL1 results using the two methods are similar to those in Tables 2-5.The numerical results obtained from both methods were similar for both short-and longterm detections (ARL0 = 370 and 500, respectively) for all shift sizes of the process mean.Once again, the computation time for the analytical IE was far less than that of the numerical IE method (less than 1 second versus approximately 33 minutes, respectively).This means that the analytical IE for determining ARL1 is a good alternative for detecting changes in the process mean on a CUSUM control chart.
A SARFIMAX process involving real data was run on both CUSUM and EWMA control charts.The performances of the control charts were compared in terms of the ARL for detecting small to moderate shifts in the process mean; the ARL constructed using the analytical IE was used on the CUSUM control chart whereas the one using the numerical IE method was used on the EWMA control chart [27] for the EWMA control chart were used with prespecified ARL0 = 370 or ARL0 = 500, and smoothing parameter  for the EWMA control chart was determined as 0.1.The results of the comparison are summarized in Table 9.For ARL0 = 370, the EWMA control chart provided a smaller ARL1 than the CUSUM control chart for shift size 0 0.3

 
whereas the CUSUM control chart provided a smaller ARL1 than the EWMA control chart for shift size 0.5 2.0.

 
However, for ARL0= 500, the CUSUM control chart performed better than the EWMA control chart for shift sizes 0.3 2.0.

 
This means that the performance of the CUSUM control chart was more powerful than that of the EWMA chart when detecting moderate shifts in the process mean for ARL0 = 370 and 500.Therefore, the CUSUM control chart was more efficacious than the EWMA control chart under these conditions, which is in accordance with the study results.We can see that the signal was first given at the 12 th and 32 nd -47 th time points for ARL0 = 370 and 500, respectively (the observations in red plotted above the upper control limit).These results confirm the sensitivity of the proposed analytical IE for the ARL on a CUSUM control chart providing 176 and 156 out-of-control signals for ARL0 = 370 and 500, respectively.This situation directly affects the stock price of PPT.Hence, tracking changes in exchange rates can be key both in terms of economy and finance.In an investor's view, the control chart above is the point in making profitable decisions for investors.

Conclusions
Performance assessment of control chart can be measured by using the ARL determined via analytical and numerical IE methods used to analyze and approximate the ARL computation, respectively.We applied both methods to determine the ARL for monitoring changes in the mean of a long-memory SARFIMAX process with underlying exponential white noise on a CUSUM control chart.Computation of the ARL by using the numerical IE method was used to check the accuracy of the analytical IE method by using the percentage accuracy benchmark and comparing their computation times.The number of division points (m) for the numerical IE method was large enough so that its ARL converged to that determined using the explicit formulas.The results show that the analytical IE was as accurate as the numerical IE method but required less computational time because its computational scheme was easier and could be completed in one stage.
The CUSUM control chart is suitable for detecting small and moderate parameter shifts for autocorrelated processes with bounded support [28,29].However, finding optimal parameters of the CUSUM control chart in these situations can be somewhat difficult.Therefore, obtaining the optimal values for the parameters of a CUSUM control chart should be studied further.The scope of the study could be extended in terms of the ARL and the methods used to compute it and could be applied to other types of processes including application to real data with other distributions of white noise such as gamma and Weibull from the exponential family.is a Cauchy sequence of ( ).  Moreover, since  is a contraction mapping, we arrive at

Appendices
for After repeating this argument n times, we obtain n n m nm dd Applying the property of the sum as a geometric series on  gives ... ... ...
Hence, we obtain where The triangular inequality is used for the supremum norm as follows: [0, ] ,.

C b
Hence, by applying Theorem 4.2, the uniqueness of the ARL based on the analytical IE such that () T is confirmed.This completes the proof.Therefore, the ARL based on the analytical IE for the CUSUM control chart for a longmemory SARFIMAX( , , , ) s P D Q X process exists and is unique.The results are expressed as percentage accuracy with computational times in parentheses for the analytical IE method (seconds) and the numerical IE method (minutes).
[1] SARFIMAX(2, 0.  The results are expressed as percentage accuracy with computational times in parentheses for the analytical IE method (seconds) and the numerical IE method (minutes).
[  The results are expressed as percentage accuracy with computational times in parentheses for the analytical IE method (seconds) and the numerical IE method (minutes).
[1] SARFIMAX(2, 0.  The results are expressed as percentage accuracy with computational times in parentheses for the analytical IE method (seconds) and the numerical IE method (minutes).
[    Model [1] where ƞ = 3.0 Model [1] where ƞ = 3.5 Model [1] where ƞ = 4.0 Model [2] where ƞ = 3.0 Model [2] where ƞ = 3.5 Model [2] where ƞ = 4.0 Model [3] where ƞ = 3.0 Model [3] where ƞ = 3.5 Funded studies Enter a statement with the following details: Initials of the authors who received each award • Grant numbers awarded to each author • The full name of each funder • URL of each funder website • Did the sponsors or funders play any role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript?• NO -Include this sentence at the end of your statement: The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.• YES -Specify the role(s) played.• * typeset Competing Interests Use the instructions below to enter a competing interest statement for this submission.On behalf of all authors, disclose any competing interests that could be perceived to bias this work-acknowledging all financial support and any other relevant financial or nonfinancial competing interests.This statement is required for submission and will appear in the published article if the submission is accepted.Please make sure it is accurate and that any funding sources listed in your Funding Information later in the submission form are also declared in your Financial Disclosure statement.View published research articles from PLOS ONE for specific examples.King Mongkut's University of Technology North Bangkok has received funding support from the National Science, Research and Innovation Fund (NSRF) (Grant No. KMUTNB-BasicR-64-20). Powered by Editorial Manager® and ProduXion Manager® from Aries Systems Corporation NO authors have competing interests Enter: The authors have declared that no competing interests exist.Authors with competing interests Enter competing interest details beginning with this statement: I have read the journal's policy and the authors of this manuscript have the following competing interests: [insert competing interests here] * typeset Ethics Statement Enter an ethics statement for this submission.This statement is required if the study involved: /A" if the submission does not require an ethics statement.General guidance is provided below.Consult the submission guidelines for detailed instructions.Make sure that all information entered here is included in the Methods section of the manuscript.Applied Statistics Powered by Editorial Manager® and ProduXion Manager® from Aries Systems Corporation Format for specific study types Human Subject Research (involving human participants and/or tissue) Give the name of the institutional review board or ethics committee that approved the study • Include the approval number and/or a statement indicating approval of this research • Indicate the form of consent obtained (written/oral) or the reason that consent was not obtained (e.g. the data were analyzed anonymously) • Animal Research (involving vertebrate animals, embryos or tissues) Provide the name of the Institutional Animal Care and Use Committee (IACUC) or other relevant ethics board that reviewed the study protocol, and indicate whether they approved this research or granted a formal waiver of ethical approval • Include an approval number if one was obtained • If the study involved non-human primates, add additional details about animal welfare and steps taken to ameliorate suffering • If anesthesia, euthanasia, or any kind of animal sacrifice is part of the study, include briefly which substances and/or methods were applied • Field Research Include the following details if this study involves the collection of plant, animal, or other materials from a natural setting: Field permit number • Name of the institution or relevant body that granted permission • Data Availability Authors are required to make all data underlying the findings described fully available, without restriction, and from the time of publication.PLOS allows rare exceptions to address legal and ethical concerns.See the PLOS Data Policy and FAQ for detailed information.Yes -all data are fully available without restriction Powered by Editorial Manager® and ProduXion Manager® from Aries Systems Corporation A Data Availability Statement describing where the data can be found is required at submission.Your answers to this question constitute the Data Availability Statement and will be published in the article, if accepted.Important: Stating 'data available on request from the author' is not sufficient.If your data are only available upon request, select 'No' for the first question and explain your exceptional situation in the text box.Do the authors confirm that all data underlying the findings described in their manuscript are fully available without restriction?Describe where the data may be found in full sentences.If you are copying our sample text, replace any instances of XXX with the appropriate details.
b C It has a unique solution, as shown in Appendix B.

Theorem 4
of integral f is chosen by applying base / bm with heights at the midpoints of the intervals of length / bm beginning at zero and () Wg is a weight function.Division points The numerical IE method   () By applying ()
seen that for  = 3, b = 4.300562 and 4.021805, respectively.Only the 0 ARL , , b and  values are involved in the design of the CUSUM control chart (see Equation (
12 SARFIMAX( , , , ) P D Q r was selected as the model.Statistical software package Eviews 10 was utilized for 12 SARFIMAX( , , , ) P D Q r filtering and estimating the model parameters.
series behaving as white noise.The suitability of the model was checked by plotting the actual ( t Y ) and fitted ( ˆt Y ) components and appears to be in close agreement, as shown in Fig 2.

Appendix A :
To partially prove the existence of the solutions for functional ( ).  Proof.Let ([0, ]) C b be a set of continuous functions ()  on interval [0, ].
Cauchy sequence, there exists unique point  M such that That is to say, a set of continuous functions of the ARL defined on [0, ], b  kz is a kernel function of the analytical IE for the ARL based on IEs obtained by using Theorem 4.1:

Table 1 .
The coefficient parameters for long-memory

Table legends Table 1 .
Calculated b for corresponding reference parameter  values for a long-memory

Table 2 .
Comparison of the ARL1 values obtained using the analytical and numerical ARL

Table 3 .
Comparison of the ARL1 values obtained using the analytical and numerical ARL

Table 4 .
Comparison of the ARL1 values obtained using the analytical and numerical ARL

Table 5 .
Comparison of the ARL1 values obtained using the analytical and numerical ARL

Table 6 .
The 12 SARFIMAX( , , , ) P D Q r coefficients for the PTT stock price dataset from January 2003 to August 2018.

Table 7 .
Testing the suitability of an exponential distribution for the exponential white noise.

Table 8 .
Comparison of the ARL1 values obtained via the analytical and numerical IE methods for detecting changes in the mean of the long-memory SARFIMAX(1, 0.496151, 1 1)12 process on a CUSUM control chart for The results are expressed as percentage accuracy with the computational times in parentheses for the analytical IE (seconds) and numerical IE (minutes) methods.

Table 9 .
Comparison of the ARL1 values obtained via the analytical IE on a CUSUM control chart and a numerical IE method on a EWMA control chart with 0.1   for detecting changes in the mean of the long-memory SARFIMAX(1, 0.496151, 1 1)12 process.