Hellenic Complex Systems Laboratory

Size: px
Start display at page:

Download "Hellenic Complex Systems Laboratory"

Transcription

1 Hellenic Complex Systems Laboratory Technical Report No IV Calculation of the confidence bounds for the fraction nonconforming of normal populations of measurements in clinical laboratory medicine Aristides T. Hatjimihail, M.D., Ph.D. Drama 2015

2 HCSLTR04 Hellenic Complex Systems Laboratory, K. Palamas 21, GR Drama, Greece Copyright Aristides T. Hatjimihail,

3 1. Abstract The fraction nonconforming is a key quality measure used in statistical quality control design in clinical laboratory medicine. The confidence bounds of normal populations of measurements for the fraction nonconforming each of the lower and upper quality specification limits when both the random and the systematic error are unknown can be calculated using the noncentral t-distribution, as it is described in detail and illustrated with examples. 3

4 2. Introduction The quality control (QC) in clinical laboratory medicine is based mainly on the control charts for variables. Although the fraction nonconforming is used in the control charts for attributes, it is also a useful quality measure related to various process capability indices (Pearn & Kotz 2006). It is defined as the ratio of the items nonconforming to the total number of items in a subgroup (Juran et al. 1999). The quality control design process (Hatjimihail 1992) is based on the quality specifications of the maximum total allowable analytical error, the upper acceptable bound for the fraction nonconforming, the lower acceptable bounds for the probabilities for the critical random and systematic error detection, and the upper acceptable bound for the false rejection (CLSI 2012; Linnet 1989). It includes the following steps: 1. Estimation of the error of the measurement procedure based on a sample of n control measurements. 2. Calculation of the critical errors of the measurement procedure based on an upper acceptable bound for the expected fraction nonconforming the maximum total allowable analytical error quality specification. This bound is usually set equal to Definition of the decision bounds for the quality control rules based on the lower acceptable bounds for the probabilities for the critical errors detection and the upper acceptable bound for false rejection. These steps are performed at typically one to three levels of the measurand. Therefore, the expected fraction nonconforming is a key measure of the quality control design process. While this measure is routinely calculated in the clinical laboratories, very little attention has been given to its confidence bounds. 4

5 In general, for the quality control estimations in laboratory medicine it is assumed that the error of a measurement procedure is normally distributed with known mean and standard deviation. Usually in practice we only know the mean x and the standard deviation s of a sample of n control measurements at each of various levels of the measurand. As both the systematic and the random error of the measurement procedure are unknown lower and upper bounds for the tail probabilities of the lower (x LSL ), and the upper (x USL ) quality specifications of the measurement x can be calculated with confidence level γ using the noncentral t-distribution (1987), as described in detail in the Materials and Methods section. Therefore, the purpose of this technical report is to introduce the noncentral t-distribution into clinical laboratory medicine and to describe its application for the calculation of the confidence bounds of the fraction nonconforming each of the lower and upper quality specification limits of populations of normal measurements in a clinical laboratory setting. 3. Materials and Methods If Z a standard normal N (0, 1) deviate and V a chi-square deviate with v degrees of freedom, the distribution of the ratio is the noncentral t-distribution with ν degrees of freedom and noncentrality parameter δ. T, Z (3.1) V The cumulative density function of a noncentral t-distribution with ν degrees of freedom and noncentrality parameter δ is defined as (Owen 1965): 5

6 2 tx v1 F, ( t) x ( x) dx 1/2 ( 2)/2 ( / 2) 2 0 v (3.2) where ( x) the gamma function and: 2 x 2 1 () x e (3.3) 2 and x ( x) ( t) dt (3.4) Scholz (Scholz 1994) presents an elegant derivation of the lower and upper bounds for the lower tail probability p(x LSL ) of a normal population with unknown mean μ and unknown standard deviation σ. Briefly, xlsl pxlsl pxlsl,, P, X xlsl Φ σ (3.5) If we denote with pˆul an upper bound for p(x LSL ) with confidence level γ then ˆ, UL LSL P p p x (3.6) Let / Z n X (3.7) a standard normal deviate, and 6

7 2 2 1 / V n S (3.8) a chi-square deviate with n-1 degrees of freedom. If n x LSL 1 nφ p xlsl (3.9) then, n x LSL s n X n x LSL x 2 n X 2 n n 1 n x LSL n1 S s Z T V 1 n1, (3.10) where the Tn 1, has by definition a noncentral t-distribution with n-1 degrees of freedom and noncentrality parameter δ. According to the probability integral transformation (Qeuesenberry 2006) the random variable U is uniformly distributed over (0,1). Therefore, n x X U F F T S LSL n1, n1, n1, (3.11) PU ( 1 ) (3.12) As x Fn 1, is a decreasing function of δ we have 7

8 F n1, n xlsl X 1 ˆ S (3.13) where ˆ is estimated by the following equation: F n x X LSL 1, ˆ 1 n S (3.14) Consequently, the upper bound with confidence level γ for the fraction non conforming the lower quality specification limit is calculated as ˆ p UL ˆ Φ n (3.15) As pˆ ˆ UL LSL pul LSL P p x P p x (3.16),, 1 the lower bound with confidence level γ for the fraction non conforming the lower quality specification limit is calculated as p pˆ ˆ 1 (3.17) LL UL In addition, if we denote with pˆuu fraction non conforming the upper quality specification limit then the upper bound with confidence level γ for the P, pˆ UU p xusl p x P ˆ, 1 puu 1 USL 1 (3.18) Therefore, the lower and upper bounds with confidence level γ for the fraction non conforming the upper quality specification limit are calculated respectively as 8

9 pˆ 1 pˆ (3.19) ' LU UU and pˆ ' UU 1 pˆ 1 (3.20) UU If we denote with c the true value of the measurand of the control material, and e max the maximum total clinically allowable analytical error (Linnet 1989), expressed as a ratio to the true value, then x c(1 e ) (3.21) LSL max and x c(1 e ) (3.22) USL max To illustrate the application of the noncentral t-distribution for the calculation of the confidence bounds of the fraction nonconforming of normal populations of measurements in a clinical laboratory setting I calculated upper and lower confidence bounds for the fraction nonconforming the upper and lower quality specification limits of populations of normal measurements, assuming unknown systematic and random errors of the measurement procedure. The calculations were performed by MATLAB, ver. 2012b, using the method described by Owen (Owen 1965). It was assumed that a reference control material with a measurand with an assigned value of units has been measured repeatedly by the measuring procedure and that the sample mean and standard deviation of the measurand were estimated from samples of 20 and 40 control measurements. It was also assumed that both the systematic and the random error of the measurement procedure were unknown and that the uncertainty of the assigned value of 9

10 the measurand of the control material was negligible. The lower and upper confidence bounds for the fraction nonconforming each of the lower and upper quality specifications of the population of the measurements were calculated by the noncentral t-distribution as described previously. The parameters of these calculations are presented in the Table 1. figures tables assigned value of the measurand μ of the control material sample mean x increments of the sample mean sample s.d. s increments of the sample s.d total allowable analytical error e max confidence level γ sample size n Table 1: The parameters of the calculations of the confidence bounds for the fraction nonconforming of the populations of the normal control measurements. 4. Results The results of the study are described in the following subsections 4.1 and 4.2. They are presented in the tables 2-5 and in the figures 1-8. They are also presented in the respective MATLAB figures included in the supporting information file MATLAB_Figures.zip so that the data can be explored interactively in the MATLAB environment. 10

