Dimension Reduction and Visualization of the Histogram Data

Size: px
Start display at page:

Download "Dimension Reduction and Visualization of the Histogram Data"

Transcription

1 The 4th Workshop n Symbolc Data Analyss (SDA 214): Tutoral Dmenson Reducton and Vsualzaton of the Hstogram Data Han-Mng Wu ( 吳漢銘 ) Department of Mathematcs Tamkang Unversty Tamsu 25137, Tawan hmwu@mal.tku.edu.tw DR for Hstogram Data June 13, / 27

2 Outlnes 1 The Hstogram Data DR for Hstogram Data June 13, / 27

3 Outlnes 1 The Hstogram Data 2 Prncpal Component Analyss for Hstogram Data () The Hstogram Arthmetc Approach () The Quantle Method () The Transformaton Approach (v) The Dstrbutonal Approach DR for Hstogram Data June 13, / 27

4 Outlnes 1 The Hstogram Data 2 Prncpal Component Analyss for Hstogram Data () The Hstogram Arthmetc Approach () The Quantle Method () The Transformaton Approach (v) The Dstrbutonal Approach 3 An Example DR for Hstogram Data June 13, / 27

5 The Hstogram data (1/2) I ( b) x ( b) x ( b+1) = (, ] p ( b) b... B DR for Hstogram Data June 13, / 27

6 The Hstogram data (1/2) I ( b) x ( b) x ( b+1) = (, ] p ( b) A hstogram: b... B h = {I, p} = {x I (b) = (x (b), x (b+1) ], p (b), b = 1,, B}. DR for Hstogram Data June 13, / 27

7 The Hstogram data (1/2) I ( b) x ( b) x ( b+1) = (, ] p ( b) A hstogram: b... B h = {I, p} = {x I (b) = (x (b), x (b+1) ], p (b), b = 1,, B}. 1 B b=1 p(b) = 1, B: the number of modaltes of h. DR for Hstogram Data June 13, / 27

8 The Hstogram data (1/2) I ( b) x ( b) x ( b+1) = (, ] p ( b) b... B A hstogram: h = {I, p} = {x I (b) = (x (b), x (b+1) ], p (b), b = 1,, B}. 1 B b=1 p(b) = 1, B: the number of modaltes of h. 2 I = b I (b) = (x (1), x (B+1) ]. DR for Hstogram Data June 13, / 27

9 The Hstogram data (2/2) Let O = {o 1, o 2,, o n } be gven objects. 1 The classcal data: x = (x 1 x p ) T, x R. DR for Hstogram Data June 13, / 27

10 The Hstogram data (2/2) Let O = {o 1, o 2,, o n } be gven objects. 1 The classcal data: x = (x 1 x p ) T, x R. 2 The hstogram data: H = (h ) n p, where h = {I, p } = {x I (b) = (x (b), x (b+1) ], p (b), b = 1,, B }. DR for Hstogram Data June 13, / 27

11 The Hstogram data (2/2) Let O = {o 1, o 2,, o n } be gven objects. 1 The classcal data: x = (x 1 x p ) T, x R. 2 The hstogram data: H = (h ) n p, where h = {I, p } = {x I (b) = (x (b), x (b+1) ], p (b), b = 1,, B }. Possble source: the result of an aggregaton, the descrpton of a populaton, or any other grouped collectve. DR for Hstogram Data June 13, / 27

12 PCA for Hstogram Data (1/2) The Hstogram Arthmetc Approach 1 Rodrguez, O., Dday, E., and Wnsberg, S., (2), Generalzaton of the prncpal component analyss to hstogram data, PKDD2, Lyon, 2. DR for Hstogram Data June 13, / 27

13 PCA for Hstogram Data (1/2) The Hstogram Arthmetc Approach 1 Rodrguez, O., Dday, E., and Wnsberg, S., (2), Generalzaton of the prncpal component analyss to hstogram data, PKDD2, Lyon, 2. The Quantle Method 1 Ichno, M., (28), Symbolc PCA for hstogram-valued data, n Proceedngs IASC. December 5-8, Yokohama, Japan. 2 Ichno, M., (211), The quantle method for symbolc prncpal component analyss, Statstcal Analyss and Data Mnng, 4(2): Verde, R. and Irpno, A., (213), Dmenson reducton technques for dstrbutonal symbolc data, SIS 213 Statstcal Conference, Advances n Latent Varables - Methods, Models and Applcatons, Unversty of Bresca, June, 19-21, 213. (Consder only one unvarate hstogram varable.) DR for Hstogram Data June 13, / 27

14 PCA for Hstogram Data (2/2) The Transformaton Approach 1 Rodrguez, O., Dday, E., and Wnsberg, S., (2), Generalzaton of the prncpal component analyss to hstogram data, PKDD2, Lyon, 2. 2 Makosso-Kallyth, S. and Dday, E., (212), Adaptaton of nterval PCA to symbolc hstogram varables, Adv Data Anal Classf, 6: DR for Hstogram Data June 13, / 27

15 PCA for Hstogram Data (2/2) The Transformaton Approach 1 Rodrguez, O., Dday, E., and Wnsberg, S., (2), Generalzaton of the prncpal component analyss to hstogram data, PKDD2, Lyon, 2. 2 Makosso-Kallyth, S. and Dday, E., (212), Adaptaton of nterval PCA to symbolc hstogram varables, Adv Data Anal Classf, 6: The Dstrbutonal Approach 1 Wang, H., Chen, M., L, N., and Wang, L., (211), Prncpal component analyss of Modal nterval-valued data wth constant numercal characterstcs, Int. Statstcal Inst.: Proc. 58th World Statstcal Congress, 211, Dubln (Sesson CPS53), pp DR for Hstogram Data June 13, / 27

16 The quantle method (1/3) Defne the cumulatve dstrbuton functon F (x) of the hstogram: 1 Assume that n bns (sub-ntervals) have unform dstrbutons. DR for Hstogram Data June 13, / 27

17 The quantle method (1/3) Defne the cumulatve dstrbuton functon F (x) of the hstogram: 1 Assume that n bns (sub-ntervals) have unform dstrbutons. F (x) = for x < x 1, F (x) = p 1 x x 1 x 2 x 1 for x 1 x < x 2, F (x) = F (x 2 ) + p 2 x x 2 x 3 x 2 for x 2 x < x 3,. F (x) = F (x n ) + p n x x n x n+1 x n for x n x < x n+1, F (x) = 1 for x n+1 x. DR for Hstogram Data June 13, / 27

18 The quantle method (2/3) Obtan quantles for a hstogram-valued feature: 1 The object o s descrbed by an (m + 1)-tuple (q 1, q 2,, q m, q m+1 ). F (x mn ) =, q 1 = x mn F (q 2 ) = 1/m, F (q 3 ) = 2/m,, F (q m ) = (m 1)/m, and F (x max ) = 1, q (m+1) = x max DR for Hstogram Data June 13, / 27

19 The quantle method (2/3) Obtan quantles for a hstogram-valued feature: 1 The object o s descrbed by an (m + 1)-tuple (q 1, q 2,, q m, q m+1 ). F (x mn ) =, q 1 = x mn F (q 2 ) = 1/m, F (q 3 ) = 2/m,, F (q m ) = (m 1)/m, and F (x max ) = 1, q (m+1) = x max 2 Express each object o under a hstogram-valued varable X j by an (m + 1)-tuple: (q 1, q 2,, q (m), q (m+1) ). DR for Hstogram Data June 13, / 27

20 The quantle method (3/4) Output: 1 PCA on the correlaton matrx of X Q [n p] = (o Q 1 o Qn ) T, where o Q = q 11 q p1 q 12 q p2... q 1m q pm q 1(m+1) q p(m+1) T, n = n (m + 1) DR for Hstogram Data June 13, / 27

21 The quantle method (3/4) Output: 1 PCA on the correlaton matrx of X Q [n p] = (o Q 1 o Qn ) T, where Vsualzaton: o Q = q 11 q p1 q 12 q p2... q 1m q pm q 1(m+1) q p(m+1) T, n = n (m + 1) 1 In the factor plane, express each object as a seres of m connected arrow lnes. DR for Hstogram Data June 13, / 27

