# A Procedure for the Selection of Single Sampling Plans by Variables Based on Pareto Distribution

Authors: Geetha Sathya Narayanan; Vijayaraghavan Rajarathinam
Journal: Journal of Quality and Reliability Engineering
 (2013)

---

## Abstract

Product control in statistical quality control is the methodology that deals with procedures for taking decisions about one or more lots of finished products manufactured by production processes. Sampling inspection by variables is one of the major classifications of product control and comprises procedures for deciding about the disposition of a lot of individual units based on sample measurements of units on a quality characteristic under study. These procedures are defined under the assumption that the quality characteristic is measurable on a continuous scale, and the functional form of the probability distribution must be known. Inspection procedures which have been developed under the implicit assumption that the quality characteristic is distributed as normal with the related properties are found in the literature of product control. The assumption of normality may not be realized often in practice, and it becomes indispensable to investigate the properties of sampling plans based on nonnormal distributions. In this paper, a single sampling plan by variables is formulated and evaluated when the quality characteristic is assumed to be distributed according to a Pareto distribution. A procedure is developed for determining the parameters of the plan for specified requirements to ensure protection to both the producer and consumer.


---

## Body

## 1. Introduction

Sampling inspection is an activity for taking decisions on one or more lots of finished products which have been submitted for inspection. The decision of either acceptance or rejection of the lots is usually taken by adopting suitable sampling inspection procedures called sampling plans. Sampling plans are generally categorized into two types, namely, lot-by-lot sampling by attributes and lot-by-lot sampling by variables. In lot-by-lot inspection by attributes, one or more samples of items are drawn from a given lot of manufactured items; each item in the sample(s) is classified as conforming or nonconforming, and the decision of acceptance or rejection of the lot is made based on a specific rule. In lot-by-lot inspection by variables, one or more samples of items are drawn from a given lot; the measurement of a quality characteristic in each sampled item is recorded, and the decision of acceptance or rejection of the lot is made as a function of such measurements. The theory of inspection by variables is applicable when the quality characteristic of sampled items is measurable on a continuous scale and the functional form of the probability distribution is assumed to be known. A variables sampling is advantageous in the sense that it generates more information from each item inspected, requires small sample, and provides same protection when compared to attributes sampling. (See Bowker and Goode [1] and Montgomery [2]).On the basis of the implicit assumption that the quality characteristic is distributed according to normal with meanμ and standard deviation σ, the concept of variables sampling inspection has been studied by many researchers. Some of the early works on variables sampling inspection are seen in Lieberman and Resnikoff [3], Schilling [4], Owen [5, 6], and Hamaker [7]. Studies relating to sampling plans when the assumption of normality of the quality characteristic fails or the functional form of the underlying distribution deviates from normal or the form of the distribution is not known are also found in the literature of acceptance sampling. Duncan [8], Srivastava [9], Das and Mitra [10], Owen [11], Takagi [12], Guenther [13, 14], Zimmer and Burr [15], Aminzadeh [16], Suresh and Ramanathan [17], Vijayaraghavan and Geetha [18–20], and Geetha and Vijayaraghavan [21, 22] are few references which deal with variables inspection using nonnormal distributions.When the quality characteristicX follows normal distribution with mean μ and standard deviation σ, the problem of designing single sampling plans by variables is addressed with ease. In the industrial situations, quite often, the assumption of normality may not be valid or the quality characteristic may be distributed according to non-normal distributions. In such cases, the selection of variable sampling plans becomes complicated. However, the literature of acceptance sampling traces several works which deal with the designing of variables plans when the quality characteristic follows a probability distribution other than normal. Owen [11] presented a detailed survey on various works related to variable sampling plans with emphasis on nonnormality.In this paper a study on single sampling plans by variables is formulated under the assumption that the quality characteristics have a Pareto distribution. A procedure for determining the parameters of the proposed plan for specified requirements in terms of producer’s and consumer’s protection is also developed.


## 2. Single Sampling Inspection Plans by Variables

A single sampling inspection plan by variables is defined under the following assumptions.(a)
The quality characteristic, denoted byX, is measurable on a continuous scale and has a known form of probability distribution, represented by FX(x;μ,σ), which is the distribution function of X with mean μ and variance σ2.(b)
Each individual unit in a submitted lot has a one-sided specification say, lower specification,L or upper specification U. If, for a unit, X>U (or X<L), the unit is classified as a nonconforming unit.The operating procedure of a variable sampling plan is as follows.Step 1.
Draw a random sample ofn units from a lot and observe the measurements x1,x2,…,xn of the quality characteristic X.Step 2.
Whenσ is known, accept the lot if x-+kσ≤U (or x--kσ≥L); otherwise, reject the lot where x- is the sample mean. When σ is unknown, accept the lot if x-+ks≤U (or x--ks≥L), where s2=(1/(n-1))∑i:1n(xi-x-)2 is an unbiased estimate of σ2.Thus, a single sampling plan by variables is designated by two parameters, namely, the sample size,n, and the acceptability constant, k. When these parameters are known, the plan could be implemented. The explicit expressions for n and k can be derived by specifying two points on the operating characteristic curve of the plan, namely, (p1,1-α) and (p2,β), where p1 and p2 are termed as acceptable quality level (AQL), and the limiting quality level (LQL) associated with the producer’s risk α and the consumer’s risk β, respectively. A sampling plan by variables is termed as a known σ or unknown σ plan according to whether σ is known or unknown.