11 parameters of the sample confidence bounds for the fraction nonconforming at a 95% confidence level size of the control measurements sample n = 20 total allowable analytical error taae = 0.10 fraction nonconforming the lower quality specification limit fraction nonconforming the upper quality specification limit x s expected lower bound upper bound expected lower bound upper bound ,62E-24 2,14E-33 2,85E-15 7,62E-24 2,14E-33 2,85E ,87E-07 8,69E-10 5,32E-05 2,87E-07 8,69E-10 5,32E ,29E-04 2,63E-05 5,36E-03 4,29E-04 2,63E-05 5,36E ,21E-03 1,09E-03 2,99E-02 6,21E-03 1,09E-03 2,99E ,28E-02 6,56E-03 7,01E-02 2,28E-02 6,56E-03 7,01E ,91E-28 5,78E-40 4,15E-18 1,13E-19 1,92E-27 1,07E ,90E-08 1,85E-11 9,75E-06 3,40E-06 2,88E-08 2,51E ,23E-04 4,53E-06 2,42E-03 1,35E-03 1,31E-04 1,12E ,98E-03 3,94E-04 1,86E-02 1,22E-02 2,79E-03 4,66E ,39E-02 3,33E-03 5,07E-02 3,59E-02 1,22E-02 9,50E ,78E-33 3,76E-47 3,30E-21 6,22E-16 4,14E-22 2,19E ,87E-10 2,78E-13 1,54E-06 3,17E-05 6,72E-07 1,02E ,17E-05 6,72E-07 1,02E-03 3,83E-03 5,59E-04 2,19E ,35E-03 1,31E-04 1,12E-02 2,28E-02 6,56E-03 7,01E ,20E-03 1,61E-03 3,58E-02 5,48E-02 2,17E-02 1,26E ,12E-39 5,81E-55 1,43E-24 1,28E-12 2,17E-17 2,48E ,02E-11 2,93E-15 2,10E-07 2,33E-04 1,11E-05 3,63E ,34E-06 8,54E-08 4,07E-04 9,82E-03 2,06E-03 4,03E ,77E-04 3,98E-05 6,47E-03 4,01E-02 1,42E-02 1,02E ,66E-03 7,34E-04 2,48E-02 8,08E-02 3,67E-02 1,64E ,79E-45 2,18E-63 3,37E-28 9,87E-10 2,78E-13 1,54E ,28E-12 2,17E-17 2,48E-08 1,35E-03 1,31E-04 1,12E ,53E-06 9,31E-09 1,52E-04 2,28E-02 6,56E-03 7,01E ,33E-04 1,11E-05 3,63E-03 6,68E-02 2,84E-02 1,44E ,56E-03 3,18E-04 1,68E-02 1,15E-01 5,88E-02 2,10E ,67E-51 1,96E-72 4,33E-32 2,87E-07 8,69E-10 5,32E ,19E-14 1,13E-19 2,51E-09 6,21E-03 1,09E-03 2,99E ,87E-07 8,69E-10 5,32E-05 4,78E-02 1,81E-02 1,15E ,84E-05 2,85E-06 1,96E-03 1,06E-01 5,25E-02 1,98E ,35E-03 1,31E-04 1,12E-02 1,59E-01 8,99E-02 2,63E-01 Table 2: Confidence bounds for the fraction nonconforming the lower and upper quality specification limits assuming a true value of the measurand of the control material equal to measurement units, a sample of 20 control measurements, and a total allowable analytical error equal to

12 parameters of the sample confidence bounds for the fraction nonconforming at a 95% confidence level size of the control measurements sample n = 40 total allowable analytical error taae = 0.10 fraction nonconforming the lower quality specification limit fraction nonconforming the upper quality specification limit x s expected lower bound upper bound expected lower bound upper bound ,62E-24 1,85E-30 1,05E-17 7,62E-24 1,85E-30 1,05E ,87E-07 5,14E-09 1,20E-05 2,87E-07 5,14E-09 1,20E ,29E-04 6,17E-05 2,61E-03 4,29E-04 6,17E-05 2,61E ,21E-03 1,86E-03 1,91E-02 6,21E-03 1,86E-03 1,91E ,28E-02 9,58E-03 5,09E-02 2,28E-02 9,58E-03 5,09E ,91E-28 2,01E-36 4,82E-21 1,13E-19 4,68E-25 1,12E ,90E-08 1,55E-10 1,65E-06 3,40E-06 1,24E-07 7,38E ,23E-04 1,25E-05 1,04E-03 1,35E-03 2,67E-04 6,12E ,98E-03 7,31E-04 1,10E-02 1,22E-02 4,38E-03 3,18E ,39E-02 5,15E-03 3,51E-02 3,59E-02 1,70E-02 7,21E ,78E-33 5,97E-43 1,08E-24 6,22E-16 3,27E-20 5,83E ,87E-10 3,41E-12 1,91E-07 3,17E-05 2,19E-06 3,81E ,17E-05 2,19E-06 3,81E-04 3,83E-03 1,01E-03 1,33E ,35E-03 2,67E-04 6,12E-03 2,28E-02 9,58E-03 5,09E ,20E-03 2,64E-03 2,36E-02 5,48E-02 2,88E-02 9,97E ,12E-39 4,90E-50 1,19E-28 1,28E-12 6,34E-16 1,50E ,02E-11 5,46E-14 1,85E-08 2,33E-04 2,82E-05 1,66E ,34E-06 3,35E-07 1,30E-04 9,82E-03 3,32E-03 2,70E ,77E-04 9,02E-05 3,25E-03 4,01E-02 1,95E-02 7,84E ,66E-03 1,29E-03 1,54E-02 8,08E-02 4,66E-02 1,34E ,79E-45 1,10E-57 6,39E-33 9,87E-10 3,41E-12 1,91E ,28E-12 6,34E-16 1,50E-09 1,35E-03 2,67E-04 6,12E ,53E-06 4,45E-08 4,11E-05 2,28E-02 9,58E-03 5,09E ,33E-04 2,82E-05 1,66E-03 6,68E-02 3,69E-02 1,16E ,56E-03 6,01E-04 9,84E-03 1,15E-01 7,22E-02 1,77E ,67E-51 6,80E-66 1,67E-37 2,87E-07 5,14E-09 1,20E ,19E-14 5,35E-18 1,02E-10 6,21E-03 1,86E-03 1,91E ,87E-07 5,14E-09 1,20E-05 4,78E-02 2,43E-02 8,97E ,84E-05 8,17E-06 8,12E-04 1,06E-01 6,50E-02 1,66E ,35E-03 2,67E-04 6,12E-03 1,59E-01 1,07E-01 2,28E-01 Table 3: Confidence bounds for the fraction nonconforming the lower and upper quality specification limits assuming a true value of the measurand of the control material equal to measurement units, a sample of 40 control measurements, and a total allowable analytical error equal to

13 parameters of the sample confidence bounds for the fraction nonconforming at a 95% confidence level size of the control measurements sample n = 20 total allowable analytical error taae = 0.20 fraction nonconforming the lower quality specification limit fraction nonconforming the upper quality specification limit x s expected lower bound upper bound expected lower bound upper bound ,75E-89 5,76E-127 1,62E-55 2,75E-89 5,76E-127 1,62E ,62E-24 2,14E-33 2,85E-15 7,62E-24 2,14E-33 2,85E ,31E-11 5,94E-16 1,05E-07 1,31E-11 5,94E-16 1,05E ,87E-07 8,69E-10 5,32E-05 2,87E-07 8,69E-10 5,32E ,17E-05 6,72E-07 1,02E-03 3,17E-05 6,72E-07 1,02E ,28E-98 9,84E-140 5,34E-61 8,53E-81 8,09E-115 2,66E ,32E-26 1,33E-36 1,17E-16 1,05E-21 2,42E-30 5,95E ,28E-12 2,17E-17 2,48E-08 1,20E-10 1,39E-14 4,16E ,60E-08 1,32E-10 2,32E-05 1,02E-06 5,22E-09 1,18E ,33E-05 1,99E-07 5,94E-04 7,23E-05 2,15E-06 1,73E ,44E-107 4,03E-153 9,58E-67 9,74E-73 2,72E-103 2,37E ,91E-28 5,78E-40 4,15E-18 1,13E-19 1,92E-27 1,07E ,12E-13 6,79E-19 5,48E-09 9,87E-10 2,78E-13 1,54E ,90E-08 1,85E-11 9,75E-06 3,40E-06 2,88E-08 2,51E ,41E-06 5,55E-08 3,36E-04 1,59E-04 6,52E-06 2,85E ,33E-117 3,95E-167 9,31E-73 4,11E-65 2,19E-92 1,15E ,60E-31 1,76E-43 1,26E-19 9,48E-18 1,06E-24 1,65E ,83E-15 1,82E-20 1,13E-09 7,28E-09 4,75E-12 5,36E ,46E-09 2,37E-12 3,95E-06 1,07E-05 1,45E-07 5,16E ,11E-06 1,47E-08 1,86E-04 3,37E-04 1,87E-05 4,59E ,39E-127 9,26E-182 4,91E-79 6,39E-58 4,23E-82 3,03E ,78E-33 3,76E-47 3,30E-21 6,22E-16 4,14E-22 2,19E ,22E-16 4,14E-22 2,19E-10 4,82E-08 6,94E-11 1,74E ,87E-10 2,78E-13 1,54E-06 3,17E-05 6,72E-07 1,02E ,93E-07 3,67E-09 1,01E-04 6,87E-04 5,08E-05 7,24E ,06E-138 5,20E-197 1,41E-85 3,67E-51 1,96E-72 4,33E ,73E-36 5,56E-51 7,40E-23 3,19E-14 1,13E-19 2,51E ,93E-17 8,08E-24 3,98E-11 2,87E-07 8,69E-10 5,32E ,05E-10 2,98E-14 5,80E-07 8,84E-05 2,85E-06 1,96E ,87E-07 8,69E-10 5,32E-05 1,35E-03 1,31E-04 1,12E-02 Table 4: Confidence bounds for the fraction nonconforming the lower and upper quality specification limits assuming a true value of the measurand of the control material equal to measurement units, a sample of 20 control measurements, and a total allowable analytical error equal to