22 The Transformaton Approach (1/2) Makosso-Kallyth, S. and Dday, E., (212), Adaptaton of nterval PCA to symbolc hstogram varables, Adv Data Anal Classf, 6: Preprocessng: 1 Assume bns of the X j varable s the same for all ndvduals. DR for Hstogram Data June 13, / 27

23 The Transformaton Approach (1/2) Makosso-Kallyth, S. and Dday, E., (212), Adaptaton of nterval PCA to symbolc hstogram varables, Adv Data Anal Classf, 6: Preprocessng: 1 Assume bns of the X j varable s the same for all ndvduals. 2 D j = (α j, β j ): Doman of all possble values of bns of X j. DR for Hstogram Data June 13, / 27

24 The Transformaton Approach (1/2) Makosso-Kallyth, S. and Dday, E., (212), Adaptaton of nterval PCA to symbolc hstogram varables, Adv Data Anal Classf, 6: Preprocessng: 1 Assume bns of the X j varable s the same for all ndvduals. 2 D j = (α j, β j ): Doman of all possble values of bns of X j. 3 δ j = nf b=1,,b {x (b+1) x (b) }: smallest length of the nterval I. DR for Hstogram Data June 13, / 27

25 The Transformaton Approach (1/2) Makosso-Kallyth, S. and Dday, E., (212), Adaptaton of nterval PCA to symbolc hstogram varables, Adv Data Anal Classf, 6: Preprocessng: 1 Assume bns of the X j varable s the same for all ndvduals. 2 D j = (α j, β j ): Doman of all possble values of bns of X j. 3 δ j = nf b=1,,b {x (b+1) 4 Determne the lower lmt and the upper lmt 1 f I (b) x (b) }: smallest length of the nterval I. = (, a j ], then I (b) = (max{a j δ j, α j }, a j ]. DR for Hstogram Data June 13, / 27

26 The Transformaton Approach (1/2) Makosso-Kallyth, S. and Dday, E., (212), Adaptaton of nterval PCA to symbolc hstogram varables, Adv Data Anal Classf, 6: Preprocessng: 1 Assume bns of the X j varable s the same for all ndvduals. 2 D j = (α j, β j ): Doman of all possble values of bns of X j. 3 δ j = nf b=1,,b {x (b+1) 4 Determne the lower lmt and the upper lmt x (b) }: smallest length of the nterval I. 1 f I (b) = (, a j ], then I (b) = (max{a j δ j, α j }, a j ]. 2 f I (b) = (b j, ), then I (b) = (b j, max{β j, b j + δ j }]. DR for Hstogram Data June 13, / 27

27 The Transformaton Approach (1/2) Makosso-Kallyth, S. and Dday, E., (212), Adaptaton of nterval PCA to symbolc hstogram varables, Adv Data Anal Classf, 6: Preprocessng: 1 Assume bns of the X j varable s the same for all ndvduals. 2 D j = (α j, β j ): Doman of all possble values of bns of X j. 3 δ j = nf b=1,,b {x (b+1) 4 Determne the lower lmt and the upper lmt x (b) }: smallest length of the nterval I. 1 f I (b) = (, a j ], then I (b) = (max{a j δ j, α j }, a j ]. 2 f I (b) = (b j, ), then I (b) = (b j, max{β j, b j + δ j }]. 5 f bns of varables don t have same unt, replace (a, b ] by an adjusted nterval (a/(b a ), b /(b a )]. DR for Hstogram Data June 13, / 27

28 The Transformaton Approach (2/2) Quantfcaton: 1 The transformed varable: s = (s (1), s (2),, s (B j ) ) 1 Codng 1: s (b) = (x (b) + x (b+1) )/2: center of the adjusted ntervals. DR for Hstogram Data June 13, / 27

29 The Transformaton Approach (2/2) Quantfcaton: 1 The transformed varable: s = (s (1), s (2),, s (B j ) ) 1 Codng 1: s (b) = (x (b) + x (b+1) )/2: center of the adjusted ntervals. 2 Codng 2: s (b) : ranks, the scores of the modaltes of the bns. DR for Hstogram Data June 13, / 27

30 The Transformaton Approach (2/2) Quantfcaton: 1 The transformed varable: s = (s (1), s (2),, s (B j ) ) Output: 1 Codng 1: s (b) = (x (b) + x (b+1) )/2: center of the adjusted ntervals. 2 Codng 2: s (b) : ranks, the scores of the modaltes of the bns. 1 PCA on G = (g ), g = B j b=1 s(b) p (b), where p (b) may be normalzed by a angular transformaton arcsn p (b). DR for Hstogram Data June 13, / 27

31 The Transformaton Approach (2/2) Quantfcaton: 1 The transformed varable: s = (s (1), s (2),, s (B j ) ) Output: 1 Codng 1: s (b) = (x (b) + x (b+1) )/2: center of the adjusted ntervals. 2 Codng 2: s (b) : ranks, the scores of the modaltes of the bns. 1 PCA on G = (g ), g = B j b=1 s(b) p (b), where p (b) may be normalzed by a angular transformaton arcsn p (b). Vsualzaton: 1 Use hypercube for representaton of concepts. 2 Use nterval lengths for representng the concepts. DR for Hstogram Data June 13, / 27

32 The Dstrbutonal Approach (1/1) Reference: 1 Bllard, L. and dday, E., (23), From the statstcs of data to the statstcs of knowledge: symbolc data analyss, Journal of the Amercan Statstcal Assocaton, 98(462), Wang, H., Chen, M., L, N., and Wang, L., (211), Prncpal component analyss of Modal nterval-valued data wth constant numercal characterstcs, Int. Statstcal Inst.: Proc. 58th World Statstcal Congress, 211, Dubln (Sesson CPS53), pp DR for Hstogram Data June 13, / 27

33 The Dstrbutonal Approach (1/1) Reference: 1 Bllard, L. and dday, E., (23), From the statstcs of data to the statstcs of knowledge: symbolc data analyss, Journal of the Amercan Statstcal Assocaton, 98(462), Wang, H., Chen, M., L, N., and Wang, L., (211), Prncpal component analyss of Modal nterval-valued data wth constant numercal characterstcs, Int. Statstcal Inst.: Proc. 58th World Statstcal Congress, 211, Dubln (Sesson CPS53), pp Assumpton: 1 Wthn each sub-nterval I (b), the random varable x s unformly dstrbuted across the sub-nterval. DR for Hstogram Data June 13, / 27

34 The Dstrbutonal Approach (2/1) The emprcal densty functon (Bllard and Dday, 23): 1 The emprcal densty functon for jth hstogram-valued varable X j : f j (u) = 1 n B n =1 b=1 x (b+1) p (b) x (b), x (b) u x (b+1), j = 1,, p. DR for Hstogram Data June 13, / 27

35 The Dstrbutonal Approach (2/1) The emprcal densty functon (Bllard and Dday, 23): 1 The emprcal densty functon for jth hstogram-valued varable X j : f j (u) = 1 n B n =1 b=1 x (b+1) p (b) x (b), x (b) u x (b+1), j = 1,, p. The symbolc sample mean and varance: 1 The symbolc sample mean for jth hstogram-valued varable X j : µ j = E(X j ) = = 1 n = 1 2n B n =1 b=1 =1 b=1 (b+1) x x (b) ξ f j (ξ) dξ ξ x (b+1) p (b) x (b) dξ B n p (b) (x (b+1) + x (b) ), j = 1,, p. DR for Hstogram Data June 13, / 27

36 The Dstrbutonal Approach (3/1) 2 The symbolc sample varance for jth hstogram-valued varable X j : Var(X j ) = E(Xj 2 ) µ 2 j B n = 1 3n =1 1 4n 2 b [ p (b) =1 b=1 (x (b+1) B n p (b) (x (b+1) + x (b) ) 2 + x (b+1) x (b) + (x (b) ) 2] ) 2, j = 1,, p. DR for Hstogram Data June 13, / 27

37 The Dstrbutonal Approach (4/1) The emprcal jont densty functon (Bllard and Dday, 23): 1 Assume that (X j, X j ) values are unformly dstrbuted across the rectangles (x (b j ), x (b j +1) ] (x (b j ), x (b j +1) ]. 2 Assume h and h are ndependent. DR for Hstogram Data June 13, / 27