## 3. Operating Characteristic Function

An important measure of the performance of a variables sampling plan is its operating characteristic function, which is a function of the proportionp of non-conforming units called incoming lot quality, and it provides the probability, Pa(p) of acceptance of a lot. The plot of Pa(p) against p results in a curve called the operating characteristic (OC) curve. For a given upper specification limit U when σ is known, p and Pa(p) are defined by
(1)p=P(X>U∣μ),(2)Pa(p)=P(x-+kσ≤U∣μ).AQL and LQL, using (1), are defined by
(3)AQL=p1=P(X>U∣μ1),LQL=p2=P(X>U∣μ2),
where μ1 and μ2 are the means of the underlying distribution which results in AQL and LQL, respectively.Assume thatX is modeled by the Pareto distribution with parameters a and b. Then, its probability density function and the cumulative distribution function are, respectively, given by
(4)f(x)=abaxa+1,x>b,a<0,(5)F(x)=1-(bx)a.The mean and the variance of the underlying distribution are given byμ=ab/(a-1) for a>1 and σ2=ab2/(a-1)2(a-2) for a>2, respectively. Further, the measures of skewness and kurtosis, which are the measures for identifying the asymmetric nature of the distribution, are expressed as α3=(2(1+a)/(a-3))(a-2)/a  for a>3 and α4=3(a-2)(3a3+a+2)/a(a-3)(a-4), respectively. From (1) and (3), the lot quality levels p,  AQL, and LQL using standard Pareto distribution are, respectively, defined by
(6)p=P(T>Kp*),(7)AQL=p1=P(T>Kp1*),(8)LQL=p2=P(T>Kp2*),
where T=(X-μ)/σ, Kp*=(U-μ)/σ, Kp1*=(U-μ1)/σ, and Kp2*=(U-μ2)/σ.The producer’s riskα and the consumer’s risk β corresponding to AQL and LQL are, respectively, defined from (2) as
(9)α=P(x-+kσ>U∣μ=μ1),(10)1-β=P(x-+kσ>U∣μ=μ2).Whenσ is unknown, the estimate s is used in the decision criterion and hence in the evaluation of α and β.


## 4. Designing Single Sampling Plans by Variables

In the industrial practice the unknown standard deviation variables plans are more realistic than the known standard deviation variables plans. If the distribution is non-normal, then the designing of unknownσ plans is much complicated. Takagi [12] attempted to provide solution to such problems and proposed a methodology for determining the parameters of variables sampling plans under the conditions of non-normal populations by introducing an expansion factor in terms of the measures of skewness and kurtosis. The procedures of designing known and unknown standard deviation sampling plans giving protection to the producer and consumer are discussed in [18–22] under the assumption that the quality characteristics under study follow gamma, Laplace, Rayleigh and logistic distributions when the measures of skewness and/or kurtosis are specified.The methodology proposed in [12] using the expansion factor will, now, be discussed for an unknown sigma plan by variables under the assumption of Pareto distribution for the quality characteristic.In the case of unknown sigma plan, the determination ofn and k is usually based on the sampling distribution of  x-+ks or x--ks. When the underlying distribution is normal, x- and s are independently distributed as normal and hence x-+ks and x--ks are also normally distributed. By this property, explicit forms of n and k could be derived for determining n and k. A linear combination of independently distributed variables tends to be normally distributed even when the individual variables depart from normality (see [15]). Further, when the underlying distribution is non-normal, according to [12], the asymptotic distributions of x-+ks and x--ks are normal with the means, respectively, given by μy=μ+kσ and μy=μ-kσ, and with the common variance expressed by
(11)σY2=σ2n[1+k24(α4-1)±kα3],
where α3 and α4 represent the measures of skewness and kurtosis of the underlying distribution.It is known thatKp*=(U-μ)/σ. Thus, for specified p1 and p2,
(12)μ1=U-Kp1*σ,μ2=U-Kp2*σ.The operating characteristic functionPa(p) for an unknown standard deviation plan is expressed as
(13)Pa(p)=Pr[x-+kUs≤U∣p]=Pr(Y≤U∣p).
Assume that p1 and p2 are specified with the risks α of rejecting the lot and β of accepting the lot. Then, corresponding to α and β, the normal deviates Kα and Kβ are, respectively, defined by
(14)Kα=U-(μ1+kσ)(σ2/n)[1+(k2/4)(α4-1)+kα3]=U-(μ1+kUσ)σeU/nU,-Kβ=U-(μ2+kσ)(σ2/n)[1+(k2/4)(α4-1)+kα3]=U-(μ2+kUσ)σeU/nU.On substituting the expressions (12) for μ1 and μ2, (14) would result in the following:
(15)Kp1*=kU+KαeUnU,Kp2*=kU-KβeUnU.The mathematical expressions forn and k of an unknown standard deviation plan for specified requirements in terms of the points (p1,1-α) and (p2,β) on the OC curve such that Pa(p1)=1-α and Pa(p2)=β, when U is specified, are derived from (15) as
(16)nU=eU[Kα+KβKp1*-Kp2*]2,(17)kU=KαKp2*+KβKp1*Kα+Kβ,
where Kα and Kβ are the standardized deviates exceeded with the probabilities α and β, respectively, and eU is an expansion factor, which is a function of α3 and α4, expressed as eU=1+(kU2/4)(α4-1)+kUα3. The expansion factor is used to obtain the known standard deviation plans with parameters n′ and k′ by using the relations nU′=nU/eU and kU′=kU. In a similar way, when L is specified, the expressions for n and k are given by
(18)nL=eL[Kα+KβK1-p1*-K1-p2*]2,kL=-KαK1-p2*+KβK1-p1*Kα+Kβ,
where eL is expressed by eL=1+(kL2/4)(α4-1)-kLα3, using which the known standard deviation plan with parameters n′ and k′ could be derived as nL′=nL/eL and kL′=kL.


