A double integration method for generating exact tolerance limit factors for normal populations

. This article introduces a new method for generating the exact one-sided and two-sided tolerance limit factors for normal populations. This method does not need to handle the noncentral t -distribution at all, but only needs to do a double integration of a joint probability density function with respect to the two independent variables “ s ” (standard deviation) and “ x ” (sample mean). The factors generated by this method are investigated through Monte Carlo simulations and compared with the existing factors. As a result, it is identi ﬁ ed that the two-sided tolerance limit factors being currently used in practical applications are inaccurate. For the right understanding, some factors generated by this method are presented in Tables along with a guidance for correct use of them. The AQL (Acceptable Quality Level) is a good, common measure about quality of a product lot which was already produced or will be produced. Therefore, when performing sampling inspection on a given lot using a tolerance limit factor, there is a necessity to know the AQL assigned to the factor. This new double integration method even makes it possible to generate the AQLs corresponding to the one-sided and two-sided tolerance limit factors.


Introduction
The one-sided and two-sided tolerance limit factors are needed for quality control by sampling inspection.In response to this necessity, for the past tens of years, a lot of methodologies for generating the factors have been developed.Wald and Wolfowitz [1] developed an approximate formula for setting the tolerance limits in 1946.Even since then, other approximate formulas were developed.Representatively, Lieberman [2] developed a formula for one-sided tolerance limit factors and it was cited widely in other popular references (e.g., Natrella [3]).Unfortunately, Lieberman's formula tends to underestimate the factors; as the sample size n decreases, the factor gets much more underestimated.For this reason, it could be used only for sample sizes larger than 50.In 1970, Guttman [4] developed an approximate formula for one-sided tolerance limit factors using the statistical properties of the noncentral t-distribution.For reference, tens of years ago, it was not easy to use directly the noncentral t-distribution.Johnson and Welch [5] encouraged applications of the noncentral t-distribution through their research results.Owen and Amos [6] suggested some computational programs needed to use the noncentral t-distribution in order to set the tolerance limits.In 1972, Abramowitz and Stegun [7] proved that the noncentral t-distribution can be approximated by the normal distribution.With the aid of Abramowitz and Stegun's theory, Link [8] also suggested an approximate formula, and then Link compared the three formulas for one-sided tolerance limit factors (i.e., Lieberman [2], Guttman [4], and Link [8]).In comparatively recent years, Janiga and Garaj [9], and Young [10] presented their respective computational methodologies.In this way, various research activities relating to the one-sided tolerance limit factors have been done for a long time.It seems that most of recent studies are focused mainly on improvement or simplification of computational procedures by statistical software rather than exact calculation of the factors.
Although related studies are still ongoing, fortunately, the sufficiently accurate one-sided tolerance limit factors which were generated by Owen [11] in 1963 are being currently used in industry.Owen created each factor by an iterative method based on the noncentral t-distribution with degrees of freedom n and noncentrality parameter d.Later in 1993, Wheeler [12] reviewed and compared two methods for generating tolerance limit factors which were 2 Method and determination of factors 2.1 Definitions of k 1 , k 2 , k 1,practical and k 2,practical The four tolerance limit factors k 1 , k 2 , k 1,practical and k 2,practical are defined.(The practical tolerance limit factors k 1,practical and k 2,practical are determined by the factors k 1 and k 2 in a proper manner.)First, imagine that an inspector needs to perform a sampling inspection on a given normal population according to the inspection parameters of sample size n, confidence level (1 À a) and proportion p.Here, 0 < p < 1 and 0 < a < 1.
Let's assume that not "at least", but "exactly" the proportion p of the given normal population lies below the upper specification limit U (or above the lower specification limit L) and let's assume that the inspector tests only this population.The factor k 1 is defined as the one-sided tolerance limit factor which will be used by the inspector to ensure with confidence level (1 À a) that the proportion p of the normal population is below U (or above L).When the inspector performs the one-sided test using the factor k 1 , the probability that the population does not pass the test is "exactly (1 À a)".
Let's assume that not "at least", but "exactly" the proportion p of the given normal population lies between L and U, and the population mean is equal to (U + L)/2.In addition, let's assume that the inspector tests only this population.The factor k 2 is defined as the two-sided tolerance limit factor which will be used by the inspector to ensure with confidence level (1 À a) that the proportion p of the normal population lies between L and U. When the inspector performs the two-sided test using the factor k 2 , the probability that the population does not pass the test is "exactly (1 À a)".
The factor k 1,practical is defined as the one-sided tolerance limit factor which can be directly used in all practical applications to ensure with confidence level (1 À a) that at least the proportion p of any given normal population is below U (or above L).In case that more than the proportion p of the normal population is below U (or above L), although k 1 is used instead of k 1,practical as the onesided tolerance limit factor, the inspector can conservatively ensure with confidence level (1 À a) that at least the proportion p of the normal population is below U (or above L).Therefore, k 1 can be selected as substitute for k 1,practical in all practical applications.
The factor k 2,practical is defined as the two-sided tolerance limit factor which can be directly used in all practical applications to ensure with confidence level (1 À a) that at least the proportion p of any given normal population lies between L and U. Here, unlike the case of k 1 and k 1,practical , k 2 cannot always be selected as substitute for k 2,practical .Whether or not k 2 can be used instead of k 2,practical as the two-sided tolerance limit factor in a practical application depends on the situation.More specifically, if k 2 is larger than its corresponding factor k 1 , k 2 can be selected as substitute for k 2,practical .However, if k 2 is smaller than k 1 , k 2 cannot be selected as substitute for k 2,practical .Instead k 1 should be selected; the detailed reason is explained in Section 3.
Throughout this article, when referring to the term "two corresponding factors k 1 and k 2 ", it includes the meaning that the inspection parameters (i.e., sample size, confidence level, and proportion) given to the factor k 1 are identical to those given to the factor k 2 .Similarly, when referring to the term "AQL corresponding to factor k 1 (or k 2 )", it includes the meaning that the inspection parameters given to the AQL are identical to those given to the factor k 1 (or k 2 ).These remarks are also effective for k 1,practical and k 2,practical .For example, in this Section 2.1, the factors k 1 , k 2 , k 1,practical and k 2,practical have been defined based on the same inspection parameters (i.e., the same sample size n, the same confidence level (1 À a), and the same proportion p).Therefore, in this case, it can be said that the four factors correspond to one another.