14 parameters of the sample confidence bounds for the fraction nonconforming at a 95% confidence level size of the sample of the measurements n = 40 total allowable analytical error taae = 0.20 fraction nonconforming the lower quality specification limit fraction nonconforming the upper quality specification limit x s expected lower bound upper bound expected lower bound upper bound ,75E-89 2,23E-115 4,21E-65 2,75E-89 2,23E-115 4,21E ,62E-24 1,85E-30 1,05E-17 7,62E-24 1,85E-30 1,05E ,31E-11 1,28E-14 8,16E-09 1,31E-11 1,28E-14 8,16E ,87E-07 5,14E-09 1,20E-05 2,87E-07 5,14E-09 1,20E ,17E-05 2,19E-06 3,81E-04 3,17E-05 2,19E-06 3,81E ,28E-98 5,78E-127 1,47E-71 8,53E-81 2,34E-104 5,88E ,32E-26 2,26E-33 2,46E-19 1,05E-21 1,09E-27 3,74E ,28E-12 6,34E-16 1,50E-09 1,20E-10 2,24E-13 4,10E ,60E-08 9,30E-10 4,56E-06 1,02E-06 2,63E-08 3,05E ,33E-05 7,21E-07 2,02E-04 7,23E-05 6,32E-06 7,01E ,44E-107 4,11E-139 2,48E-78 9,74E-73 6,75E-94 4,00E ,91E-28 2,01E-36 4,82E-21 1,13E-19 4,68E-25 1,12E ,12E-13 2,72E-17 2,55E-10 9,87E-10 3,41E-12 1,91E ,90E-08 1,55E-10 1,65E-06 3,40E-06 1,24E-07 7,38E ,41E-06 2,26E-07 1,04E-04 1,59E-04 1,73E-05 1,25E ,33E-117 7,98E-152 2,04E-85 4,11E-65 5,33E-84 1,32E ,60E-31 1,28E-39 7,90E-23 9,48E-18 1,45E-22 2,78E ,83E-15 1,01E-18 4,01E-11 7,28E-09 4,51E-11 8,21E ,46E-09 2,40E-11 5,74E-07 1,07E-05 5,42E-07 1,71E ,11E-06 6,74E-08 5,21E-05 3,37E-04 4,53E-05 2,18E ,39E-127 4,24E-165 8,18E-93 6,39E-58 1,15E-74 2,13E ,78E-33 5,97E-43 1,08E-24 6,22E-16 3,27E-20 5,83E ,22E-16 3,27E-20 5,83E-12 4,82E-08 5,17E-10 3,27E ,87E-10 3,41E-12 1,91E-07 3,17E-05 2,19E-06 3,81E ,93E-07 1,91E-08 2,54E-05 6,87E-04 1,13E-04 3,70E ,06E-138 6,16E-179 1,59E-100 3,67E-51 6,80E-66 1,67E ,73E-36 2,01E-46 1,24E-26 3,19E-14 5,35E-18 1,02E ,93E-17 9,16E-22 7,82E-13 2,87E-07 5,14E-09 1,20E ,05E-10 4,49E-13 6,07E-08 8,84E-05 8,17E-06 8,12E ,87E-07 5,14E-09 1,20E-05 1,35E-03 2,67E-04 6,12E-03 Table 5: Confidence bounds for the fraction nonconforming the lower and upper quality specification limits assuming a true value of the measurand of the control material equal to measurement units, a sample of 40 control measurements, and a total allowable analytical error equal to Fraction nonconforming the lower quality specification limit Figures 1, 3, 5, and 7 present the expected fraction nonconforming the lower quality specification limit (surface II) and the lower (surface III) and upper (surface I) confidence bounds for the fraction nonconforming at a 95% confidence level versus the sample mean and the sample standard deviation. The confidence bounds were calculated using the noncentral t- distribution. It was assumed a known true value of the measurand of the control material equal to measurement units. For the figures 1 and 3 it was assumed a total allowable analytical error equal to 0.10 while for the figures 5 and 7 it was assumed a total allowable 14

15 analytical error equal to The respective samples were samples of 20 control measurements for the figures 1 and 5, and samples of 40 control measurements for the figures 3 and Fraction nonconforming the upper quality specification limit Figures 2, 4, 6, and 8 present the expected fraction nonconforming the upper quality specification limit (surface II) and the lower (surface III) and upper (surface I) confidence bounds for the fraction nonconforming at a 95% confidence level versus the sample mean and the sample standard deviation. The confidence bounds were calculated using the noncentral t- distribution. It was assumed a known true value of the measurand of the control material equal to measurement units. For the figures 2 and 4 it was assumed a total allowable analytical error equal to 0.10 while for the figures 6 and 8 it was assumed a total allowable analytical error equal to The respective samples were samples of 20 control measurements for the figures 2 and 6, and samples of 40 control measurements for the figures 4 and 8. The results are summarized in the tables 2 and 3, assuming a total allowable analytical error equal to 0.10, and in the tables 4 and 5, assuming a total allowable analytical error equal to The respective samples were samples of 20 control measurements for the tables 2 and 4, and samples of 40 control measurements for the tables 3 and 5. The confidence bounds increase with the estimated systematic and random error of the measurements and decrease with the sample size and the total allowable analytical error. Although there have been assumed samples with a positive systematic error, it can be easily shown that the bounds for the fraction nonconforming the upper quality specification limit for a systematic error se are equal to the respective bounds for the fraction nonconforming the lower quality specification limit for a systematic error se. 15

16 Figure 1: Confidence bounds for the fraction nonconforming the lower quality specification limit. The expected fraction nonconforming the lower quality specification limit (surface II) and the lower (surface III) and upper (surface I) confidence bounds for the fraction nonconforming at a 95% confidence level versus the sample mean and the sample standard deviation. It is assumed a known true value of the measurand of the control material equal to measurement units, a sample of 20 control measurements, and a total allowable analytical error equal to

17 Figure 2: Confidence bounds for the fraction nonconforming the upper quality specification limit. The expected fraction nonconforming the upper quality specification limit (surface II) and the lower (surface III) and upper (surface I) confidence bounds for the fraction nonconforming at a 95% confidence level versus the sample mean and the sample standard deviation. It is assumed a known true value of the measurand of the control material equal to measurement units, a sample of 20 control measurements, and a total allowable analytical error equal to

18 Figure 3: Confidence bounds for the fraction nonconforming the lower quality specification limit. The expected fraction nonconforming the lower quality specification limit (surface II) and the lower (surface III) and upper (surface I) confidence bounds for the fraction nonconforming at a 95% confidence level versus the sample mean and the sample standard deviation. It is assumed a known true value of the measurand of the control material equal to measurement units, a sample of 40 control measurements, and a total allowable analytical error equal to

19 Figure 4: Confidence bounds for the fraction nonconforming the upper quality specification limit. The expected fraction nonconforming the upper quality specification limit (surface II) and the lower (surface III) and upper (surface I) confidence bounds for the fraction nonconforming at a 95% confidence level versus the sample mean and the sample standard deviation. It is assumed a known true value of the measurand of the control material equal to measurement units, a sample of 40 control measurements, and a total allowable analytical error equal to

20 Figure 5: Confidence bounds for the fraction nonconforming the lower quality specification limit. The expected fraction nonconforming the lower quality specification limit (surface II) and the lower (surface III) and upper (surface I) confidence bounds for the fraction nonconforming at a 95% confidence level versus the sample mean and the sample standard deviation. It is assumed a known true value of the measurand of the control material equal to measurement units, a sample of 20 control measurements, and a total allowable analytical error equal to

21 Figure 6: Confidence bounds for the fraction nonconforming the upper quality specification limit. The expected fraction nonconforming the upper quality specification limit (surface II) and the lower (surface III) and upper (surface I) confidence bounds for the fraction nonconforming at a 95% confidence level versus the sample mean and the sample standard deviation. It is assumed a known true value of the measurand of the control material equal to measurement units, a sample of 20 control measurements, and a total allowable analytical error equal to