## 5. Numerical Example

Suppose that a set of measurements yieldsα3=2.2671 and α4=11.3626. It is desired to determine a variables sampling plan giving protection to the producer and the consumer in terms of (p1=0.01,α=0.05) and (p2=0.06,β=0.10). For the given requirements, the value of a is calculated from the expression of α3 of Pareto distribution as 25.  Corresponding to p1 and p2, the values of Kp1* and Kp2* are determined, respectively, from (7) and (8) as 3.6970 and 1.7828. The normal deviates Kα and Kβ are obtained as 1.64485 and 1.28155 by satisfying (9) and (10) for specified values of α and β. Substituting these values in (16) and (17), the parameters of the desired plan are determined as n=57.8256 which when rounded becomes 58 and k=2.6211. The value of eU is obtained as 24.7401. Thus, the parameters of a known standard deviation plan are computed as nU′=nU/eU=2.3373≈3 and kU′=2.6211.


## 6. A Note on Deriving the Plans Based on Exact Pareto Distribution

It is known that Pareto distribution is a very heavy-tailed distribution for which therth moment is finite when the shape parameter a is greater than r. Thus, large observations can be expected in samples from this distribution. Further, the distribution function and the maximum likelihood estimators of the parameters can be written in closed form. As the distribution is not in general symmetric and the results given in the previous sections are based only on the assumption of a symmetric distribution, it would be appropriate if one could derive variables sampling plans using the exact Pareto distribution. Whereas it would be reasonable to assume that the mean is finite, no assumption about finite moments has to be made if one uses the median in lieu of mean. It is to be noted that the median could be used as a measure of location when a distribution is skewed  (say, for instance, a Pareto distribution) or when end values are not known or when one requires reduced importance to be attached to outliers. When sample median is used in the operating procedure of a sampling plan for making a decision about the lot, the explicit expressions for the plan parameters are not plausible to derive. In such a case, the procedure for finding the parameters based on exact Pareto distribution will be based on the sampling distribution of sample median and will be much complicated.


## 7. Conclusion

Sampling inspection by variables is the methodology used for deciding about the disposition of a lot based on sample measurements of units on a quality characteristic under study.  Such a procedure is defined under the assumption that the quality characteristic is measurable on a continuous scale and the functional form of the probability distribution is known. Though sampling plans by variables under the implicit assumption that the quality characteristic is distributed as normal have been developed with the related properties, it should be realized that the assumption of normality is seldom hold in practice. Hence, it is inevitable to develop sampling plans by variables when the underlying quality characteristics follow non-normal distributions, which may be skewed or symmetric. In this paper a single sampling plan by variables is formulated and evaluated when the quality characteristic is distributed according to a Pareto distribution. The approximation of normal distribution for a linear combination of independent variables even when the individual variables depart from normality is effectively used in developing the procedure for determining the parameters of the plan when requirements for producer’s and the consumer’s protection are specified.



---

*Source: 808741-2013-09-22.xml*