Properties of s-distribution
A distribution named s-distribution is defined.Let's assume that a normal population of the random variable x is given and n samples (i.e., sample size is n) are randomly selected from it.If the same sampling process is repeated numerously, we can envision the probability distribution of the standard deviation s (s = { Here, such probability distribution is defined as s-distribution.To proceed further, some basic information is reminded.If x is the normal random variable and its variance is s 2 , then the sample mean x is normally distributed independently of s and the statistic (n À 1)s 2 / s 2 follows the chi-square distribution with degrees of freedom n (= n À 1).Based on such definition of sdistribution, we can proceed as follows.
Let s 0 = {(n À 1) 1/2 /s}s = as, a = (n À 1) 1/2 /s, and h = s 02 .In addition, let ' (s) be the probability density function of the continuous variable s, let c (s 0 )be the probability density function of s 0 , and let f (h) be the probability density function of h (= s 02 ).Then the following can be found.
This equality holds for all real numbers t (t > 0). Similarly,

Determination of k 1 by double integration
According to the definition of k 1 in Section 2.1, let's assume that not "at least", but "exactly" the proportion p of the given normal population lies below the upper specification limit U and let's assume that an inspector tests only this population.In addition, suppose that the inspector needs to ensure with confidence level (1 À a) by a sampling inspection that the proportion p of the normal population is below U, where 0 < p < 1 and 0 < a < 1.If the mean and the variance of the population are m and s 2 respectively and the sample size is n, then the one-sided tolerance limit factor k 1 to be used for the sampling inspection should satisfy the following equation.
The factor K p is the critical value corresponding to the proportion p in one-tailed test of normal population, and F x ð Þ is the probability density function of the normal variable x with the mean m and the variance s 2 /n.Using equations ( 1) and (2), equation ( 3) can be rewritten in the following form.
It should be reminded that f (s '2 ) is the probability density function of the variable s 02 with degrees of freedom n (= n À 1).Equation (4) does not need to handle the noncentral t-distribution.Only the probability density functions of normal distribution and chi-square distribution are needed to solve it.However, this equation can be solved numerically, not analytically.
In real situation, to determine numerically the factor k 1 , the probability density function in equation ( 4) should be integrated with respect to s '.In principle, its integration range should be from 0 to ∞.However, the range of 0 to 10 (or the range of 0 to 20 for large sample size, i.e., n > 40) is sufficient, because f(s '2 ) and/or the definite integral of F x ð Þ is nearly 0 for s 0 greater than 10 (or 20) and hence, does not affect the integration results.The definite integral with respect to s 0 is evaluated numerically by "division quadrature" with the aid of a statistical software (e.g., program "MINITAB 15").Each unit interval of s ' from 0 to 10 (or 20) is divided into 100,000 equal parts for the division quadrature.
64485 and a = (15−1) 1/2 /5 = 0.74833 have been substituted into equation ( 4).Both by "division quadrature" and by "trial and error", we find that k 1 is 2.566.(It will take a few minutes for an expert accustomed to this process to generate a factor needed.)In this way, all k 1 factors can be numerically determined.Table 1 shows some k 1 factors thus obtained.For reference, the one-sided tolerance limit factors being currently used in industry (Owen [11], Wheeler [12] and ISO [24]) coincide with the k 1 factors which are determined by the above method.

