{'fraction_non_alphanumeric': 0.04004042542952019, 'fraction_numerical': 0.019875836180759422, 'mean_word_length': 4.013268998793728, 'pattern_counts': {'":': 0, '<': 2, '<?xml version=': 0, '>': 3, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 8, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
One-sided Shewhart control charts for monitoring the ratio of two normal variables in Short Production Runs

Monitoring the ratio of two normal random variables plays an important role in several manufacturing environments. For short production runs, however, the control charts assumed infinite processes cannot function effectively to detect anomalies. In this paper, we tackle this problem by proposing two one-sided Shewhart-type charts to monitor the ratio of two normal random variables for an infinite horizon production. The statistical performance of the proposed charts is investigated using the truncated average run length as a performance measure in short production runs. In order to help the quality practitioner to implement these control charts, we have provided ready-to-use tables of the control limit parameters. An illustrative example from the food industry is given for illustration.


Introduction
Quality control (QC) enables the manufacturer to produce high-quality products according to the needs of the customer. Its tool kit helps to improve product quality by minimizing product cost, increase the efficiency of the process by reducing product waste. Control charts are famous tools for monitoring the assignable causes, they also tell us when corrective action must be taken or timely notify us when corrective action must be taken to improve the process behavior (see Shewhart 20 ). With growing competition in customer markets, manufacturers extremely depend on the quality of their products and services for their survival. The short-run manufacturing production process is become quite common for achieving the satisfaction of the customers like job-shops, which are categorized by a high amount of flexibility and manufacture diversity. Therefore, the life cycle of the products is decreasing rapidly. Nowadays, the production lines in several manufacturings and engineering processes have limited. In short runs production, some sources are fixed or the time span of the product is too short, maybe one hour or day. For example, any warehouse which may be on the lease. This warehouse is operating in the short run because it has a limited place. The owner cannot extend their business or shifted it to another place. When the agreement expires, he will have in position to enlarge business or shift to a large place. The same is the case, in manufacturing industries, like robotics manufacturing industry incorporates the limited production runs of automatic parts within the flexible production cells and semiconductor industry assembly of electronic boards and in beverages industries where the high volume of production and filling of soft drinks in every 24 to 48 hour required the 20 to 30 inspection between consecutive session.
There is a myth among some manufacturing organizations, as they mostly feel that control charts are less applicable for the short-run process as the duration of their production cycle is too short. Recent studies signify control charts for the production process draw the attention of quality practitioners. Quality professionals are more concerned about the quality of the product. They are doing continuous improvement in product quality by reducing variability. The repute of any industry depends on the quality of goods or services that they offer to the customer. Quality is a major strategy that increases the productivity of any industry.
There are various SPC (statistical process control) short run control charts presented in literature to serve the purpose. Short run control charts are more effective and useful for the small lot manufacturing runs with limited production data. Ladany 11 first introduced optimized-p chart for short production run. Later on, numerous authors introduced effective designs for such economic process, see for instance Ladany 9 , Naseri et al. 13 and Khoo et al. 10 .
This study is planned to design the Shewhart control chart for monitoring the ratio of two variables for a short production run. There are several production or manufacturing processes, where the quality characteristics of interest formed the ratio of two normal variables. Production strategies, where numerous components need to be blended together to get a product composition can require monitoring the ratios of random variables when quality experts are frequently interested in the relative evaluation of the same property for two-components. In fact, guaranteeing the stable ratio between different components permits the product specifications to be encountered. Automotive, aerospace, electronics, pharmaceutical, materials production, food preparation, and packaging industries are typical applications of these manufacturing environments.
Celano and Castagliola 4 designed the ratio type Shewhart control chart by using the data from the food processing company. They took the ratio of two seeds (pumpkin and flex) to meet the nutrition facts. Tran et al. 22 investigated the performance of the ratio chart using the run-rules scheme. Tran et al. 22 designed EWMA type control charts for efficient monitoring of the ratio. Tran et al. 23 assessed RZ chart efficiency and performance in presence of measurement error. Tran and Knoth 25 did an analysis to evaluate the steady-state ARL performance of the RZ chart. Tran et al. 24 designed the CUSUM control chart which was more efficient than its RZ setup. Nguyen et al. 16 and Nguyen, Tran and Heuchenne 17 introduced EWMA and CUSUM design for ratio type control chart at several sampling intervals. Nguyen, Tran and Goh 19 applied the variable sampling interval idea to RZ chart to evaluate its efficacy and performance. Nguyen and Tran 18 presented the one-sided RZ chart which was more efficient in case of measurement error present in the process.
The articles cited above are planned to observe the ratio over a production horizon considered as infinite. However, there are numerous situations, where the production run no longer is infinite. As far as we know, no research has been done regarding the observing of the ratio (of two normal variables) infinite horizon framework. The aim of this study is to fill this gap by introducing two one-sided RZ charts aiming to observe the decrease or increase in the finite runs. In the remainder of the study, they are denoted by "Sh + RZ Chart" for the upper sided chart and by "Sh − RZ Chart" for the lower sided one.
The rest of the article organized as follows, Section 2 denotes the main properties of ratio distribution. Section 3 deals with the design structure of one-sided RZ (Sh + RZ and Sh − RZ ) Shewhart charts for a limited production run. The truncated run length (TRL), as the performance measure of the proposed short run charts, is briefly described in Section 4. To assess the behavior of the new charts, a numerical analysis is done in Section 5. Real-life illustration of the proposed technique presented in Section 6. Finally, conclusion and recommendation for future research are given in Section 7.

Derivation of properties of the ratio distribution
Suppose that X and Y be the two normal random variables such that W = (X, Y ) T ∼ N (µ µ µ W , Σ Σ Σ W ), i.e. W is a bivariate normal random vector with mean vector and variance-covariance matrix as follows: where µ X and µ Y are the means of two variables and ρ is the correlation coefficient between them. Coefficients of variation (γ X , γ Y ) and standard-deviation ratio (ω) of X and Y are denoted by γ Celano and Castagliola 4 derived an adequate approximation for the c.d.f (cumulative distribution function) of Z as a function of γ X , γ Y , ω, and ρ as: where and Φ is the c.d.f of standard normal distribution (Sdn). After some tedious derivations, approximated p.d.f (probability density function) of the ratio Z is where φ(·) is the p.d.f. of Sdn. Solving the equation F Z (z |γ X , γ Y , ω, ρ) = p allows obtaining an approximate expression for the i.d.f (inverse distribution function) F −1 Z (p|γ X , γ Y , ω, ρ). We have where C 1 , C 2 , and C 3 are functions of ω, ρ, γ X , γ Y , p and they are: and where Φ −1 (·) is the i.d.f. of Sdn.
3 Design of two one sided Shewhart RZ charts for short production Runs The production run is planned to produce small size lot having N parts after a fixed rolling length H. Let I be the number of planned inspections of rolling horizon H and assume that no inspection takes place at the end of the run. By these settings, the sampling frequency (the time interval between two consecutive inspections) will be S h = H (I+1) hours. The observed values of random variables X and Y are used to calculate the ratio for two one sided Shewhart charts for short run. In order to monitor the ratio Z, samples of size n will be taken at every sampling interval from the process and the quality characteristic W will be measured for each item. Let [W i,1 , W i,2 , . . . , W i,n ] be the collected sample in which the couples To design the charts, let γ X and γ Y be the known and constant coefficients of variation and let z 0 = µ X,i µ Y,i and ρ 0 be the known in-control values of the ratio and the coefficient of correlation that will ensure the stability of process. We further assume that linear relationships are held between the sample standard deviation and the sample mean for both variables X and Y , i.e., σ X,i = γ X µ X,i and σ Y,i = γ Y µ Y,i for every i ≥ 1. This implies that the standard deviation of each sample can change proportionally to its mean such that their ratio remains constant. There are several quality characteristics in practice (such as weights, tensile strengths and linear dimensions) that can have a dispersion proportional to the population mean. Finally, the sample units are considered free to vary from sample to sample and thus it is possible to have µ µ µ W,i = µ µ µ W,k and Σ Σ Σ W,i = Σ Σ Σ W,k , for i = k. The monitoring statistic of the proposed chart is the ratio of sample means that should be calculated at the inspections i = 1, 2, , 3, ... as: To obtain the c.d.f and i.d.f of the statisticẐ i , one needs to some distributional proportions of the sample meansX i andȲ i . It is easy to see that By the definition and the mentioned assumptions, the standard deviation ratio ω i in each inspection can be calculated as: where ω 0 is the in-control standard deviations ratio. We are now in a position to obtain the c.d.f and i.d.f ofẐ i based on the c.d.f and i.d.f of Z in (2) and (4) as (Nguyen,Tran and Goh 19 ): The control limits of two separate one sided ratio charts are as follows The downward ratio uses to monitor the decrease in ratio of two variables. The lower and the upper control limits LCL − and U CL − of said chart are: where α 0 is the probability of type I error.
2. One sided Shewhart-RZ + chart (Sh + RZ Chart) The upward ratio uses to monitor the increase in ratio of two variables. The lower and the upper control limits LCL + and U CL + of said chart are: The ratio of two normal variablesẐ i will be plotted against the limits. We will count all points which fall outside of the limits of both charts. The process is consider to be out-of-control (occ) ifẐ i fulfills either of these conditions, Z i < LCL − orẐ i > U CL + and assignable causes will be removed. As the production run is too small, the traditional run length will not be applicable to determine the run length properties of the control charts. To serve the purpose truncated run-length TRL and its average (TARL) will be used to determine the performance of the charts. The details and computations of TRL and TARL are provided in what follows.

Properties of Truncated Run Length TRL
In certain manufacturing processes, it could not be practical to gather a sufficient amount of data at the beginning of the manufacturing process for trial limits in phase-I. For example, in the aerospace industry, the production rate of huge components can be finite. It takes too much time to gather enough samples for the applying control charts. On the other hand, it is often needed that the quality control procedure starts as early as possible since the rate of each product is high. In such circumstances, traditional quality charts are not as effective. It may need more samples to create control limits until they are precise adequate to screen the whole process. Traditional run-length does not give the reliable monitoring of assignable cause for such a process in which a finite (or limited) number I of samples is available in the production horizon having finite length H. Therefore, TRL as a modified version of traditional run length is used to determine the behavior of the process short production runs. It is define as "the number of samples until a shift is identified or until completion of process either happen first". TRL is a discreet random variable with support T RL ∈ {1, 2, ..., I, I + 1}. The event T RL = l for l ∈ {1, 2, ..., I} means that an ooc signal is declared by the chart before the termination of the production run and T RL = I + 1 expresses that the run is completed without detecting any shift in I inspections. The p.m.f (probability mass function) f T RL (l) and the c.d.f F T RL (l) of TRL can be derived as follows: and where p is the probability of an occ signal. Let TARL be the expected value of TRL that can be used in place of the traditional ARL to evaluate the performance of the control charts in finite horizon production. Moreover, let TARL 0 and TARL 1 be the in-control and out-of-control values of TARL, respectively. Substituting p by α (α is the probability of type I error) or 1 − β (β is the probability of type II error) in (14) and taking expectation over the TRL support, TARL 0 and TARL 1 will be derived as (Nenes and Tagaras 15 ): The error probabilities α and β can be calculated for the proposed charts by: and where z 1 = τ z 0 (τ > 0) and ρ 1 are the out-of-control values of the ratio and the coefficient of correlation, respectively. While the values of shift size τ ∈ (0, 1) correspond to decrease of nominal ratio z 0 , the values of τ > 1 correspond to an increase of nominal ratio z 0 .

Numerical Performance Analysis
In this section, we evaluate the statistical performance of the two one-sided Sh − RZ and Sh + RZ charts for the short production run. The performance of the chart is evaluated by the startup of the short run when a particular shift (τ ) happens. To serve the purpose, the TARL metric is used as the performance measure instead of using traditional ARL. Initially, let T ARL 0 = I when the process expected to be in control. Statistical performance has been calculated by considering the following parameters values.
The findings of proposed design in Tables 1-3 are as  follows: 1. Influence of sample size (n) to chart performance The width of control limits decreases as sample size (n) increases for particular values of γ X , γ Y , ρ 0 and I.  Table  1). Similar pattern can be observed for other values.     Table 3: 2. Effect of γ X and γ Y values to chart performance Mostly, the control limits LCL − and U CL + do not hold symmetry around 1. However, symmetry approximated attained by increasing sample size (n) for strong correlation and smaller values of (γ X , γ Y ). In general,  Table 1). Similar pattern may be observe for I = 30&50(γ X = γ Y ) and which is accordance to Celano and Castagliola 4 .  The following results can be drawn from Tables 4-9.
1. When shift size τ decrease or increase at fixed values of I, n and ρ 0 , the T ARL will decrease or increase accordingly. It can be observe that T ARL 1 < T ARL 0 in all cases. For instance, at fixed values of I = 10, Since, the ratio distribution is skewed but one sided Shewhart charts attain symmetry.