38 The Dstrbutonal Approach (4/1) The emprcal jont densty functon (Bllard and Dday, 23): 1 Assume that (X j, X j ) values are unformly dstrbuted across the rectangles (x (b j ), x (b j +1) ] (x (b j ), x (b j +1) ]. 2 Assume h and h are ndependent. 3 The emprcal jont densty functon for jth and j th hstogram-valued varables (X j, X j ): f jj (u, v) = 1 n B n =1 B b j =1 b j =1 p (b j ) p (b j ) Z (b j,b j ), where Z (b j,b j ) s the area of the rectangle Z (b j,b j ) = (x (b j ), x (b j +1) ] (x (b j ), x (b j +1) ]. (u, v) Z (b j,b j ), DR for Hstogram Data June 13, / 27

39 The Dstrbutonal Approach (5/1) The symbolc sample covarance: 1 The symbolc sample covarance functon s Cov(X j, X j ) = 1 n B B 4n =1 b j =1 b j =1 1 n B 4n 2 =1 b j =1 B n =1 b j =1 p (b j ) p (b j ) (x (b j ) + x (b j +1) ) (x (b j ) + x (b j +1) p (b j ) (x (b j ) + x (b j +1) ) p (b j ) (x (b j ) + x (b j +1) ) ), j, j = 1,, p, j j. DR for Hstogram Data June 13, / 27

40 The Dstrbutonal Approach (6/1) Lnear combnaton of p ntervals Gven p nterval-valued varables X 1,, X p, all wth n observatons and real numbers α j, a lnear combnaton of p nterval-valued varables s defned as Y = p α j X j = ([y 1, y 1 ], [y 2, y 2 ],, [y n, y n ]) T, j=1 where y = p α j (τx + (1 τ)x ), and y = j=1 p α j ((1 τ)x + τx ), j=1 wth τ = I(α j > ). DR for Hstogram Data June 13, / 27

41 The Dstrbutonal Approach (7/1) Lnear combnaton of p hstograms Gven p hstogram varables X 1,, X p, and real numbers α j, defne a hstogram varable Y as a lnear combnaton of X 1,, X p Y = p α j X j = (y 1,, y n ) T, j=1 where y = {I, p } = {y I (b) where B = max j {B }. = (y (b), y (b+1) ], p (b), b = 1,, B }, DR for Hstogram Data June 13, / 27

42 The Dstrbutonal Approach (8/1) 1 For the th hstogram h, wth the number of modaltes B, each sub-nterval has ts densty functon p (b j ), b j = 1,, B. DR for Hstogram Data June 13, / 27

43 The Dstrbutonal Approach (8/1) 1 For the th hstogram h, wth the number of modaltes B, each sub-nterval has ts densty functon p (b j ), b j = 1,, B. 2 Hence, the th observaton, y, contans B p j=1 B hypercubes n p dmensonal space. DR for Hstogram Data June 13, / 27

44 The Dstrbutonal Approach (8/1) 1 For the th hstogram h, wth the number of modaltes B, each sub-nterval has ts densty functon p (b j ), b j = 1,, B. 2 Hence, the th observaton, y, contans B p j=1 B hypercubes n p dmensonal space. 3 Let p u denote the densty of each hypercube of the th observaton, that s the product of the densty of correspondng subntervals (I (b 1) 1,, I (bp) p ), u = 1,, B, p u p u(i (b 1) 1,, I (bp) p ) = p j=1 p (b j ), b j {1,, B }, j = 1,, p. DR for Hstogram Data June 13, / 27

45 The Dstrbutonal Approach (9/1) 4 Calculatng lnear combnaton to all the hypercubes, (I (b 1) 1,, I (bp) p ), we obtan I u, u = 1,, B : I u I u(i (b 1) 1,, I (bp) p ) = p j=1 α j I (b j ), b j {1,, B }, j = 1,, p. DR for Hstogram Data June 13, / 27

46 The Dstrbutonal Approach (9/1) 4 Calculatng lnear combnaton to all the hypercubes, (I (b 1) 1,, I (bp) p ), we obtan I u, u = 1,, B : I u I u(i (b 1) 1,, I (bp) p ) = p j=1 α j I (b j ), b j {1,, B }, j = 1,, p. 5 Let I = (y (1) mn u {I u }, y (B +1) max u {I u }]. The subnterval I (b) = (y (b), y (b+1) ], b = 1,, B, can be obtaned by dvdng I nto B parts: y (b) = y (1) + y (B +1) y (1) (b 1), b = 1,, B. B DR for Hstogram Data June 13, / 27

47 The Dstrbutonal Approach (1/1) 6 Projectng I u to I (b), we obtan p (b) = B u=1 p I u I (b) u I u, b = 1,, B, = 1,, n. DR for Hstogram Data June 13, / 27

48 The Dstrbutonal Approach (1/1) 6 Projectng I u to I (b), we obtan p (b) = B u=1 p I u I (b) u I u, b = 1,, B, = 1,, n. Vsualzaton: 1 none. (Wang, Chen, L, and Wang, 211) DR for Hstogram Data June 13, / 27

49 Vsualzaton: A hstogram of hstograms A hstogram of all of the observed hstograms (Bllard and Dday, 23): Let {(I (g), p (g) ), g = 1,, r} represent the relatve frequency hstogram for the combned set of observed hstograms. DR for Hstogram Data June 13, / 27

50 Vsualzaton: A hstogram of hstograms A hstogram of all of the observed hstograms (Bllard and Dday, 23): Let {(I (g), p (g) ), g = 1,, r} represent the relatve frequency hstogram for the combned set of observed hstograms. 1 Let I = (mn,b {x (b) }, max,b {x (b+1) }]: the nterval that spans all of the observed values of X. DR for Hstogram Data June 13, / 27

51 Vsualzaton: A hstogram of hstograms A hstogram of all of the observed hstograms (Bllard and Dday, 23): Let {(I (g), p (g) ), g = 1,, r} represent the relatve frequency hstogram for the combned set of observed hstograms. 1 Let I = (mn,b {x (b) }, max,b {x (b+1) }]: the nterval that spans all of the observed values of X. 2 Partton I nto r subntervals, I (g) = (ξ (g), ξ (g+1) ], g = 1,, r. DR for Hstogram Data June 13, / 27

52 Vsualzaton: A hstogram of hstograms A hstogram of all of the observed hstograms (Bllard and Dday, 23): Let {(I (g), p (g) ), g = 1,, r} represent the relatve frequency hstogram for the combned set of observed hstograms. 1 Let I = (mn,b {x (b) }, max,b {x (b+1) }]: the nterval that spans all of the observed values of X. 2 Partton I nto r subntervals, I (g) = (ξ (g), ξ (g+1) ], g = 1,, r. 3 The relatve observed frequency for the nterval I (g) s p (g) = 1 n n B =1 b=1 I (b) I (g) I (b) p (b), g = 1,, r, where I (b) (x (b), x (b+1) ]. DR for Hstogram Data June 13, / 27

53 Vsualzaton: a jont hstogram for Z 1 and Z 2 Let {R (g 1,g 2 ), p (g 1,g 2 ) } be the jont hstogram for two hstogram-valued varables, Z 1 and Z 2 (Bllard and Dday, 23): DR for Hstogram Data June 13, / 27

54 Vsualzaton: a jont hstogram for Z 1 and Z 2 Let {R (g 1,g 2 ), p (g 1,g 2 ) } be the jont hstogram for two hstogram-valued varables, Z 1 and Z 2 (Bllard and Dday, 23): 1 R (g 1,g 2 ) (ξ (g 1) 1, ξ (g 1+1) 1 ] (ξ (g 2) 2, ξ (g 2+1) 2 ], g 1 = 1,, r 1, g 2 = 1,, r 2. DR for Hstogram Data June 13, / 27

55 Vsualzaton: a jont hstogram for Z 1 and Z 2 Let {R (g 1,g 2 ), p (g 1,g 2 ) } be the jont hstogram for two hstogram-valued varables, Z 1 and Z 2 (Bllard and Dday, 23): 1 R (g 1,g 2 ) (ξ (g 1) 1, ξ (g 1+1) 1 ] (ξ (g 2) 2, ξ (g 2+1) 2 ], g 1 = 1,, r 1, g 2 = 1,, r 2. 2 p (g 1,g 2 ) = 1 n n B 1 B 2 =1 b 1 =1 b 2 =1 Π (b 1,b 2 ) R (g 1,g 2 ) Π (b 1,b 2 ) where Π (b 1,b 2 ) (x (b 1) 1, x (b 1+1) 1 ] (x (b 2) 2, x (b 2+1) 2 ]. p (b 1) 1 p (b 2) 2, DR for Hstogram Data June 13, / 27