22 Figure 7: Confidence bounds for the fraction nonconforming the lower quality specification limit. The expected fraction nonconforming the lower quality specification limit (surface II) and the lower (surface III) and upper (surface I) confidence bounds for the fraction nonconforming at a 95% confidence level versus the sample mean and the sample standard deviation. It is assumed a known true value of the measurand of the control material equal to measurement units, a sample of 40 control measurements, and a total allowable analytical error equal to

23 Figure 8: Confidence bounds for the fraction nonconforming the upper quality specification limit. The expected fraction nonconforming the upper quality specification limit (surface II) and the lower (surface III) and upper (surface I) confidence bounds for the fraction nonconforming at a 95% confidence level versus the sample mean and the sample standard deviation. It is assumed a known true value of the measurand of the control material equal to measurement units, a sample of 40 control measurements, and a total allowable analytical error equal to

24 5. Discussion In clinical laboratory medicine the fraction nonconforming besides being a quality measure related to various process capability indices, it is used for the definition of the critical errors during the QC design process. To estimate the fraction nonconforming the upper and lower quality specification bounds for the measurement error are needed. These are defined by the maximum total clinically allowable analytical error (Linnet 1989). As both the systematic and the random error of any measurement procedure are actually unknown, we can calculate the respective confidence bounds for the fraction nonconforming each of the upper and lower quality specification limits, using the noncentral t-distribution. The calculation is mathematically complicated and computationally intensive. It is not possible to calculate the confidence bounds for the total fraction nonconforming either the upper or the lower quality specification limits, without additional assumptions (Owen 1965). The results of the study show that the difference of the confidence bounds for the fraction nonconforming from the expected fraction nonconforming is considerable for greater random and systematic errors and for smaller sample sizes, even when we assume an insignificant uncertainty of the assigned value of the measurand of the control material. Consequently it is considerable when the expected fraction nonconforming equals its acceptable upper bound that is usually set equal to Therefore, the uncertainty of the calculation of the critical errors is considerable as well, as Gernand and Wojciechit have pointed out previously (Gernand 2006). Accordingly, the uncertainty of the whole quality control design process is significant, even assuming a negligible uncertainty of the assigned value of the measurand of the control material. Although it was assumed a reference control material with a measurand with an assigned value of , the tables and figures of this study can be used for the estimation of the 24

25 confidence bounds for the fraction nonconforming for any measurand with any assigned value. 6. Conclusion The confidence bounds for the fraction nonconforming are useful quality measures, as the expected fraction nonconforming is a key measure of the quality design process. The confidence bounds for the fraction nonconforming each of the upper and lower quality specification limits can be calculated using the noncentral t-distribution. 7. Acknowledgements I am deeply grateful to Dr. Fritz Scholz for providing me with the technical report on the applicatios of the noncentral t- distribution (Scholz 1994), that was the basis of this technical report, as well as for his suggestions and corrections. 25

26 8. References Sampling distributions. In: Ito K, ed. The encyclopedic dictionary of mathematics. 2nd ed. Boston, Massachusetts: The MIT Press, CLSI Expression of Uncertainty in Measurement in Laboratory Medicine. Wayne, Pennsylvania: Clinical and Laboratory Standards Institute. p Gernand W How Reliable Are Critical Error Calculations? Clinical Chemistry 52: Hatjimihail AT A tool for the design and evaluation of alternative quality-control procedures [published erratum appears in Clin Chem 1992 Nov;38(11):2360] [see comments]. Clin Chem 38: Juran JM, Godfrey AB, Hoogstoel RE, and Schilling EG Juran's quality handbook: McGraw Hill. Linnet K Choosing quality-control systems to detect maximum clinically allowable analytical errors. Clin Chem 35: Owen DB A special case of a bivariate non-central t-distribution. Biometrika 52: Pearn W, and Kotz S Encyclopedia and handbook of process capability indices: World Scientific Publishing. Qeuesenberry CP Probability integral transformations. In: Kotz S, Balakrishnan N, Read C, and Vidakovic B, eds. Encyclopedia of statistical sciences. 2nd ed. Hoboken, New Jersey: John Wiley and Sons, Scholz FW Applications of the Noncentral t-distribution. Mathematics and Engineering Analysis Library Boeing Computer Services. p

Acceptance Sampling by Variables

Acceptance Sampling by Variables Acceptance Sampling by Variables Advantages of Variables Sampling o Smaller sample sizes are required o Measurement data usually provide more information about the manufacturing process o When AQLs are

More information

Part One of this article (1) introduced the concept

Part One of this article (1) introduced the concept Establishing Acceptance Limits for Uniformity of Dosage Units: Part Two Pramote Cholayudth The concept of sampling distribution of acceptance value (AV) was introduced in Part One of this article series.

More information

What is Process Capability?

What is Process Capability? 6. Process or Product Monitoring and Control 6.1. Introduction 6.1.6. What is Process Capability? Process capability compares the output of an in-control process to the specification limits by using capability

More information

Bootstrap Confidence Interval of the Difference Between Two Process Capability Indices

Bootstrap Confidence Interval of the Difference Between Two Process Capability Indices Int J Adv Manuf Technol (2003) 21:249 256 Ownership and Copyright 2003 Springer-Verlag London Limited Bootstrap Confidence Interval of the Difference Between Two Process Capability Indices J.-P. Chen 1

More information

Computing Optimal Strata Bounds Using Dynamic Programming

Computing Optimal Strata Bounds Using Dynamic Programming Computing Optimal Strata Bounds Using Dynamic Programming Eric Miller Summit Consulting, LLC 7/27/2012 1 / 19 Motivation Sampling can be costly. Sample size is often chosen so that point estimates achieve

More information

Cpk: What is its Capability? By: Rick Haynes, Master Black Belt Smarter Solutions, Inc.

Cpk: What is its Capability? By: Rick Haynes, Master Black Belt Smarter Solutions, Inc. C: What is its Capability? By: Rick Haynes, Master Black Belt Smarter Solutions, Inc. C is one of many capability metrics that are available. When capability metrics are used, organizations typically provide

More information

Modified S-Control Chart for Specified value of Cp

Modified S-Control Chart for Specified value of Cp American International Journal of Research in Science, Technology, Engineering & Mathematics Available online at http://www.iasir.net ISSN (Print): 38-349, ISSN (Online): 38-358, ISSN (CD-ROM): 38-369

More information

Econ 3790: Business and Economics Statistics. Instructor: Yogesh Uppal

Econ 3790: Business and Economics Statistics. Instructor: Yogesh Uppal Econ 3790: Business and Economics Statistics Instructor: Yogesh Uppal Email: yuppal@ysu.edu Chapter 8: Interval Estimation Population Mean: Known Population Mean: Unknown Margin of Error and the Interval

More information

Pre-control and Some Simple Alternatives

Pre-control and Some Simple Alternatives Pre-control and Some Simple Alternatives Stefan H. Steiner Dept. of Statistics and Actuarial Sciences University of Waterloo Waterloo, N2L 3G1 Canada Pre-control, also called Stoplight control, is a quality

More information

nquery Sample Size & Power Calculation Software Validation Guidelines

nquery Sample Size & Power Calculation Software Validation Guidelines nquery Sample Size & Power Calculation Software Validation Guidelines Every nquery sample size table, distribution function table, standard deviation table, and tablespecific side table has been tested

More information

= = P. IE 434 Homework 2 Process Capability. Kate Gilland 10/2/13. Figure 1: Capability Analysis

= = P. IE 434 Homework 2 Process Capability. Kate Gilland 10/2/13. Figure 1: Capability Analysis Kate Gilland 10/2/13 IE 434 Homework 2 Process Capability 1. Figure 1: Capability Analysis σ = R = 4.642857 = 1.996069 P d 2 2.326 p = 1.80 C p = 2.17 These results are according to Method 2 in Minitab.

More information

CONFIDENCE INTERVALS FROM A SINGLE SIMULATION USING THE DELTA METHOD 1

CONFIDENCE INTERVALS FROM A SINGLE SIMULATION USING THE DELTA METHOD 1 CONFIDENCE INTERVALS FROM A SINGLE SIMULATION USING THE DELTA METHOD Christoph Roser Masaru Nakano Toyota Central Research and Development Laboratories Nagakute, Aichi, 480-9, JAPAN croser@robotics.tytlabs.co.jp