Real Life Illustration
The following example has been introduced in Celano et al. 5 for the implementation of a two-sided Shewhart control chart for monitoring the ratio of two    10.0 9.8 9.8 9.7 9.6 9.9 9.7 9.7 9.6 9.5 1.02 9.9 9.6 9.4 9.3 9.0 9.7 9.3 9. 9.7 9.5 9.5 9.4 9.3 9.9 9.7 9.6 9.6 9.4 1.00 9.9 9.9 9.9 9.9 9.9 9.9 9.9 9.9 9.9 9.9 1.01 9.9 9.7 9.6 9.6 9.4 9.7 9.5 9.5 9.4 9.3 1.02 9.8 9.4 9.3 9.1 8.9 9.5 9. 9.7 9.6 9.5 9.4 9.3 9.9 9.7 9.7 9.6 9.5 1.00 10.0 10.0 10.0 10.0 10.0 9.9 9.9 9.9 9.9 9.9 1.01 9.9 9.7 9.7 9.6 9.5 9.7 9.6 9.5 9.4 9.3 1.02 9.8 9.5 9.             normal variables in a long production run context and will be adapted, in this paper, to a short production run context. This example considers a real quality control problem from the food industry but with the simulation data. A muesli brand recipe is produced by using a mixture of several ingredients including sunflower oil, wildflower honey, seeds (pumpkin, flaxseeds, sesame, poppy), coconut milk powder and rolled oats. To meet the food's nutrient composition requirements declared in the brand packaging label and to preserve the mixture taste, the recipe calls for equal weights of "pumpkin seeds" and "flaxseeds". Furthermore, their nominal proportions to the total weight of box content are both fixed at p p = p f = 0.1. To satisfy the needs of customers, the brand boxes are manufactured in different sizes. According to Celano et al. 5 , whichever is the package dimension, the quality practitioner wants to perform on-line SPC monitoring at regular intervals i = 1, 2, . . . to check deviations from the in-control ratio z 0 = µp,i µ f,i = 0.1 0.1 = 1 due to problems occurring at the dosing machine. Here, µ p,i and µ f,i are the mean weights for "pumpkin seeds" and "flaxseeds", respectively. The quality practitioner collects a sample of n = 5 boxes every 60 minutes. Because the box size is allowed to change from one sample to another, it is possible to have µ p,i = µ p,k and µ f,i = µ f,k , ∀i = k. In the quality control laboratory, a mechanical procedure separates the "pumpkin seeds" and "flaxseeds" from the muesli mixture filling each box and the sample average weightsW p,i = 1 n n j=1 W p,i,j andW f,i = 1 n n j=1 W f,i,j are recorded. Finally, the ratioẐ i =W p,ī W f,i is computed. For this example, like in Celano et al. 5 , for i = 1, 2, . . . , and j = 1, 2, . . . , n both W p,i,j and W f,i,j can be well approximated as normal variables with constant coefficients of variation γ p = 0.02 and γ f = 0.01, i.e. W p,i,j ∼ N (µ p,i , 0.02 × µ p,i ) and W f,i,j ∼ N (µ f,i , 0.01 × µ f,i ). Moreover, the in-control correlation coefficient between these two variables is ρ 0 = 0.8. The process engineer has decided to implement the Sh + RZ chart in order to monitor the ratio for a short run production of H = 16 hours calling for I = 15 inspections, i.e. an inspection every hour.
For n = 5 and ρ 0 = 0.8, the control limit of the Sh + RZ control chart to potentially detect unexpected increase in the ratio for a short run production is U CL = 1.01421. Table -16 shows the set of simulated sample data collected from the process (based on Celano et al. 5 ), the corresponding box sizes 250-500 gr, and theẐ i statistic. The process is assumed to run in-control up to sample #10. Then, between samples #10 and #11, Celano et al. 5 have simulated the occurrence of an assignable cause shifting z 0 = 1 to z 1 = 1.01 × z 0 , i.e. a ratio percentage increase equal to 1%. Figure 2 shows the Sh + RZ control chart, which signals the occurrence of the out-of-control condition by plotting point #11 above the control limit U CL + = 1.01421 (see also bold values in Table  -16). The process is allowed to continue, while corrective actions are started by the repair crew who find and eliminate the assignable cause after sample #14 and restore the process back to the in-control condition.

Conclusion and Recommendations
In this article we examined the statistical performance of two one-sided Shewhart charts (Sh − RZ and Sh + RZ ) in monitoring the ratio of two normal variables using truncated average run length as a performance measure in short production runs. For practical consideration of quality experts, we have presented the probability limits and TARL values of both charts by varying the number of inspections. The sample size(n) and the number of inspections(M) affect the width of the control limits. Quality practitioners maybe adjust control limits according to the requirement of the production process. We observed the TARL performance of the both designs at fixed values n, ρ 0 , ρ 1 , γ X, γ Y , τ , Z 0 at prespecified T ARL 0 = 10, 30 & 50 . The main conclusions maybe drive from results. The performance of chart one-sided Shewhart charts (Sh − RZ and Sh + RZ ) is depend on the values n, γ X , γ Y and ρ 0 . When shift size increase or decrease the values of TRL increase or decrease accordingly. There are several extensions for future research like one-sided EWMA type and CUSUM type ratio chart for short production runs.