56 Example: rs data (15 4, y=(5, 5, 5)) PCA of rs data PCA Setosa Verscolor Vrgnca PCA1 DR for Hstogram Data June 13, / 27

57 Example: rs data (15 4, y=(5, 5, 5)) PCA of rs data PCA Setosa Verscolor Vrgnca PCA1 Repeat (Wthn each class, 3 observatons were randomly sampled to generate hstograms for four varables) 5 tmes. DR for Hstogram Data June 13, / 27

58 Hstogram-PCA..3 DR for Hstogram Data June 13, / 27

59 Jont hstogram of hstogram-pca Jont Hstogram of Hstograms: ( DR1, DR2) [1.9,2.3) [2.3,2.8) [2.8,3.3) [3.3,3.8) [3.8,4.2) [4.2,4.7) [4.7,5.2) [5.2,5.7) [5.7,6.2) [1.9,2.7) [3.6,4.4) [5.3,6.1) [6.9,7.8) [8.6,9.5) DR for Hstogram Data June 13, / 27

60 Jont hstogram for each speces Jont Hstogram of Hstograms: ( DR1, DR2) [1.9,2.3) [2.3,2.8) [2.8,3.3) [3.3,3.8) [3.8,4.2) [4.2,4.7) [4.7,5.2) [5.2,5.7) [5.7,6.2) [1.9,2.7) [3.6,4.4) [5.3,6.1) [6.9,7.8) [8.6,9.5) DR for Hstogram Data June 13, / 27

61 Jont hstogram for each speces Jont Hstogram of Hstograms: ( DR1, DR2) [1.9,2.3) [2.3,2.8) [2.8,3.3) [3.3,3.8) [3.8,4.2) [4.2,4.7) [4.7,5.2) [5.2,5.7) [5.7,6.2) [1.9,2.7) [3.6,4.4) [5.3,6.1) [6.9,7.8) [8.6,9.5) Jont Hstogram of Hstograms: ( DR1, DR2), cluster=1 [1.9,2.3) [2.3,2.8) [2.8,3.3) [3.3,3.8).8.4 [3.8,4.2) [4.2,4.7) [4.7,5.2) [5.2,5.7).2.1 [5.7,6.2) [1.9,2.7) [3.6,4.4) [5.3,6.1) [6.9,7.8) [8.6,9.5) DR for Hstogram Data June 13, / 27

62 Jont hstogram for each speces Jont Hstogram of Hstograms: ( DR1, DR2) [1.9,2.3) [2.3,2.8) [2.8,3.3) [3.3,3.8) [3.8,4.2) [4.2,4.7) [4.7,5.2) [5.2,5.7) [5.7,6.2) [1.9,2.7) [3.6,4.4) [5.3,6.1) [6.9,7.8) [8.6,9.5) Jont Hstogram of Hstograms: ( DR1, DR2), cluster=1 [1.9,2.3) [2.3,2.8) [2.8,3.3) [3.3,3.8).8.4 [3.8,4.2) [4.2,4.7) [4.7,5.2) [5.2,5.7).2.1 [5.7,6.2) [1.9,2.7) [3.6,4.4) [5.3,6.1) [6.9,7.8) [8.6,9.5) Jont Hstogram of Hstograms: ( DR1, DR2), cluster=2 [1.9,2.3) [2.3,2.8).1.1 [2.8,3.3) [3.3,3.8) [3.8,4.2) [4.2,4.7) [4.7,5.2) [5.2,5.7) [5.7,6.2) [1.9,2.7) [3.6,4.4) [5.3,6.1) [6.9,7.8) [8.6,9.5) DR for Hstogram Data June 13, / 27

63 Jont hstogram for each speces Jont Hstogram of Hstograms: ( DR1, DR2) [1.9,2.3) [2.3,2.8) [2.8,3.3) [3.3,3.8) [3.8,4.2) [4.2,4.7) [4.7,5.2) [5.2,5.7) [5.7,6.2) [1.9,2.7) [3.6,4.4) [5.3,6.1) [6.9,7.8) [8.6,9.5) Jont Hstogram of Hstograms: ( DR1, DR2), cluster=1 [1.9,2.3) [2.3,2.8) [2.8,3.3) [3.3,3.8).8.4 [3.8,4.2) [4.2,4.7) [4.7,5.2) [5.2,5.7).2.1 [5.7,6.2) [1.9,2.7) [3.6,4.4) [5.3,6.1) [6.9,7.8) [8.6,9.5) Jont Hstogram of Hstograms: ( DR1, DR2), cluster=2 Jont Hstogram of Hstograms: ( DR1, DR2), cluster=3 [1.9,2.3) [1.9,2.3) [2.3,2.8).1.1 [2.8,3.3) [3.3,3.8) [3.8,4.2) [4.2,4.7) [4.7,5.2) [5.2,5.7) [5.7,6.2) [1.9,2.7) [3.6,4.4) [5.3,6.1) [6.9,7.8) [8.6,9.5) [2.3,2.8).1.1 [2.8,3.3) [3.3,3.8) [3.8,4.2) [4.2,4.7) [4.7,5.2) [5.2,5.7) [5.7,6.2) [1.9,2.7) [3.6,4.4) [5.3,6.1) [6.9,7.8) [8.6,9.5) DR for Hstogram Data June 13, /

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family IOSR Journal of Mathematcs IOSR-JM) ISSN: 2278-5728. Volume 3, Issue 3 Sep-Oct. 202), PP 44-48 www.osrjournals.org Usng T.O.M to Estmate Parameter of dstrbutons that have not Sngle Exponental Famly Jubran

More information

Statistical analysis using matlab. HY 439 Presented by: George Fortetsanakis

Statistical analysis using matlab. HY 439 Presented by: George Fortetsanakis Statstcal analyss usng matlab HY 439 Presented by: George Fortetsanaks Roadmap Probablty dstrbutons Statstcal estmaton Fttng data to probablty dstrbutons Contnuous dstrbutons Contnuous random varable X

More information

x = , so that calculated

x = , so that calculated Stat 4, secton Sngle Factor ANOVA notes by Tm Plachowsk n chapter 8 we conducted hypothess tests n whch we compared a sngle sample s mean or proporton to some hypotheszed value Chapter 9 expanded ths to

More information

Statistical Inference. 2.3 Summary Statistics Measures of Center and Spread. parameters ( population characteristics )

Statistical Inference. 2.3 Summary Statistics Measures of Center and Spread. parameters ( population characteristics ) Ismor Fscher, 8//008 Stat 54 / -8.3 Summary Statstcs Measures of Center and Spread Dstrbuton of dscrete contnuous POPULATION Random Varable, numercal True center =??? True spread =???? parameters ( populaton

More information

Maximum Likelihood Estimation of Binary Dependent Variables Models: Probit and Logit. 1. General Formulation of Binary Dependent Variables Models

Maximum Likelihood Estimation of Binary Dependent Variables Models: Probit and Logit. 1. General Formulation of Binary Dependent Variables Models ECO 452 -- OE 4: Probt and Logt Models ECO 452 -- OE 4 Maxmum Lkelhood Estmaton of Bnary Dependent Varables Models: Probt and Logt hs note demonstrates how to formulate bnary dependent varables models

More information

Statistics for Economics & Business

Statistics for Economics & Business Statstcs for Economcs & Busness Smple Lnear Regresson Learnng Objectves In ths chapter, you learn: How to use regresson analyss to predct the value of a dependent varable based on an ndependent varable

More information

Statistics for Business and Economics

Statistics for Business and Economics Statstcs for Busness and Economcs Chapter 11 Smple Regresson Copyrght 010 Pearson Educaton, Inc. Publshng as Prentce Hall Ch. 11-1 11.1 Overvew of Lnear Models n An equaton can be ft to show the best lnear

More information

Chat eld, C. and A.J.Collins, Introduction to multivariate analysis. Chapman & Hall, 1980