Determination of k 2 by double integration
According to the definition of k 2 in Section 2.1, let's assume that not "at least", but "exactly" the proportion p of the given normal population lies between the lower specification limit L and the upper specification limit U, and the population mean is equal to (U + L)/2.In addition, let's assume that an inspector tests only this population.Also suppose that the inspector needs to ensure with confidence level (1 À a) by a sampling inspection that the proportion p of the normal population lies between L and U, where 0 < p < 1 and 0 < a < 1.If the mean and the variance of the normal population are m and s 2 respectively, and the sample size is n, then the two-sided tolerance limit factor k 2 to be used for the sampling inspection should satisfy the following equation.
The factor K ' p is the critical value corresponding to the proportion p in two-tailed test of normal population.On the right-hand side of equation ( 5), the first term includes the definite integral with respect to the normal variable x greater than m, and the second term includes the definite integral with respect to x smaller than m.In the same manner as used in the case of one-sided, equation ( 6) is changed into the integral expression with respect to the variable s ' as follows Table 1.Numerically determined factor k 1 .

Sample size n
1 À a = 0.75 p = 0.75 Only the probability density functions of normal distribution and chi-square distribution are needed to solve equation (7).However, this equation cannot be solved analytically, but can be solved numerically like equation (4).
Similar to the case of k 1 , let's determine numerically a factor k 2 by equation (7).First, let m = 500 and s 2 = 5 2 , and let p = 0.95, (1 À a) = 0.95, and n = 15.Then equation ( 7) is rewritten in the following form to determine its corresponding k 2 .

See equation below.
Both by "division quadrature" and by "trial and error", we obtain k 2 = 2.616.At this time, K ' p = K ' 0.95 = 1.95996 has been used as the critical value corresponding to the proportion 0.95 in two-tailed test of normal population.In this way, all k 2 factors can be numerically determined.Table 2 shows some k 2 factors thus obtained.From Tables 1  and 2, when the sample size n is small, some k 2 factors are smaller than their respective corresponding k 1 factors.As the proportion p and the confidence level (1 À a) come closer to 1, this phenomenon becomes more apparent.For reference, the two-sided tolerance limit factors being currently used in industry (Wheeler [12] and ISO [24]) are different from the k 2 factors which are determined by the above method.

