Statistical Process Control Methods from the Viewpoint of Industrial Application
|
|
- Sheryl Morgan
- 5 years ago
- Views:
Transcription
1 c Heldermann Verlag Economic Quality Control ISSN Vol 16 (2001), No. 1, Statistical Process Control Methods from the Viewpoint of Industrial Application Constantin Anghel Abstract: Statistical process control is a major part of industrial statistics and consists not only of control charting, but also of capability analysis, of design of experiments and other statistical techniques. This paper reviews and comments some available techniques for the univariate as well as for the multivariate case from the viewpoint of industrial application. 1 Distribution of a Quality Characteristic Many statistical process control activities assume that the quality characteristic under study exhibits a stable performance or in other words that the process is stable and capable. The stability of a process can be defined as the stability of the underlying probability distribution over time and very often this can be described as the stability of the distribution parameters over time. Only if the stability assumption is met by the process, the calculation of capability indices is meaningful and may be used in practice for evaluating process performance. Mathematically, we can, of course, calculate always capability indices, but for an unstable process these indices have no real significance, because there are not identified assignable causes in the process. Thus, a correct identification of the type of probability distribution is not sufficient without the assurance that it is stable over time. In the case that the process is not stable, the probability distribution of the quality characteristic may vary from time-point to time-point. In such a case one could use a mixture of probability distributions as model, which, therefore, can be looked upon as an indication for an unstable process. Each quality characteristic has a proper nature defined by physical and technical conditions. For instance a meaningful flatness-characteristic is defined as a parameter naturally bounded by null. In such a case the normal distribution cannot describe reality sufficiently well and a different model has to be selected. One of many possibilities would be to select the folded normal distribution [1] and [8]. Of course, any model should always be justified by experiments, i.e. by samples. If a sample leads to reject the hypothesis of the assumed model, then the situation must be analyzed again, with respect to various features as the resolution of the measurement gauge, mixture appearance, etc. Table 1 gives some recommendations for the selection of an appropriate model for some frequently relevant quality characteristics (see also [1]). A decision for a specific model to be used to describe the random variations exhibited by a quality characteristic, should be based on any available information about the physical and
2 50 Constantin Anghel technical conditions on the one hand and on a sufficiently large sample on the other hand, collected under ideal conditions (one operator, one supplier, optimal resolution of the gauge). Subsequent samples of the quality characteristic should be utilized for verifying the hypothesis that the initial distribution has not changed. In case of a rejection of the null hypothesis, an explanation has to be found by an additional analysis. Table 1: Recommended Models for Various Quality Characteristics Characteristic Model Characteristic Model geometric dimensions: straightness folded normal parallelism folded normal evenness folded normal rectangularity folded normal round form folded normal inclination (angularity) folded normal cylinder form folded normal position eccentric Rayleigh line form folded normal coaxiality (concentricity) eccentric Rayleigh surface form folded normal symmetry folded normal circularity flatness folded normal a) form folded normal b) position eccentric Rayleigh life length Weibull, Hjorth resistance normal voltage normal capacitance normal pressure normal viscosity normal These models should be looked upon as a first and preliminary attempt which, in any case, must be justified or discarded by a subsequent analysis of the situation at hand. 2 Data Collection and the Resolution of the Gauge Any data collection is the result of a measurement process and, therefore, the gauge has a high importance beginning with its resolution and its capability. A resolution below 2% of the tolerance is considered as consistent with industrial practice. Above this level the measurement values are assigned to classes by the measurement procedure and the random character of the sample will partly be lost. If we take samples of size two, then the expectation of range is times the process standard deviation. This means that the expected distance between any two randomly selected measurements is about 1.1sσ, where σ denotes the process standard deviation. Whenever the difference between two successive measurements is less than the resolution of the gauge, the difference is set to zero. Thus, a larger resolution yields a larger number of zeros, and consequently results in a too small value of the average range (process variability) leading to action limits for an X-chart as well as for an R-chart, which are too tight (Figure 1). The consequences are too many false alarms, and often the wrong decision that the process is not stable [13].
3 Statistical Process Control Methods 51 1 Example Consider a quality characteristic with tolerance interval [L, U] =[28.7, 31.3], which is evaluated by a gauge with resolution r =0.01 ( 4% of the length of the tolerance interval). 50 samples of size n = 2 are taken the results of which are represented in Table 1 and Figures 1. Table 2: Observed values: 50 samples of size n =2. observed value x frequency Figure 1: X-chart and R-chart for the 50 samples of size n =2. The uncertainty about the measurement result determines the gauge capability, which often is specified by four times the standard deviation 4σ. Clearly, the gauge capability affects the determination of the actual value of the process capability given by C p. If the ratio of the measurement uncertainty (4σ) andthetolerance U L, withl lower tolerance limit and U upper tolerance limit, is 25%, then a measured process capability (with measurement variability added to process variability) of 1.33 means for the actual index (without measurement variability) a value of 1.55 and a measured value of 1.67 corresponds to a real value of 2.2 (Figure 2).
4 52 Constantin Anghel Figure 2: The effect of measurement variability (uncertainty) on the determination of process capability C p. Conversely, for a given actual process capability (without measurement variability) of 1.33, the capability with added measurement variability ( 4σ = 25%) reaches only 1.2 and T an actual value of 1.67 is decreased to 1.41 (Figure 3). Figure 3: The effect of measurement variability (uncertainty) on the verification of a given process capability C p. Thus, information about gauge capability is very important to avoid false decisions with respect to process capability and with respect to releasing false alarms.
5 Statistical Process Control Methods 53 3 Capability Indices As mentioned above, it makes sense to define capability indices for a stable process. Of course, it is possible to use the definitions and to calculate these indices, even if the distribution of the relevant quality characteristics change in time. However, only a stable process assures the predictions based on the stated capability index and, therefore, using the knowledge of what happened today for predicting what we expect tomorrow makes sense only if the situation does not change. From the viewpoint of application, it is desirable to define capability indices being independent of the special type of distribution allowing the comparison of two or more processes with different underlying distributions. The most widely used capability indices for the normal distribution are the following: C p = U L (1) 6σ C pk = 1 ( U µ 3 min, µ L ) (2) σ σ where µ is the process mean and σ 2 the process variance. Introducing the probability p U of exceeding the upper specification limit U and p L of not reaching the lower specification limit L, another representation of the capability indices is obtained. ( ) ( ) U µ L µ p U =1 Φ, p L =Φ (3) σ σ and therefore C p = 1 ( ) Φ 1 (1 p U ) Φ 1 (p L ) (4) 6 C pk = 1 ) (Φ 3 min 1 (1 p U ), Φ 1 (p L ) = 1 ) (Φ 3 min 1 (1 p U ), Φ 1 (1 p L ) = 1 3 Φ 1 (1 p) (5) where Φ 1 denotes the inverse distribution function of the standard normal distribution and p =max(p U,p L ). The representation (4) and (5) also suggests how to define process capability indices for the general case given by an arbitrary distribution function F (x). In the general case the probability of exceeding the upper specification limit denoted by p (F ) U and the probability of not exceeding the lower specification limit denoted by p (F ) L are given by the following expressions: p (F ) U =1 F (U), p(f ) L = F (L) (6) Replacing p U and p L in (4) and (5) by the expressions (6) yields the following formulas for C p and C pk in the general case.
6 54 Constantin Anghel C p = 1 ( ) Φ 1 (1 p (F ) U 6 ) Φ 1 (p (F ) L ) (7) C pk = 1 ( ) 3 min Φ 1 (1 p (F ) U ), Φ 1 (p (F ) L ) = 1 3 Φ 1 (1 p) (8) We will call (8) the classical definition of C pk, as it defines process capability by means of the maximum nonconformance probability max(p U,p L ) on the one or the other side of the specification interval. In the relevant literature (e.g. [5]) a different proposal for defining C p and C pk for the general case is made. The definition is based on the following representation of the denominator of (1) which is valid under the normal model. 6σ =(3σ + µ) ( 3σ + µ) =q q (9) where q = µ+3σ is the quantile and q = µ 3σ is the quantile of the normal distribution N(µ, σ 2 ). Replacing the denominators in (1) and (2) by the expressions obtained in (9), we arrive at the following formula for the general case. Cp U L = q (F ) q (F ) ( Cpk U µ = min q (F ) µ, µ L µ q (F ) ) where q (F ) and q (F ) , respectively, are quantiles of an arbitrary distribution function F (x). Clearly in the normal case the two capability indices (11) and (8) are equivalent. However, this is not true in the general case as illustrated by the following example. 2 Example Consider a quality characteristic X with specification limits L = 1.0 and U = Moreover, let X follow a log-normal distribution [3] given by its density function f X (x): 1 f X (x) = σ(x θ) 1 2( ln(x θ) µ 2π e σ ) 2 (12) with θ = 7.5, µ =2.56, σ =0.27. The expectation of X is given by σ2 µ+ E[X] =θ + e 2 (13) and utilizing the fact that ln(x θ) N(µ, σ 2 ), the p-quantile of the log-normal distribution q p (F ) may be determined by the following formula. q p (F ) = θ + e qpσ+µ (14) where q p is the p-quantile of the standardized normal distribution. Thus, we obtain: (10) (11)
7 Statistical Process Control Methods 55 Table 3: Distribution characteristic of X characteristic actual value process mean E[X] = quantile q = quantile q = With the values given in Table 3, we obtain the capability index C pk U E[X] q E[X] =0.592, E[X] L =0.642 (15) E[X] q and thus: Cpk =0.592 (16) where (16) refers to process quality with respect to exceeding the upper specification limit U. Next, the classical definition of C pk, as given by (8), is used for calculating the process capability. To this end the probability p U of exceeding the upper specification limit and the probability p L of not reaching the lower specification limit have to be determined. ( ) ln(u θ) µ p U = 1 F X (U) =1 Φ σ = 1 Φ(2.08) = (17) ( ) ln(l θ) µ p L = F X (L) =Φ σ = Φ( 1.56) = (18) The probability of not reaching the lower specification limit is more than three times the probability of exceeding the upper specification limit. Thus, process capability with respect to meeting the upper specification is much better than process capability with respect to meeting the lower specification and this fact should be reflected by any meaningful capability index. Unfortunately, as shown by (15), Cpk violates this requirement. For the classical capability index C pk, we obtain: Φ 1 (1 p U ) = Φ 1 ( ) = 2.08 (19) Φ 1 (1 p L ) = Φ 1 ( ) = 1.56 (20) and therefore C pk =0.520 (21) which refers to the process capability of meeting the lower specification limit L.
8 56 Constantin Anghel It is generally accepted in industry that the requirement for the classical C pk in the normal case should be C pk 1.33 with the meaning that on an average there are not more than 31.7 ppm of nonconforming items produced on either side of the specification interval. Taking Cpk with C pk 1.33 in the case of an exponentially distributed quality characteristic may lead to an out of specification proportion of 212 ppm on one side of the specification interval, which is about seven times the tolerated number. Therefore, Cpk should not be used in the non-normal case. 4 Median Versus Expectation For non-normal distributions it may be more appropriate to use the median rather than the expectation for controlling the process location, since the expectation of a non-normal distribution does not necessarily correspond to a quantile of fixed order as in the case of the normal distribution (50%). Figure 5 shows the action lines and the operation of 500 simulated samples of size n = 3 of an X-charts for an exponential distribution with =0.3 and action limit (99%) determined according to a Shewhart-chart and determined based on the distribution of the sample mean. For either chart a number of false alarms occur which is considerably larger for the Shewhart chart. Figure 4: Performance of X-charts for controlling the process mean in the exponential case. Figure 5 shows the operation of a median chart for the same simulated samples. As can be seen, there is no false alarm for the median chart.
9 Statistical Process Control Methods 57 Figure 5: Performance of median chart for controlling the process mean in the exponential case. As can be seen from Figure 5, there is no false alarm for the median chart. Schneider et al.[11] pointed out that for the traditional three-sigma control limit (Shewhart) the probability for a false alarm is 0.27% for the normal distribution but 1.8% for an exponential distribution with µ = 1 (Figure 6). Figure 6: False alarm probability using three-sigma control limit in the normal and the exponential case. During process control with control charts it is not necessary to check the distribution and to calculate capability indices as long as the control charts do not indicate a change. The distribution as well as the capability indices as determined at the latest process capability study of the output characteristic are valid until an alarm is released and the process is analyzed and a new process capability study is carried out.
10 58 Constantin Anghel 5 Multivariate Process Analysis In industrial practice product quality is generally determined not by only one characteristic but by a vector of several characteristics which may be correlated. For such situations it is generally not appropriate to decide about each variable, or to use quality charts, or to calculate capability indices for each characteristic separately by means of the marginal distributions and thus neglecting the correlations. Therefore, multivariate techniques are necessary taking into account the correlations between quality characteristics for calculating quality charts, tolerances (see e.g. [6]), and capability indices. Figure 8 and 9 illustrate independent control of two variables with respect to location and variability compared with a joint control (Figure 9) using quality charts for the multivariate mean variation. The better performance of the joint charts is evident. Figure 7: Independent control of locations of the two quality characteristics X, Y. Figure 8: Independent control of variability of two quality characteristics X, Y. Figure 9: Joint control of location of (X, Y ) and variability of (X, Y ). 5.1 Reduction of a Multivariate to a Univariate Case For some problems is possible to reduce the multivariate case to a one-dimensional one, so that the established one-dimensional capability indices can be used or the corresponding one-dimensional control charts can be implemented for process control.
11 Statistical Process Control Methods 59 3 Example Consider the location of the center of a hole defined by its 2-dimensional coordinates (X, Y ). Assume that X and Y are independent and normally distributed random variables with the same variance σ 2. Let the deviation of X and of Y from given targets x 0 and y 0 be the quality characteristics. Setting x 0 =0andy 0 = 0 leaves the two-dimensional quality characteristic (X, Y ). The pair (X, Y ) has a two-dimensional normal distribution with density function f (X,Y ). f (X,Y ) = 1 1 σ 2 (x2 +y 2 ) (22) 2πσ 2 e The deviation from target determines the quality of the product. Thus, instead of (X, Y ) one can introduce the deviation Z as single quality characteristic defined by: Z = (X) 2 +(Y ) 2 (23) Of course, the random variable Z is not normally distributed but has a so-called Raleigh distribution with density function f Z,[2]. f Z (z) = z z2 e 2σ σ2 2 (24) By means of Z the two-dimensional problem is transferred to an equivalent one-dimensional one. A process capability analysis based on Z by means of the normal distribution yields inconsistent results [9]. Moreover, because of the equivalence one can take the two or the one dimensional characteristic for a capability analysis, without priority [7]. A graphical illustration of the problem is presented in Figure 10. Figure 10: Illustration of the quality characteristics (X, Y ), f (X,Y ) and f Z. 5.2 Multivariate Quality Analysis Let X =(X 1,...,X m )beam-dimensional quality characteristic. The multi-dimensional tolerance region for X is often given by a hypercube or an ellipsoid. The first step of the quality analysis consists of calculating the non-conformance probability p (see for instance [10]). p =1... f X (x 1,...,x m )dx 1...dx m (25) tol.region
12 60 Constantin Anghel Next, a tolerance ellipsoid E T is introduced, i.e. an ellipsoid with center at the target value t =(t 1,...,t m ) representing an event which occurs with probability of (1 p). E T is defined as solution of the following quadratic equation: ( x t) T Σ 1 ( x t) =χ 2 1 p;m 1 (26) where χ 2 1 p;m 1 is the (1 p)-quantile of the χ 2 -distribution with (m 1) degrees of freedom. For m = 1, i.e. the one-dimensional case, and the normal distribution the above introduced tolerance ellipsoid reduces to the given tolerance interval. [ ] Φ 1 (p L ), Φ 1 (1 p U ) σ =[L, U] (27) The relation of the tolerance region and the tolerance ellipsoid is illustrated for two dimensions in Figure 11. Figure 11: Tolerance region and tolerance ellipsoid in a two-dimensional case. The tolerance ellipsoid enables a straightforward generalization of capability indices to cover the multivariate case. volume of tolerance ellipsoid C p = (28) volume of 99.73%-ellipsoid where the 99.73%-ellipsoid is defined as solution of ( x t) T Σ 1 ( x t) =χ ;m 1 (29) Taam et.al. use instead of (28) a different definition. They introduce instead of the tolerance ellipsoid a modified tolerance ellipsoid, which is an ellipsoid within the tolerance region, centered at target and having maximum volume. Clearly, the modified tolerance ellipsoid is independent of the underlying multivariate distribution and, thus, in particular independent of the correlation coefficients. C p (T aam) = volume of modified tolerance ellipsoid volume of 99.73%-ellipsoid (30)
13 Statistical Process Control Methods 61 If C p (T aam) = 1 then the modified tolerance ellipsoid contains 99.73% of the possible process space and the tolerance itself in general more than 99.73% which is not in line with the definition of a C p -index. Let U be the vector of the upper specification limits, then the second index taking into account the actual value E[ X] of a multi-dimensional location parameter is denoted by C pk and C pm (T aam), respectively. ( E[ C pk = (1 k)c p where k = ) T ( X] t E[ ) X] t ( ) T ( (31) U t U t) C pm = 1+ C p (T aam) ( E[ X] T ( t) Σ 1 E[ X] ) t (32) where it is assumed here that the target value t is the center of the tolerance region. 4 Example Consider a two-dimensional quality characteristic distributed according to a bivariate normal distribution. The input parameters are the following: ( ) ( ) ( ) 0.0 t =, L 3.2 =, U 3.2 = (33) E[ X]= ( ), V[ X]= ( ) (34) The joint density function is given by: 1 e (x ) 2 2ρ( x )( y )+( y ) 2 f (X,Y ) (x, y) = 2(1 ρ 2 ) (35) 2πσ X σ Y 1 ρ 2 In Figure 12 the capability indices C p and C pk and C p (T aam) andc pm are displayed as functions of the correlation coefficient. 6 Conclusions In many industrial instances product quality depends on a multitude of dependent characteristics and, therefore, there is an urgent need for appropriate multivariate models and methods. One problem is to define suitable tolerance regions taking into account the correlation structure among the variables. The problem of tolerance regions is closely connected to the problem of deriving process capability indices. Another problem refers to process monitoring for detecting changes in process distribution. Assuming the normal model might lead to inferior methods and wrong decisions.
14 62 Constantin Anghel Figure 12: C p and C pk and C p (T aam) andc pm as functions of ρ. References [1] C. Anghel, H. Hausberger, W. Streinz (1992): Unsymmetriegrößen erster und zweiter Art richtig auswerten (Teil 1). QZ 37, [2] C. Anghel, H. Hausberger, W. Streinz (1993): Unsymmetriegrößen erster und zweiter Art richtig auswerten (Teil 2). QZ 38, ) [3] C. Anghel (1998): Die Prozesszentrierung richtig beurteilen. QZ 43, [4] Castagliola, P. (1996): Evaluation of non-normal process capability indices using Burrs distribution. Quality Engineering 4, pp [5] Clement, J.A. (1989): Process capability calculation for non-normal distributions. Quality Progress Sept., [6] Collani, E.v. and Killmann, F. (2001): A note on the convolution of uniform and related distributions and their use in quality control. EQC, 16,... [7] Dietrich, E., Schulze, A. (1996): Fähigkeitsbeurteilung bei Positionstoleranzen. QZ 41, [8] Leone, F.C., Nelson, L.S. and Nottingham, R.B. (1996): The folded normal distribution. Technometrics 3, [9] Littig, S.J. and Pollock, S.M. (1992): Capability measurements for multivariate processes: Definitions and an example for a gear carrier. Technical Report 92-42, Dept. Industrial and Operational Engineering, University of Michigan, Ann Arbor.
15 Statistical Process Control Methods 63 [10] Hamilton, D.C. and Lesperance, M.L. (1995): A comparison of methods for univariate and multivariate acceptance sampling by variables. Technometrics 37, [11] H. Schneider, W.J. Kasperski, T. Ledford and W. Kraushaar (1995,1996): Control charts for skewed and censored data. Quality Engineering 8, [12] W. Taam, P. Subbaiah and J.W. Liddy (1993): A note on multivariate capability indices. Journal of Applied Statistics 20, [13] D.J. Wheeler, D.S. Chambers (1989): Understanding Statistical Process Control. Statistical Process Controls, Inc. Dr. C. Anghel Quality Department BMW AG D Dingolfing Germany
STT 843 Key to Homework 1 Spring 2018
STT 843 Key to Homework Spring 208 Due date: Feb 4, 208 42 (a Because σ = 2, σ 22 = and ρ 2 = 05, we have σ 2 = ρ 2 σ σ22 = 2/2 Then, the mean and covariance of the bivariate normal is µ = ( 0 2 and Σ
More informationDiskussionsbeiträge des Fachgebietes Unternehmensforschung
Diskussionsbeiträge des Fachgebietes Unternehmensforschung New Methods in Multivariate Statistical Process Control MSPC) Lorenz Braun November 2001 Universität Hohenheim 510 B) Institut für Betriebswirtschaftslehre
More informationPractice Problems Section Problems
Practice Problems Section 4-4-3 4-4 4-5 4-6 4-7 4-8 4-10 Supplemental Problems 4-1 to 4-9 4-13, 14, 15, 17, 19, 0 4-3, 34, 36, 38 4-47, 49, 5, 54, 55 4-59, 60, 63 4-66, 68, 69, 70, 74 4-79, 81, 84 4-85,
More informationJRF (Quality, Reliability and Operations Research): 2013 INDIAN STATISTICAL INSTITUTE INSTRUCTIONS
JRF (Quality, Reliability and Operations Research): 2013 INDIAN STATISTICAL INSTITUTE INSTRUCTIONS The test is divided into two sessions (i) Forenoon session and (ii) Afternoon session. Each session is
More informationObjective Experiments Glossary of Statistical Terms
Objective Experiments Glossary of Statistical Terms This glossary is intended to provide friendly definitions for terms used commonly in engineering and science. It is not intended to be absolutely precise.
More informationTUTORIAL 2 STA437 WINTER 2015
TUTORIAL 2 STA437 WINTER 2015 AL NOSEDAL Contents 1. Assessing the assumption of Normality 1 1.1. Evaluating the Normality of the Univariate Marginal Distributions 1 1.2. Evaluating Bivariate Normality
More informationA Bivariate Weibull Regression Model
c Heldermann Verlag Economic Quality Control ISSN 0940-5151 Vol 20 (2005), No. 1, 1 A Bivariate Weibull Regression Model David D. Hanagal Abstract: In this paper, we propose a new bivariate Weibull regression
More informationM.Sc. (Final) DEGREE EXAMINATION, MAY Final Year. Statistics. Paper I STATISTICAL QUALITY CONTROL. Answer any FIVE questions.
(DMSTT ) M.Sc. (Final) DEGREE EXAMINATION, MAY 0. Final Year Statistics Paper I STATISTICAL QUALITY CONTROL Time : Three hours Maximum : 00 marks Answer any FIVE questions. All questions carry equal marks..
More informationConfirmation Sample Control Charts
Confirmation Sample Control Charts Stefan H. Steiner Dept. of Statistics and Actuarial Sciences University of Waterloo Waterloo, NL 3G1 Canada Control charts such as X and R charts are widely used in industry
More informationHow Measurement Error Affects the Four Ways We Use Data
Measurement error is generally considered to be a bad thing, and yet there is very little written about how measurement error affects the way we use our measurements. This column will consider these effects
More information15 Discrete Distributions
Lecture Note 6 Special Distributions (Discrete and Continuous) MIT 4.30 Spring 006 Herman Bennett 5 Discrete Distributions We have already seen the binomial distribution and the uniform distribution. 5.
More informationSTAT 135 Lab 6 Duality of Hypothesis Testing and Confidence Intervals, GLRT, Pearson χ 2 Tests and Q-Q plots. March 8, 2015
STAT 135 Lab 6 Duality of Hypothesis Testing and Confidence Intervals, GLRT, Pearson χ 2 Tests and Q-Q plots March 8, 2015 The duality between CI and hypothesis testing The duality between CI and hypothesis
More informationMultivariate Distribution Models
Multivariate Distribution Models Model Description While the probability distribution for an individual random variable is called marginal, the probability distribution for multiple random variables is
More informationIntroductory Econometrics. Review of statistics (Part II: Inference)
Introductory Econometrics Review of statistics (Part II: Inference) Jun Ma School of Economics Renmin University of China October 1, 2018 1/16 Null and alternative hypotheses Usually, we have two competing
More informationOptimal SPRT and CUSUM Procedures using Compressed Limit Gauges
Optimal SPRT and CUSUM Procedures using Compressed Limit Gauges P. Lee Geyer Stefan H. Steiner 1 Faculty of Business McMaster University Hamilton, Ontario L8S 4M4 Canada Dept. of Statistics and Actuarial
More informationCH.9 Tests of Hypotheses for a Single Sample
CH.9 Tests of Hypotheses for a Single Sample Hypotheses testing Tests on the mean of a normal distributionvariance known Tests on the mean of a normal distributionvariance unknown Tests on the variance
More informationSummarizing Measured Data
Summarizing Measured Data 12-1 Overview Basic Probability and Statistics Concepts: CDF, PDF, PMF, Mean, Variance, CoV, Normal Distribution Summarizing Data by a Single Number: Mean, Median, and Mode, Arithmetic,
More informationSequential Procedure for Testing Hypothesis about Mean of Latent Gaussian Process
Applied Mathematical Sciences, Vol. 4, 2010, no. 62, 3083-3093 Sequential Procedure for Testing Hypothesis about Mean of Latent Gaussian Process Julia Bondarenko Helmut-Schmidt University Hamburg University
More informationNotes on the Multivariate Normal and Related Topics
Version: July 10, 2013 Notes on the Multivariate Normal and Related Topics Let me refresh your memory about the distinctions between population and sample; parameters and statistics; population distributions
More informationEco517 Fall 2004 C. Sims MIDTERM EXAM
Eco517 Fall 2004 C. Sims MIDTERM EXAM Answer all four questions. Each is worth 23 points. Do not devote disproportionate time to any one question unless you have answered all the others. (1) We are considering
More informationDirectionally Sensitive Multivariate Statistical Process Control Methods
Directionally Sensitive Multivariate Statistical Process Control Methods Ronald D. Fricker, Jr. Naval Postgraduate School October 5, 2005 Abstract In this paper we develop two directionally sensitive statistical
More information7. Estimation and hypothesis testing. Objective. Recommended reading
7. Estimation and hypothesis testing Objective In this chapter, we show how the election of estimators can be represented as a decision problem. Secondly, we consider the problem of hypothesis testing
More informationEE 302: Probabilistic Methods in Electrical Engineering
EE : Probabilistic Methods in Electrical Engineering Print Name: Solution (//6 --sk) Test II : Chapters.5 4 //98, : PM Write down your name on each paper. Read every question carefully and solve each problem
More informationCourse: ESO-209 Home Work: 1 Instructor: Debasis Kundu
Home Work: 1 1. Describe the sample space when a coin is tossed (a) once, (b) three times, (c) n times, (d) an infinite number of times. 2. A coin is tossed until for the first time the same result appear
More informationProbability, Statistics, and Reliability for Engineers and Scientists FUNDAMENTALS OF STATISTICAL ANALYSIS
CHAPTER Probability, Statistics, and Reliability for Engineers and Scientists FUNDAMENTALS OF STATISTICAL ANALYSIS Second Edition A. J. Clark School of Engineering Department of Civil and Environmental
More informationInstitute of Actuaries of India
Institute of Actuaries of India Subject CT3 Probability & Mathematical Statistics May 2011 Examinations INDICATIVE SOLUTION Introduction The indicative solution has been written by the Examiners with the
More informationSolutions to Final STAT 421, Fall 2008
Solutions to Final STAT 421, Fall 2008 Fritz Scholz 1. (8) Two treatments A and B were randomly assigned to 8 subjects (4 subjects to each treatment) with the following responses: 0, 1, 3, 6 and 5, 7,
More informationSTAT 6350 Analysis of Lifetime Data. Probability Plotting
STAT 6350 Analysis of Lifetime Data Probability Plotting Purpose of Probability Plots Probability plots are an important tool for analyzing data and have been particular popular in the analysis of life
More informationarxiv: v1 [stat.me] 14 Jan 2019
arxiv:1901.04443v1 [stat.me] 14 Jan 2019 An Approach to Statistical Process Control that is New, Nonparametric, Simple, and Powerful W.J. Conover, Texas Tech University, Lubbock, Texas V. G. Tercero-Gómez,Tecnológico
More informationSubject CS1 Actuarial Statistics 1 Core Principles
Institute of Actuaries of India Subject CS1 Actuarial Statistics 1 Core Principles For 2019 Examinations Aim The aim of the Actuarial Statistics 1 subject is to provide a grounding in mathematical and
More informationCorrelation and Regression Theory 1) Multivariate Statistics
Correlation and Regression Theory 1) Multivariate Statistics What is a multivariate data set? How to statistically analyze this data set? Is there any kind of relationship between different variables in
More information18.05 Practice Final Exam
No calculators. 18.05 Practice Final Exam Number of problems 16 concept questions, 16 problems. Simplifying expressions Unless asked to explicitly, you don t need to simplify complicated expressions. For
More informationThe Delta Method and Applications
Chapter 5 The Delta Method and Applications 5.1 Local linear approximations Suppose that a particular random sequence converges in distribution to a particular constant. The idea of using a first-order
More informationThe Model Building Process Part I: Checking Model Assumptions Best Practice
The Model Building Process Part I: Checking Model Assumptions Best Practice Authored by: Sarah Burke, PhD 31 July 2017 The goal of the STAT T&E COE is to assist in developing rigorous, defensible test
More information1. Density and properties Brief outline 2. Sampling from multivariate normal and MLE 3. Sampling distribution and large sample behavior of X and S 4.
Multivariate normal distribution Reading: AMSA: pages 149-200 Multivariate Analysis, Spring 2016 Institute of Statistics, National Chiao Tung University March 1, 2016 1. Density and properties Brief outline
More informationStatistical Methods in Particle Physics
Statistical Methods in Particle Physics Lecture 3 October 29, 2012 Silvia Masciocchi, GSI Darmstadt s.masciocchi@gsi.de Winter Semester 2012 / 13 Outline Reminder: Probability density function Cumulative
More informationProbability Theory and Statistics. Peter Jochumzen
Probability Theory and Statistics Peter Jochumzen April 18, 2016 Contents 1 Probability Theory And Statistics 3 1.1 Experiment, Outcome and Event................................ 3 1.2 Probability............................................
More informationx. Figure 1: Examples of univariate Gaussian pdfs N (x; µ, σ 2 ).
.8.6 µ =, σ = 1 µ = 1, σ = 1 / µ =, σ =.. 3 1 1 3 x Figure 1: Examples of univariate Gaussian pdfs N (x; µ, σ ). The Gaussian distribution Probably the most-important distribution in all of statistics
More informationLecture 2: Basic Concepts and Simple Comparative Experiments Montgomery: Chapter 2
Lecture 2: Basic Concepts and Simple Comparative Experiments Montgomery: Chapter 2 Fall, 2013 Page 1 Random Variable and Probability Distribution Discrete random variable Y : Finite possible values {y
More informationTail dependence in bivariate skew-normal and skew-t distributions
Tail dependence in bivariate skew-normal and skew-t distributions Paola Bortot Department of Statistical Sciences - University of Bologna paola.bortot@unibo.it Abstract: Quantifying dependence between
More informationSTAT 461/561- Assignments, Year 2015
STAT 461/561- Assignments, Year 2015 This is the second set of assignment problems. When you hand in any problem, include the problem itself and its number. pdf are welcome. If so, use large fonts and
More informationThe Model Building Process Part I: Checking Model Assumptions Best Practice (Version 1.1)
The Model Building Process Part I: Checking Model Assumptions Best Practice (Version 1.1) Authored by: Sarah Burke, PhD Version 1: 31 July 2017 Version 1.1: 24 October 2017 The goal of the STAT T&E COE
More informationCUMULATIVE SUM CHARTS FOR HIGH YIELD PROCESSES
Statistica Sinica 11(2001), 791-805 CUMULATIVE SUM CHARTS FOR HIGH YIELD PROCESSES T. C. Chang and F. F. Gan Infineon Technologies Melaka and National University of Singapore Abstract: The cumulative sum
More informationUnit 10: Planning Life Tests
Unit 10: Planning Life Tests Ramón V. León Notes largely based on Statistical Methods for Reliability Data by W.Q. Meeker and L. A. Escobar, Wiley, 1998 and on their class notes. 11/2/2004 Unit 10 - Stat
More informationProblem 1 (20) Log-normal. f(x) Cauchy
ORF 245. Rigollet Date: 11/21/2008 Problem 1 (20) f(x) f(x) 0.0 0.1 0.2 0.3 0.4 0.0 0.2 0.4 0.6 0.8 4 2 0 2 4 Normal (with mean -1) 4 2 0 2 4 Negative-exponential x x f(x) f(x) 0.0 0.1 0.2 0.3 0.4 0.5
More information401 Review. 6. Power analysis for one/two-sample hypothesis tests and for correlation analysis.
401 Review Major topics of the course 1. Univariate analysis 2. Bivariate analysis 3. Simple linear regression 4. Linear algebra 5. Multiple regression analysis Major analysis methods 1. Graphical analysis
More information2.830J / 6.780J / ESD.63J Control of Manufacturing Processes (SMA 6303) Spring 2008
MIT OpenCourseWare http://ocw.mit.edu.830j / 6.780J / ESD.63J Control of Manufacturing Processes (SMA 6303) Spring 008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationContents. Preface to Second Edition Preface to First Edition Abbreviations PART I PRINCIPLES OF STATISTICAL THINKING AND ANALYSIS 1
Contents Preface to Second Edition Preface to First Edition Abbreviations xv xvii xix PART I PRINCIPLES OF STATISTICAL THINKING AND ANALYSIS 1 1 The Role of Statistical Methods in Modern Industry and Services
More informationReview of Statistics 101
Review of Statistics 101 We review some important themes from the course 1. Introduction Statistics- Set of methods for collecting/analyzing data (the art and science of learning from data). Provides methods
More information01 Probability Theory and Statistics Review
NAVARCH/EECS 568, ROB 530 - Winter 2018 01 Probability Theory and Statistics Review Maani Ghaffari January 08, 2018 Last Time: Bayes Filters Given: Stream of observations z 1:t and action data u 1:t Sensor/measurement
More informationFall 2017 STAT 532 Homework Peter Hoff. 1. Let P be a probability measure on a collection of sets A.
1. Let P be a probability measure on a collection of sets A. (a) For each n N, let H n be a set in A such that H n H n+1. Show that P (H n ) monotonically converges to P ( k=1 H k) as n. (b) For each n
More informationEconomics 520. Lecture Note 19: Hypothesis Testing via the Neyman-Pearson Lemma CB 8.1,
Economics 520 Lecture Note 9: Hypothesis Testing via the Neyman-Pearson Lemma CB 8., 8.3.-8.3.3 Uniformly Most Powerful Tests and the Neyman-Pearson Lemma Let s return to the hypothesis testing problem
More informationOne-Way Repeated Measures Contrasts
Chapter 44 One-Way Repeated easures Contrasts Introduction This module calculates the power of a test of a contrast among the means in a one-way repeated measures design using either the multivariate test
More informationBasic Statistics. 1. Gross error analyst makes a gross mistake (misread balance or entered wrong value into calculation).
Basic Statistics There are three types of error: 1. Gross error analyst makes a gross mistake (misread balance or entered wrong value into calculation). 2. Systematic error - always too high or too low
More informationChapter 4. One-sided Process Capability Assessment in the Presence of Gauge Measurement Errors
hapter One-sided Process apability Assessment in the Presence of Gauge Measurement Errors n the manufacturing industry, many product characteristics are of one-sided specifications The process capability
More informationEE/CpE 345. Modeling and Simulation. Fall Class 9
EE/CpE 345 Modeling and Simulation Class 9 208 Input Modeling Inputs(t) Actual System Outputs(t) Parameters? Simulated System Outputs(t) The input data is the driving force for the simulation - the behavior
More informationReview of Statistics
Review of Statistics Topics Descriptive Statistics Mean, Variance Probability Union event, joint event Random Variables Discrete and Continuous Distributions, Moments Two Random Variables Covariance and
More informationChapter 3 sections. SKIP: 3.10 Markov Chains. SKIP: pages Chapter 3 - continued
Chapter 3 sections Chapter 3 - continued 3.1 Random Variables and Discrete Distributions 3.2 Continuous Distributions 3.3 The Cumulative Distribution Function 3.4 Bivariate Distributions 3.5 Marginal Distributions
More informationMultivariate probability distributions and linear regression
Multivariate probability distributions and linear regression Patrik Hoyer 1 Contents: Random variable, probability distribution Joint distribution Marginal distribution Conditional distribution Independence,
More information2. What are the tradeoffs among different measures of error (e.g. probability of false alarm, probability of miss, etc.)?
ECE 830 / CS 76 Spring 06 Instructors: R. Willett & R. Nowak Lecture 3: Likelihood ratio tests, Neyman-Pearson detectors, ROC curves, and sufficient statistics Executive summary In the last lecture we
More informationLectures on Simple Linear Regression Stat 431, Summer 2012
Lectures on Simple Linear Regression Stat 43, Summer 0 Hyunseung Kang July 6-8, 0 Last Updated: July 8, 0 :59PM Introduction Previously, we have been investigating various properties of the population
More informationA Multivariate EWMA Control Chart for Skewed Populations using Weighted Variance Method
OPEN ACCESS Int. Res. J. of Science & Engineering, 04; Vol. (6): 9-0 ISSN: 3-005 RESEARCH ARTICLE A Multivariate EMA Control Chart for Skewed Populations using eighted Variance Method Atta AMA *, Shoraim
More informationECE 313: Conflict Final Exam Tuesday, May 13, 2014, 7:00 p.m. 10:00 p.m. Room 241 Everitt Lab
University of Illinois Spring 1 ECE 313: Conflict Final Exam Tuesday, May 13, 1, 7: p.m. 1: p.m. Room 1 Everitt Lab 1. [18 points] Consider an experiment in which a fair coin is repeatedly tossed every
More informationSTA 2201/442 Assignment 2
STA 2201/442 Assignment 2 1. This is about how to simulate from a continuous univariate distribution. Let the random variable X have a continuous distribution with density f X (x) and cumulative distribution
More informationparameter space Θ, depending only on X, such that Note: it is not θ that is random, but the set C(X).
4. Interval estimation The goal for interval estimation is to specify the accurary of an estimate. A 1 α confidence set for a parameter θ is a set C(X) in the parameter space Θ, depending only on X, such
More informationWithin Cases. The Humble t-test
Within Cases The Humble t-test 1 / 21 Overview The Issue Analysis Simulation Multivariate 2 / 21 Independent Observations Most statistical models assume independent observations. Sometimes the assumption
More informationQuestion. Hypothesis testing. Example. Answer: hypothesis. Test: true or not? Question. Average is not the mean! μ average. Random deviation or not?
Hypothesis testing Question Very frequently: what is the possible value of μ? Sample: we know only the average! μ average. Random deviation or not? Standard error: the measure of the random deviation.
More informationMultivariate Simulations
Multivariate Simulations Katarína Starinská starinskak@gmail.com Charles University Faculty of Mathematics and Physics Prague, Czech Republic 21.11.2011 Katarína Starinská Multivariate Simulations 1 /
More informationApplied Econometrics (QEM)
Applied Econometrics (QEM) based on Prinicples of Econometrics Jakub Mućk Department of Quantitative Economics Jakub Mućk Applied Econometrics (QEM) Meeting #3 1 / 42 Outline 1 2 3 t-test P-value Linear
More informationMathematical and Computer Modelling. Economic design of EWMA control charts based on loss function
Mathematical and Computer Modelling 49 (2009) 745 759 Contents lists available at ScienceDirect Mathematical and Computer Modelling journal homepage: www.elsevier.com/locate/mcm Economic design of EWMA
More informationControl charts are used for monitoring the performance of a quality characteristic. They assist process
QUALITY AND RELIABILITY ENGINEERING INTERNATIONAL Qual. Reliab. Engng. Int. 2009; 25:875 883 Published online 3 March 2009 in Wiley InterScience (www.interscience.wiley.com)..1007 Research Identifying
More informationPreliminary Statistics Lecture 5: Hypothesis Testing (Outline)
1 School of Oriental and African Studies September 2015 Department of Economics Preliminary Statistics Lecture 5: Hypothesis Testing (Outline) Gujarati D. Basic Econometrics, Appendix A.8 Barrow M. Statistics
More informationJRF (Quality, Reliability & Operations Research): 2015 INDIAN STATISTICAL INSTITUTE INSTRUCTIONS
JRF (Quality, Reliability & Operations Research): 2015 INDIAN STATISTICAL INSTITUTE INSTRUCTIONS The test is divided into two sessions (i) Forenoon session and (ii) Afternoon session. Each session is of
More information2.830J / 6.780J / ESD.63J Control of Manufacturing Processes (SMA 6303) Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 2.830J / 6.780J / ESD.63J Control of Manufacturing Processes (SMA 6303) Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationProject Report for STAT571 Statistical Methods Instructor: Dr. Ramon V. Leon. Wage Data Analysis. Yuanlei Zhang
Project Report for STAT7 Statistical Methods Instructor: Dr. Ramon V. Leon Wage Data Analysis Yuanlei Zhang 77--7 November, Part : Introduction Data Set The data set contains a random sample of observations
More informationTable of z values and probabilities for the standard normal distribution. z is the first column plus the top row. Each cell shows P(X z).
Table of z values and probabilities for the standard normal distribution. z is the first column plus the top row. Each cell shows P(X z). For example P(X.04) =.8508. For z < 0 subtract the value from,
More informationUnit 14: Nonparametric Statistical Methods
Unit 14: Nonparametric Statistical Methods Statistics 571: Statistical Methods Ramón V. León 8/8/2003 Unit 14 - Stat 571 - Ramón V. León 1 Introductory Remarks Most methods studied so far have been based
More informationSPRING 2007 EXAM C SOLUTIONS
SPRING 007 EXAM C SOLUTIONS Question #1 The data are already shifted (have had the policy limit and the deductible of 50 applied). The two 350 payments are censored. Thus the likelihood function is L =
More informationParameter Estimation and Fitting to Data
Parameter Estimation and Fitting to Data Parameter estimation Maximum likelihood Least squares Goodness-of-fit Examples Elton S. Smith, Jefferson Lab 1 Parameter estimation Properties of estimators 3 An
More information2. Variance and Covariance: We will now derive some classic properties of variance and covariance. Assume real-valued random variables X and Y.
CS450 Final Review Problems Fall 08 Solutions or worked answers provided Problems -6 are based on the midterm review Identical problems are marked recap] Please consult previous recitations and textbook
More informationHypothesis Testing One Sample Tests
STATISTICS Lecture no. 13 Department of Econometrics FEM UO Brno office 69a, tel. 973 442029 email:jiri.neubauer@unob.cz 12. 1. 2010 Tests on Mean of a Normal distribution Tests on Variance of a Normal
More informationStatistical techniques for data analysis in Cosmology
Statistical techniques for data analysis in Cosmology arxiv:0712.3028; arxiv:0911.3105 Numerical recipes (the bible ) Licia Verde ICREA & ICC UB-IEEC http://icc.ub.edu/~liciaverde outline Lecture 1: Introduction
More informationPolitical Science 236 Hypothesis Testing: Review and Bootstrapping
Political Science 236 Hypothesis Testing: Review and Bootstrapping Rocío Titiunik Fall 2007 1 Hypothesis Testing Definition 1.1 Hypothesis. A hypothesis is a statement about a population parameter The
More informationDistribution Fitting (Censored Data)
Distribution Fitting (Censored Data) Summary... 1 Data Input... 2 Analysis Summary... 3 Analysis Options... 4 Goodness-of-Fit Tests... 6 Frequency Histogram... 8 Comparison of Alternative Distributions...
More information2.830J / 6.780J / ESD.63J Control of Manufacturing Processes (SMA 6303) Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 2.830J / 6.780J / ESD.63J Control of Processes (SMA 6303) Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationLecture Note 1: Probability Theory and Statistics
Univ. of Michigan - NAME 568/EECS 568/ROB 530 Winter 2018 Lecture Note 1: Probability Theory and Statistics Lecturer: Maani Ghaffari Jadidi Date: April 6, 2018 For this and all future notes, if you would
More informationChapter 23: Inferences About Means
Chapter 3: Inferences About Means Sample of Means: number of observations in one sample the population mean (theoretical mean) sample mean (observed mean) is the theoretical standard deviation of the population
More informationStatistical Hypothesis Testing
Statistical Hypothesis Testing Dr. Phillip YAM 2012/2013 Spring Semester Reference: Chapter 7 of Tests of Statistical Hypotheses by Hogg and Tanis. Section 7.1 Tests about Proportions A statistical hypothesis
More informationHomework 10 (due December 2, 2009)
Homework (due December, 9) Problem. Let X and Y be independent binomial random variables with parameters (n, p) and (n, p) respectively. Prove that X + Y is a binomial random variable with parameters (n
More informationCOMPARISON OF FIVE TESTS FOR THE COMMON MEAN OF SEVERAL MULTIVARIATE NORMAL POPULATIONS
Communications in Statistics - Simulation and Computation 33 (2004) 431-446 COMPARISON OF FIVE TESTS FOR THE COMMON MEAN OF SEVERAL MULTIVARIATE NORMAL POPULATIONS K. Krishnamoorthy and Yong Lu Department
More informationTesting Simple Hypotheses R.L. Wolpert Institute of Statistics and Decision Sciences Duke University, Box Durham, NC 27708, USA
Testing Simple Hypotheses R.L. Wolpert Institute of Statistics and Decision Sciences Duke University, Box 90251 Durham, NC 27708, USA Summary: Pre-experimental Frequentist error probabilities do not summarize
More informationA Comparison of Equivalence Criteria and Basis Values for HEXCEL 8552 IM7 Unidirectional Tape computed from the NCAMP shared database
Report No: Report Date: A Comparison of Equivalence Criteria and Basis Values for HEXCEL 8552 IM7 Unidirectional Tape computed from the NCAMP shared database NCAMP Report Number: Report Date: Elizabeth
More informationLecture 13: p-values and union intersection tests
Lecture 13: p-values and union intersection tests p-values After a hypothesis test is done, one method of reporting the result is to report the size α of the test used to reject H 0 or accept H 0. If α
More informationYou can compute the maximum likelihood estimate for the correlation
Stat 50 Solutions Comments on Assignment Spring 005. (a) _ 37.6 X = 6.5 5.8 97.84 Σ = 9.70 4.9 9.70 75.05 7.80 4.9 7.80 4.96 (b) 08.7 0 S = Σ = 03 9 6.58 03 305.6 30.89 6.58 30.89 5.5 (c) You can compute
More informationEconometrics Summary Algebraic and Statistical Preliminaries
Econometrics Summary Algebraic and Statistical Preliminaries Elasticity: The point elasticity of Y with respect to L is given by α = ( Y/ L)/(Y/L). The arc elasticity is given by ( Y/ L)/(Y/L), when L
More informationStatistical Inference
Statistical Inference Bernhard Klingenberg Institute of Statistics Graz University of Technology Steyrergasse 17/IV, 8010 Graz www.statistics.tugraz.at February 12, 2008 Outline Estimation: Review of concepts
More informationInferences about a Mean Vector
Inferences about a Mean Vector Edps/Soc 584, Psych 594 Carolyn J. Anderson Department of Educational Psychology I L L I N O I S university of illinois at urbana-champaign c Board of Trustees, University
More informationQualifying Exam CS 661: System Simulation Summer 2013 Prof. Marvin K. Nakayama
Qualifying Exam CS 661: System Simulation Summer 2013 Prof. Marvin K. Nakayama Instructions This exam has 7 pages in total, numbered 1 to 7. Make sure your exam has all the pages. This exam will be 2 hours
More informationGlossary. The ISI glossary of statistical terms provides definitions in a number of different languages:
Glossary The ISI glossary of statistical terms provides definitions in a number of different languages: http://isi.cbs.nl/glossary/index.htm Adjusted r 2 Adjusted R squared measures the proportion of the
More informationLab 6. Current Balance
Lab 6. Current Balance Goals To explore and verify the right-hand rule governing the force on a current-carrying wire immersed in a magnetic field. To determine how the force on a current-carrying wire
More information