Chat eld, C. and A.J.Collins, Introduction to multivariate analysis. Chapman & Hall, 1980 MT07: Multvarate Statstcal Methods Mke Tso: emal mke.tso@manchester.ac.uk Webpage for notes: http://www.maths.manchester.ac.uk/~mkt/new_teachng.htm. Introducton to multvarate data. Books Chat eld, C. and

More information

Statistics for Managers Using Microsoft Excel/SPSS Chapter 13 The Simple Linear Regression Model and Correlation

Statistics for Managers Using Microsoft Excel/SPSS Chapter 13 The Simple Linear Regression Model and Correlation Statstcs for Managers Usng Mcrosoft Excel/SPSS Chapter 13 The Smple Lnear Regresson Model and Correlaton 1999 Prentce-Hall, Inc. Chap. 13-1 Chapter Topcs Types of Regresson Models Determnng the Smple Lnear

More information

j) = 1 (note sigma notation) ii. Continuous random variable (e.g. Normal distribution) 1. density function: f ( x) 0 and f ( x) dx = 1

j) = 1 (note sigma notation) ii. Continuous random variable (e.g. Normal distribution) 1. density function: f ( x) 0 and f ( x) dx = 1 Random varables Measure of central tendences and varablty (means and varances) Jont densty functons and ndependence Measures of assocaton (covarance and correlaton) Interestng result Condtonal dstrbutons

More information

Chapter 3 Describing Data Using Numerical Measures

Chapter 3 Describing Data Using Numerical Measures Chapter 3 Student Lecture Notes 3-1 Chapter 3 Descrbng Data Usng Numercal Measures Fall 2006 Fundamentals of Busness Statstcs 1 Chapter Goals To establsh the usefulness of summary measures of data. The

More information

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction ECONOMICS 5* -- NOTE (Summary) ECON 5* -- NOTE The Multple Classcal Lnear Regresson Model (CLRM): Specfcaton and Assumptons. Introducton CLRM stands for the Classcal Lnear Regresson Model. The CLRM s also

More information

The Ordinary Least Squares (OLS) Estimator

The Ordinary Least Squares (OLS) Estimator The Ordnary Least Squares (OLS) Estmator 1 Regresson Analyss Regresson Analyss: a statstcal technque for nvestgatng and modelng the relatonshp between varables. Applcatons: Engneerng, the physcal and chemcal

More information

Uncertainty as the Overlap of Alternate Conditional Distributions

Uncertainty as the Overlap of Alternate Conditional Distributions Uncertanty as the Overlap of Alternate Condtonal Dstrbutons Olena Babak and Clayton V. Deutsch Centre for Computatonal Geostatstcs Department of Cvl & Envronmental Engneerng Unversty of Alberta An mportant

More information

RELIABILITY ASSESSMENT

RELIABILITY ASSESSMENT CHAPTER Rsk Analyss n Engneerng and Economcs RELIABILITY ASSESSMENT A. J. Clark School of Engneerng Department of Cvl and Envronmental Engneerng 4a CHAPMAN HALL/CRC Rsk Analyss for Engneerng Department

More information

Int. Statistical Inst.: Proc. 58th World Statistical Congress, 2011, Dublin (Session CPS053) p.5185

Int. Statistical Inst.: Proc. 58th World Statistical Congress, 2011, Dublin (Session CPS053) p.5185 Int. Statstcal Inst.: Proc. 58th World Statstcal Congress, 0, Dubln (Sesson CPS05).585 Prncal Coonent Analyss of Modal nterval-valued Data wth Constant Nuercal Characterstcs WANG Huwen,, CHEN Melng,, LI

More information

LINEAR REGRESSION ANALYSIS. MODULE VIII Lecture Indicator Variables

LINEAR REGRESSION ANALYSIS. MODULE VIII Lecture Indicator Variables LINEAR REGRESSION ANALYSIS MODULE VIII Lecture - 7 Indcator Varables Dr. Shalabh Department of Maematcs and Statstcs Indan Insttute of Technology Kanpur Indcator varables versus quanttatve explanatory

More information

Maximum Likelihood Estimation of Binary Dependent Variables Models: Probit and Logit. 1. General Formulation of Binary Dependent Variables Models

Maximum Likelihood Estimation of Binary Dependent Variables Models: Probit and Logit. 1. General Formulation of Binary Dependent Variables Models ECO 452 -- OE 4: Probt and Logt Models ECO 452 -- OE 4 Mamum Lkelhood Estmaton of Bnary Dependent Varables Models: Probt and Logt hs note demonstrates how to formulate bnary dependent varables models for

More information

Department of Quantitative Methods & Information Systems. Time Series and Their Components QMIS 320. Chapter 6

Department of Quantitative Methods & Information Systems. Time Series and Their Components QMIS 320. Chapter 6 Department of Quanttatve Methods & Informaton Systems Tme Seres and Ther Components QMIS 30 Chapter 6 Fall 00 Dr. Mohammad Zanal These sldes were modfed from ther orgnal source for educatonal purpose only.

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

Composite Hypotheses testing