More information

SPSS Basics for Probability Distributions

SPSS Basics for Probability Distributions Built-in Statistical Functions in SPSS Begin by defining some variables in the Variable View of a data file, save this file as Probability_Distributions.sav and save the corresponding output file as Probability_Distributions.spo.

More information

Process capability for a complete electronic product assembly

Process capability for a complete electronic product assembly Rochester Institute of Technology RIT Scholar Works Theses Thesis/Dissertation Collections 6-206 Process capability for a complete electronic product assembly Flavia Carvalho Resende fcr4827@rit.edu Follow

More information

Generating random samples from user-defined distributions

Generating random samples from user-defined distributions The Stata Journal (2011) 11, Number 2, pp. 299 304 Generating random samples from user-defined distributions Katarína Lukácsy Central European University Budapest, Hungary lukacsy katarina@phd.ceu.hu Abstract.

More information

Chapter 3. Bootstrap. 3.1 Introduction. 3.2 The general idea

Chapter 3. Bootstrap. 3.1 Introduction. 3.2 The general idea Chapter 3 Bootstrap 3.1 Introduction The estimation of parameters in probability distributions is a basic problem in statistics that one tends to encounter already during the very first course on the subject.

More information

AC : DETERMINING PROCESS CAPABILITY OF AN INDUSTRIAL PROCESS IN LABORATORY USING COMPUTER AIDED HARDWARE AND SOFTWARE TOOLS

AC : DETERMINING PROCESS CAPABILITY OF AN INDUSTRIAL PROCESS IN LABORATORY USING COMPUTER AIDED HARDWARE AND SOFTWARE TOOLS AC 007-150: DETERMINING PROCESS CAPABILITY OF AN INDUSTRIAL PROCESS IN LABORATORY USING COMPUTER AIDED HARDWARE AND SOFTWARE TOOLS Akram Hossain, Purdue University-Calumet Akram Hossain is a professor

More information

Study of t Resolution Function

Study of t Resolution Function Belle-note 383 Study of t Resolution Function Takeo Higuchi and Hiroyasu Tajima Department of Physics, University of Tokyo (January 6, 200) Abstract t resolution function is studied in detail. It is used

More information

John A. Conte, P.E. 2/22/2012 1

John A. Conte, P.E. 2/22/2012 1 John A. Conte, P.E. 2/22/2012 1 Objectives Excited to be here! Students, faculty, engineers Share my engineering career Some thoughts on Six Sigma Some thoughts on Process Capability Cp, Cpk, Pp and Ppk

More information

SIMULATION AND MONTE CARLO

SIMULATION AND MONTE CARLO JHU course no. 550.790, Modeling, Simulation, and Monte Carlo SIMULATION AND MONTE CARLO Some General Principles James C. Spall Johns Hopkins University Applied Physics Laboratory September 2004 Overview

More information

We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors

We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists 4,000 116,000 120M Open access books available International authors and editors Downloads Our

More information

QUALITY AND RELIABILITY CORNER IJQRM 21,1

QUALITY AND RELIABILITY CORNER IJQRM 21,1 The Emerald Research Register for this journal is available at www.em eraldinsight.com/res earchregister The current issue and full text archive of this journal is available at www.em eraldinsight.com/0265-671x.htm

More information

Sharif University of Technology. Scientia Iranica Transactions E: Industrial Engineering

