Two applications of max stable distributions: Random sketches and Heavy tail exponent estimation
|
|
- Cecil Blair
- 6 years ago
- Views:
Transcription
1 Two applications of max stable distributions: Random sketches and Heavy tail exponent estimation Stilian Stoev Department of Statistics University of Michigan Boston University Probability Seminar March 1, 2007 Based on joint works with Marios Hadjieleftheriou, George Kollios, George Michailidis and Murad Taqqu 1
2 Outline: Max stable distributions: background Max stable sketches: norm, distance and dominance norms Max spectrum: heavy tail exponent estimation 2
3 Max stable distributions Th (Fisher-Tippett-Gnedenko) Let X 1,..., X n,... be iid. If 1 a n (X 1 X n ) b n d Z, (n ) (1) then the non degenerate Z has one of the three types of distributions: (i) Frechet: G α (x) := exp{ x α }, x > 0, with α > 0. (ii) Negative Frechet: Ψ α (x) := exp{ ( x) α }, x < 0, with α > 0. (iii) Gumbel: Λ(x) := exp{ e x }, x R. Analog of the CLT when replaces +. Roughly: the tails of X i determine the type of limit. Def A rv X is max stable if, a, b > 0, c > 0, d: where X, X, X iid. ax bx d = cx + d, Th The limit Z in (1) are max stable. Conversely any max stable Z appears in a limit of iid maxima as in (1). 3
4 More on Frechet distributions Def Z is α Fréchet (α > 0) if P{Z x} = exp{ cx α }, x > 0, (c > 0). Properties scale family: For all a > 0, P{aZ x} = exp{ (a α c)x α }, x > 0 Notation: For the scale c 1/α : Z α := c 1/α, so az α = a Z α, a > 0. heavy tails: As x, P{Z > x} = 1 e cx α cx α. moments: EZ p < iff p < α and EZ p = 0 x p de cx α = c p Γ(1 p/α). Max stability: For iid Z, Z 1,..., Z n : Indeed: Z 1 Z n d = n 1/α Z. P{Z i x, i = 1,..., n} = P{Z x} n = e ncx α, x > 0. 4
5 Application I: Random Sketches First, mathematical concepts, then explanation & use. 5
6 Max stable sketches Signal : f(i) 0, i = 1,..., n, with large n. Def For α > 0, and k N, let n E j (f) := f(i)z j (i), 1 j k, i=1 for iid standard α Fréchet Z j (i) s P{Z j (i) x} = exp{ x α }, x > 0 i.e. Z j (i) α = 1. Max stable sketch of f: E j (f), 1 j k. Idea: use {E j (f)} with k << n as a proxy to f! More on this later... Properties E j (f) are α Fréchet with n E j (f) α = ( f(i) α ) 1/α =: f lα. Indeed, for x > 0: i=1 P{E j (f) x} = i P{Z x/f(i)} = exp{ ( i f(i) α )x α }. max linearity: For a, b > 0 and another signal g(i) E j (af bg) = ae j (f) be j (g), 1 j k, where (af bg)(i) = af(i) bg(i). 6
7 Max stable sketches: estimation Observations The E j (f) s are iid in j = 1,..., k The scale E j (f) α equals the norm f lα. Estimators: (moment) For 0 < p < α, define L p (f) := c k p,α E j (f) p k with c p,α := Γ(1 p/α) 1. j=1 Motivation: 1 E j (f) p EE j (f) p = Γ(1 p/α) E j (f) p k α. j (median) Define, Motivation: 1/p M(f) := ln(2) 1/α med{e j (f), 1 j k}. e xα = 1/2 x = ln(2) 1/α., 7
8 Norms and distances Given the max stable sketch {E j (f)} of a signal f. For large k s, we can estimate the norm: f lα L p (f) or f lα M(f). Performance guarantees later... Given another signal g, define the distance: ρ α (f, g) := n f(i) α g(i) α. i=1 Note that a b = 2(a b) a b, a, b 0, so: ρ α (f, g) = 2 i f(i) α g(i) α i f(i) α i g(i) α, and hence: ρ α (f, g) = 2 f g α l α f α l α g α l α. Thus, given the sketches E j (f) and E j (g), we can estimate ρ α (f, g) by or by ρ α (f, g) 2L p (f g) α L p (f) α L p (g) α ρ α (f, g) 2M(f g) α M(f) α M(g) α The key is: E j (f g) = E j (f) E j (g) so E j (f g) is available from E j (f) and E j (g). 8
9 Dominance norms Given are several signals f t (i), 1 i n, indexed by t = 1,..., T. Def The dominant signal is: f (i) := max 1 t T f t(i), 1 i n. The dominance norm of the f t s is f lα. Problem: Estimate f lα from the sketches {E j (f t )}, t = 1,..., T. Without accessing the signals f t directly! Solution: Max linearity implies E j (f ) = E j (f 1 ) E j (f T ), 1 j k. The sketch E j (f ) is readily available! Proceed as above and estimate f lα by: f lα L p (f ) or by f lα M(f ). 9
10 Motivation: Data streams Signals which cannot be stored and processed with conventional methods are emerging in: Telephone call monitoring (unmanageably large data bases). Internet traffic monitoring (vast, rapidly growing amounts of data). Sensing applications (restrictions on power, bandwidth, etc.). Online transactions (Banking, Stock market data, E-commerce, etc.). The streaming model a general framework. Let f(i), 1 i n be a signal with large domain size n. Starting with f(i) = 0, update f sequentially in time: cash register Observe (i, a(i)), and update f: f(i) := f(i) + a(i) aggregated Observe (i, f(i)) directly. 10
11 Examples: Call detail records & IP monitoring Phone calls: Each time a phone call ends, a CDR packet is sent to a main station. Think of: CDR = (i, a(i), ), with i = (phone no.) a(i) = time used. Given a CDR (i, a(i), update f(i) := f(i) + a(i). So the signal f(i) contains info on all phone numbers! Multiple signals: f t (i), t = 1,...,7 capture the calling patterns for different days of the week. IP monitoring: Monitor a fast Internet link during a given day (or hour). How many distinct IPs used it? How many bytes were transmitted from a given range of IPs? Can you give me estimates in real time??? Signal: f(i) = Bytes transmitted by IP i. Stream items: A packed from IP i contains info: (i, b(i)), i = (IP number), b(i) = (Bytes) On observing (i, b(i)), update: f(i) := f(i) + b(i). Impossible to store and process all data in real time. 11
12 Classical (additive) random sketches Instead of storing the signal f(i), 1 i n, keep: n S j (f) = f(i)x j (i), 1 j k, i=1 where k << n and where X j (i) are iid random variables. The S j (f) s 1 j k is called the random sketch of f. Important: Can generate X j (i) s on the fly from small O(klog 2 (n)) set of seeds in main memory. Can update the sketch sequentially in the most general, cash register model. Can approximate norms and inner products! If Var(X j (i)) = σ 2, then Var(S j (f)) = f 2 l 2 σ 2, f 2 l 2 = i f(i)2. So f 2 l 2 1 K S j (f) 2. For two signals f and g, by linearity, j f g 2 l 2 1 k k S j (f) S j (g) 2, j=1 12
13 Back to max stable sketches Strengths 1. Max stable sketches are max linear. Natural use in: sensor data, Internet traffic monitoring, Transaction record tracking, max sequential updates. 2. Work for any α > 0 and not only for 0 < α 2 (unlike the sum stable sketches). 3. Tailor made for: dominance norm estimation. Illustration: Signals: f t (i) = Bytes of IP number i during t th hour of monitoring. Large domain: 1 i n where n = Problem: Estimate the norm of the heaviest load scenario during the day, i.e. for f (i) = max 1 t 24 f t(i) estimate the dominance norm f lα? Problem: Compare several network traffic scenarios f, g, h,... without accessing the enormous data directly. Compute pair wise distances, norms and dominance norms. Max stable sketches provide efficient solutions. Limitation: Max stable sketches are not linear cannot be updated in the linear cash register format! 13
14 Max stable sketches: Performance guarantees Norms Given a max stable sketch E j (f), j = 1,..., k of a signal f(i) 0. Recall that M(f) is the median based estimator of f lα. Th 1(SHKT 07) Let ǫ, δ (0,1). Then, P{ M(f)/ f lα 1 ǫ} 1 δ, if k C log(1/δ)/ǫ 2, for some C > 0. For optimal value of the constant, as ǫ, δ 0, we have k = k(ǫ, δ) 2 α 2 ln 2 (2) ln(1/δ) ǫ 2. Distances Let now g be another signal and let D(f, g) := 2M(f g) α M(f) M(g) ρ α (f, g) = f α g α l1. Th 4 (SHKT 07) Let ǫ, δ, η (0,1). If ρ α (f) η f g α l α, then P{ D(f, g)/ρ α (f, g) 1 ǫ} 1 δ, for k C log(1/δ)/ǫ 2, with some C > 0. 14
15 Sketches: summary & credits Some seminal papers & algorithmic applications Alon, Matias & Szegedy (1996) Feingenbaum, Kannanm Strauss & Viswanathan (1999) Fong & Strauss (2000) Indyk (2000) Gilbert, Kotidis, Muthukrishnan & Strauss (2001, 2002) Cormode & Muthukrishnan (2003) Datar, Immorlica, Indyk & Mirrokni (2004) Review paper: Muthukrishnan (2003) More on max stable sketches: Stoev, Hadjieleftheriou, Kollios & Taqqu (2007) Code: By Marios Hadjieleftheriou (AT&T Research) A lesson: (For probabilists and computer scientists) Study each other s fields! 15
16 Application II: Heavy tail exponents 16
17 Heavy tailed data A random variable X is said to be heavy tailed if P{ X x} L(x)x α, as x, for some α > 0 and a slowly varying function L. Here we focus on the simpler but important context: X 0, a.s. and P{X > x} Cx α, as x. X (infinite moments) For p > 0, EX p < if and only if p < α. In particular, and 0 < α 2 Var(X) = 0 < α 1 E X =. The estimation of the heavy tail exponent α is an important problem with rich history. Why do we need heavy tail models? Every finite sample X 1,..., X n has finite sample mean, variance and all sample moments! Why consider heavy tailed models in practice?! 17
18 Why use heavy tailed models? All models are wrong, but some are useful. George Box Let F and G be any two distributions with positive densities on (0, ). Let ǫ > 0 and x 1,..., x n (0, ) be arbitrary, then both: and P F {X i (x i ǫ, x i + ǫ), i = 1,..., n} > 0 are positive! P G {X i (x i ǫ, x i + ǫ), i = 1,..., n} > 0 For a given sample, very many models apply. The ones that continue to work as the sample grows are most suitable. We next present real data sets of Financial, Insurance and Internet data. They can be very heavy tailed. 18
19 Traded volumes on the Intel stock 10 x Traded Volumes No. Stocks, INTC, Nov 1, x x
20 Insurance claims due to fire loss 250 Danish Fire Loss Data: Hill plot: α (k) = H order statistics Max Spectrum H= ( ), α = Scales j 20
21 TCP flow sizes (in number of packets) x 10 4 TCP Flow Sizes (packets): UNC link 2001 (~ 36 min) time x 10 4 The first minute
22 History Hill (1975) worked out the MLE in the Pareto model P{X > x} = x α, x 1 and introduced the Hill plot: α H (k) := ( 1 k k log(x i,n ) log(x k+1,n )) 1, i=1 where X 1,n X 2,n X k,n are the top K order statistics of the sample. How to choose k? pick k where the plot of α H (k) vs. k stabilizes. serious problems in practice: volatile & hard to interpret: Hill horror plot confidence intervals robustness Consistency and asymptotic normality resolved: Weissman (1978), Hall (1982) in semi parametric setting. Many other estimators: kernel based Csörgő, Deheuvels and Mason (1985), moment Dekkers, Einmahl and de Haan (1989), among many others. Most estimators exploit the top order statistics. 22
23 Hill horror plots: TCP flow sizes 2 Hill plot 1.5 α Order statistics k x 10 4 α H (300) = α H (12000) = α 1 α Order statistics k Order statistics k x
24 Fréchet max stable laws Consider i.i.d. X i s with P{X 1 > x} Cx α, x. By the Fisher Tippett Gnedenko theorem: 1 max n X 1/α i 1 n d X 1 i n n 1/α i Z, as n, i=1 where Z is α Fréchet extreme value distribution: P{Z x} = exp{ Cx α }, x > 0. The extreme value distributions are max stable. In particular, for i.i.d. α Fréchet Z,& Z i s: Z 1 Z n d = n 1/α Z. A time series of i.i.d. α Fréchet {Z k } is 1/α max-selfsimilar: m N : { 1 i m Z m(k 1)+i } k N d = m 1/α {Z k } k N. Block maxima of size m have the same distribution as the original data, rescaled by the factor m 1/α. Any heavy tailed {X k } (i.i.d.) data set is asymptotically max self similar. 24
25 Max spectrum Given a positive sample X 1,..., X n, consider the dyadic block maxima: D(j, k) := max 1 i 2 j X 2j (k 1)+i with k = 1,..., n j = [n/2 j ]. In view of the asymptotic scaling: observe that Y j := 1 n j 1 d 2j/αD(j, k) Z, n j j=1 where c := Elog 2 Z. 2 j i=1 (j ), X 2j (k 1)+i, log 2 D(j, k) j/α + c, (j, n j ) The last asymptotics follow from the LLN since D(j, k) s are independent in k. The max spectrum of the data is defined as the statistics: Y j, j = 1,...,[log 2 (n)]. Can identify α from the slope of the Y j s vs. j, for large j s. 25
26 Max spectrum estimators of α Given the max spectrum Y j, j = 1,...,[log 2 (n)], define j 2 α(j 1, j 2 ) := ( w j Y j ) 1, j=j 1 where j jw j = 1 and j w j = 0. That is, use linear regression to estimate the slope 1/α. Issues: weighted or generalized least squares must be used, since Var(Y j ) 1/n j 2 j, The choice of the scales j 1 and j 2 is critical. These issues and confidence intervals, are addressed in Stoev, Michailidis and Taqqu (2006). 26
27 Examples of max spectra Frechet(α = 1.5), Estimated α = Intel Volume: α = Max Spectrum 10 5 Max Spectrum Scales j 1.3 stable(β=1), α = Scales j T dist (df = 3), α = Max Spectrum Max Spectrum Scales j Scales j 27
28 Robustness of the max spectrum 20 Max self similarity: α(5,10) = , α(10,12) = Max Spectrum Scales j Compare with the Hill horror plot of the Internet flow size data (next slide). 28
29 Hill horror plots: TCP flow sizes 2 Hill plot 1.5 α Order statistics k x 10 4 α H (300) = α H (12000) = α 1 α Order statistics k Order statistics k x 10 4 Compare with the max spectrum plot of the Internet flow size data (previous slide). 29
30 Asymptotic normality of the max spectrum Let P{X k x} =: F(x), where 1 F(x) = Cx α (1 + O(x β )), (x ), (2) with some β > 0. Consider the max spectrum {Y j } for the range of scales: r(n) + j 1 j r(n) + j 2, with fixed j 1 < j 2 and r(n), n. Theorem Under (2) and another technical condition, P{ n j2 +r(( θ, Y r ) ( θ, µ r )) x} Φ(x/σ θ ) sup x R C θ (1/2 rβ/α + r2 r/2 / n), (3) where Φ is the standard Normal c.d.f. Here Y r = (Y r+j ) j 2 j=j 1 and and ( θ, Y r ) = j 2 j=j 1 θ j Y r+j, σ 2 θ = ( θ,σ j 2 αθ) := µ r (j) = (r + j)/α + c i,j=j 1 θ i Σ α (i, j)θ j > 0. 30
31 The covariance structure of the max spectrum The asymptotic covariance matrix Σ is the same as if the data were i.i.d. α Fréchet! The intuition is that the block maxima {D(r + j, k), k = 1,..., n r+j } behave like i.i.d. α Fréchet variables, as r = r(n). The covariance entries are given by: with Σ α (i, j) = α 2 2 i j j 2ψ( i j ), ψ(a) := Cov(log 2 (Z 1 ),log 2 (Z 1 (2 a 1)Z 2 ), a 0, for i.i.d. 1 Fréchet Z 1 & Z 2. Note that α appears only as a factor in Σ α : Σ α = α 2 Σ 1. It does not affect the correlation structure of the max spectrum. GLS estimators for α use the matrix Σ α. Asymptotic normality for α(r + j 1, r + j 2 ) follows from the last theorem. Stoev, Michailidis and Taqqu (2006) has details on confidence intervals for α and automatic selection of j 1 &j 2. 31
32 Automatic selection of scales Null hypothesis : Y = {Y j } J j=1 is multivariate Normal with covariance matrix α 2 Σ 1 and means EY j = j/α + const,1 j J. Given a range of scales j 1 j j 2, set Ĥ(j 1, j 2 ) 1/ α(j 1, j 2 ) := w j 1,j 2 Y, for the GLS estimator of the slope H = 1/α of Y j vs j. Algorithm: 1. Set j 2 J = [log 2 (n)] and j 1 := j Compute Ĥ +1 := Ĥ(j 1 1, j 2 ) and Ĥ 0 := Ĥ(j 1, j 2 ). 3. Compute a conf int for Ĥ +1 Ĥ 0 (based on the null hypothesis). 4. If the conf int does not contain 0, then stop. Else, set j 1 := j 1 1 and repeat Step 2. Important details: Plug in for α used in the covariance matrix Level of conf int to be set by the user Multiple testing issue not addressed The method works very well in practice! 32
33 Automatic selection of scales: an illustration Simulated: a mixture of 1 Fréchet (10%) and Exponential with mean 5 (90%). Sample size: n = 2 17 = 131,072. Root MSE Automatic choices of j 1 : p= j 1 Square root MSE j Automatic j 1 : α estimates α MSE best j 1 =10: α estimates α 33
34 Rates for moment functionals of maxima Let X 1,..., X n be i.i.d. from F and recall that M n := 1 d X n 1/α i Z, (n ). 1 i n Are there results on the rate of convergence for reasonable f s? Ef(M n ) Ef(Z), (n ) Pickands (1975) shows only the convergence of moments (no rates). Our approach: consider with F(x) = exp{ σ α (x)x α }, x R, (4) σ α (x) C, (x ). Note: 1 F(x) Cx α, x is equivalent to (4). Extra assumption: σ α (x) C Dx β, for large x. 34
35 Rates (cont d) But (4) is very convenient to handle rates! Note that Thus, Ef(M n ) = and also P{M n x} Ef(Z) = 0 0 = P{X 1 n 1/α x} n = F(n 1/α x) n = exp{ σ α (n 1/α x)x α }. (5) f(x)df n (x) = f(x)dg(x) = 0 Now, integration by parts yields: E(f(M n ) f(z)) = 0 0 f(x)dexp{ σ α (n 1/α x)x α }, f(x)dexp{ Cx α }. (G(x) F n (x))f (x)dx. However G(x) and F n (x) are of the same exponential form! 35
36 Rates (cont d) By the mean value theorem: F n (x) G(x) = exp{ σ α (n 1/α x)x α } exp{ Cx α } σ α (n 1/α x) C x α exp{ cx α } Dn β/α x (α+β) exp{ cx α }, as x, Thus, for any ǫ > 0, as n, ǫ F n (x) G(x) f (x) dx Dn β/α f (x) x (α+β) e cx α dx = O(n β/α ). 0 By taking ǫ = ǫ n 0, and using a mild technical condition we can also bound the integral near zero ǫn 0 F n (x) G(x) f (x) dx. Th If σ α (x) Dx β, as x, then as n, n β/α (Ef(M n ) Ef(Z)) D 0 x (α+β) f (x)e Cx α dx, provided mild technical conditions on σ(x) at 0 and on f at 0 and hold. 36
37 Rates (cont d) We have thus obtained exact rates for moment functionals Ef( 1 max n X 1/α i) 1 i n They are valid for a large class of absolutely continuous f, including: f(x) = log 2 (x), and f(x) = x p, p (0, α), for example. More details can be found in Stoev, Michailidis and Taqqu (2006). As a corollary, we also get rates of convergence for Cov(log 2 D(r+j 1, k),log 2 D(r+j 2, k)), as r = r(n). These tools and the Berry Esseen Theorem, yield the uniform asymptotic normality of the max spectrum. Thank you! 37
38 References Alon, N., Matias, Y. & Szegedy, M. (1996), The space complexity of approximating the frequency moments, in Proceedings of the 28th ACM Symposium on Theory of Computing, pp Cormode, G. & Muthukrishnan, S. (2003), Algorithms - ESA 2003, Vol of Lecture Notes in Computer Science, Springer, chapter Estimating Dominance Norms of Multiple Data Streams, pp Csörgő, S., Deheuvels, P. & Mason, D. (1985), Kernel estimates of the tail index of a distribution, Annals of Statistics 13(3), Datar, M., Immorlica, N., Indyk, P. & Mirrokni, V. S. (2004), Locality sensitive hashing scheme based on p stable distributions, in SCG 04: Proceedings of the twentieth annual symposium on Computational geometry, ACM Press, New York, NY, USA, pp Dekkers, A., Einmahl, J. & de Haan, L. (1989), A moment estimator for the index of an extreme value distribution, Ann. Statist. 17(4), Feigenbaum, J., Kannan, S., Strauss, M. & Viswanathan, M. (1999), An approximate l1-difference algorithm for massive data streams, in FOCS 99: Proceedings of the 40th Annual Symposium on Foundations of Computer Science, IEEE Computer Society, Washington, DC, USA, p Fong, J. H. & Strauss, M. (2000), An approximate lp-difference algorithm for massive data streams, in STACS 00: Proceedings of the 17th Annual Symposium on Theoretical Aspects of Computer Science, Springer-Verlag, London, UK, pp Gilbert, A., Kotidis, Y., Muthukrishnan, S. & Strauss, M. (2001), Surfing wavelets on streams: one-pass summaries for approximate aggregate queries, in Procedings of VLDB, Rome, Italy.
39 Gilbert, A., Kotidis, Y., Muthukrishnan, S. & Strauss, M. (2002), How to summarize the universe: Dynamic maintenance of quantiles. Hall, P. (1982), On some simple estimates of an exponent of regular variation, J. Roy. Stat. Assoc. 44, Series B. Hill, B. M. (1975), A simple general approach to inference about the tail of a distribution, The Annals of Statistics 3, Indyk, P. (2000), Stable distributions, pseudorandom generators, embeddings and data stream computation, in IEEE Symposium on Foundations of Computer Science, pp Muthukrishnan, S. (2003), Data streams: algorithms and applications, Preprint. muthu/. Pickands, J. (1975), Statistical inference using extreme order statistics, Ann. Statist. 3, Stoev, S., Hadjieleftheriou, M., Kollios, G. & Taqqu, M. (2007), Norm, point and distance estimation over multiple signals using max stable distributions, in 23rd International Conference on Data Engineering, ICDE, IEEE. To appear. Stoev, S., Michailidis, G. & Taqqu, M. (2006), Estimating heavy tail exponents through max self similarity, Preprint. Weissman, I. (1978), Estimation of parameters and large quantiles based on the k largest observations, Journal of the American Statistical Association 73,
On the estimation of the heavy tail exponent in time series using the max spectrum. Stilian A. Stoev
On the estimation of the heavy tail exponent in time series using the max spectrum Stilian A. Stoev (sstoev@umich.edu) University of Michigan, Ann Arbor, U.S.A. JSM, Salt Lake City, 007 joint work with:
More informationarxiv:math/ v1 [math.st] 6 Sep 2006
Estimating heavy tail exponents through max self similarity arxiv:math/69163v1 [math.st] 6 Sep 26 1. Introduction Stilian A. Stoev University of Michigan, Ann Arbor George Michailidis University of Michigan,
More informationarxiv: v1 [cs.ds] 24 May 2010
arxiv:1005.4344v1 [cs.ds] 24 May 2010 Max stable sketches: estimation of l α norms, dominance norms and point queries for non negative signals Stilian A. Stoev Department of Statistics University of Michigan,
More informationWavelet decomposition of data streams. by Dragana Veljkovic
Wavelet decomposition of data streams by Dragana Veljkovic Motivation Continuous data streams arise naturally in: telecommunication and internet traffic retail and banking transactions web server log records
More informationAnalysis methods of heavy-tailed data
Institute of Control Sciences Russian Academy of Sciences, Moscow, Russia February, 13-18, 2006, Bamberg, Germany June, 19-23, 2006, Brest, France May, 14-19, 2007, Trondheim, Norway PhD course Chapter
More informationMax stable Processes & Random Fields: Representations, Models, and Prediction
Max stable Processes & Random Fields: Representations, Models, and Prediction Stilian Stoev University of Michigan, Ann Arbor March 2, 2011 Based on joint works with Yizao Wang and Murad S. Taqqu. 1 Preliminaries
More informationIntroduction to Algorithmic Trading Strategies Lecture 10
Introduction to Algorithmic Trading Strategies Lecture 10 Risk Management Haksun Li haksun.li@numericalmethod.com www.numericalmethod.com Outline Value at Risk (VaR) Extreme Value Theory (EVT) References
More informationData Stream Methods. Graham Cormode S. Muthukrishnan
Data Stream Methods Graham Cormode graham@dimacs.rutgers.edu S. Muthukrishnan muthu@cs.rutgers.edu Plan of attack Frequent Items / Heavy Hitters Counting Distinct Elements Clustering items in Streams Motivating
More informationExtreme Value Analysis and Spatial Extremes
Extreme Value Analysis and Department of Statistics Purdue University 11/07/2013 Outline Motivation 1 Motivation 2 Extreme Value Theorem and 3 Bayesian Hierarchical Models Copula Models Max-stable Models
More informationMath 576: Quantitative Risk Management
Math 576: Quantitative Risk Management Haijun Li lih@math.wsu.edu Department of Mathematics Washington State University Week 11 Haijun Li Math 576: Quantitative Risk Management Week 11 1 / 21 Outline 1
More informationEstimating Dominance Norms of Multiple Data Streams Graham Cormode Joint work with S. Muthukrishnan
Estimating Dominance Norms of Multiple Data Streams Graham Cormode graham@dimacs.rutgers.edu Joint work with S. Muthukrishnan Data Stream Phenomenon Data is being produced faster than our ability to process
More informationLinear Sketches A Useful Tool in Streaming and Compressive Sensing
Linear Sketches A Useful Tool in Streaming and Compressive Sensing Qin Zhang 1-1 Linear sketch Random linear projection M : R n R k that preserves properties of any v R n with high prob. where k n. M =
More informationMFM Practitioner Module: Quantitiative Risk Management. John Dodson. October 14, 2015
MFM Practitioner Module: Quantitiative Risk Management October 14, 2015 The n-block maxima 1 is a random variable defined as M n max (X 1,..., X n ) for i.i.d. random variables X i with distribution function
More informationMultivariate Heavy Tails, Asymptotic Independence and Beyond
Multivariate Heavy Tails, endence and Beyond Sidney Resnick School of Operations Research and Industrial Engineering Rhodes Hall Cornell University Ithaca NY 14853 USA http://www.orie.cornell.edu/ sid
More informationEstimating Dominance Norms of Multiple Data Streams
Estimating Dominance Norms of Multiple Data Streams Graham Cormode and S. Muthukrishnan Abstract. There is much focus in the algorithms and database communities on designing tools to manage and mine data
More informationWe introduce methods that are useful in:
Instructor: Shengyu Zhang Content Derived Distributions Covariance and Correlation Conditional Expectation and Variance Revisited Transforms Sum of a Random Number of Independent Random Variables more
More informationLecture 2. Frequency problems
1 / 43 Lecture 2. Frequency problems Ricard Gavaldà MIRI Seminar on Data Streams, Spring 2015 Contents 2 / 43 1 Frequency problems in data streams 2 Approximating inner product 3 Computing frequency moments
More informationA MODIFICATION OF HILL S TAIL INDEX ESTIMATOR
L. GLAVAŠ 1 J. JOCKOVIĆ 2 A MODIFICATION OF HILL S TAIL INDEX ESTIMATOR P. MLADENOVIĆ 3 1, 2, 3 University of Belgrade, Faculty of Mathematics, Belgrade, Serbia Abstract: In this paper, we study a class
More informationSketching and Streaming for Distributions
Sketching and Streaming for Distributions Piotr Indyk Andrew McGregor Massachusetts Institute of Technology University of California, San Diego Main Material: Stable distributions, pseudo-random generators,
More informationSTAT 200C: High-dimensional Statistics
STAT 200C: High-dimensional Statistics Arash A. Amini May 30, 2018 1 / 59 Classical case: n d. Asymptotic assumption: d is fixed and n. Basic tools: LLN and CLT. High-dimensional setting: n d, e.g. n/d
More informationRandom Variables. Random variables. A numerically valued map X of an outcome ω from a sample space Ω to the real line R
In probabilistic models, a random variable is a variable whose possible values are numerical outcomes of a random phenomenon. As a function or a map, it maps from an element (or an outcome) of a sample
More informationChapter 5. Chapter 5 sections
1 / 43 sections Discrete univariate distributions: 5.2 Bernoulli and Binomial distributions Just skim 5.3 Hypergeometric distributions 5.4 Poisson distributions Just skim 5.5 Negative Binomial distributions
More informationMax stable processes: representations, ergodic properties and some statistical applications
Max stable processes: representations, ergodic properties and some statistical applications Stilian Stoev University of Michigan, Ann Arbor Oberwolfach March 21, 2008 The focus Let X = {X t } t R be a
More informationOverview of Extreme Value Theory. Dr. Sawsan Hilal space
Overview of Extreme Value Theory Dr. Sawsan Hilal space Maths Department - University of Bahrain space November 2010 Outline Part-1: Univariate Extremes Motivation Threshold Exceedances Part-2: Bivariate
More informationQuantile-quantile plots and the method of peaksover-threshold
Problems in SF2980 2009-11-09 12 6 4 2 0 2 4 6 0.15 0.10 0.05 0.00 0.05 0.10 0.15 Figure 2: qqplot of log-returns (x-axis) against quantiles of a standard t-distribution with 4 degrees of freedom (y-axis).
More informationThe Convergence Rate for the Normal Approximation of Extreme Sums
The Convergence Rate for the Normal Approximation of Extreme Sums Yongcheng Qi University of Minnesota Duluth WCNA 2008, Orlando, July 2-9, 2008 This talk is based on a joint work with Professor Shihong
More informationDoes k-th Moment Exist?
Does k-th Moment Exist? Hitomi, K. 1 and Y. Nishiyama 2 1 Kyoto Institute of Technology, Japan 2 Institute of Economic Research, Kyoto University, Japan Email: hitomi@kit.ac.jp Keywords: Existence of moments,
More informationExtreme Value Theory and Applications
Extreme Value Theory and Deauville - 04/10/2013 Extreme Value Theory and Introduction Asymptotic behavior of the Sum Extreme (from Latin exter, exterus, being on the outside) : Exceeding the ordinary,
More informationSTAT 512 sp 2018 Summary Sheet
STAT 5 sp 08 Summary Sheet Karl B. Gregory Spring 08. Transformations of a random variable Let X be a rv with support X and let g be a function mapping X to Y with inverse mapping g (A = {x X : g(x A}
More informationShape of the return probability density function and extreme value statistics
Shape of the return probability density function and extreme value statistics 13/09/03 Int. Workshop on Risk and Regulation, Budapest Overview I aim to elucidate a relation between one field of research
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 informationModel Fitting. Jean Yves Le Boudec
Model Fitting Jean Yves Le Boudec 0 Contents 1. What is model fitting? 2. Linear Regression 3. Linear regression with norm minimization 4. Choosing a distribution 5. Heavy Tail 1 Virus Infection Data We
More informationB669 Sublinear Algorithms for Big Data
B669 Sublinear Algorithms for Big Data Qin Zhang 1-1 2-1 Part 1: Sublinear in Space The model and challenge The data stream model (Alon, Matias and Szegedy 1996) a n a 2 a 1 RAM CPU Why hard? Cannot store
More informationA Few Notes on Fisher Information (WIP)
A Few Notes on Fisher Information (WIP) David Meyer dmm@{-4-5.net,uoregon.edu} Last update: April 30, 208 Definitions There are so many interesting things about Fisher Information and its theoretical properties
More informationData Streams & Communication Complexity
Data Streams & Communication Complexity Lecture 1: Simple Stream Statistics in Small Space Andrew McGregor, UMass Amherst 1/25 Data Stream Model Stream: m elements from universe of size n, e.g., x 1, x
More informationFormulas for probability theory and linear models SF2941
Formulas for probability theory and linear models SF2941 These pages + Appendix 2 of Gut) are permitted as assistance at the exam. 11 maj 2008 Selected formulae of probability Bivariate probability Transforms
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 informationf (1 0.5)/n Z =
Math 466/566 - Homework 4. We want to test a hypothesis involving a population proportion. The unknown population proportion is p. The null hypothesis is p = / and the alternative hypothesis is p > /.
More informationA New Estimator for a Tail Index
Acta Applicandae Mathematicae 00: 3, 2003. 2003 Kluwer Academic Publishers. Printed in the Netherlands. A New Estimator for a Tail Index V. PAULAUSKAS Department of Mathematics and Informatics, Vilnius
More informationMaintaining Significant Stream Statistics over Sliding Windows
Maintaining Significant Stream Statistics over Sliding Windows L.K. Lee H.F. Ting Abstract In this paper, we introduce the Significant One Counting problem. Let ε and θ be respectively some user-specified
More informationSTAT 430/510: Lecture 16
STAT 430/510: Lecture 16 James Piette June 24, 2010 Updates HW4 is up on my website. It is due next Mon. (June 28th). Starting today back at section 6.7 and will begin Ch. 7. Joint Distribution of Functions
More informationChapter 2. Continuous random variables
Chapter 2 Continuous random variables Outline Review of probability: events and probability Random variable Probability and Cumulative distribution function Review of discrete random variable Introduction
More informationLecture Lecture 25 November 25, 2014
CS 224: Advanced Algorithms Fall 2014 Lecture Lecture 25 November 25, 2014 Prof. Jelani Nelson Scribe: Keno Fischer 1 Today Finish faster exponential time algorithms (Inclusion-Exclusion/Zeta Transform,
More informationEstimation of risk measures for extreme pluviometrical measurements
Estimation of risk measures for extreme pluviometrical measurements by Jonathan EL METHNI in collaboration with Laurent GARDES & Stéphane GIRARD 26th Annual Conference of The International Environmetrics
More informationTail negative dependence and its applications for aggregate loss modeling
Tail negative dependence and its applications for aggregate loss modeling Lei Hua Division of Statistics Oct 20, 2014, ISU L. Hua (NIU) 1/35 1 Motivation 2 Tail order Elliptical copula Extreme value copula
More informationA Deterministic Algorithm for Summarizing Asynchronous Streams over a Sliding Window
A Deterministic Algorithm for Summarizing Asynchronous Streams over a Sliding Window Costas Busch 1 and Srikanta Tirthapura 2 1 Department of Computer Science Rensselaer Polytechnic Institute, Troy, NY
More informationSTAT Chapter 5 Continuous Distributions
STAT 270 - Chapter 5 Continuous Distributions June 27, 2012 Shirin Golchi () STAT270 June 27, 2012 1 / 59 Continuous rv s Definition: X is a continuous rv if it takes values in an interval, i.e., range
More informationSYSM 6303: Quantitative Introduction to Risk and Uncertainty in Business Lecture 4: Fitting Data to Distributions
SYSM 6303: Quantitative Introduction to Risk and Uncertainty in Business Lecture 4: Fitting Data to Distributions M. Vidyasagar Cecil & Ida Green Chair The University of Texas at Dallas Email: M.Vidyasagar@utdallas.edu
More informationTitle: Count-Min Sketch Name: Graham Cormode 1 Affil./Addr. Department of Computer Science, University of Warwick,
Title: Count-Min Sketch Name: Graham Cormode 1 Affil./Addr. Department of Computer Science, University of Warwick, Coventry, UK Keywords: streaming algorithms; frequent items; approximate counting, sketch
More information8 Laws of large numbers
8 Laws of large numbers 8.1 Introduction We first start with the idea of standardizing a random variable. Let X be a random variable with mean µ and variance σ 2. Then Z = (X µ)/σ will be a random variable
More informationEC212: Introduction to Econometrics Review Materials (Wooldridge, Appendix)
1 EC212: Introduction to Econometrics Review Materials (Wooldridge, Appendix) Taisuke Otsu London School of Economics Summer 2018 A.1. Summation operator (Wooldridge, App. A.1) 2 3 Summation operator For
More information4. Distributions of Functions of Random Variables
4. Distributions of Functions of Random Variables Setup: Consider as given the joint distribution of X 1,..., X n (i.e. consider as given f X1,...,X n and F X1,...,X n ) Consider k functions g 1 : R n
More informationConditional Sampling for Max Stable Random Fields
Conditional Sampling for Max Stable Random Fields Yizao Wang Department of Statistics, the University of Michigan April 30th, 0 th GSPC at Duke University, Durham, North Carolina Joint work with Stilian
More informationAsymptotic Statistics-III. Changliang Zou
Asymptotic Statistics-III Changliang Zou The multivariate central limit theorem Theorem (Multivariate CLT for iid case) Let X i be iid random p-vectors with mean µ and and covariance matrix Σ. Then n (
More informationSome functional (Hölderian) limit theorems and their applications (II)
Some functional (Hölderian) limit theorems and their applications (II) Alfredas Račkauskas Vilnius University Outils Statistiques et Probabilistes pour la Finance Université de Rouen June 1 5, Rouen (Rouen
More informationStatistical inference on Lévy processes
Alberto Coca Cabrero University of Cambridge - CCA Supervisors: Dr. Richard Nickl and Professor L.C.G.Rogers Funded by Fundación Mutua Madrileña and EPSRC MASDOC/CCA student workshop 2013 26th March Outline
More informationNetwork Traffic Characteristic
Network Traffic Characteristic Hojun Lee hlee02@purros.poly.edu 5/24/2002 EL938-Project 1 Outline Motivation What is self-similarity? Behavior of Ethernet traffic Behavior of WAN traffic Behavior of WWW
More informationFinancial Econometrics and Volatility Models Extreme Value Theory
Financial Econometrics and Volatility Models Extreme Value Theory Eric Zivot May 3, 2010 1 Lecture Outline Modeling Maxima and Worst Cases The Generalized Extreme Value Distribution Modeling Extremes Over
More informationChallenges in implementing worst-case analysis
Challenges in implementing worst-case analysis Jon Danielsson Systemic Risk Centre, lse,houghton Street, London WC2A 2AE, UK Lerby M. Ergun Systemic Risk Centre, lse,houghton Street, London WC2A 2AE, UK
More informationChapter 2: Fundamentals of Statistics Lecture 15: Models and statistics
Chapter 2: Fundamentals of Statistics Lecture 15: Models and statistics Data from one or a series of random experiments are collected. Planning experiments and collecting data (not discussed here). Analysis:
More informationPerhaps the simplest way of modeling two (discrete) random variables is by means of a joint PMF, defined as follows.
Chapter 5 Two Random Variables In a practical engineering problem, there is almost always causal relationship between different events. Some relationships are determined by physical laws, e.g., voltage
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 informationClassical Extreme Value Theory - An Introduction
Chapter 1 Classical Extreme Value Theory - An Introduction 1.1 Introduction Asymptotic theory of functions of random variables plays a very important role in modern statistics. The objective of the asymptotic
More informationEstimation of the long Memory parameter using an Infinite Source Poisson model applied to transmission rate measurements
of the long Memory parameter using an Infinite Source Poisson model applied to transmission rate measurements François Roueff Ecole Nat. Sup. des Télécommunications 46 rue Barrault, 75634 Paris cedex 13,
More informationBasic concepts of probability theory
Basic concepts of probability theory Random variable discrete/continuous random variable Transform Z transform, Laplace transform Distribution Geometric, mixed-geometric, Binomial, Poisson, exponential,
More informationLecture 6 September 13, 2016
CS 395T: Sublinear Algorithms Fall 206 Prof. Eric Price Lecture 6 September 3, 206 Scribe: Shanshan Wu, Yitao Chen Overview Recap of last lecture. We talked about Johnson-Lindenstrauss (JL) lemma [JL84]
More information1: PROBABILITY REVIEW
1: PROBABILITY REVIEW Marek Rutkowski School of Mathematics and Statistics University of Sydney Semester 2, 2016 M. Rutkowski (USydney) Slides 1: Probability Review 1 / 56 Outline We will review the following
More informationLecture 2: Repetition of probability theory and statistics
Algorithms for Uncertainty Quantification SS8, IN2345 Tobias Neckel Scientific Computing in Computer Science TUM Lecture 2: Repetition of probability theory and statistics Concept of Building Block: Prerequisites:
More informationWavelet and SiZer analyses of Internet Traffic Data
Wavelet and SiZer analyses of Internet Traffic Data Cheolwoo Park Statistical and Applied Mathematical Sciences Institute Fred Godtliebsen Department of Mathematics and Statistics, University of Tromsø
More informationNon-parametric Inference and Resampling
Non-parametric Inference and Resampling Exercises by David Wozabal (Last update. Juni 010) 1 Basic Facts about Rank and Order Statistics 1.1 10 students were asked about the amount of time they spend surfing
More informationLecture 1: August 28
36-705: Intermediate Statistics Fall 2017 Lecturer: Siva Balakrishnan Lecture 1: August 28 Our broad goal for the first few lectures is to try to understand the behaviour of sums of independent random
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 informationStat 710: Mathematical Statistics Lecture 31
Stat 710: Mathematical Statistics Lecture 31 Jun Shao Department of Statistics University of Wisconsin Madison, WI 53706, USA Jun Shao (UW-Madison) Stat 710, Lecture 31 April 13, 2009 1 / 13 Lecture 31:
More informationMath 494: Mathematical Statistics
Math 494: Mathematical Statistics Instructor: Jimin Ding jmding@wustl.edu Department of Mathematics Washington University in St. Louis Class materials are available on course website (www.math.wustl.edu/
More informationLecture Notes 5 Convergence and Limit Theorems. Convergence with Probability 1. Convergence in Mean Square. Convergence in Probability, WLLN
Lecture Notes 5 Convergence and Limit Theorems Motivation Convergence with Probability Convergence in Mean Square Convergence in Probability, WLLN Convergence in Distribution, CLT EE 278: Convergence and
More informationThe Count-Min-Sketch and its Applications
The Count-Min-Sketch and its Applications Jannik Sundermeier Abstract In this thesis, we want to reveal how to get rid of a huge amount of data which is at least dicult or even impossible to store in local
More informationAsymptotic distribution of the sample average value-at-risk
Asymptotic distribution of the sample average value-at-risk Stoyan V. Stoyanov Svetlozar T. Rachev September 3, 7 Abstract In this paper, we prove a result for the asymptotic distribution of the sample
More informationNonlinear Time Series Modeling
Nonlinear Time Series Modeling Part II: Time Series Models in Finance Richard A. Davis Colorado State University (http://www.stat.colostate.edu/~rdavis/lectures) MaPhySto Workshop Copenhagen September
More informationLecture 32: Asymptotic confidence sets and likelihoods
Lecture 32: Asymptotic confidence sets and likelihoods Asymptotic criterion In some problems, especially in nonparametric problems, it is difficult to find a reasonable confidence set with a given confidence
More informationThis does not cover everything on the final. Look at the posted practice problems for other topics.
Class 7: Review Problems for Final Exam 8.5 Spring 7 This does not cover everything on the final. Look at the posted practice problems for other topics. To save time in class: set up, but do not carry
More informationThe Behavior of Multivariate Maxima of Moving Maxima Processes
The Behavior of Multivariate Maxima of Moving Maxima Processes Zhengjun Zhang Department of Mathematics Washington University Saint Louis, MO 6313-4899 USA Richard L. Smith Department of Statistics University
More informationBivariate distributions
Bivariate distributions 3 th October 017 lecture based on Hogg Tanis Zimmerman: Probability and Statistical Inference (9th ed.) Bivariate Distributions of the Discrete Type The Correlation Coefficient
More informationUC Berkeley Department of Electrical Engineering and Computer Sciences. EECS 126: Probability and Random Processes
UC Berkeley Department of Electrical Engineering and Computer Sciences EECS 6: Probability and Random Processes Problem Set 3 Spring 9 Self-Graded Scores Due: February 8, 9 Submit your self-graded scores
More informationRefining the Central Limit Theorem Approximation via Extreme Value Theory
Refining the Central Limit Theorem Approximation via Extreme Value Theory Ulrich K. Müller Economics Department Princeton University February 2018 Abstract We suggest approximating the distribution of
More informationNorm, Point, and Distance Estimation Over Multiple Signals Using Max Stable Distributions
Norm, Point, and Distance Estimation Over Multiple Signals Using Max Stable Distributions Stilian A. Stoev Marios Hadjieleftheriou George Kollios Murad S. Taqqu University of Michigan, Ann Arbor, AT&T
More informationBasic concepts of probability theory
Basic concepts of probability theory Random variable discrete/continuous random variable Transform Z transform, Laplace transform Distribution Geometric, mixed-geometric, Binomial, Poisson, exponential,
More informationLecture 01 August 31, 2017
Sketching Algorithms for Big Data Fall 2017 Prof. Jelani Nelson Lecture 01 August 31, 2017 Scribe: Vinh-Kha Le 1 Overview In this lecture, we overviewed the six main topics covered in the course, reviewed
More informationA NOTE ON SECOND ORDER CONDITIONS IN EXTREME VALUE THEORY: LINKING GENERAL AND HEAVY TAIL CONDITIONS
REVSTAT Statistical Journal Volume 5, Number 3, November 2007, 285 304 A NOTE ON SECOND ORDER CONDITIONS IN EXTREME VALUE THEORY: LINKING GENERAL AND HEAVY TAIL CONDITIONS Authors: M. Isabel Fraga Alves
More informationNonparametric Bayesian Methods (Gaussian Processes)
[70240413 Statistical Machine Learning, Spring, 2015] Nonparametric Bayesian Methods (Gaussian Processes) Jun Zhu dcszj@mail.tsinghua.edu.cn http://bigml.cs.tsinghua.edu.cn/~jun State Key Lab of Intelligent
More informationA Conditional Approach to Modeling Multivariate Extremes
A Approach to ing Multivariate Extremes By Heffernan & Tawn Department of Statistics Purdue University s April 30, 2014 Outline s s Multivariate Extremes s A central aim of multivariate extremes is trying
More informationStable Process. 2. Multivariate Stable Distributions. July, 2006
Stable Process 2. Multivariate Stable Distributions July, 2006 1. Stable random vectors. 2. Characteristic functions. 3. Strictly stable and symmetric stable random vectors. 4. Sub-Gaussian random vectors.
More informationRegression and Statistical Inference
Regression and Statistical Inference Walid Mnif wmnif@uwo.ca Department of Applied Mathematics The University of Western Ontario, London, Canada 1 Elements of Probability 2 Elements of Probability CDF&PDF
More informationThe Moment Method; Convex Duality; and Large/Medium/Small Deviations
Stat 928: Statistical Learning Theory Lecture: 5 The Moment Method; Convex Duality; and Large/Medium/Small Deviations Instructor: Sham Kakade The Exponential Inequality and Convex Duality The exponential
More informationInference and Regression
Inference and Regression Assignment 3 Department of IOMS Professor William Greene Phone:.998.0876 Office: KMC 7-90 Home page:www.stern.nyu.edu/~wgreene Email: wgreene@stern.nyu.edu Course web page: www.stern.nyu.edu/~wgreene/econometrics/econometrics.htm.
More informationSTAT 430/510: Lecture 10
STAT 430/510: Lecture 10 James Piette June 9, 2010 Updates HW2 is due today! Pick up your HW1 s up in stat dept. There is a box located right when you enter that is labeled "Stat 430 HW1". It ll be out
More informationHIERARCHICAL MODELS IN EXTREME VALUE THEORY
HIERARCHICAL MODELS IN EXTREME VALUE THEORY Richard L. Smith Department of Statistics and Operations Research, University of North Carolina, Chapel Hill and Statistical and Applied Mathematical Sciences
More informationNonparametric Inference via Bootstrapping the Debiased Estimator
Nonparametric Inference via Bootstrapping the Debiased Estimator Yen-Chi Chen Department of Statistics, University of Washington ICSA-Canada Chapter Symposium 2017 1 / 21 Problem Setup Let X 1,, X n be
More informationESTIMATING BIVARIATE TAIL
Elena DI BERNARDINO b joint work with Clémentine PRIEUR a and Véronique MAUME-DESCHAMPS b a LJK, Université Joseph Fourier, Grenoble 1 b Laboratoire SAF, ISFA, Université Lyon 1 Framework Goal: estimating
More informationProbabilities & Statistics Revision
Probabilities & Statistics Revision Christopher Ting Christopher Ting http://www.mysmu.edu/faculty/christophert/ : christopherting@smu.edu.sg : 6828 0364 : LKCSB 5036 January 6, 2017 Christopher Ting QF
More informationBasic concepts of probability theory
Basic concepts of probability theory Random variable discrete/continuous random variable Transform Z transform, Laplace transform Distribution Geometric, mixed-geometric, Binomial, Poisson, exponential,
More information