Composite Hypotheses testing Composte ypotheses testng In many hypothess testng problems there are many possble dstrbutons that can occur under each of the hypotheses. The output of the source s a set of parameters (ponts n a parameter

More information

U-Pb Geochronology Practical: Background

U-Pb Geochronology Practical: Background U-Pb Geochronology Practcal: Background Basc Concepts: accuracy: measure of the dfference between an expermental measurement and the true value precson: measure of the reproducblty of the expermental result

More information

A be a probability space. A random vector

A be a probability space. A random vector Statstcs 1: Probablty Theory II 8 1 JOINT AND MARGINAL DISTRIBUTIONS In Probablty Theory I we formulate the concept of a (real) random varable and descrbe the probablstc behavor of ths random varable by

More information

Modeling and Simulation NETW 707

Modeling and Simulation NETW 707 Modelng and Smulaton NETW 707 Lecture 5 Tests for Random Numbers Course Instructor: Dr.-Ing. Magge Mashaly magge.ezzat@guc.edu.eg C3.220 1 Propertes of Random Numbers Random Number Generators (RNGs) must

More information

Multi-dimensional Central Limit Theorem

Multi-dimensional Central Limit Theorem Mult-dmensonal Central Lmt heorem Outlne ( ( ( t as ( + ( + + ( ( ( Consder a sequence of ndependent random proceses t, t, dentcal to some ( t. Assume t = 0. Defne the sum process t t t t = ( t = (; t

More information

Lecture 3 Stat102, Spring 2007

Lecture 3 Stat102, Spring 2007 Lecture 3 Stat0, Sprng 007 Chapter 3. 3.: Introducton to regresson analyss Lnear regresson as a descrptve technque The least-squares equatons Chapter 3.3 Samplng dstrbuton of b 0, b. Contnued n net lecture

More information

Probability Theory. The nth coefficient of the Taylor series of f(k), expanded around k = 0, gives the nth moment of x as ( ik) n n!

Probability Theory. The nth coefficient of the Taylor series of f(k), expanded around k = 0, gives the nth moment of x as ( ik) n n! 8333: Statstcal Mechancs I Problem Set # 3 Solutons Fall 3 Characterstc Functons: Probablty Theory The characterstc functon s defned by fk ep k = ep kpd The nth coeffcent of the Taylor seres of fk epanded

More information

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:

More information

x i1 =1 for all i (the constant ).

x i1 =1 for all i (the constant ). Chapter 5 The Multple Regresson Model Consder an economc model where the dependent varable s a functon of K explanatory varables. The economc model has the form: y = f ( x,x,..., ) xk Approxmate ths by

More information

Strong Markov property: Same assertion holds for stopping times τ.

Strong Markov property: Same assertion holds for stopping times τ. Brownan moton Let X ={X t : t R + } be a real-valued stochastc process: a famlty of real random varables all defned on the same probablty space. Defne F t = nformaton avalable by observng the process up

More information

Computation of Higher Order Moments from Two Multinomial Overdispersion Likelihood Models

Computation of Higher Order Moments from Two Multinomial Overdispersion Likelihood Models Computaton of Hgher Order Moments from Two Multnomal Overdsperson Lkelhood Models BY J. T. NEWCOMER, N. K. NEERCHAL Department of Mathematcs and Statstcs, Unversty of Maryland, Baltmore County, Baltmore,

More information

UNIVERSITY OF TORONTO Faculty of Arts and Science. December 2005 Examinations STA437H1F/STA1005HF. Duration - 3 hours

UNIVERSITY OF TORONTO Faculty of Arts and Science. December 2005 Examinations STA437H1F/STA1005HF. Duration - 3 hours UNIVERSITY OF TORONTO Faculty of Arts and Scence December 005 Examnatons STA47HF/STA005HF Duraton - hours AIDS ALLOWED: (to be suppled by the student) Non-programmable calculator One handwrtten 8.5'' x

More information

/ n ) are compared. The logic is: if the two

/ n ) are compared. The logic is: if the two STAT C141, Sprng 2005 Lecture 13 Two sample tests One sample tests: examples of goodness of ft tests, where we are testng whether our data supports predctons. Two sample tests: called as tests of ndependence

More information

Probability Theory (revisited)

Probability Theory (revisited) Probablty Theory (revsted) Summary Probablty v.s. plausblty Random varables Smulaton of Random Experments Challenge The alarm of a shop rang. Soon afterwards, a man was seen runnng n the street, persecuted

More information

Lecture 6 More on Complete Randomized Block Design (RBD)

Lecture 6 More on Complete Randomized Block Design (RBD) Lecture 6 More on Complete Randomzed Block Desgn (RBD) Multple test Multple test The multple comparsons or multple testng problem occurs when one consders a set of statstcal nferences smultaneously. For

More information

Lecture 16 Statistical Analysis in Biomaterials Research (Part II)

Lecture 16 Statistical Analysis in Biomaterials Research (Part II) 3.051J/0.340J 1 Lecture 16 Statstcal Analyss n Bomaterals Research (Part II) C. F Dstrbuton Allows comparson of varablty of behavor between populatons usng test of hypothess: σ x = σ x amed for Brtsh statstcan

More information

4 Analysis of Variance (ANOVA) 5 ANOVA. 5.1 Introduction. 5.2 Fixed Effects ANOVA

4 Analysis of Variance (ANOVA) 5 ANOVA. 5.1 Introduction. 5.2 Fixed Effects ANOVA 4 Analyss of Varance (ANOVA) 5 ANOVA 51 Introducton ANOVA ANOVA s a way to estmate and test the means of multple populatons We wll start wth one-way ANOVA If the populatons ncluded n the study are selected

More information

LECTURE 9 CANONICAL CORRELATION ANALYSIS

LECTURE 9 CANONICAL CORRELATION ANALYSIS LECURE 9 CANONICAL CORRELAION ANALYSIS Introducton he concept of canoncal correlaton arses when we want to quantfy the assocatons between two sets of varables. For example, suppose that the frst set of

More information

First Year Examination Department of Statistics, University of Florida

First Year Examination Department of Statistics, University of Florida Frst Year Examnaton Department of Statstcs, Unversty of Florda May 7, 010, 8:00 am - 1:00 noon Instructons: 1. You have four hours to answer questons n ths examnaton.. You must show your work to receve

More information

Basic Business Statistics, 10/e

Basic Business Statistics, 10/e Chapter 13 13-1 Basc Busness Statstcs 11 th Edton Chapter 13 Smple Lnear Regresson Basc Busness Statstcs, 11e 009 Prentce-Hall, Inc. Chap 13-1 Learnng Objectves In ths chapter, you learn: How to use regresson

More information

Credit Card Pricing and Impact of Adverse Selection

Credit Card Pricing and Impact of Adverse Selection Credt Card Prcng and Impact of Adverse Selecton Bo Huang and Lyn C. Thomas Unversty of Southampton Contents Background Aucton model of credt card solctaton - Errors n probablty of beng Good - Errors n

More information

Chapter 8 Indicator Variables

Chapter 8 Indicator Variables Chapter 8 Indcator Varables In general, e explanatory varables n any regresson analyss are assumed to be quanttatve n nature. For example, e varables lke temperature, dstance, age etc. are quanttatve n

More information

Primer on High-Order Moment Estimators

Primer on High-Order Moment Estimators Prmer on Hgh-Order Moment Estmators Ton M. Whted July 2007 The Errors-n-Varables Model We wll start wth the classcal EIV for one msmeasured regressor. The general case s n Erckson and Whted Econometrc

More information

Feb 14: Spatial analysis of data fields

Feb 14: Spatial analysis of data fields Feb 4: Spatal analyss of data felds Mappng rregularly sampled data onto a regular grd Many analyss technques for geophyscal data requre the data be located at regular ntervals n space and/or tme. hs s

More information

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4) I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes

More information

MULTISPECTRAL IMAGE CLASSIFICATION USING BACK-PROPAGATION NEURAL NETWORK IN PCA DOMAIN

MULTISPECTRAL IMAGE CLASSIFICATION USING BACK-PROPAGATION NEURAL NETWORK IN PCA DOMAIN MULTISPECTRAL IMAGE CLASSIFICATION USING BACK-PROPAGATION NEURAL NETWORK IN PCA DOMAIN S. Chtwong, S. Wtthayapradt, S. Intajag, and F. Cheevasuvt Faculty of Engneerng, Kng Mongkut s Insttute of Technology

More information

Regularized Discriminant Analysis for Face Recognition

Regularized Discriminant Analysis for Face Recognition 1 Regularzed Dscrmnant Analyss for Face Recognton Itz Pma, Mayer Aladem Department of Electrcal and Computer Engneerng, Ben-Guron Unversty of the Negev P.O.Box 653, Beer-Sheva, 845, Israel. Abstract Ths

More information

: 5: ) A

: 5: ) A Revew 1 004.11.11 Chapter 1: 1. Elements, Varable, and Observatons:. Type o Data: Qualtatve Data and Quanttatve Data (a) Qualtatve data may be nonnumerc or numerc. (b) Quanttatve data are always numerc.

More information

Predictive Analytics : QM901.1x Prof U Dinesh Kumar, IIMB. All Rights Reserved, Indian Institute of Management Bangalore

Predictive Analytics : QM901.1x Prof U Dinesh Kumar, IIMB. All Rights Reserved, Indian Institute of Management Bangalore Sesson Outlne Introducton to classfcaton problems and dscrete choce models. Introducton to Logstcs Regresson. Logstc functon and Logt functon. Maxmum Lkelhood Estmator (MLE) for estmaton of LR parameters.

More information

P R. Lecture 4. Theory and Applications of Pattern Recognition. Dept. of Electrical and Computer Engineering /

P R. Lecture 4. Theory and Applications of Pattern Recognition. Dept. of Electrical and Computer Engineering / Theory and Applcatons of Pattern Recognton 003, Rob Polkar, Rowan Unversty, Glassboro, NJ Lecture 4 Bayes Classfcaton Rule Dept. of Electrcal and Computer Engneerng 0909.40.0 / 0909.504.04 Theory & Applcatons

More information

8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS

8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS SECTION 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS 493 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS All the vector spaces you have studed thus far n the text are real vector spaces because the scalars

More information

Statistics Chapter 4

Statistics Chapter 4 Statstcs Chapter 4 "There are three knds of les: les, damned les, and statstcs." Benjamn Dsrael, 1895 (Brtsh statesman) Gaussan Dstrbuton, 4-1 If a measurement s repeated many tmes a statstcal treatment

More information

Here is the rationale: If X and y have a strong positive relationship to one another, then ( x x) will tend to be positive when ( y y)

Here is the rationale: If X and y have a strong positive relationship to one another, then ( x x) will tend to be positive when ( y y) Secton 1.5 Correlaton In the prevous sectons, we looked at regresson and the value r was a measurement of how much of the varaton n y can be attrbuted to the lnear relatonshp between y and x. In ths secton,

More information

2016 Wiley. Study Session 2: Ethical and Professional Standards Application

2016 Wiley. Study Session 2: Ethical and Professional Standards Application 6 Wley Study Sesson : Ethcal and Professonal Standards Applcaton LESSON : CORRECTION ANALYSIS Readng 9: Correlaton and Regresson LOS 9a: Calculate and nterpret a sample covarance and a sample correlaton

More information

Chapter 11: Simple Linear Regression and Correlation

Chapter 11: Simple Linear Regression and Correlation Chapter 11: Smple Lnear Regresson and Correlaton 11-1 Emprcal Models 11-2 Smple Lnear Regresson 11-3 Propertes of the Least Squares Estmators 11-4 Hypothess Test n Smple Lnear Regresson 11-4.1 Use of t-tests

More information

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur Analyss of Varance and Desgn of Experment-I MODULE VII LECTURE - 3 ANALYSIS OF COVARIANCE Dr Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur Any scentfc experment s performed

More information

PHYS 450 Spring semester Lecture 02: Dealing with Experimental Uncertainties. Ron Reifenberger Birck Nanotechnology Center Purdue University

PHYS 450 Spring semester Lecture 02: Dealing with Experimental Uncertainties. Ron Reifenberger Birck Nanotechnology Center Purdue University PHYS 45 Sprng semester 7 Lecture : Dealng wth Expermental Uncertantes Ron Refenberger Brck anotechnology Center Purdue Unversty Lecture Introductory Comments Expermental errors (really expermental uncertantes)

More information

Why Monte Carlo Integration? Introduction to Monte Carlo Method. Continuous Probability. Continuous Probability

Why Monte Carlo Integration? Introduction to Monte Carlo Method. Continuous Probability. Continuous Probability Introducton to Monte Carlo Method Kad Bouatouch IRISA Emal: kad@rsa.fr Wh Monte Carlo Integraton? To generate realstc lookng mages, we need to solve ntegrals of or hgher dmenson Pel flterng and lens smulaton

More information

Bivariate copulas on the exponentially weighted moving average control chart

Bivariate copulas on the exponentially weighted moving average control chart Songklanakarn J. Sc. Technol. 38 (5), 569-574, Sep. - Oct. 6 http://www.sjst.psu.ac.th Orgnal Artcle Bvarate copulas on the exponentally weghted movng average control chart Sasgarn Kuvattana, Pyapatr Busababodhn,

More information

Fall 2012 Analysis of Experimental Measurements B. Eisenstein/rev. S. Errede. . For P such independent random variables (aka degrees of freedom): 1 =

Fall 2012 Analysis of Experimental Measurements B. Eisenstein/rev. S. Errede. . For P such independent random variables (aka degrees of freedom): 1 = Fall Analss of Epermental Measurements B. Esensten/rev. S. Errede More on : The dstrbuton s the.d.f. for a (normalzed sum of squares of ndependent random varables, each one of whch s dstrbuted as N (,.

More information

Spatial Statistics and Analysis Methods (for GEOG 104 class).

Spatial Statistics and Analysis Methods (for GEOG 104 class). Spatal Statstcs and Analyss Methods (for GEOG 104 class). Provded by Dr. An L, San Dego State Unversty. 1 Ponts Types of spatal data Pont pattern analyss (PPA; such as nearest neghbor dstance, quadrat

More information

Statistical tools to perform Sensitivity Analysis in the Context of the Evaluation of Measurement Uncertainty

Statistical tools to perform Sensitivity Analysis in the Context of the Evaluation of Measurement Uncertainty Statstcal tools to perform Senstvty Analyss n the Contet of the Evaluaton of Measurement Uncertanty N. Fscher, A. Allard Laboratore natonal de métrologe et d essas (LNE) MATHMET PTB Berln nd June Outlne

More information

Homogenization of reaction-diffusion processes in a two-component porous medium with a non-linear flux-condition on the interface

Homogenization of reaction-diffusion processes in a two-component porous medium with a non-linear flux-condition on the interface Homogenzaton of reacton-dffuson processes n a two-component porous medum wth a non-lnear flux-condton on the nterface Internatonal Conference on Numercal and Mathematcal Modelng of Flow and Transport n

More information

Chapter 9: Statistical Inference and the Relationship between Two Variables

Chapter 9: Statistical Inference and the Relationship between Two Variables Chapter 9: Statstcal Inference and the Relatonshp between Two Varables Key Words The Regresson Model The Sample Regresson Equaton The Pearson Correlaton Coeffcent Learnng Outcomes After studyng ths chapter,

More information

Chapter 7 Channel Capacity and Coding

Chapter 7 Channel Capacity and Coding Chapter 7 Channel Capacty and Codng Contents 7. Channel models and channel capacty 7.. Channel models Bnary symmetrc channel Dscrete memoryless channels Dscrete-nput, contnuous-output channel Waveform

More information

Statistical inference for generalized Pareto distribution based on progressive Type-II censored data with random removals

Statistical inference for generalized Pareto distribution based on progressive Type-II censored data with random removals Internatonal Journal of Scentfc World, 2 1) 2014) 1-9 c Scence Publshng Corporaton www.scencepubco.com/ndex.php/ijsw do: 10.14419/jsw.v21.1780 Research Paper Statstcal nference for generalzed Pareto dstrbuton

More information

Non-Mixture Cure Model for Interval Censored Data: Simulation Study ABSTRACT

Non-Mixture Cure Model for Interval Censored Data: Simulation Study ABSTRACT Malaysan Journal of Mathematcal Scences 8(S): 37-44 (2014) Specal Issue: Internatonal Conference on Mathematcal Scences and Statstcs 2013 (ICMSS2013) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES Journal

More information

Statistical tables are provided Two Hours UNIVERSITY OF MANCHESTER. Date: Wednesday 4 th June 2008 Time: 1400 to 1600

Statistical tables are provided Two Hours UNIVERSITY OF MANCHESTER. Date: Wednesday 4 th June 2008 Time: 1400 to 1600 Statstcal tables are provded Two Hours UNIVERSITY OF MNCHESTER Medcal Statstcs Date: Wednesday 4 th June 008 Tme: 1400 to 1600 MT3807 Electronc calculators may be used provded that they conform to Unversty

More information

Some basic statistics and curve fitting techniques

Some basic statistics and curve fitting techniques Some basc statstcs and curve fttng technques Statstcs s the dscplne concerned wth the study of varablty, wth the study of uncertanty, and wth the study of decsonmakng n the face of uncertanty (Lndsay et

More information

Conjugacy and the Exponential Family

Conjugacy and the Exponential Family CS281B/Stat241B: Advanced Topcs n Learnng & Decson Makng Conjugacy and the Exponental Famly Lecturer: Mchael I. Jordan Scrbes: Bran Mlch 1 Conjugacy In the prevous lecture, we saw conjugate prors for the

More information

Testing for seasonal unit roots in heterogeneous panels

Testing for seasonal unit roots in heterogeneous panels Testng for seasonal unt roots n heterogeneous panels Jesus Otero * Facultad de Economía Unversdad del Rosaro, Colomba Jeremy Smth Department of Economcs Unversty of arwck Monca Gulett Aston Busness School

More information

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 30 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 2 Remedes for multcollnearty Varous technques have

More information

Association for the Chi-square Test

Association for the Chi-square Test Assocaton for the Ch-square Test Davd J Olve Southern Illnos Unversty February 8, 2012 Abstract A problem wth measures of assocaton for the ch-square test s that the measures depend on the number of observatons

More information

A Hybrid Variational Iteration Method for Blasius Equation

A Hybrid Variational Iteration Method for Blasius Equation Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method

More information

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,

More information

1. Inference on Regression Parameters a. Finding Mean, s.d and covariance amongst estimates. 2. Confidence Intervals and Working Hotelling Bands

1. Inference on Regression Parameters a. Finding Mean, s.d and covariance amongst estimates. 2. Confidence Intervals and Working Hotelling Bands Content. Inference on Regresson Parameters a. Fndng Mean, s.d and covarance amongst estmates.. Confdence Intervals and Workng Hotellng Bands 3. Cochran s Theorem 4. General Lnear Testng 5. Measures of

More information

Mining Data Streams-Estimating Frequency Moment

Mining Data Streams-Estimating Frequency Moment Mnng Data Streams-Estmatng Frequency Moment Barna Saha October 26, 2017 Frequency Moment Computng moments nvolves dstrbuton of frequences of dfferent elements n the stream. Frequency Moment Computng moments

More information

IRO0140 Advanced space time-frequency signal processing

IRO0140 Advanced space time-frequency signal processing IRO4 Advanced space tme-frequency sgnal processng Lecture Toomas Ruuben Takng nto account propertes of the sgnals, we can group these as followng: Regular and random sgnals (are all sgnal parameters determned

More information

AS-Level Maths: Statistics 1 for Edexcel

AS-Level Maths: Statistics 1 for Edexcel 1 of 6 AS-Level Maths: Statstcs 1 for Edecel S1. Calculatng means and standard devatons Ths con ndcates the slde contans actvtes created n Flash. These actvtes are not edtable. For more detaled nstructons,

More information

Tracking with Kalman Filter

Tracking with Kalman Filter Trackng wth Kalman Flter Scott T. Acton Vrgna Image and Vdeo Analyss (VIVA), Charles L. Brown Department of Electrcal and Computer Engneerng Department of Bomedcal Engneerng Unversty of Vrgna, Charlottesvlle,

More information

Statistics and Probability Theory in Civil, Surveying and Environmental Engineering

Statistics and Probability Theory in Civil, Surveying and Environmental Engineering Statstcs and Probablty Theory n Cvl, Surveyng and Envronmental Engneerng Pro. Dr. Mchael Havbro Faber ETH Zurch, Swtzerland Contents o Todays Lecture Overvew o Uncertanty Modelng Random Varables - propertes

More information

EIGENVALUE DISTRIBUTION OF THE CORRELATION MATRIX IN L`-FILTERS. Constantine Kotropoulos and Ioannis Pitas

EIGENVALUE DISTRIBUTION OF THE CORRELATION MATRIX IN L`-FILTERS. Constantine Kotropoulos and Ioannis Pitas EIGEVALUE DISTRIBUTIO OF THE CORRELATIO MATRI I L`-FILTERS Constantne Kotropoulos and Ioanns Ptas Dept of Informatcs, Arstotle Unversty of Thessalonk Box, Thessalonk 0 0, GREECE fcostas, ptas g@zeuscsdauthgr

More information

Comparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method

Comparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method Appled Mathematcal Scences, Vol. 7, 0, no. 47, 07-0 HIARI Ltd, www.m-hkar.com Comparson of the Populaton Varance Estmators of -Parameter Exponental Dstrbuton Based on Multple Crtera Decson Makng Method

More information

be a second-order and mean-value-zero vector-valued process, i.e., for t E

be a second-order and mean-value-zero vector-valued process, i.e., for t E CONFERENCE REPORT 617 DSCUSSON OF TWO PROCEDURES FOR EXPANDNG A VECTOR-VALUED STOCHASTC PROCESS N AN ORTHONORMAL WAY by R. GUTkRREZ and M. J. VALDERRAMA 1. ntroducton Snce K. Karhunen [l] and M. Lo&e [2]

More information

Econometrics: What's It All About, Alfie?

Econometrics: What's It All About, Alfie? ECON 351* -- Introducton (Page 1) Econometrcs: What's It All About, Ale? Usng sample data on observable varables to learn about economc relatonshps, the unctonal relatonshps among economc varables. Econometrcs

More information

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests Smulated of the Cramér-von Mses Goodness-of-Ft Tests Steele, M., Chaselng, J. and 3 Hurst, C. School of Mathematcal and Physcal Scences, James Cook Unversty, Australan School of Envronmental Studes, Grffth

More information

MDL-Based Unsupervised Attribute Ranking

MDL-Based Unsupervised Attribute Ranking MDL-Based Unsupervsed Attrbute Rankng Zdravko Markov Computer Scence Department Central Connectcut State Unversty New Brtan, CT 06050, USA http://www.cs.ccsu.edu/~markov/ markovz@ccsu.edu MDL-Based Unsupervsed

More information

PubH 7405: REGRESSION ANALYSIS. SLR: INFERENCES, Part II

PubH 7405: REGRESSION ANALYSIS. SLR: INFERENCES, Part II PubH 7405: REGRESSION ANALSIS SLR: INFERENCES, Part II We cover te topc of nference n two sessons; te frst sesson focused on nferences concernng te slope and te ntercept; ts s a contnuaton on estmatng

More information

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle

More information

Feature Selection & Dynamic Tracking F&P Textbook New: Ch 11, Old: Ch 17 Guido Gerig CS 6320, Spring 2013

Feature Selection & Dynamic Tracking F&P Textbook New: Ch 11, Old: Ch 17 Guido Gerig CS 6320, Spring 2013 Feature Selecton & Dynamc Trackng F&P Textbook New: Ch 11, Old: Ch 17 Gudo Gerg CS 6320, Sprng 2013 Credts: Materal Greg Welch & Gary Bshop, UNC Chapel Hll, some sldes modfed from J.M. Frahm/ M. Pollefeys,

More information

SUPPLEMENTARY INFORMATION

SUPPLEMENTARY INFORMATION do: 0.08/nature09 I. Resonant absorpton of XUV pulses n Kr + usng the reduced densty matrx approach The quantum beats nvestgated n ths paper are the result of nterference between two exctaton paths of

More information

Chapter 7 Channel Capacity and Coding

Chapter 7 Channel Capacity and Coding Wreless Informaton Transmsson System Lab. Chapter 7 Channel Capacty and Codng Insttute of Communcatons Engneerng atonal Sun Yat-sen Unversty Contents 7. Channel models and channel capacty 7.. Channel models

More information

= = = (a) Use the MATLAB command rref to solve the system. (b) Let A be the coefficient matrix and B be the right-hand side of the system.

= = = (a) Use the MATLAB command rref to solve the system. (b) Let A be the coefficient matrix and B be the right-hand side of the system. Chapter Matlab Exercses Chapter Matlab Exercses. Consder the lnear system of Example n Secton.. x x x y z y y z (a) Use the MATLAB command rref to solve the system. (b) Let A be the coeffcent matrx and

More information

2.3 Nilpotent endomorphisms

2.3 Nilpotent endomorphisms s a block dagonal matrx, wth A Mat dm U (C) In fact, we can assume that B = B 1 B k, wth B an ordered bass of U, and that A = [f U ] B, where f U : U U s the restrcton of f to U 40 23 Nlpotent endomorphsms

More information

Color Rendering Uncertainty

Color Rendering Uncertainty Australan Journal of Basc and Appled Scences 4(10): 4601-4608 010 ISSN 1991-8178 Color Renderng Uncertanty 1 A.el Bally M.M. El-Ganany 3 A. Al-amel 1 Physcs Department Photometry department- NIS Abstract:

More information

Fuzzy Boundaries of Sample Selection Model

Fuzzy Boundaries of Sample Selection Model Proceedngs of the 9th WSES Internatonal Conference on ppled Mathematcs, Istanbul, Turkey, May 7-9, 006 (pp309-34) Fuzzy Boundares of Sample Selecton Model L. MUHMD SFIIH, NTON BDULBSH KMIL, M. T. BU OSMN

More information

The Prncpal Component Transform The Prncpal Component Transform s also called Karhunen-Loeve Transform (KLT, Hotellng Transform, oregenvector Transfor

The Prncpal Component Transform The Prncpal Component Transform s also called Karhunen-Loeve Transform (KLT, Hotellng Transform, oregenvector Transfor Prncpal Component Transform Multvarate Random Sgnals A real tme sgnal x(t can be consdered as a random process and ts samples x m (m =0; ;N, 1 a random vector: The mean vector of X s X =[x0; ;x N,1] T

More information

An (almost) unbiased estimator for the S-Gini index

An (almost) unbiased estimator for the S-Gini index An (almost unbased estmator for the S-Gn ndex Thomas Demuynck February 25, 2009 Abstract Ths note provdes an unbased estmator for the absolute S-Gn and an almost unbased estmator for the relatve S-Gn for

More information

Efficient, General Point Cloud Registration with Kernel Feature Maps

Efficient, General Point Cloud Registration with Kernel Feature Maps Effcent, General Pont Cloud Regstraton wth Kernel Feature Maps Hanchen Xong, Sandor Szedmak, Justus Pater Insttute of Computer Scence Unversty of Innsbruck 30 May 2013 Hanchen Xong (Un.Innsbruck) 3D Regstraton

More information

Hidden Markov Models & The Multivariate Gaussian (10/26/04)

Hidden Markov Models & The Multivariate Gaussian (10/26/04) CS281A/Stat241A: Statstcal Learnng Theory Hdden Markov Models & The Multvarate Gaussan (10/26/04) Lecturer: Mchael I. Jordan Scrbes: Jonathan W. Hu 1 Hdden Markov Models As a bref revew, hdden Markov models

More information