Control Charts Based on Alternative Hypotheses
|
|
- Agnes Underwood
- 6 years ago
- Views:
Transcription
1 Control Charts Based on Alternative Hypotheses A. Di Bucchianico, M. Hušková (Prague), P. Klášterecky (Prague), W.R. van Zwet (Leiden) Dortmund, January 11, /48
2 Goals of this talk introduce hypothesis testing framework for control charts in SPC develop monitoring procedures for practical out-of-control situations 2/48
3 Contents of talk Background: statistical process control: control charts change-point problems sequential analysis Testing and control charts Some alternative hypotheses Some thoughts on performance measures of monitoring procedures Likelihood ratio tests Asymptotics for critical values Simulations Future work 3/48
4 Background 1. SPC (Statistical Process Control) background in industry (Shewhart 1924) uses control charts as monitoring tools (detection of out-of-control situations) emphasis on on-line monitoring (Phase II) 2. changepoint analysis (cf. Lai, J. Roy. Stat. Soc. B 57 (1995)) background in mathematical statistics aims at estimation of changepoint main emphasis on retrospective analysis 3. sequential analysis (cf. Lai, Stat. Sinica 11 (2001)) developed in military context (Wald, Wolfowitz 1940 s) initial emphasis on hypothesis testing 4/48
5 SPC: Shewhart charts introduced by Shewhart in 1924 practical tool without theoretical background specific terminology: in-control (common causes), out-of-control (special causes), rational subgroups,... chart signals if summary statistic of ith group is above or below 3σ T variants: X, R, S, MR, attribute control charts,... additions: VSR (Variable Sampling Rate), runs rules (Western Electric 1956), warning zones,... 5/48
6 Shewhart X-chart with control lines. µ + 3σ/ n µ UCL CL µ 3σ/ n LCL 6/48
7 SPC: CUSUM charts introduced in 1954 by Page cumulative sums enable to detect small changes of the mean (< 1.5σ ) recursive practical form: threshold on cumulative sums: Q i = max{0, Q i 1 + X i k)} optimality with respect to ARL proved by Moustakides and Ritov performance quickly deteriorates away from optimal alternative (finetuning of k and decision threshold) additions: FIR (Fast Initial Response) by Lucas monograph Hawkins and Olwell 7/48
8 SPC: EWMA chart based on ideas from Girshick, Rubin, Roberts and Shiryaev (early 1960 s) inspired by Bayesian analysis (prior distribution on changepoint) V i = λ(x i ) + (1 λ)v i 1, V 0 = 0 λ 0: CUSUM, λ = 1: Shewhart practical choice for λ: 0.1 < λ < 0.3 performance nearly as good as CUSUM, but less sensitive to nonnormality decision threshold changes with i 8/48
9 SPC: other charts control charts based on robust statistics combined control charts (e.g., Shewhart-CUSUM) cuscore charts (Box; based on Fisher s efficient score statistics) nonparametric control charts linear rank statistics (Wilcoxon,... ) precedence statistics (Chakraborti and V.d. Laan) regression-type control charts (monitoring linear profiles)... 9/48
10 SPC: Phase I and II Phase I retrospective (usually pilot study of new production process) determination of in-control parameter values Phase II on-line (full scale process) uses in-control parameter values from Phase I Detection performance of control charts in Phase II may heavily deteriorate when using estimated parameters in the control limits (see e.g., Chakraborti, Comm. Stat. Simul. 29 (2000)). Robust estimation of parameters is required in noisy environments (see e.g., Gather et al., Estadistica 53 (2001)). Hawkins et al. (J. Qual. Techn. 35 (2003)) argue that application likelihood ratio methods from changepoint analysis in SPC context makes Phase I/II distinction superfluous. / department of mathematics and computer science 10/48
11 sample X 1,..., X n H 0 : µ = µ 0 12 Neyman-Pearson Lemma H 1 : µ = µ 1 reject H 0 if P H 1 (x 1,..., x n ) P H0 (x 1,..., x n ) > c. This test has maximal power under all tests with the same type I error. Frisén and De Maré, Biometrika 78 (1991), have version of Neyman-Pearson for detection of critical events at given time point. Tests for composite hypotheses H 0 : µ µ 0 against H 1 : µ > µ 0 may be treated similarly if likelihood ratio is monotone. 11/48
12 Wald SPRT (Sequential Probability Ratio Test) sample X 1, X 2,... H 0 : µ = µ 0 against H 1 : µ = µ 1 accept H 0 accept H 1 if if P H1 (x 1,..., x n ) P H0 (x 1,..., x n ) < a P H1 (x 1,..., x n ) P H0 (x 1,..., x n ) > b continue testing if a P H 1 (x 1,..., x n ) P H0 (x 1,..., x n ) b This procedure simultaneously minimizes the ASN s (Average Sample Number) under H 0 and H 1, given P H0 (reject H 0 ) α and P H1 (reject H 0 ) β. 12/48
13 Summary of existing procedures Existing procedures either are not based on optimality criteria or are optimal w.r.t. to ARL (Average Run Length). However, run length distributions are highly skewed (often close to geometric distribution). Performance should be judged on other features of run length distributions (e.g., quantiles). Cf. PSD (Probability of Successful Detection) in medical applications (Frisén, Stat. Medicine 11 (1992) and Frisé and Wessmann, Comm. Stat. Simul. 28 (1999)). In case of active surveillance: predictive value of alarm. It would be nice to reconcile classical hypothesis testing with sequential detection. This would allow to develop optimal monitoring procedures for specific alternative hypotheses. 13/48
14 Examples of alternative hypotheses Usually theoretical studies of control charts focus on persistent changes of the mean (mathematical convenience?!). Examples of other alternative hypotheses (cf. Gitlow et al., Quality Management: Tools and Methods for Improvement, Chapter 8): persistent threshold crossing (sea and river levels) persistent monotone threshold crossing (tool wear, e.g., chisels) persistent shift or drift in variance (wear of bearing) epidemic alternatives (e.g., joint SPC APC scheme: feedback controller removes special cause )... One may also think of monitoring cycles (e.g., business cycles or medical cycles (Frisén)) or profiles (cf. Kim et al., J. Qual. Technology 35 (2003) ). 14/48
15 Our Model Sequential observations X 1,..., X k (k n) n fixed beforehand (batch industry, medical event,... ) X i = µ i + e i, i = 1,..., n, e 1,..., e n are i.i.d. error terms with density f symmetric around 0 µ 1,..., µ n are unknown parameters. Test statistic based on likelihood ratios (details later). Critical value of test statistics based on false alarm rate (details later). 15/48
16 Possible alternative hypotheses Persistent change of mean (µ 0 and δ known) H 0 : µ 1 =... = µ n { µ i = µ 0, i = 1,..., m, H 1 : µ i = µ 0 + δ, i = m + 1,..., n Epidemic alternative (µ 0 and δ known) H 0 : µ 1 =... = µ n µ i = µ 0, i = 1,..., l, H 1 : µ i = µ 0 + δ, i = l + 1,..., m, µ i = µ 0, i = m,..., n 16/48
17 Possible alternative hypotheses: continued Persistent non-monotone threshold crossing (δ known) H 0 : µ i δ, i = 1,..., n, { µ i δ, i = 1,..., m, H 1 : µ i > δ, i = m + 1,..., n, Persistent monotone threshold crossing (δ known) H 0 : µ 1... µ n δ { µ 1... µ m < δ i = 1,..., m H 1 : δ < µ m+1 µ m+2... µ n i = m + 1,..., n Similar alternative hypotheses for the variance are also important (cf. philosophy SPC). 17/48
18 Toy example: Shewhart chart + persistent threshold crossing alternative X i independent with densities f (x µ i ), where f is a symmetric density that is nonincreasing for x > 0. Hypotheses with known δ 0 : against H 0i : E X i = µ i δ 0 H 1i : E X i = µ i > δ 0 The ML-estimators µ i0 and µ i1 based on X i under H 0i and H 1i, respectively, are given by: { { X i if X i δ 0 δ 0 if X i δ 0 µ i0 = µ i1 = δ 0 if X i > δ 0 X i if X i > δ 0 18/48
19 The loglikelihood ratio based on X i thus equals log ( f (X i µ i1 )/f (X i µ 0i ) ) Under normality with common known variance σ 2 this reduces up to a multiplicative constant 1/2 to: Z i (δ 0 ) = (X i δ 0 ) 2 σ 2 sign (X i δ 0 ). Since (t δ 0 ) 2 sign (δ 0 t) is a monotone function of t, we restrict ourselves to the least favourable situation of the null hypothesis H 0 : the restricted null hypothesis H 0. H 0 : µ 1 = µ 2 =... = δ 0. N S = { min{1 k n; Z k (δ 0 ) c} if max 1 k n Z k c if max 1 k n Z k < c. 19/48
20 Critical values Z i Standard choice of 3σ control limits yields in-control ARL of 1/ E H 0 N S = 1 1 H(c), where H is distribution function of Z i. Note that E H 0 N S is not meaningful because there is no good way to incorporate the event {max 1 k n Z k < c}. Alternative 1: jointly considering H(c) n, the probability of no alarm at all, and the conditional ARL: E H 0 (N S N S n) = 1 1 H(c) n(h(c))n 1 H(c) n. Alternative 2: 100α% quantile of the distribution of N S. ( ) 1 H(c) ns(α) = P H 0 max Z i(δ 0 ) c = P H 1 i n S (α) 0 (N S n S (α)) α 20/48
21 Critical values Z i : continued Alternative 2: 100α% quantile of the distribution of N S. The lower bound n S (α) for the 100α% quantile of the distribution of N S satisfies ( ) P H 0 max Z i(δ 0 ) c = P H 1 i n S (α) 0 (N S n S (α)) α In other words, the test H 0 against H 1 based on n S (α) observations has significance level α. Hence, this choice of the critical value can be easily connected to the Neyman-Pearson framework. Thinking in terms of hypothesis testing with level α it will be typically chosen small, while thinking in terms of the median of the run length distribution we think more about α = 1/2. The next slides present tables with numerical values. 21/48
22 quantile of N S n E H 0 (N S N S n) P H 0 (N S > n) 1% 5% 50% 80% 90% Influence of n on the run length distribution: c = 9. 22/48
23 β = 0.8 β = 0.9 β = 1 n 90% 95% 99% 90% 95% 99% P H0 (N S > n) = Critical values for Shewhart charts with prescribed lower bound for quantiles of run length distribution and for prescribed false alarm rate: n s = β n. 23/48
24 Detection performance Z i distribution function H of Z i (δ 0 ) if E X i = δ has the form: ( ) x δ δ 0 σ, x < 0 H(x) = ( x ) δ δ 0 σ, x 0 In the next table, we assume that no false alarm takes place before the changepoint, which gives a good impression of the performance, because the Shewhart chart is memoryless. 24/48
25 H 1 -quantile of N S m N S > m n m c P H1 (N S > n N s > m) 1% 5% 50% 80% 90% Detection performance for Shewhart charts with prescribed false alarm probability 0.05 under jump alternative with δ 1 = δ 0 + σ : E X i = δ 0 + σ for i > m. 25/48
26 Detection performance Z i : uniform alternative Another interesting alternative is that the mean of X i after the changepoint is uniformly distributed on [δ 0, δ 0 + γ σ ] with γ > 0. We then have P H1 (Z i c) = 1 γ σ γ σ 0 P(Z i c µ i = δ 0 +x) dx = 1 γ σ γ σ 0 ( c x ) σ Writing as an integral and interchanging the order of integration, we obtain that b a (u) du = b (b) a (a) + ϕ(b) ϕ(a), where ϕ is the density of the standard normal distribution. Combining everything, we obtain P H1 (Z i c) = c ( c) ( c γ ) ( c γ ) + ϕ( c) ϕ( c γ ). dx 26/48
27 General approach: likelihood ratio tests max 1 m k max θ mi=1 f (x i ; θ) max η ki=m+1 f i (x i ; η) max θ ki=1 f (x i ; θ) Persistent non-monotone threshold crossing and Persistent monotone threshold crossing for normally distributed data with known variance: 1 k Q n = max m<k 2σ 2 (X i δ) 2 sign (X i δ), k = 2,..., n. i=m+1 Standardization necessary to obtain distributional results (cf. critical values). Q kn = 1 n Q k Q kk = 1 k Q k. Q kn may perform poorly for early changes. 27/48
28 Likelihood ratio: windowed versions (Willsky - Jones) Q k = max m<k 1 2σ 2 k i=m+1 (X i δ) 2 sign (X i δ), k = 2,..., n. The above expression (standardized or not) needs to be maximized once or twice with respect to m. Problems: maximization may be time consuming asymptotical distributional results may be hard to obtain Q k,g = 1 G max k G m<k 1 2σ 2 Q simp k,g = 1 1 k G 2σ 2 (X i δ) 2 sign (X i δ). i=k G+1 28/48 k i=m+1 (X i δ) 2 sign (X i δ),
29 Isotonic regression Chang and Fricker (J. Qual. Technology 31 (1999)) : isotonic regression method for Persistent monotone threshold crossing alternative M k = = k k (X i Z i ) 2 (X i Y i ) 2 i=1 i=1 0 for Y k δ k (X i δ) 2 k (X i Y i ) 2 for Y k > δ, i=j i=j where Z 1,..., Z k denotes an isotonic regression restricted to increase up to at most δ Y 1,..., Y k is the corresponding unrestricted isotonic regression J = min{i : Y i > δ}. 29/48
30 Intermezzo: scale alternatives X i = µ + σ i e i, i = 1,..., n, where µ is a known constant, σ i, i = 1,..., n are unknown positive parameters, e 1,..., e n are i.i.d. random variables with N(0, 1) distribution. Hypotheses: H 0 : σ i δ, i = 1,..., against H 1 : there is an m such that σ i δ, i = 1,..., σ i > δ i = m+1,..., n, The estimator σ i0 2 (H 0) of σi 2 under the restriction that σi 2 of the minimization problem: { min log σ σi 2 i δ2 σi 2 (X i µ) 2}. δ 2 is the solution 30/48
31 Scale alternatives: continued σ i0 2 = min((x i µ) 2, δ 2 ) The estimator σ i1 2 of σ i 2 under the restriction that σi 2 σ 2 i1 = max((x i µ) 2, δ 2 ). Hence, the loglikelihood ratio for the i-th observation is ( ( log δ 2 ) (X i µ) (X i µ) 2 δ 2 ) sign δ 2 is ( (Xi µ) 2 ) δ 2 1. q(t) = ( log(t 2 ) 1 + t 2 )sign (t 1), t > 0 is nondecreasing Q σ 2 k = max 0 j k k i= j+1 ( log ( (X i µ) 2 ) δ 2 1+ (X i µ) 2 ) δ 2 sign ((X i µ) 2 /δ 2 1). Clearly, the test statistics Q σ 2 k are similar to the Q k s. 31/48
32 Critical values The most common way to control false alarm in SPC is to use run lengths, in particular the in-control ARL (average run length). Sometimes also the SRL (standard deviation of the run length) is taken into account, or even better one uses quantiles. We reject H 0 as soon as for some k n T k > c n,α, where T k denotes any of the statistics presented in this talk and c N,α is chosen in such a way the significance level is α, i.e. ( ) P H0 max T k > c n,α = α. 1 k n Since the distribution of our statistics do not admit closed-form expression, these critical values have to be obtained by simulations or by asymptotics. 32/48
33 Asymptotics of critical values P H 0 (N C (c) > n C (α)) = ( )1 α for a prechosen α and a prechosen integer n C (α). We write c = c C (α) and rewrite the above relations as ( ) P H 0 max Q k < c C (α) 1 α, 1 k n C (α) so the probability of false alarm based on n C (α) observations is α. ( lim P H n 0 max Q k < c ) ( ) n(var H 1<k n 0 Z 1 (δ 0 )) 1/2 = P sup W (t) W (s) < c 0<s<t<1 where {W (t), t (0, 1)} is a Wiener process. This implies that lim n c C(α) n = c α ( VarH 0 Z 1 (δ 0 ) ) 1/2,, where c α is such that P ( ) sup 0<s<t<1 W (t) W (s) < c α = 1 α. 33/48
34 Asymptotics of critical values continued Hence, c C (α) can be approximated by c α n(varh 0 Z 1 (δ 0 )) 1/2 and as soon as n is large so does c C (α) and the procedure will perform poorly in detecting early changes. Alternative: weighted CUSUM N C,β (c n (β), n) = min{k 1; (k/n) β Q k nc n (β)}, β [0, 1/2], where c n (β) is determined in such a way that ( P H 0 max 1 k n (k/n) β Q k < ) nc n (β) 1 α. Hence, we reject H 0 as soon as there is a k such that Q k n c n (β)(k/n) β. Obviously, choosing β = 0 we have back the original procedure but for β (0, 1/2] larger weights are assigned to smaller k s. / department of mathematics and computer science 34/48
35 Asymptotics of critical values: weighted version lim c n(β) = d α (β) ( Var H n 0 Z 1 (δ 0 ) ) 1/2, β [0, 1/2), where d α (β) is determined as ( P sup 0<s<t<1 1 t β (W (t) W (s)) d α(β) ) = 1 α, where {W (t), t (0, 1)} is a Wiener process. It implies that c n (β)(var H 0 Z 1 (δ 0 )/n) 1/2 can be approximated by d β (α). For β = 1/2 we have for c n (1/2) the approximation d α,n (1/2)(Var H 0 Z 1 (δ 0 )) 1/2, where P ( sup 1 s t n ) 1 (W (t) W (s)) d α,n (1/2) = 1 α. t We can show that lim n c n (β) = O(1), β [0, 1/2), while c n (1/2) = O( log log n). / department of mathematics and computer science 35/48
36 Background of asymptotics I Komlós et al. (Z. Wahr. Verw. Gebiete 32 (1975) and 34 (1976)) Y 1,..., Y n i.i.d. with zero mean, unit variance and E X i 2+ < for some > 0. Then there exists a sequence of Wiener processes W n = {W n (t), t 0}, n n 0, such that, as n, ( ) and P max 1 k n k Y i W n (k) > x i=1 max 1 k n max 1 k n 1 k 1/(2+ ) C 1 nx (2+ ), x > 0 k Y i W n (k) = O P(1) i=1 k X i W n (k) = o P(n 1/(2+ ) ). i=1 36/48
37 Background of asymptotics II Erdös-Darling theorems and P(d 1 (log n) max 1<k n 1 k k Y i y + d 2 (log n)) exp{ exp{ y}}. i=1 P(d 1 (log n) max 1<k n 1 k k Y i y + d 2 (log n)) exp{ 2 exp{ y}}, i=1 where d 1 (t) = 2 log t, t > 1, d 2 (t) = 2 log t + 1 log log t log(π), log t > /48
38 Simulation of critical values (Ba- Sample quantiles ξ α are asymptotically normal with variance α(1 α) f 2 (ξ α ) hadur) We estimated f 2 (ξ α ) by applying a kernel density estimator to a pre-run of reasonable size. Example: for the 90% quantile the standard deviation equals ( ) /2, n because the the density at the quantile approximately equals Hence, if n = 10000, then the standard deviation of this quantile approximately equals /48
39 Performance simulation Simulation performed using R (see version on a Unix platform. Its new built-in procedure isoreg() is much faster in computing the isotonic regression (needed for simulations of the M n statistics) than the PAVA algorithm. Source codes of all procedures we used are available at klaster/compstat04 39/48
40 Simulation improvement Since Var N S E (N S ) 2, so increase in ARL leads to increase in variance of estimates. Variance reduction techniques have been developed for these simulations. Jun and Choi, Comm. Stat. A 22 (1993) present two techniques for CUSUM simulations: 1. hazard controlled estimator (following a suggestion of Ross to use the total hazard of a Markov chain) 2. ratio estimators using a cycle. 40/48
41 Simulation setup n = 100, α = 0.05, G = 0.1n = 10, σ = 1, 1000 repetitions. Two types of changes, each of which has a small and a large version: small version corresponds to change of mean of σ large version corresponds to change of mean of 3σ general change mean is uniformly distributed on [δ r, δ] before changepoint, and uniformly distributed on [δ, δ + r] after the change-point, where r = 1 (small change) or r = 3 (large change), respectively (so the expected mean increases by 1 and 3, respectively). gradual change mean increases on 5 equally spaced intervals before change point from µ to δ and on 5 equally spaced intervals after the change-point from δ to 2δ µ, where δ = 1 and µ = 0.5 (small change) or δ = 3 and µ = 1.5 (large change). 41/48
42 Simulation interpretation Detection delays must be interpreted with care. Suppose change at x and number of observations equals n, then in our simulations the number n x indicates either a detection at time point n or no detection at all Summary statistics of simulation include (run lengths are skewed!): mean standard deviation quantiles 42/48
43 Simulation results: general small change Q kn Q kk Q n,g Q simp n,g M n m mean; sd mean; sd mean; sd mean; sd mean; sd ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; /48
44 Gradual small change Q kn Q kk Q k,g Q simp k,g M k m mean; sd mean; sd mean; sd mean; sd mean; sd ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; /48
45 General large change Q kn Q kk Q k,g Q simp k,g M k m mean; sd mean; sd mean; sd mean; sd mean; sd ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; /48
46 Gradual large change Q kn Q kk Q k,g Q simp k,g M k m mean; sd mean; sd mean; sd mean; sd mean; sd ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; /48
47 Conclusions Chang and Fricker s M k has the smallest mean (both general and gradual case) Q kn has a considerably smaller standard deviation. M k has better low quantiles Q kn has better median and higher quantiles Summary: In most cases Q kn performs better, but M k sometimes detects extremely fast. 47/48
48 Future work other alternative hypotheses (in particular, changes of the variance) composite alternative hypotheses (asymptotical) distributional results combination of procedures application to real data... 48/48
The Shiryaev-Roberts Changepoint Detection Procedure in Retrospect - Theory and Practice
The Shiryaev-Roberts Changepoint Detection Procedure in Retrospect - Theory and Practice Department of Statistics The Hebrew University of Jerusalem Mount Scopus 91905 Jerusalem, Israel msmp@mscc.huji.ac.il
More informationSequential Detection. Changes: an overview. George V. Moustakides
Sequential Detection of Changes: an overview George V. Moustakides Outline Sequential hypothesis testing and Sequential detection of changes The Sequential Probability Ratio Test (SPRT) for optimum hypothesis
More informationSurveillance of Infectious Disease Data using Cumulative Sum Methods
Surveillance of Infectious Disease Data using Cumulative Sum Methods 1 Michael Höhle 2 Leonhard Held 1 1 Institute of Social and Preventive Medicine University of Zurich 2 Department of Statistics University
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 informationChange Detection Algorithms
5 Change Detection Algorithms In this chapter, we describe the simplest change detection algorithms. We consider a sequence of independent random variables (y k ) k with a probability density p (y) depending
More informationDesign and Implementation of CUSUM Exceedance Control Charts for Unknown Location
Design and Implementation of CUSUM Exceedance Control Charts for Unknown Location MARIEN A. GRAHAM Department of Statistics University of Pretoria South Africa marien.graham@up.ac.za S. CHAKRABORTI Department
More informationChapter 4. Theory of Tests. 4.1 Introduction
Chapter 4 Theory of Tests 4.1 Introduction Parametric model: (X, B X, P θ ), P θ P = {P θ θ Θ} where Θ = H 0 +H 1 X = K +A : K: critical region = rejection region / A: acceptance region A decision rule
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 informationThe University of Hong Kong Department of Statistics and Actuarial Science STAT2802 Statistical Models Tutorial Solutions Solutions to Problems 71-80
The University of Hong Kong Department of Statistics and Actuarial Science STAT2802 Statistical Models Tutorial Solutions Solutions to Problems 71-80 71. Decide in each case whether the hypothesis is simple
More informationSTAT 135 Lab 5 Bootstrapping and Hypothesis Testing
STAT 135 Lab 5 Bootstrapping and Hypothesis Testing Rebecca Barter March 2, 2015 The Bootstrap Bootstrap Suppose that we are interested in estimating a parameter θ from some population with members x 1,...,
More informationH 2 : otherwise. that is simply the proportion of the sample points below level x. For any fixed point x the law of large numbers gives that
Lecture 28 28.1 Kolmogorov-Smirnov test. Suppose that we have an i.i.d. sample X 1,..., X n with some unknown distribution and we would like to test the hypothesis that is equal to a particular distribution
More informationDetection and Diagnosis of Unknown Abrupt Changes Using CUSUM Multi-Chart Schemes
Sequential Analysis, 26: 225 249, 2007 Copyright Taylor & Francis Group, LLC ISSN: 0747-4946 print/532-476 online DOI: 0.080/0747494070404765 Detection Diagnosis of Unknown Abrupt Changes Using CUSUM Multi-Chart
More informationTesting Hypothesis. Maura Mezzetti. Department of Economics and Finance Università Tor Vergata
Maura Department of Economics and Finance Università Tor Vergata Hypothesis Testing Outline It is a mistake to confound strangeness with mystery Sherlock Holmes A Study in Scarlet Outline 1 The Power Function
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 informationAn Adaptive Exponentially Weighted Moving Average Control Chart for Monitoring Process Variances
An Adaptive Exponentially Weighted Moving Average Control Chart for Monitoring Process Variances Lianjie Shu Faculty of Business Administration University of Macau Taipa, Macau (ljshu@umac.mo) Abstract
More informationFirst Semester Dr. Abed Schokry SQC Chapter 9: Cumulative Sum and Exponential Weighted Moving Average Control Charts
Department of Industrial Engineering First Semester 2014-2015 Dr. Abed Schokry SQC Chapter 9: Cumulative Sum and Exponential Weighted Moving Average Control Charts Learning Outcomes After completing this
More informationStatistical Models and Algorithms for Real-Time Anomaly Detection Using Multi-Modal Data
Statistical Models and Algorithms for Real-Time Anomaly Detection Using Multi-Modal Data Taposh Banerjee University of Texas at San Antonio Joint work with Gene Whipps (US Army Research Laboratory) Prudhvi
More information40.530: Statistics. Professor Chen Zehua. Singapore University of Design and Technology
Singapore University of Design and Technology Lecture 9: Hypothesis testing, uniformly most powerful tests. The Neyman-Pearson framework Let P be the family of distributions of concern. The Neyman-Pearson
More informationHomework 7: Solutions. P3.1 from Lehmann, Romano, Testing Statistical Hypotheses.
Stat 300A Theory of Statistics Homework 7: Solutions Nikos Ignatiadis Due on November 28, 208 Solutions should be complete and concisely written. Please, use a separate sheet or set of sheets for each
More informationMaster s Written Examination - Solution
Master s Written Examination - Solution Spring 204 Problem Stat 40 Suppose X and X 2 have the joint pdf f X,X 2 (x, x 2 ) = 2e (x +x 2 ), 0 < x < x 2
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 informationDecentralized Sequential Hypothesis Testing. Change Detection
Decentralized Sequential Hypothesis Testing & Change Detection Giorgos Fellouris, Columbia University, NY, USA George V. Moustakides, University of Patras, Greece Outline Sequential hypothesis testing
More informationMaster s Written Examination
Master s Written Examination Option: Statistics and Probability Spring 016 Full points may be obtained for correct answers to eight questions. Each numbered question which may have several parts is worth
More informationStatistical Inference
Statistical Inference Robert L. Wolpert Institute of Statistics and Decision Sciences Duke University, Durham, NC, USA Spring, 2006 1. DeGroot 1973 In (DeGroot 1973), Morrie DeGroot considers testing the
More informationsimple if it completely specifies the density of x
3. Hypothesis Testing Pure significance tests Data x = (x 1,..., x n ) from f(x, θ) Hypothesis H 0 : restricts f(x, θ) Are the data consistent with H 0? H 0 is called the null hypothesis simple if it completely
More informationUncertainty. Jayakrishnan Unnikrishnan. CSL June PhD Defense ECE Department
Decision-Making under Statistical Uncertainty Jayakrishnan Unnikrishnan PhD Defense ECE Department University of Illinois at Urbana-Champaign CSL 141 12 June 2010 Statistical Decision-Making Relevant in
More informationAN ANALYSIS OF SHEWHART QUALITY CONTROL CHARTS TO MONITOR BOTH THE MEAN AND VARIABILITY KEITH JACOB BARRS. A Research Project Submitted to the Faculty
AN ANALYSIS OF SHEWHART QUALITY CONTROL CHARTS TO MONITOR BOTH THE MEAN AND VARIABILITY by KEITH JACOB BARRS A Research Project Submitted to the Faculty of the College of Graduate Studies at Georgia Southern
More informationGOTEBORG UNIVERSITY. Department of Statistics ON PERFORMANCE OF METHODS FOR STATISTICAL SURVEILLANCE
GOTEBORG UNIVERSITY Department of Statistics RESEARCH REPORT 1994:7 ISSN 0349-8034 ON PERFORMANCE OF METHODS FOR STATISTICAL SURVEILLANCE by o Goran Akermo Statistiska institutionen Goteborgs Universitet
More informationQuality Control & Statistical Process Control (SPC)
Quality Control & Statistical Process Control (SPC) DR. RON FRICKER PROFESSOR & HEAD, DEPARTMENT OF STATISTICS DATAWORKS CONFERENCE, MARCH 22, 2018 Agenda Some Terminology & Background SPC Methods & Philosophy
More informationHypothesis Testing. 1 Definitions of test statistics. CB: chapter 8; section 10.3
Hypothesis Testing CB: chapter 8; section 0.3 Hypothesis: statement about an unknown population parameter Examples: The average age of males in Sweden is 7. (statement about population mean) The lowest
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 informationDistribution-Free Monitoring of Univariate Processes. Peihua Qiu 1 and Zhonghua Li 1,2. Abstract
Distribution-Free Monitoring of Univariate Processes Peihua Qiu 1 and Zhonghua Li 1,2 1 School of Statistics, University of Minnesota, USA 2 LPMC and Department of Statistics, Nankai University, China
More informationSTA 732: Inference. Notes 2. Neyman-Pearsonian Classical Hypothesis Testing B&D 4
STA 73: Inference Notes. Neyman-Pearsonian Classical Hypothesis Testing B&D 4 1 Testing as a rule Fisher s quantification of extremeness of observed evidence clearly lacked rigorous mathematical interpretation.
More informationApplication of Variance Homogeneity Tests Under Violation of Normality Assumption
Application of Variance Homogeneity Tests Under Violation of Normality Assumption Alisa A. Gorbunova, Boris Yu. Lemeshko Novosibirsk State Technical University Novosibirsk, Russia e-mail: gorbunova.alisa@gmail.com
More informationSequential Change-Point Approach for Online Community Detection
Sequential Change-Point Approach for Online Community Detection Yao Xie Joint work with David Marangoni-Simonsen H. Milton Stewart School of Industrial and Systems Engineering Georgia Institute of Technology
More informationDirection: This test is worth 250 points and each problem worth points. DO ANY SIX
Term Test 3 December 5, 2003 Name Math 52 Student Number Direction: This test is worth 250 points and each problem worth 4 points DO ANY SIX PROBLEMS You are required to complete this test within 50 minutes
More informationLECTURE 10: NEYMAN-PEARSON LEMMA AND ASYMPTOTIC TESTING. The last equality is provided so this can look like a more familiar parametric test.
Economics 52 Econometrics Professor N.M. Kiefer LECTURE 1: NEYMAN-PEARSON LEMMA AND ASYMPTOTIC TESTING NEYMAN-PEARSON LEMMA: Lesson: Good tests are based on the likelihood ratio. The proof is easy in the
More informationSome General Types of Tests
Some General Types of Tests We may not be able to find a UMP or UMPU test in a given situation. In that case, we may use test of some general class of tests that often have good asymptotic properties.
More informationAn Investigation of Combinations of Multivariate Shewhart and MEWMA Control Charts for Monitoring the Mean Vector and Covariance Matrix
Technical Report Number 08-1 Department of Statistics Virginia Polytechnic Institute and State University, Blacksburg, Virginia January, 008 An Investigation of Combinations of Multivariate Shewhart and
More informationChapter 7. Hypothesis Testing
Chapter 7. Hypothesis Testing Joonpyo Kim June 24, 2017 Joonpyo Kim Ch7 June 24, 2017 1 / 63 Basic Concepts of Testing Suppose that our interest centers on a random variable X which has density function
More informationLecture 8: Information Theory and Statistics
Lecture 8: Information Theory and Statistics Part II: Hypothesis Testing and I-Hsiang Wang Department of Electrical Engineering National Taiwan University ihwang@ntu.edu.tw December 23, 2015 1 / 50 I-Hsiang
More informationLet us first identify some classes of hypotheses. simple versus simple. H 0 : θ = θ 0 versus H 1 : θ = θ 1. (1) one-sided
Let us first identify some classes of hypotheses. simple versus simple H 0 : θ = θ 0 versus H 1 : θ = θ 1. (1) one-sided H 0 : θ θ 0 versus H 1 : θ > θ 0. (2) two-sided; null on extremes H 0 : θ θ 1 or
More informationA Modified Poisson Exponentially Weighted Moving Average Chart Based on Improved Square Root Transformation
Thailand Statistician July 216; 14(2): 197-22 http://statassoc.or.th Contributed paper A Modified Poisson Exponentially Weighted Moving Average Chart Based on Improved Square Root Transformation Saowanit
More informationOptimum CUSUM Tests for Detecting Changes in Continuous Time Processes
Optimum CUSUM Tests for Detecting Changes in Continuous Time Processes George V. Moustakides INRIA, Rennes, France Outline The change detection problem Overview of existing results Lorden s criterion and
More informationDefinition 3.1 A statistical hypothesis is a statement about the unknown values of the parameters of the population distribution.
Hypothesis Testing Definition 3.1 A statistical hypothesis is a statement about the unknown values of the parameters of the population distribution. Suppose the family of population distributions is indexed
More informationSpring 2012 Math 541B Exam 1
Spring 2012 Math 541B Exam 1 1. A sample of size n is drawn without replacement from an urn containing N balls, m of which are red and N m are black; the balls are otherwise indistinguishable. Let X denote
More informationarxiv:math/ v2 [math.st] 15 May 2006
The Annals of Statistics 2006, Vol. 34, No. 1, 92 122 DOI: 10.1214/009053605000000859 c Institute of Mathematical Statistics, 2006 arxiv:math/0605322v2 [math.st] 15 May 2006 SEQUENTIAL CHANGE-POINT DETECTION
More informationOn detection of unit roots generalizing the classic Dickey-Fuller approach
On detection of unit roots generalizing the classic Dickey-Fuller approach A. Steland Ruhr-Universität Bochum Fakultät für Mathematik Building NA 3/71 D-4478 Bochum, Germany February 18, 25 1 Abstract
More informationFaculty of Science and Technology MASTER S THESIS
Faculty of Science and Technology MASTER S THESIS Study program/ Specialization: Spring semester, 20... Open / Restricted access Writer: Faculty supervisor: (Writer s signature) External supervisor(s):
More informationChapter 6. Hypothesis Tests Lecture 20: UMP tests and Neyman-Pearson lemma
Chapter 6. Hypothesis Tests Lecture 20: UMP tests and Neyman-Pearson lemma Theory of testing hypotheses X: a sample from a population P in P, a family of populations. Based on the observed X, we test a
More informationA Very Brief Summary of Statistical Inference, and Examples
A Very Brief Summary of Statistical Inference, and Examples Trinity Term 2008 Prof. Gesine Reinert 1 Data x = x 1, x 2,..., x n, realisations of random variables X 1, X 2,..., X n with distribution (model)
More informationSelf-Starting Control Chart for Simultaneously Monitoring Process Mean and Variance
International Journal of Production Research Vol. 00, No. 00, 15 March 2008, 1 14 Self-Starting Control Chart for Simultaneously Monitoring Process Mean and Variance Zhonghua Li a, Jiujun Zhang a,b and
More informationSEQUENTIAL CHANGE-POINT DETECTION WHEN THE PRE- AND POST-CHANGE PARAMETERS ARE UNKNOWN. Tze Leung Lai Haipeng Xing
SEQUENTIAL CHANGE-POINT DETECTION WHEN THE PRE- AND POST-CHANGE PARAMETERS ARE UNKNOWN By Tze Leung Lai Haipeng Xing Technical Report No. 2009-5 April 2009 Department of Statistics STANFORD UNIVERSITY
More informationOn ARL-unbiased c-charts for i.i.d. and INAR(1) Poisson counts
On ARL-unbiased c-charts for iid and INAR(1) Poisson counts Manuel Cabral Morais (1) with Sofia Paulino (2) and Sven Knoth (3) (1) Department of Mathematics & CEMAT IST, ULisboa, Portugal (2) IST, ULisboa,
More informationCHAPTER 8. Test Procedures is a rule, based on sample data, for deciding whether to reject H 0 and contains:
CHAPTER 8 Test of Hypotheses Based on a Single Sample Hypothesis testing is the method that decide which of two contradictory claims about the parameter is correct. Here the parameters of interest are
More informationTesting Statistical Hypotheses
E.L. Lehmann Joseph P. Romano Testing Statistical Hypotheses Third Edition 4y Springer Preface vii I Small-Sample Theory 1 1 The General Decision Problem 3 1.1 Statistical Inference and Statistical Decisions
More informationSummary of Chapters 7-9
Summary of Chapters 7-9 Chapter 7. Interval Estimation 7.2. Confidence Intervals for Difference of Two Means Let X 1,, X n and Y 1, Y 2,, Y m be two independent random samples of sizes n and m from two
More informationLecture 7 Introduction to Statistical Decision Theory
Lecture 7 Introduction to Statistical Decision Theory I-Hsiang Wang Department of Electrical Engineering National Taiwan University ihwang@ntu.edu.tw December 20, 2016 1 / 55 I-Hsiang Wang IT Lecture 7
More informationCentral Limit Theorem ( 5.3)
Central Limit Theorem ( 5.3) Let X 1, X 2,... be a sequence of independent random variables, each having n mean µ and variance σ 2. Then the distribution of the partial sum S n = X i i=1 becomes approximately
More informationA Note on Data-Adaptive Bandwidth Selection for Sequential Kernel Smoothers
6th St.Petersburg Workshop on Simulation (2009) 1-3 A Note on Data-Adaptive Bandwidth Selection for Sequential Kernel Smoothers Ansgar Steland 1 Abstract Sequential kernel smoothers form a class of procedures
More informationEFFICIENT CHANGE DETECTION METHODS FOR BIO AND HEALTHCARE SURVEILLANCE
EFFICIENT CHANGE DETECTION METHODS FOR BIO AND HEALTHCARE SURVEILLANCE A Thesis Presented to The Academic Faculty by Sung Won Han In Partial Fulfillment of the Requirements for the Degree Doctor of Philosophy
More informationLecture Notes in Statistics 180 Edited by P. Bickel, P. Diggle, S. Fienberg, U. Gather, I. Olkin, S. Zeger
Lecture Notes in Statistics 180 Edited by P. Bickel, P. Diggle, S. Fienberg, U. Gather, I. Olkin, S. Zeger Yanhong Wu Inference for Change-Point and Post-Change Means After a CUSUM Test Yanhong Wu Department
More informationLecture 5: Likelihood ratio tests, Neyman-Pearson detectors, ROC curves, and sufficient statistics. 1 Executive summary
ECE 830 Spring 207 Instructor: R. Willett Lecture 5: Likelihood ratio tests, Neyman-Pearson detectors, ROC curves, and sufficient statistics Executive summary In the last lecture we saw that the likelihood
More informationAsymptotically Optimal Quickest Change Detection in Distributed Sensor Systems
This article was downloaded by: [University of Illinois at Urbana-Champaign] On: 23 May 2012, At: 16:03 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954
More informationSTAT 830 Hypothesis Testing
STAT 830 Hypothesis Testing Hypothesis testing is a statistical problem where you must choose, on the basis of data X, between two alternatives. We formalize this as the problem of choosing between two
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 informationLeast Favorable Distributions for Robust Quickest Change Detection
Least Favorable Distributions for Robust Quickest hange Detection Jayakrishnan Unnikrishnan, Venugopal V. Veeravalli, Sean Meyn Department of Electrical and omputer Engineering, and oordinated Science
More informationLecture 8: Information Theory and Statistics
Lecture 8: Information Theory and Statistics Part II: Hypothesis Testing and Estimation I-Hsiang Wang Department of Electrical Engineering National Taiwan University ihwang@ntu.edu.tw December 22, 2015
More informationThree Non-parametric Control Charts Based on Ranked Set Sampling Waliporn Tapang, Adisak Pongpullponsak* and Sukuman Sarikavanij
Chiang Mai J. Sci. 16; 43(4) : 914-99 http://epg.science.cmu.ac.th/ejournal/ Contributed Paper Three Non-parametric s Based on Ranked Set Sampling Waliporn Tapang, Adisak Pongpullponsak* and Sukuman Sarikavanij
More informationThe Multinomial Model
The Multinomial Model STA 312: Fall 2012 Contents 1 Multinomial Coefficients 1 2 Multinomial Distribution 2 3 Estimation 4 4 Hypothesis tests 8 5 Power 17 1 Multinomial Coefficients Multinomial coefficient
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 informationSection II: Assessing Chart Performance. (Jim Benneyan)
Section II: Assessing Chart Performance (Jim Benneyan) 1 Learning Objectives Understand concepts of chart performance Two types of errors o Type 1: Call an in-control process out-of-control o Type 2: Call
More informationAsymptotic Statistics-VI. Changliang Zou
Asymptotic Statistics-VI Changliang Zou Kolmogorov-Smirnov distance Example (Kolmogorov-Smirnov confidence intervals) We know given α (0, 1), there is a well-defined d = d α,n such that, for any continuous
More informationarxiv: v2 [eess.sp] 20 Nov 2017
Distributed Change Detection Based on Average Consensus Qinghua Liu and Yao Xie November, 2017 arxiv:1710.10378v2 [eess.sp] 20 Nov 2017 Abstract Distributed change-point detection has been a fundamental
More informationDA Freedman Notes on the MLE Fall 2003
DA Freedman Notes on the MLE Fall 2003 The object here is to provide a sketch of the theory of the MLE. Rigorous presentations can be found in the references cited below. Calculus. Let f be a smooth, scalar
More informationMonitoring binary outcomes using risk-adjusted charts: a comparative study
Research Article Received 25 March 20, Accepted May 20 Published online 22 July 20 in Wiley Online Library (wileyonlinelibrary.com) DOI: 0.002/sim.4305 Monitoring binary outcomes using risk-adjusted charts:
More informationRegenerative Likelihood Ratio control schemes
Regenerative Likelihood Ratio control schemes Emmanuel Yashchin IBM Research, Yorktown Heights, NY XIth Intl. Workshop on Intelligent Statistical Quality Control 2013, Sydney, Australia Outline Motivation
More informationConfidence Regions For The Ratio Of Two Percentiles
Confidence Regions For The Ratio Of Two Percentiles Richard Johnson Joint work with Li-Fei Huang and Songyong Sim January 28, 2009 OUTLINE Introduction Exact sampling results normal linear model case Other
More informationA BINARY CONTROL CHART TO DETECT SMALL JUMPS
A BINARY CONTROL CHART TO DETECT SMALL JUMPS Ewaryst Rafaj lowicz Institute of Computer Engineering, Control and Robotics Wroc law University of Technology, Poland Ansgar Steland 1 Institute of Statistics
More informationEXTENSIONS TO THE EXPONENTIALLY WEIGHTED MOVING AVERAGE PROCEDURES
EXTENSIONS TO THE EXPONENTIALLY WEIGHTED MOVING AVERAGE PROCEDURES CHAO AN-KUO NATIONAL UNIVERSITY OF SINGAPORE 2016 EXTENSIONS TO THE EXPONENTIALLY WEIGHTED MOVING AVERAGE PROCEDURES CHAO AN-KUO (MS,
More informationBTRY 4090: Spring 2009 Theory of Statistics
BTRY 4090: Spring 2009 Theory of Statistics Guozhang Wang September 25, 2010 1 Review of Probability We begin with a real example of using probability to solve computationally intensive (or infeasible)
More informationLecture 12 November 3
STATS 300A: Theory of Statistics Fall 2015 Lecture 12 November 3 Lecturer: Lester Mackey Scribe: Jae Hyuck Park, Christian Fong Warning: These notes may contain factual and/or typographic errors. 12.1
More informationA Study on the Power Functions of the Shewhart X Chart via Monte Carlo Simulation
A Study on the Power Functions of the Shewhart X Chart via Monte Carlo Simulation M.B.C. Khoo Abstract The Shewhart X control chart is used to monitor shifts in the process mean. However, it is less sensitive
More informationChapter 9: Hypothesis Testing Sections
Chapter 9: Hypothesis Testing Sections 9.1 Problems of Testing Hypotheses 9.2 Testing Simple Hypotheses 9.3 Uniformly Most Powerful Tests Skip: 9.4 Two-Sided Alternatives 9.6 Comparing the Means of Two
More information557: MATHEMATICAL STATISTICS II HYPOTHESIS TESTING: EXAMPLES
557: MATHEMATICAL STATISTICS II HYPOTHESIS TESTING: EXAMPLES Example Suppose that X,..., X n N, ). To test H 0 : 0 H : the most powerful test at level α is based on the statistic λx) f π) X x ) n/ exp
More informationModified Simes Critical Values Under Positive Dependence
Modified Simes Critical Values Under Positive Dependence Gengqian Cai, Sanat K. Sarkar Clinical Pharmacology Statistics & Programming, BDS, GlaxoSmithKline Statistics Department, Temple University, Philadelphia
More informationEfficient scalable schemes for monitoring a large number of data streams
Biometrika (2010, 97, 2,pp. 419 433 C 2010 Biometrika Trust Printed in Great Britain doi: 10.1093/biomet/asq010 Advance Access publication 16 April 2010 Efficient scalable schemes for monitoring a large
More informationImproving Biosurveillance System Performance. Ronald D. Fricker, Jr. Virginia Tech July 7, 2015
Improving Biosurveillance System Performance Ronald D. Fricker, Jr. Virginia Tech July 7, 2015 What is Biosurveillance? The term biosurveillance means the process of active datagathering of biosphere data
More informationApplying the Benjamini Hochberg procedure to a set of generalized p-values
U.U.D.M. Report 20:22 Applying the Benjamini Hochberg procedure to a set of generalized p-values Fredrik Jonsson Department of Mathematics Uppsala University Applying the Benjamini Hochberg procedure
More informationOptimal Cusum Control Chart for Censored Reliability Data with Log-logistic Distribution
CMST 21(4) 221-227 (2015) DOI:10.12921/cmst.2015.21.04.006 Optimal Cusum Control Chart for Censored Reliability Data with Log-logistic Distribution B. Sadeghpour Gildeh, M. Taghizadeh Ashkavaey Department
More informationData-Efficient Quickest Change Detection
Data-Efficient Quickest Change Detection Venu Veeravalli ECE Department & Coordinated Science Lab University of Illinois at Urbana-Champaign http://www.ifp.illinois.edu/~vvv (joint work with Taposh Banerjee)
More informationOn robust stopping times for detecting changes in distribution
On robust stopping times for detecting changes in distribution by Yuri Golubev, Mher Safarian No. 116 MAY 2018 WORKING PAPER SERIES IN ECONOMICS KIT Die Forschungsuniversität in der Helmholtz-Gemeinschaft
More informationTHE CUSUM MEDIAN CHART FOR KNOWN AND ESTIMATED PARAMETERS
THE CUSUM MEDIAN CHART FOR KNOWN AND ESTIMATED PARAMETERS Authors: Philippe Castagliola Université de Nantes & LS2N UMR CNRS 6004, Nantes, France (philippe.castagliola@univ-nantes.fr) Fernanda Otilia Figueiredo
More informationA CUSUM approach for online change-point detection on curve sequences
ESANN 22 proceedings, European Symposium on Artificial Neural Networks, Computational Intelligence and Machine Learning. Bruges Belgium, 25-27 April 22, i6doc.com publ., ISBN 978-2-8749-49-. Available
More informationCorrection factors for Shewhart and control charts to achieve desired unconditional ARL
International Journal of Production Research ISSN: 0020-7543 (Print) 1366-588X (Online) Journal homepage: http://www.tandfonline.com/loi/tprs20 Correction factors for Shewhart and control charts to achieve
More informationThe Changepoint Model for Statistical Process Control
The Changepoint Model for Statistical Process Control DOUGLAS M. HAWKINS and PEIHUA QIU University of Minnesota, Minneapolis, MN 55455 CHANG WOOK KANG Hanyang University, Seoul, Korea Statistical process
More informationMathematical Statistics
Mathematical Statistics MAS 713 Chapter 8 Previous lecture: 1 Bayesian Inference 2 Decision theory 3 Bayesian Vs. Frequentist 4 Loss functions 5 Conjugate priors Any questions? Mathematical Statistics
More informationMathematics Ph.D. Qualifying Examination Stat Probability, January 2018
Mathematics Ph.D. Qualifying Examination Stat 52800 Probability, January 2018 NOTE: Answers all questions completely. Justify every step. Time allowed: 3 hours. 1. Let X 1,..., X n be a random sample from
More informationHypothesis testing (cont d)
Hypothesis testing (cont d) Ulrich Heintz Brown University 4/12/2016 Ulrich Heintz - PHYS 1560 Lecture 11 1 Hypothesis testing Is our hypothesis about the fundamental physics correct? We will not be able
More informationLecture 17: The Exponential and Some Related Distributions
Lecture 7: The Exponential and Some Related Distributions. Definition Definition: A continuous random variable X is said to have the exponential distribution with parameter if the density of X is e x if
More information