Sharif University of Technology. Scientia Iranica Transactions E: Industrial Engineering Scientia Iranica E (214) 21(3), 144{15 Sharif University of Technology Scientia Iranica Transactions E: Industrial Engineering www.scientiairanica.com Developing functional process capability indices for

More information

Evaluating Robot Systems

Evaluating Robot Systems Evaluating Robot Systems November 6, 2008 There are two ways of constructing a software design. One way is to make it so simple that there are obviously no deficiencies. And the other way is to make it

More information

Z-TEST / Z-STATISTIC: used to test hypotheses about. µ when the population standard deviation is unknown

Z-TEST / Z-STATISTIC: used to test hypotheses about. µ when the population standard deviation is unknown Z-TEST / Z-STATISTIC: used to test hypotheses about µ when the population standard deviation is known and population distribution is normal or sample size is large T-TEST / T-STATISTIC: used to test hypotheses

More information

LOCATION AND DISPERSION EFFECTS IN SINGLE-RESPONSE SYSTEM DATA FROM TAGUCHI ORTHOGONAL EXPERIMENTATION

LOCATION AND DISPERSION EFFECTS IN SINGLE-RESPONSE SYSTEM DATA FROM TAGUCHI ORTHOGONAL EXPERIMENTATION Proceedings of the International Conference on Manufacturing Systems ICMaS Vol. 4, 009, ISSN 184-3183 University POLITEHNICA of Bucharest, Machine and Manufacturing Systems Department Bucharest, Romania

More information

StatsMate. User Guide

StatsMate. User Guide StatsMate User Guide Overview StatsMate is an easy-to-use powerful statistical calculator. It has been featured by Apple on Apps For Learning Math in the App Stores around the world. StatsMate comes with

More information

BESTFIT, DISTRIBUTION FITTING SOFTWARE BY PALISADE CORPORATION

BESTFIT, DISTRIBUTION FITTING SOFTWARE BY PALISADE CORPORATION Proceedings of the 1996 Winter Simulation Conference ed. J. M. Charnes, D. J. Morrice, D. T. Brunner, and J. J. S\vain BESTFIT, DISTRIBUTION FITTING SOFTWARE BY PALISADE CORPORATION Linda lankauskas Sam

More information

ECONOMIC DESIGN OF STATISTICAL PROCESS CONTROL USING PRINCIPAL COMPONENTS ANALYSIS AND THE SIMPLICIAL DEPTH RANK CONTROL CHART

ECONOMIC DESIGN OF STATISTICAL PROCESS CONTROL USING PRINCIPAL COMPONENTS ANALYSIS AND THE SIMPLICIAL DEPTH RANK CONTROL CHART ECONOMIC DESIGN OF STATISTICAL PROCESS CONTROL USING PRINCIPAL COMPONENTS ANALYSIS AND THE SIMPLICIAL DEPTH RANK CONTROL CHART Vadhana Jayathavaj Rangsit University, Thailand vadhana.j@rsu.ac.th Adisak

More information

Simulation Studies of the Basic Packet Routing Problem

Simulation Studies of the Basic Packet Routing Problem Simulation Studies of the Basic Packet Routing Problem Author: Elena Sirén 48314u Supervisor: Pasi Lassila February 6, 2001 1 Abstract In this paper the simulation of a basic packet routing problem using

More information

Lab 5 - Risk Analysis, Robustness, and Power

Lab 5 - Risk Analysis, Robustness, and Power Type equation here.biology 458 Biometry Lab 5 - Risk Analysis, Robustness, and Power I. Risk Analysis The process of statistical hypothesis testing involves estimating the probability of making errors

More information

The Bootstrap and Jackknife

The Bootstrap and Jackknife The Bootstrap and Jackknife Summer 2017 Summer Institutes 249 Bootstrap & Jackknife Motivation In scientific research Interest often focuses upon the estimation of some unknown parameter, θ. The parameter

More information

Snell s Law. Introduction

Snell s Law. Introduction Snell s Law Introduction According to Snell s Law [1] when light is incident on an interface separating two media, as depicted in Figure 1, the angles of incidence and refraction, θ 1 and θ 2, are related,

More information

ACCURACY AND EFFICIENCY OF MONTE CARLO METHOD. Julius Goodman. Bechtel Power Corporation E. Imperial Hwy. Norwalk, CA 90650, U.S.A.

ACCURACY AND EFFICIENCY OF MONTE CARLO METHOD. Julius Goodman. Bechtel Power Corporation E. Imperial Hwy. Norwalk, CA 90650, U.S.A. - 430 - ACCURACY AND EFFICIENCY OF MONTE CARLO METHOD Julius Goodman Bechtel Power Corporation 12400 E. Imperial Hwy. Norwalk, CA 90650, U.S.A. ABSTRACT The accuracy of Monte Carlo method of simulating

More information

A New Statistical Procedure for Validation of Simulation and Stochastic Models

A New Statistical Procedure for Validation of Simulation and Stochastic Models Syracuse University SURFACE Electrical Engineering and Computer Science L.C. Smith College of Engineering and Computer Science 11-18-2010 A New Statistical Procedure for Validation of Simulation and Stochastic

More information

Equivalence Tests for Two Means in a 2x2 Cross-Over Design using Differences

Equivalence Tests for Two Means in a 2x2 Cross-Over Design using Differences Chapter 520 Equivalence Tests for Two Means in a 2x2 Cross-Over Design using Differences Introduction This procedure calculates power and sample size of statistical tests of equivalence of the means of

More information

Unit 1 Review of BIOSTATS 540 Practice Problems SOLUTIONS - Stata Users

Unit 1 Review of BIOSTATS 540 Practice Problems SOLUTIONS - Stata Users BIOSTATS 640 Spring 2018 Review of Introductory Biostatistics STATA solutions Page 1 of 13 Key Comments begin with an * Commands are in bold black I edited the output so that it appears here in blue Unit

More information

Quality Improvement Tools

Quality Improvement Tools CHAPTER SIX SUPPLEMENT Quality Improvement Tools McGraw-Hill/Irwin Copyright 2011 by the McGraw-Hill Companies, Inc. All rights reserved. Learning Objectives 1. Apply quality management tools for problem

More information

Lecture 14. December 19, Department of Biostatistics Johns Hopkins Bloomberg School of Public Health Johns Hopkins University.

Lecture 14. December 19, Department of Biostatistics Johns Hopkins Bloomberg School of Public Health Johns Hopkins University. geometric Lecture 14 Department of Biostatistics Johns Hopkins Bloomberg School of Public Health Johns Hopkins University December 19, 2007 geometric 1 2 3 4 geometric 5 6 7 geometric 1 Review about logs

More information

Determination of the Process Mean and Production Run Using Simulation *

Determination of the Process Mean and Production Run Using Simulation * International Conference on Industrial Engineering and Systems Management IESM 27 May 3 - June 2, 27 BEIJING - CHINA Determination of the Process Mean and Production Run Using Simulation * A. M. Al-Ahmari

More information

THIS IS NOT REPRESNTATIVE OF CURRENT CLASS MATERIAL. STOR 455 Midterm 1 September 28, 2010

THIS IS NOT REPRESNTATIVE OF CURRENT CLASS MATERIAL. STOR 455 Midterm 1 September 28, 2010 THIS IS NOT REPRESNTATIVE OF CURRENT CLASS MATERIAL STOR 455 Midterm September 8, INSTRUCTIONS: BOTH THE EXAM AND THE BUBBLE SHEET WILL BE COLLECTED. YOU MUST PRINT YOUR NAME AND SIGN THE HONOR PLEDGE

More information

MBI REU Matlab Tutorial

MBI REU Matlab Tutorial MBI REU Matlab Tutorial Lecturer: Reginald L. McGee II, Ph.D. June 8, 2017 MATLAB MATrix LABoratory MATLAB is a tool for numerical computation and visualization which allows Real & Complex Arithmetics

More information

Process Capability Analysis in Case Study of Specimens for Rice Polished Cylinder

Process Capability Analysis in Case Study of Specimens for Rice Polished Cylinder International Science Index Vol: 8 No: Part V Process Capability Analysis in Case Study of Specimens for ice Polished Cylinder T. Boonkang, S. Bangphan, P. Bangphan, T. Pothom Abstract Process capability

More information

1. Estimation equations for strip transect sampling, using notation consistent with that used to

1. Estimation equations for strip transect sampling, using notation consistent with that used to Web-based Supplementary Materials for Line Transect Methods for Plant Surveys by S.T. Buckland, D.L. Borchers, A. Johnston, P.A. Henrys and T.A. Marques Web Appendix A. Introduction In this on-line appendix,

More information

Practical Design of Experiments: Considerations for Iterative Developmental Testing

Practical Design of Experiments: Considerations for Iterative Developmental Testing Practical Design of Experiments: Considerations for Iterative Developmental Testing Best Practice Authored by: Michael Harman 29 January 2018 The goal of the STAT COE is to assist in developing rigorous,

More information

IQR = number. summary: largest. = 2. Upper half: Q3 =

IQR = number. summary: largest. = 2. Upper half: Q3 = Step by step box plot Height in centimeters of players on the 003 Women s Worldd Cup soccer team. 157 1611 163 163 164 165 165 165 168 168 168 170 170 170 171 173 173 175 180 180 Determine the 5 number

More information

Developing Effect Sizes for Non-Normal Data in Two-Sample Comparison Studies

Developing Effect Sizes for Non-Normal Data in Two-Sample Comparison Studies Developing Effect Sizes for Non-Normal Data in Two-Sample Comparison Studies with an Application in E-commerce Durham University Apr 13, 2010 Outline 1 Introduction Effect Size, Complementory for Hypothesis

More information

Interpolation by Spline Functions

Interpolation by Spline Functions Interpolation by Spline Functions Com S 477/577 Sep 0 007 High-degree polynomials tend to have large oscillations which are not the characteristics of the original data. To yield smooth interpolating curves

More information

One Factor Experiments

One Factor Experiments One Factor Experiments 20-1 Overview Computation of Effects Estimating Experimental Errors Allocation of Variation ANOVA Table and F-Test Visual Diagnostic Tests Confidence Intervals For Effects Unequal

More information

Department of Industrial Engineering. Chap. 8: Process Capability Presented by Dr. Eng. Abed Schokry

Department of Industrial Engineering. Chap. 8: Process Capability Presented by Dr. Eng. Abed Schokry Department of Industrial Engineering Chap. 8: Process Capability Presented by Dr. Eng. Abed Schokry Learning Outcomes: After careful study of this chapter, you should be able to do the following: Investigate

More information

Multivariate Normal Random Numbers

Multivariate Normal Random Numbers Multivariate Normal Random Numbers Revised: 10/11/2017 Summary... 1 Data Input... 3 Analysis Options... 4 Analysis Summary... 5 Matrix Plot... 6 Save Results... 8 Calculations... 9 Summary This procedure

More information

ITT Technical Institute. NT2740 Advanced Networking Devices Onsite and Online Course SYLLABUS

ITT Technical Institute. NT2740 Advanced Networking Devices Onsite and Online Course SYLLABUS ITT Technical Institute NT2740 Advanced Networking Devices Onsite and Online Course SYLLABUS Credit hours: 3 Contact/Instructional hours: 40 (20 Theory Hours, 20 Lab Hours) Prerequisite(s) and/or Corequisite(s):

More information

Philip E. Plantz. Application Note. SL-AN-08 Revision B. Provided By: Microtrac, Inc. Particle Size Measuring Instrumentation

Philip E. Plantz. Application Note. SL-AN-08 Revision B. Provided By: Microtrac, Inc. Particle Size Measuring Instrumentation A Conceptual, Non Mathematical Explanation on the Use of Refractive Index in Laser Particle Size Measurement: (Understanding the concept of refractive index and Mie Scattering in Microtrac Instruments

More information

If your problem requires a series of related binomial coefficients, a good idea is to use recurrence relations, for example. = n k

If your problem requires a series of related binomial coefficients, a good idea is to use recurrence relations, for example. = n k 6. Incomplete Gamma Function 9 If your problem requires a series of related binomial coefficients, a good idea is to use recurrence relations, for example ( n +1 = n +1 ( ( ( n n n = + k ( n k +1 n k +1

More information

Statistical Methods. Instructor: Lingsong Zhang. Any questions, ask me during the office hour, or me, I will answer promptly.

Statistical Methods. Instructor: Lingsong Zhang. Any questions, ask me during the office hour, or  me, I will answer promptly. Statistical Methods Instructor: Lingsong Zhang 1 Issues before Class Statistical Methods Lingsong Zhang Office: Math 544 Email: lingsong@purdue.edu Phone: 765-494-7913 Office Hour: Monday 1:00 pm - 2:00

More information

Truss structural configuration optimization using the linear extended interior penalty function method

Truss structural configuration optimization using the linear extended interior penalty function method ANZIAM J. 46 (E) pp.c1311 C1326, 2006 C1311 Truss structural configuration optimization using the linear extended interior penalty function method Wahyu Kuntjoro Jamaluddin Mahmud (Received 25 October

More information

Approximations and Errors

Approximations and Errors The Islamic University of Gaza Faculty of Engineering Civil Engineering Department Numerical Analysis ECIV 3306 Chapter 3 Approximations and Errors Associate Prof. Mazen Abualtayef Civil Engineering Department,

More information

Data Analysis & Probability

Data Analysis & Probability Unit 5 Probability Distributions Name: Date: Hour: Section 7.2: The Standard Normal Distribution (Area under the curve) Notes By the end of this lesson, you will be able to Find the area under the standard

More information

MHPE 494: Data Analysis. Welcome! The Analytic Process

MHPE 494: Data Analysis. Welcome! The Analytic Process MHPE 494: Data Analysis Alan Schwartz, PhD Department of Medical Education Memoona Hasnain,, MD, PhD, MHPE Department of Family Medicine College of Medicine University of Illinois at Chicago Welcome! Your

More information

Multivariate distribution theory by holonomic gradient method. Akimichi Takemura, Shiga University June 28, 2016, IMS-APRM 2016

Multivariate distribution theory by holonomic gradient method. Akimichi Takemura, Shiga University June 28, 2016, IMS-APRM 2016 Multivariate distribution theory by holonomic gradient method Akimichi Takemura, Shiga University June 28, 2016, IMS-APRM 2016 Outline 1. An advertisement and overview of my research on holonomic gradient

More information

Capability Calculations: Are AIAG SPC Appendix F Conclusions Wrong?

Capability Calculations: Are AIAG SPC Appendix F Conclusions Wrong? WHITE PAPER Capability Calculations: Are AIAG SPC Appendix F Conclusions Wrong? Bob Doering CorrectSPC Page 0 Appendix 7 of the AIAG SPC book contains sample data set and calculations for capability. They

More information

Chapters 5-6: Statistical Inference Methods

Chapters 5-6: Statistical Inference Methods Chapters 5-6: Statistical Inference Methods Chapter 5: Estimation (of population parameters) Ex. Based on GSS data, we re 95% confident that the population mean of the variable LONELY (no. of days in past

More information

Bootstrap confidence intervals Class 24, Jeremy Orloff and Jonathan Bloom

Bootstrap confidence intervals Class 24, Jeremy Orloff and Jonathan Bloom 1 Learning Goals Bootstrap confidence intervals Class 24, 18.05 Jeremy Orloff and Jonathan Bloom 1. Be able to construct and sample from the empirical distribution of data. 2. Be able to explain the bootstrap

More information

UPPER CONFIDENCE LIMIT FOR NON-NORMAL INCAPABILITY INDEX: A CASE STUDY

UPPER CONFIDENCE LIMIT FOR NON-NORMAL INCAPABILITY INDEX: A CASE STUDY مو تمر الا زهر الهندسي الدولي العاشر AL-AZHAR ENGINEERING TENTH INTERNATIONAL ONFERENE December 4-6, 008 ode: M 01 UPPER ONFIDENE LIMIT FOR NON-NORMAL INAPABILITY INDEX: A ASE STUDY A. Rotondo 1 and Amir

More information

STATGRAPHICS PLUS for WINDOWS

STATGRAPHICS PLUS for WINDOWS TUTORIALS FOR Quality Control Analyses STATGRAPHICS PLUS for WINDOWS SEPTEMBER 1999 MANUGISTICS, INC 2115 East Jefferson Street Rockville, Maryland 20852 Introduction This manual contains tutorials for

More information

Statistical Process Control: A Case-Study on Haleeb Foods Ltd., Lahore

Statistical Process Control: A Case-Study on Haleeb Foods Ltd., Lahore 11 ISSN 1684 8403 Journal of Statistics Vol: 12, No.1 (2005) Statistical Process Control: A Case-Study on Haleeb Foods Ltd., Lahore Sarwat Zahara Khan *, Muhammad Khalid Pervaiz * and Mueen-ud-Din Azad

More information

Descriptive Statistics, Standard Deviation and Standard Error

Descriptive Statistics, Standard Deviation and Standard Error AP Biology Calculations: Descriptive Statistics, Standard Deviation and Standard Error SBI4UP The Scientific Method & Experimental Design Scientific method is used to explore observations and answer questions.

More information

Using Modified Euler Method (MEM) for the Solution of some First Order Differential Equations with Initial Value Problems (IVPs).

Using Modified Euler Method (MEM) for the Solution of some First Order Differential Equations with Initial Value Problems (IVPs). Using Modified Euler Method (MEM) for the Solution of some First Order Differential Equations with Initial Value Problems (IVPs). D.I. Lanlege, Ph.D. * ; U.M. Garba, B.Sc.; and A. Aluebho, B.Sc. Department

More information

A web app for guided and interactive generation of multimarker panels (www.combiroc.eu)

A web app for guided and interactive generation of multimarker panels (www.combiroc.eu) COMBIROC S TUTORIAL. A web app for guided and interactive generation of multimarker panels (www.combiroc.eu) Overview of the CombiROC workflow. CombiROC delivers a simple workflow to help researchers in

More information

Lecture 8: Jointly distributed random variables

Lecture 8: Jointly distributed random variables Lecture : Jointly distributed random variables Random Vectors and Joint Probability Distributions Definition: Random Vector. An n-dimensional random vector, denoted as Z = (Z, Z,, Z n ), is a function

More information

Primary Vertex Reconstruction at LHCb

Primary Vertex Reconstruction at LHCb LHCb-PUB-214-44 October 21, 214 Primary Vertex Reconstruction at LHCb M. Kucharczyk 1,2, P. Morawski 3, M. Witek 1. 1 Henryk Niewodniczanski Institute of Nuclear Physics PAN, Krakow, Poland 2 Sezione INFN

More information

3D Nearest-Nodes Finite Element Method for Solid Continuum Analysis

3D Nearest-Nodes Finite Element Method for Solid Continuum Analysis Adv. Theor. Appl. Mech., Vol. 1, 2008, no. 3, 131-139 3D Nearest-Nodes Finite Element Method for Solid Continuum Analysis Yunhua Luo Department of Mechanical & Manufacturing Engineering, University of

More information

Three Different Algorithms for Generating Uniformly Distributed Random Points on the N-Sphere

Three Different Algorithms for Generating Uniformly Distributed Random Points on the N-Sphere Three Different Algorithms for Generating Uniformly Distributed Random Points on the N-Sphere Jan Poland Oct 4, 000 Abstract We present and compare three different approaches to generate random points

More information

Introduction to the Practice of Statistics Fifth Edition Moore, McCabe

Introduction to the Practice of Statistics Fifth Edition Moore, McCabe Introduction to the Practice of Statistics Fifth Edition Moore, McCabe Section 1.3 Homework Answers Assignment 5 1.80 If you ask a computer to generate "random numbers between 0 and 1, you uniform will

More information

PARAMETRIC MODEL SELECTION TECHNIQUES

PARAMETRIC MODEL SELECTION TECHNIQUES PARAMETRIC MODEL SELECTION TECHNIQUES GARY L. BECK, STEVE FROM, PH.D. Abstract. Many parametric statistical models are available for modelling lifetime data. Given a data set of lifetimes, which may or

More information

PASS Sample Size Software. Randomization Lists

PASS Sample Size Software. Randomization Lists Chapter 880 Introduction This module is used to create a randomization list for assigning subjects to one of up to eight groups or treatments. Six randomization algorithms are available. Four of the algorithms

More information

A New Efficient and Useful Robust Optimization Approach Design for Multi-Objective Six Sigma

A New Efficient and Useful Robust Optimization Approach Design for Multi-Objective Six Sigma A New Efficient and Useful Robust Optimization Approach Design for Multi-Objective Six Sigma Koji Shimoyama Department of Aeronautics and Astronautics University of Tokyo 3-1-1 Yoshinodai Sagamihara, Kanagawa,

More information

Continuous Improvement Toolkit. Normal Distribution. Continuous Improvement Toolkit.

Continuous Improvement Toolkit. Normal Distribution. Continuous Improvement Toolkit. Continuous Improvement Toolkit Normal Distribution The Continuous Improvement Map Managing Risk FMEA Understanding Performance** Check Sheets Data Collection PDPC RAID Log* Risk Analysis* Benchmarking***

More information

Unit 5: Estimating with Confidence

Unit 5: Estimating with Confidence Unit 5: Estimating with Confidence Section 8.3 The Practice of Statistics, 4 th edition For AP* STARNES, YATES, MOORE Unit 5 Estimating with Confidence 8.1 8.2 8.3 Confidence Intervals: The Basics Estimating

More information

Xbar/R Chart for x1-x3

Xbar/R Chart for x1-x3 Chapter 6 Selected roblem Solutios Sectio 6-5 6- a) X-bar ad Rage - Iitial Study Chartig roblem 6- X-bar Rage ----- ----- UCL:. sigma 7.4 UCL:. sigma 5.79 Ceterlie 5.9 Ceterlie.5 LCL: -. sigma.79 LCL:

More information

Mean Tests & X 2 Parametric vs Nonparametric Errors Selection of a Statistical Test SW242

Mean Tests & X 2 Parametric vs Nonparametric Errors Selection of a Statistical Test SW242 Mean Tests & X 2 Parametric vs Nonparametric Errors Selection of a Statistical Test SW242 Creation & Description of a Data Set * 4 Levels of Measurement * Nominal, ordinal, interval, ratio * Variable Types

More information

Introduction to hypothesis testing

Introduction to hypothesis testing Introduction to hypothesis testing Mark Johnson Macquarie University Sydney, Australia February 27, 2017 1 / 38 Outline Introduction Hypothesis tests and confidence intervals Classical hypothesis tests

More information

APPROACHES TO THE PROCESS CAPABILITY ANALYSIS IN THE CASE OF NON- NORMALLY DISTRIBUTED PRODUCT QUALITY CHARACTERISTIC

APPROACHES TO THE PROCESS CAPABILITY ANALYSIS IN THE CASE OF NON- NORMALLY DISTRIBUTED PRODUCT QUALITY CHARACTERISTIC APPROACHES TO THE PROCESS CAPABILITY ANALYSIS IN THE CASE OF NON- NORMALLY DISTRIBUTED PRODUCT QUALITY CHARACTERISTIC Jiří PLURA, Milan ZEMEK, Pavel KLAPUT VŠB-Technical University of Ostrava, Faculty

More information

Philip E. Plantz. Application Note. SL-AN-08 Revision C. Provided By: Microtrac, Inc. Particle Size Measuring Instrumentation

Philip E. Plantz. Application Note. SL-AN-08 Revision C. Provided By: Microtrac, Inc. Particle Size Measuring Instrumentation A Conceptual, Non-Mathematical Explanation on the Use of Refractive Index in Laser Particle Size Measurement (Understanding the concept of refractive index and Mie Scattering in Microtrac Instruments and

More information

Harris: Quantitative Chemical Analysis, Eight Edition CHAPTER 05: QUALITY ASSURANCE AND CALIBRATION METHODS

Harris: Quantitative Chemical Analysis, Eight Edition CHAPTER 05: QUALITY ASSURANCE AND CALIBRATION METHODS Harris: Quantitative Chemical Analysis, Eight Edition CHAPTER 05: QUALITY ASSURANCE AND CALIBRATION METHODS 5-0. International Measurement Evaluation Program Sample: Pb in river water (blind sample) :

More information

GRAPmCAL EVALUATION OF PRODUCT CHARACTERISTICS

GRAPmCAL EVALUATION OF PRODUCT CHARACTERISTICS GRAPmCAL EVALUATION OF PRODUCT CHARACTERISTICS Melissa A. Durfee, Wyman-Gordon Company Abstract Utilizing SAS/GRAPH, many aspects of product characteristics may be displayed and analyzed. Besid.es using

More information

Six Sigma Green Belt Part 5

Six Sigma Green Belt Part 5 Six Sigma Green Belt Part 5 Process Capability 2013 IIE and Aft Systems, Inc. 5-1 Process Capability Is the measured, inherent reproducibility of the product turned out by the process. It can be quantified

More information

IJQRM 21,6. Received March 2004 Revised March 2004 Accepted March 2004

IJQRM 21,6. Received March 2004 Revised March 2004 Accepted March 2004 The Emerald Research Register for this journal is available at wwwem eraldinsightcom/res earchregister The current issue and full text archive of this journal is available at wwwem eraldinsight com/0265-671x

More information

Physics 736. Experimental Methods in Nuclear-, Particle-, and Astrophysics. - Statistical Methods -

Physics 736. Experimental Methods in Nuclear-, Particle-, and Astrophysics. - Statistical Methods - Physics 736 Experimental Methods in Nuclear-, Particle-, and Astrophysics - Statistical Methods - Karsten Heeger heeger@wisc.edu Course Schedule and Reading course website http://neutrino.physics.wisc.edu/teaching/phys736/

More information

COMPARISON OF THE ASSUMED SHAPE AND FINITE ELEMENT METHODS FOR A VIBRATING BRACE-PLATE SYSTEM

COMPARISON OF THE ASSUMED SHAPE AND FINITE ELEMENT METHODS FOR A VIBRATING BRACE-PLATE SYSTEM Proceedings of The Canadian Society for Mechanical Engineering Forum CSME FORUM June 7-9,, Victoria, British Columbia, Canada COMPARISON OF THE ASSUMED SHAPE AND FINITE ELEMENT METHODS FOR A VIBRATING

More information

SD 372 Pattern Recognition

SD 372 Pattern Recognition SD 372 Pattern Recognition Lab 2: Model Estimation and Discriminant Functions 1 Purpose This lab examines the areas of statistical model estimation and classifier aggregation. Model estimation will be

More information

Test Oracles and Randomness

Test Oracles and Randomness Test Oracles and Randomness UNIVERSITÄ T ULM DOCENDO CURANDO SCIENDO Ralph Guderlei and Johannes Mayer rjg@mathematik.uni-ulm.de, jmayer@mathematik.uni-ulm.de University of Ulm, Department of Stochastics

More information

Supplemental Material for Spatially-Balanced Designs that Incorporate Legacy Sites

Supplemental Material for Spatially-Balanced Designs that Incorporate Legacy Sites Supplemental Material for Spatially-Balanced Designs that Incorporate Legacy Sites S.D. Foster, G.R. Hosack, E. Lawrence, R. Przeslawski, P. Hedge, M.J. Caley, N.S. Barrett, A. Williams, J. Li, T. Lynch,

More information

VARIANCE REDUCTION TECHNIQUES IN MONTE CARLO SIMULATIONS K. Ming Leung

VARIANCE REDUCTION TECHNIQUES IN MONTE CARLO SIMULATIONS K. Ming Leung POLYTECHNIC UNIVERSITY Department of Computer and Information Science VARIANCE REDUCTION TECHNIQUES IN MONTE CARLO SIMULATIONS K. Ming Leung Abstract: Techniques for reducing the variance in Monte Carlo

More information

SAR SURFACE ICE COVER DISCRIMINATION USING DISTRIBUTION MATCHING

SAR SURFACE ICE COVER DISCRIMINATION USING DISTRIBUTION MATCHING SAR SURFACE ICE COVER DISCRIMINATION USING DISTRIBUTION MATCHING Rashpal S. Gill Danish Meteorological Institute, Ice charting and Remote Sensing Division, Lyngbyvej 100, DK 2100, Copenhagen Ø, Denmark.Tel.

More information

IQC monitoring in laboratory networks

IQC monitoring in laboratory networks IQC for Networked Analysers Background and instructions for use IQC monitoring in laboratory networks Modern Laboratories continue to produce large quantities of internal quality control data (IQC) despite

More information

AME 20213: Fundamentals of Measurements and Data Analysis. LEX-3: Wind Tunnel Pressure Measurement Uncertainty 9 October 2006

AME 20213: Fundamentals of Measurements and Data Analysis. LEX-3: Wind Tunnel Pressure Measurement Uncertainty 9 October 2006 AME 20213: Fundamentals of Measurements and Data Analysis LEX-3: Wind Tunnel Pressure Measurement Uncertainty 9 October 2006 TA: John Schmitz Office: B028 Hessert Phone: 1-2317 E-mail : jschmitz@nd.edu

More information

MODIFIED VERSION OF: An introduction to Matlab for dynamic modeling ***PART 3 ***

MODIFIED VERSION OF: An introduction to Matlab for dynamic modeling ***PART 3 *** MODIFIED VERSION OF: An introduction to Matlab for dynamic modeling ***PART 3 *** Stephen P. Ellner 1 and John Guckenheimer 2 1 Department of Ecology and Evolutionary Biology, and 2 Department of Mathematics

More information

α-pareto optimal solutions for fuzzy multiple objective optimization problems using MATLAB

α-pareto optimal solutions for fuzzy multiple objective optimization problems using MATLAB Advances in Modelling and Analysis C Vol. 73, No., June, 18, pp. 53-59 Journal homepage:http://iieta.org/journals/ama/ama_c α-pareto optimal solutions for fuzzy multiple objective optimization problems

More information

Exact Calculation of Beta Inequalities

Exact Calculation of Beta Inequalities Exact Calculation of Beta Inequalities John Cook Department of Biostatistics, Box 447 The University of Texas, M. D. Anderson Cancer Center 1515 Holcombe Blvd., Houston, Texas 77030, USA cook@mdanderson.org

More information