Simulation study
The factors k 1 and k 2 presented in Tables 1 and 2 can be investigated by Monte Carlo simulations.As an example, one two-sided tolerance limit factor is investigated to identify its correctness.In the case of n = 10, p = 0.90, and (1a) = 0.90, the two-sided tolerance limit factor k 2 numerically determined is 2.112 from Table 2. Before simulating with this factor, suppose that the mean m and the variance s 2 of the normal population to be two-sided tested are 500 and 5 2 respectively.After then, 491.77575 (= m À K 0 p s = 500-1.64485Â 5) and 508.22425 (= m + K ' p s = 500 + 1.64485 Â 5) are assigned as the lower specification limit L and the upper specification limit U respectively so that exactly the proportion 0.90 of the population lies between L and U.As the first stage of the simulation, the sample mean x and the standard deviation s of the ten independent random numbers generated from the given normal population (m = 500, s 2 = 5 2 ) by a statistical software (e.g., program "MINITAB 15") are calculated.Then it is observed whether or not the lower tolerance limit ðx À k 2 s) is larger than L, and whether or not the upper tolerance limit ðx + k 2 s) is smaller than U.This observation is repeated 1,000,000 times.Only in cases where (x À k 2 s) is larger than L and simultaneously ðx + k 2 s) is smaller than U, the two-sided test confirms that the population is acceptable.a Each of these two-sided tolerance limit factors is smaller than its corresponding one-sided tolerance limit factor.For example, if n = 5, (1 À a) = 0.99, and p = 0.99, then k 2 is 8.291 and k 1 is 8.939 (see Table 1 In the simulation, only 99,959 observations showed that ðx À k 2 s) was larger than L (= 491.77575) and simultaneously ðx + k 2 s) was smaller than U (= 508.22425).From this, the observed a value is 99,959/1,000,000 = 0.099959; in other words, the confidence level identified by the simulation is (1 À a) = 0.900041.Therefore, it can be said that the twosided tolerance limit factor k 2 = 2.112 determined numerically by the new integration method is correct.In this way, the correctness of all the one-sided and two-sided tolerance limit factors can be investigated by Monte Carlo simulations.Table 3 shows the results from the simulations that were conducted to investigate some k 1 and k 2 factors.
In Table 4, for n = 10, p = 0.90, and (1-a) = 0.90, the two-sided tolerance limit factor being used in industry is 2.546.This factor is significantly different from the above factor 2.112; the factor 2.546 is the one that was unduly overestimated.Table 4 additionally shows that regardless of sample size, confidence level, and proportion, the k 2,existing factors are always larger than the k 2 factors.Therefore, from Tables 3 and 4, it can be said that the twosided tolerance limit factors being used in industry are inaccurate.However, we can observe that as the sample size n increases, the difference between k 2 and k 2,existing gets smaller.

Practical factors and guidance for correct use
According to the definition of k 1,practical in Section 2.1, even in case that more than the proportion p of any given normal population is below U (or above L), the use of k 1,practical should enable inspectors to ensure with confidence level (1 À a) that at least the proportion p of the normal population is below U (or above L).Fortunately, although its corresponding factor k 1 is used as the one-sided tolerance limit factor instead of k 1,practical , the inspectors can conservatively ensure the requirement with confidence level (1 À a).Therefore, k 1 can be selected as substitute for k 1,practical in all practical applications.Thus, we can say According to the definition of k 2,practical , even in case that more than the proportion p of any given normal population lies between L and U, and moreover, the population mean is not equal to (U + L)/2, the use of k 2,practical should enable inspectors to ensure with confidence level (1 À a) that at least the proportion p of the normal population lies between L and U. Here, unlike the case of k 1 and k 1,practical , k 2 cannot always be selected as substitute for k 2,practical .Whether or not k 2 can be used instead of k 2,practical depends on the situation.Imagine the two corresponding factors k 1 and k 2 which were determined according to the same inspection parameters (i.e., the same sample size n, the same confidence level (1 À a), and the same proportion p).Here, if k 2 is larger than k 1 , k 2 can be selected as substitute for k 2,practical .On the contrary, if k 2 is smaller than k 1 (see the k 2 factors marked with superscript "a" in Table 2), k 2 cannot be selected as substitute for k 2,practical .Instead k 1 should be selected.To summarize, In case that the population mean m is not equal to (U + L)/2, i.e., the population mean does not coincide with the center of the specification range, the factor k 2,practical will be affected by the deviation magnitude of m from the center.If the population mean comes closer to U (or L) and the variance gets smaller, the sampling inspection will become more similar to the one-sided test.Moreover, the inspector performs the test, not knowing the mean and variance of the given normal population.Therefore, if k 2 is not larger than its corresponding factor k 1 , k 1 (not k 2 ) should be selected as substitute for k 2,practical for the twosided test.This is the reason that the determination of k 2,practical depends on whether k 2 is larger than k 1 or not.For reference, the two-sided tolerance limit factors being currently used in industry (Wheeler [12] and ISO [24]) are always larger than their respective corresponding onesided tolerance limit factors.
Although this assumption is not realistic, let's assume that the inspector knows the mean of a given normal population to be tested.In this case, k 2,practical needed for the two-sided test can be directly and numerically determined, not relying on k 1 and k 2 .As an example, suppose that the inspector needs to perform a two-sided test on the normal population according to the inspection parameters of n = 5, (1 À a) = 0.95, p = 0.95, L = 490, and U = 510.In addition, let's assume that the population mean m is known to be 503.4206.Then, k 2,practical (here, denoted by k' 2 for convenience during calculation) can be determined by the following equation.(Before calculating, we assign 4 2 as the variance s 2 of the population so that nearly the proportion p = 0.95 of the population lies between L and U.) (503.4206= 510-1.64485Â 4) Notably, in the last term on the right-hand side of equation (10), U À k 0 2 s 0 /a is smaller than m.If we substitute related parameters, i.e., n = 5, a = 0.05, L = 490, U = 510, m = 503.4206,s 2 = 4 2 and a = (n À 1) 1/2 /s = 0.5, into equation ( 10) and then solve it using both "division quadrature" and "trial and error" in a similar manner to determination of k 1 and k 2 , we obtain k 0 2 (= k 2,practical ) = 4.156, which is needed for the two-sided test.However, if n is 15, a different result is obtained.At this time, if we substitute n =15, a = 0.05, L = 490, U = 510, m = 503.4206,s 2 = 4 2 and a = (n À 1) 1/2 /s = 0.93541 into equation ( 10) and then solve it, we obtain k 0 2 (= k 2,practical ) ≒ 2.566.The above two calculation results are summarized as follows.when n = 5; Thus, irrespective of sample size n, the factor k 2,practical calculated directly is always placed somewhere between the two corresponding factors k 1 and k 2 .However, as previously mentioned, inspectors usually perform two-sided tests (including one-sided tests), not knowing the mean and variance of normal populations, and hence cannot calculate the k 2,practical factors.Therefore, k 2,practical should be selected from the two corresponding factors k 1 and k 2 by equation ( 9).Tables 5 and 6 show some k 1,practical and k 2,practical factors that should be used in practical applications.

AQLs corresponding to factors
According to Owen [11], AQL is defined as "percentage defective in a lot corresponding to a 0.95 chance of accepting the lot".The AQL is a good, common measure about quality of a product lot as a whole, which was already produced or will be produced.Therefore, when testing a given lot using a tolerance limit factor, there is a necessity to know the AQL which is assigned to the sampling inspection.In this regard, the method for generating the AQLs corresponding to the tolerance limit factors is introduced.
Suppose that the mean m and the variance s 2 of the normal population to be two-sided tested are 500 and 5 2 respectively and that the inspection parameters are n = 5, (1 À a) = 0.90, and p = 0.90.Also assume that the population Table 5. Practical factor k 1,practical (= k 1 ).
When rewriting, a = (n À 1) 1/2 /s = (5-1) 1/2 /5 = 0.4, m = 500, and k 2 = 2.597 were substituted into equation (7).Now, both by "division quadrature" and by "trial and error", we obtain K ' p = 1.38253.This two-tailed critical value is matched to the proportion p = 0.8332.Therefore, the AQL corresponding to k 2 = 2.597 is (1 À p) Â 100 = 16.68%.In this way, the AQLs corresponding to the one-sided and two-sided tolerance limit factors can be obtained.Table 7 shows some AQLs thus obtained.The unit "%" is omitted in practical uses of AQLs, however, it is not omitted for clear discernment of the AQLs from the tolerance limit factors in this article.
In Table 7, although two corresponding factors k 1 and k 2 are distinctly different from each other, the difference between the two AQLs matched to the two factors k 1 and k 2 is very small.For example, the AQL in the one-sided test of n = 5, (1 À a) = p = 0.90, and k 1 = 2.742 is 16.66% and it is nearly equal to the AQL 16.68% in the two-sided test of n = 5, (1 À a) = p = 0.90, and k 2 = 2.597.As another particular observation, in the case of (1 À a) = p = 0.99, as the sample size n increases, the AQL also increases.In contrast, in the case of (1 À a) = p = 0.90, the opposite phenomenon is seen.
The AQLs shown in Table 7 are based on the factors k 1 and k 2 .However, the practical AQLs should be used in practical applications in the same manner as k 1,practical and k 2,practical are used in practical applications.Table 8 presents some practical AQLs corresponding to k 1,practical and k 2,practical .It is common that people estimate the quality level of a lot from the AQL assigned to the sampling inspection.When establishing a sampling inspection plan, it is recommended that the fixed confidence level "(1 À a) =0.95" be adopted, because, in this case, simply the value "(1À p)" itself becomes the AQL corresponding to the planned sampling inspection.

Results and discussion
The factors generated by the new double integration method were investigated through Monte Carlo simulations and it was identified that they are correct and exact.During calculation and simulation, the one-sided tolerance limit factors were compared with their respective corresponding two-sided tolerance limit factors generated by the new method.As a result, a new fact was found.When the sample size n is small, some k 2 factors are smaller than their respective corresponding k 1 factors.As the proportion p and the confidence level (1 À a) come closer to 1, this phenomenon becomes more apparent.Particularly, in Table 4, through direct comparison of the two-sided tolerance limit factors determined by the new method with a The unit "%" was not omitted for clear discernment of the AQLs from the tolerance limit factors.b These two-sided tolerance limit factors and AQLs are the ones derived based on the assumption that the population mean coincides with the center of the specification range.
P. Kang: Int.J. Metrol.Qual.Eng. 13, 16 (2022) the ones generated by the existing methods, we can see that the existing two-sided tolerance limit factors being used in industry are inaccurate.As an example, for n = 10 and (1a) = p = 0.90, the two-sided tolerance limit factor k 2 determined by the new method is 2.112.However, according to the existing methods, k 2 is 2.546.Here, it was already demonstrated by the simulation results in Table 3 that the factor k 2 (= 2.112) is correct and exact.Therefore, we can say that the factor k 2 (= 2.546) being currently used is inaccurate.
If the population mean m does not coincide with the center of the specification range, the practical two-sided tolerance limit factor k 2,practical is affected by the deviation magnitude of m from the center.Even in this case, if the population mean is known, the factor k 2,practical can be calculated by the new method.In Section 3, it was illustrated that the factor k 2,practical thus calculated is always placed somewhere between its two corresponding factors k 1 and k 2 .As an example, for the inspection parameters of n = 5, (1 À a) = p = 0.95, L = 490, and U = 510, if the population mean is known to be 503.4206, in this case, we can calculate the practical factor k 2,practical (= 4.156) by the new method, and hence, we can identify that the practical factor lies between its two corresponding factors k 1 (= 4.203) and k 2 (= 3.917).However, inspectors have to conduct the sampling inspection, not knowing the mean and variance of the normal population to be tested, and hence cannot calculate k 2,practical directly.Therefore, k 2,practical should be selected from the two corresponding factors k 1 and k 2 in a proper manner.If k 2 is smaller than k 1 , k 1 (not k 2 ) should be used as substitute for k 2,practical .Unlike k 2 , k 1 can always be used as substitute for k 1,practical .
The AQL is regarded as a common, familiar index representing the expected quality level of a product lot from the aspect of statistical quality control by sampling inspection.Therefore, it is necessary to match the tolerance limit factors to their respective corresponding AQLs.Fortunately, the AQL can also be derived by double integration of a joint probability density function of normal distribution and chi-square distribution in the same manner as used to generate the tolerance limit factor itself.In the AQL Tables produced by the new method, although two corresponding factors k 1 and k 2 are distinctly different from each other, the AQL matched to k 1 is nearly (but not exactly) equal to that matched to k 2 .As an example, in Table 7, for n = 5 and (1 À a) = p = 0.90, the factor k 1 and its corresponding AQL determined by the new integration method are 2.742 and 16.66% respectively.In addition, the factor k 2 and its corresponding AQL determined in the same way are 2.597 and 16.68% respectively.Therefore, we can see that although the two factors k 1 and k 2 are distinctly different from each other, the two AQLs are nearly equal.As another interesting observation, in the case of p = (1 À a) = 0.99, as the sample size n increases, the AQL also increases.In contrast, in the case of p = (1 À a) = 0.90, the opposite phenomenon is seen.Especially, in case that (1 À a) is 0.95, the value "(1-p)" itself can be regarded as the AQL corresponding to the tolerance limit factor.Therefore, it is recommended that when establishing a sampling inspection plan, the fixed confidence level "(1 À a) = 0.95" be reflected in the plan.

Conclusion
The tolerance limit factors are important to conducting the statistical quality control by sampling inspection.The methodology for generating the one-sided tolerance limit factors was established in early 1960s.The factors  [11] and Wheeler [12] are correct and exact, and are being fully utilized in industry.However, in the case of two-sided, any correct methodology has not been established yet, and hence the inaccurate two-sided tolerance limit factors are being currently used.In addition, it seems that conceptual understanding needed for applications of the two-sided tolerance limit factors is still insufficient.
To resolve such problems, this article introduced a new double integration method for generating correctly the onesided and two-sided tolerance limit factors for normal populations along with a guidance for correct use of them.The new method generates the factors, not relying on the noncentral t-distribution, by double integration of a joint probability density function with respect to "s" (standard deviation) and "x" (sample mean).Even AQLs corresponding to the factors can be calculated in the same way.In this article, the factors thus generated were investigated through Monte Carlo simulations and it was identified that they are correct and exact.
In the nuclear fuel fabrication plants, great numbers of fuel pellets are manufactured in the pellet production process every day.Each cylindrical pellet is approximately 5 grams in weight and approximately 1 cm in length.After the production of each lot, from 5 to 50 pellets as per quality inspection item are randomly sampled from the lot, and then pellet's diameter, length, density, etc. are precisely measured for quality control and inspection.However, in order to decide "conformance" or "non-conformance" for each of the pellet's quality inspection items (i.e., diameter, length, density, etc.), the inaccurate two-sided tolerance limit factors have been used till now.Now we know that those factors being currently used in the fabrication plants are the ones unduly overestimated by the existing methods.
In addition, it is not easy for the inspectors to use other special one-sided and two-sided tolerance limit factors than the ones that are frequently used in practical applications and can be obtained directly from the statistical handbooks.As an example, let's assume that the inspectors need the one-sided or two-sided tolerance limit factor to be used specifically for "n = 53 and p = (1 À a) = 0.88".In this special case, the inspectors cannot easily obtain the needed factor.The reason is that any statistical handbook cannot contain all the factors, therefore, the books contain only the factors to be frequently used.Moreover, in most cases, the inspectors cannot calculate directly the needed factor, because the existing methods are difficult and complicated to handle.Additionally, even in cases where the needed factors are found in the handbooks, the inspectors should be satisfied with the simple use of the factors only for inspection, not knowing the AQLs corresponding to those factors; because there have been no means of matching the factors to the AQLs till now.
However, from now on, if necessary, the inspectors or engineers in industry can directly generate the exact onesided and two-sided tolerance limit factors, including the AQLs corresponding to the factors, by using the new double integration method introduced in this article.If anyone becomes accustomed to the use of this new method, it will presumably take a few minutes to generate a factor needed.Finally, we hope that this new method is widely disseminated and hence numerous inspectors performing the sampling inspection practically use the tolerance limit factors and AQLs generated by the new method.
a Each of the (1 À a) values was estimated from 1,000,000 observations.

Table 4 .
Comparison of k 2 and k 2,existing.
[24]determined by new integration method and k 2,existing being currently used in industry a These factors are the ones that were provided by Wheeler[12].bThesefactorsare the ones that were provided by ISO[24].cN/Ameans"not applicable".Either Wheeler or ISO did not provide the two-sided tolerance limit factors for some sampling inspection parameters.6P.Kang: Int.J. Metrol.Qual.Eng.13,16 (2022)

Table 7 .
AQLs corresponding to k 1 and k 2