Statististics of the SK estimator
|
|
- Stewart Neil Harrison
- 5 years ago
- Views:
Transcription
1 New Jersey Institute of Technology E-mil: Dle E. Gry New Jersey Institute of Technology E-mil: Spectrl Kurtosis (SK is sttisticl pproch for detecting nd removing rdio frequency interference (RFI in rdio stronomy dt. In this study, the sttisticl properties of the SK estimtor re investigted nd ll moments of its probbility density function re nlyticlly determined. These moments provide mens to determine the til probbilities of the estimtor tht re essentil to defining the thresholds for RFI discrimintion. It is shown tht, for number of ccumulted spectr M 24, the first SK stndrd moments stisfy the conditions required by Person Type IV probbility density function (PDF, which is shown to ccurtely reproduce the observed distributions. The cumultive function (CF of the Person Type IV is then found, in both nlyticl nd numericl forms, suitble for ccurte estimtion of the til probbilities of the SK estimtor. This sme frmework is lso shown to be pplicble to the relted Time Domin Kurtosis (TDK estimtor, whose PDF corresponds to Person Type IV when the number of time-domin smples is M 46. The PDF nd CF re determined for this cse lso. RFI mitigtion workshop - RFI2010, Mrch 29-31, 2010 Groningen, the Netherlnds Speker. This work ws supported by NSF grnt AST to NJIT c Copyright owned by the uthor(s under the terms of the Cretive Commons Attribution-NonCommercil-ShreAlike Licence.
2 1. Introduction The Spectrl Kurtosis estimtor (ŜK ws originlly proposed by Nit et l. [7] s sttisticl tool for rel-time rdio frequency interference (RFI detection nd excision in Fst Fourier Trnsform (FFT rdio spectrogrph. The first spectrogrph designed for ŜK, the Koren Solr Rdio Burst Loctor [KSRBL; 1], demonstrted the effectiveness of the SK lgorithm, but lso reveled the need for more ccurte clcultion of the theoreticl RFI detection thresholds thn initilly proposed. Consequently, Nit & Gry [8] derived the exct nlyticl expressions for the sttisticl moments of ŜK nd, bsed on its first four stndrd moments, ssigned to it Person Type IV probbility curve [9], which ws shown to be in very good greement with the Monte Crlo simulted ŜK probbility distribution function (PDF, s well s with the distribution derived from direct experimentl observtions mde with the KSRBL instrument [2]. As extensively described in the previous ppers, wht mkes n SK spectrogrph with N spectrl chnnels distinct from trditionl one is the fct tht it ccumultes not only set of M instntneous power spectrl density (PSD estimtes, denoted S 1, but lso the squred spectrl power denoted S 2. These sums, which hve n implicit dependence on frequency chnnel f k, re used to compute the verged power spectrum P = S 1 /M, s well s the quntity ŜK = M + 1 ( MS2 M 1 S1 2 1, (1.1 which is cumulnt-bsed estimtor of the spectrl vribility corresponding to the signl prent popultion, Vk 2 = σ k 2 µ k 2, (1.2 where µ k nd σk 2 re the frequency-dependent PSD popultion mens nd vrinces, respectively. For normlly distributed time domin signl, i.e. n RFI-free signl, Nit & Gry [8] showed tht the estimtor given by eqution (1.1 is unbised, i.e. E(ŜK = Vk 2 = 1. This presenttion compiles the min results of Nit & Gry [8] leding to the definition of ŜK given by eqution (1.1, nd provides prcticl guide for computing the RFI thresholds with predefined flse lrm probbility (PFA level. 2. The Unbised Spectrl Kurtosis nd Time Domin Kurtosis Estimtors The unbised estimtor ŜK my be derived bsed on the properties of the generlized gmm distribution [GGD, 10] defined s f (x,,d, p = pxd 1 e ( x p d Γ(d/p, (2.1 where Γ(z = 0 tz 1 e t dt is the well known Euler s Gmm function. The expected smple moments bout origin corresponding to GGD function re given by E(x n = Γ( d+n p Γ(d/p n, (2.2 2
3 which reduces to the known result E(x n = n!µ n in the prticulr cse of n exponentil distribution, f (x, µ,1,1. As shown by Nit & Gry [8], if one considers set of M independent rndom vribles tht re re individully distributed ccording with p = 1 GGD function, f (x,, d, 1, (which is nothing else thn stndrd gmm distribution, the probbility distribution of the men x = Σ M i=1 x i is the GGD function f ( x,/m,md,1, (2.3 nd the probbility distribution of the squred men is the GGD function [ ( 2, p( x 2 = f x 2 Md, M 2, 1 ], (2.4 2 which, mking use of eqution (2.2, provides the expecttions of n rbitrry power of the squred men, E( x 2n = Γ(Md + 2n ( 2n. (2.5 Γ(Md M Although closed form for the PDF of the men of squres x 2 = Σ M i=1 x2 i hs not been found, Nit & Gry [8] proved tht the sttisticl moments of the men of squres cn be exctly computed ccording with the formul E( x 2 n = (/ M 2n [Γ(d] M [ n 1 t n r=0 r! Γ(2r + dtr] M t=0. (2.6 n For the prticulr cses of the exponentil nd χ 2 distributions, f (x, µ,1,1 nd f (x, µ,1/2,1 respectively, Nit & Gry [8] proved tht x 2 / x 2 nd x 2 re uncorrelted rndom vribles, i.e. tht they hve null covrince, which immeditely led to the exct moments of the rtio between the men of squres nd squre of men, which cn be expressed s [( x 2 n ] E x 2 = E( x2 n E( x 2n. (2.7 However, s we will show here, the null covrince of x 2 / x 2 nd x 2 is n intrinsic property of ny p = 1 GGD function f (x,,d,1, which mkes eqution (2.7 hold for ny vlue of d. To prove this key property, we mke use of the expression ( x 2 cov x 2, x2 = 2 [ E(x 3 M E(x 2 ] 2E(x2 E(x 2 + E(x2 tht is vlid for ny prticulr PDF [Eq.19, 8], in which we enter the explicit moments provided by eqution (2.2 to obtin (2.8 ( x 2 cov x 2, x2 = 0. (2.9 Therefore, provided tht the observble x is distributed ccording to GGD function f (x,,d,1, nd tking in considertion equtions (2.5 nd (2.6, the sttisticl moments of the rtio between 3
4 the men of squres nd squre of men given by eqution (2.7 my be written in the compct form [( x 2 n ] M n Γ(Md n [ n 1 E x 2 = Γ(d M Γ(Md + 2n t n r=0 r! Γ(2r + dtr] M t=0, (2.10 which is independent of the scling prmeter. Nit & Gry [8] used the prticulr form of eqution (2.10 corresponding to d = 1 to derive the SK estimtor given by eqution (1.1, nd the prticulr form corresponding to d = 1/2 to derive time domin kurtosis (TDK estimtor T DK = M + 2 ( MS2 1, (2.11 M 1 both of them being unbised estimtors of the spectrl vribility corresponding to the underlying probbility distribution f (x,,d,1, which ccording to eqution (1.2 is S 2 1 V 2 = E(x2 E(x 2 E(x 2 = 1 d, (2.12 n expression tht is 1 for the d = 1 (SK cse, nd 2 for the d = 1/2 (TDK cse. 3. The Person Type IV Approximtion of the Spectrl nd Time Domin Kurtosis Estimtors Using eqution (2.10, the first stndrd moments of the ŜK estimtor given by eqution (1.1 my be written s µ 1 = 1; µ 2 = β 1 = 4M 2 (M 1(M + 2(M + 3 4(M + 2(M + 3(5M 72 (M 1(M (M ; β 2 = 3(M + 2(M + 3(M3 + 98M 2 185M + 78, (M 1(M + 4(M + 5(M + 6(M + 7 where β 1 = µ 3 2/µ3 2 nd β 2 = µ 4 /µ 2 2 re directly relted to the more commonly used skewness, γ 1 = β 1 nd kurtosis excess, γ 2 = β 2 3. Similrly, the first stndrd moments of the T DK estimtor re µ 1 = 2; µ 2 = 24M 2 (M 1(M + 4(M + 6 β 1 = 216(M 22 (M + 4(M + 6 (M 1(M (M ; β 2 = 3(M + 4(M + 6(M M 2 474M (M 1(M + 8(M + 10(M + 12(M The pproprite PDF pproximtion for the ŜK nd T DK estimtors my be found by investigting the Person s criterion defined s [5, p. 151]: κ = (3.1 (3.2 β 1 (β (4β 2 3β 1 (2β 2 3β 1 6, (3.3 where the exct vlues of the prmeters β 1 nd β 2 re provided by equtions (3.1 nd (3.2, respectively. Person s criterion indictes tht for M 24 the ŜK estimtor my be pproximted by Person Type IV probbility curve. The sme PDF pproximtion my be used for the T DK estimtor for M 46. 4
5 Figure 1: Comprison between the SK distributions (blck lines obtined by numericl simultion for different ccumultion lengths M nd their corresponding Person Type IV pproximtions (red lines. The four Person Type IV prmeters m, ν, λ, nd re displyed on ech plot. 3.1 Person Type IV Probbility Distribution Function The most generl nlyticl form of the Person Type IV PDF originlly introduced by Person [9], including its non-trivil normliztion fctor, ws given by Nghr [6] s p(x = 1 Γ(m + i ν 2 Γ(m i ν 2 [ ( x λ π Γ(m 1 2 Γ(m ] mexp [ νarctn( x λ ], (3.4 where the four prmeters m, µ,, nd λ my be expressed in terms of the centrl moments of the distribution s [4]: r = 6(β 2 β 1 1 2β 2 3β 1 6 ; m = r ; ν = r(r 2 β1 (3.5 16(r 1 β1 (r 2 2 = 1 µ 2 (6(r 1 β 1 (r ; λ = µ 1 4 (r 2 µ 2 β 1. Figure 1 compres the Person IV pproximtions of the ŜK PDF with the Monte Crlo simulted distributions for M = 32, 1024, 4096, nd By visul inspection, we my conclude tht the Person IV pproximtions ccurtely reproduce the shpes of the numericlly simulted histogrms for different orders of mgnitude of the ccumultion length. To compute the til probbilities of the Person Type IV PDF, one hs to compute the cumultive function P(x, (CF, nd the complementry cumultive function, 1 P(x, (CCF, P(x = x p(xdx; 1 P(x = for which Heinrich [4] provided the following closed form: 1 + P 1 (m,ν,,λ,x, x < λ 3 P(x = P 2 (m,ν,,λ,x x λ < 3 1 P 1 (m, ν,, λ, x, x > λ + 3, x p(xdx, (3.6 (3.7 where ( P 1 (m,ν,,λ,x = i x λ ( F 1,m + i ν 2m 1 2,2m, 2 1 i x λ p(x 5
6 Figure 2: ŜK threshold computtion for M = The lower nd upper thresholds displyed by the two verticl lines hve been estimted s the intersection points of the horizontl nd numericl integrtion lines. Their vlues of = / 6104 nd = / 6104, respectively, hve to be compred with the symmetric thresholds of 1 ± 6/ 6104 (verticl dotted lines. Figure 3: T DK threshold computtion for M = The two solid verticl lines, hving the ordintes nd , represent the RFI detection thresholds corresponding to symmetric stndrd flse lrm probbility level of %. These thresholds hve to be compred with with the less ccurte symmetric thresholds of 2 ± 3 24/12208 = 2 ± (verticl dotted lines. nd P 2 (m,ν,,λ,x = 1 1 e (ν+i2mπ i [ 1 + iν 2m + 2 F (1,2 2m,2 m + i ν 2, 1 + i x λ 2 F(α,β,δ,z = 1 + αβ α(α + 1β(β + 1 z + 1!δ 2!δ(δ + 1 z = ( x λ 2 ] p(x, k=0 α (k β (k k!δ (k z k (3.8 is the Guss hypergeometric series. Alterntively, insted of eqution (3.7, one my use the formul provided by Willink [11] where P(m,ν,,λ,x = e [λ i(2 2m]Φ R 1 e [λ i(2 2m]π 1, (3.9 Φ = π ( x λ 2 + rctn ; u = m i 2 ν; R = F(2 2m,u,u+1,eiΦ F(2 2m,u,u+1,1. However, if one wnts to void the numericl difficulties relted to the evlution of the hypergeometric series, one my choose to perform direct numericl integrtion (eqution [3.6] of eqution (3.4, which my chieve resonble ccurcy with fr less computtionl effort, especilly if tilored integrtion methods, [e.g. 6], re employed. Figure 2 displys the numericl results for M = 6104, computed ccording to eqution (3.7, (tringulr symbols, eqution (3.9, (squre symbols, nd by direct integrtion of equtions (3.6, 6
7 (solid line. The hypergeometric series where computed using the hypergeom function in Mple 11 (MpleSoft, nd the numericl integrtion ws performed using the int_tbulted function in IDL 6.4 (ITT. The plots disply both CF (rising nd CCF (descending needed to evlute the RFI thresholds equivlent to norml distribution s ±3σ level (probbility %, horizontl solid line. It my be concluded tht, in the region of interest, ll three methods provide similr numericl results. However, it ws found tht, for SK vlues well before the distribution pek, the numericl ccurcy of eqution (3.9 is better thn tht of eqution (3.7, while the direct numericl integrtion of CF gives similr results s eqution (3.9. After the pek of the ŜK distribution, the numericl ccurcy of eqution (3.7 is better thn tht of eqution (3.9, while the numericl integrtion of CCF gives similr results s eqution (3.7. Therefore we conclude tht the numericl evlution of eqution (3.9 gives more ccurte estimtion of the CF nd the numericl evlution of eqution (3.7 gives more ccurte estimtion of the CCF, while the direct numericl integrtion of equtions (3.6 gives results of comprble ccurcy t both sides of the ŜK distribution. The lower nd higher thresholds displyed by the two verticl lines hve been estimted s the intersection points of the horizontl nd numericl integrtion lines. Their vlues of = / 6104 nd = / 6104, respectively, re compred with the symmetric thresholds of 1 ± 6/ 6104 (verticl dotted lines originlly proposed by Nit & Gry [8]. Although this correction seems smll in bsolute vlue for the lrge-m cse, e.g. M = 6104 illustrted in Figure 2, we clculte tht, compred with the symmetric thresholds, the new thresholds ccount for 67% less rejection of vlid dt s flse RFI occurrences t the upper bound of the distribution, nd provide better rejection of true RFI signls of low signl to noise rtio t the lower bound. In combintion, the result is n overll better performnce of the ŜK bsed RFI rejection lgorithm. The correction becomes more importnt for lower M. Figure 3 displys the PDF of the T DK estimtor corresponding to n ccumultion length of M = 12208, chosen to mtch the sme frequency nd time resolution of DFT-bsed spectrogrph with M = 6104 (the exmple used in Fig. 2; see Nit et l. [7] for more detiled motivtion of this choice. Despite its lrge ccumultion length, the estimtor K still hs noticeble skewness, which needs to be properly considered in order to obtin the flse lrm probbility levels equivlent to ±3σ for norml distribution. Compred with the symmetric thresholds of 2±3 24/12208, the new thresholds would reject 71% less vlid dt t the higher end of the distribution, while the shifted lower threshold would improve the sensitivity of RFI detection t the lower end of the distribution. 4. Conclusion We hve investigted the sttisticl properties of the SK estimtor nd determined nlyticl expressions for its PDF nd CF with the gol of improving the selection of thresholds for RFI discrimintion. We lso improved the definition of ŜK (eqution [1.1] reltive to its originl definition [7] to form n unbised estimtor, nd introduced TDK unbised estimtor (eqution [2.11] to be used for RFI detection t the DC nd Nyquist frequency bins of DFT-bsed spectrogrph, or t ny frequency bin of FIR bsed spectrogrph. We hve derived closed form nlyticl expressions for the complete set of the centrl moments of the SK nd TDK estimtors, nd estblished common frmework tht llows ccurte estimtion of the RFI thresholds bsed on the first four stndrd moments of their probbility distributions (equtions [3.1] nd [3.2], which, for ny 7
8 ccumultion length M 24 nd M 46, respectively, re used to compute the four prmeters (eqution [3.5] tht completely determine the Person IV pproximtions (eqution [3.4] of their true PDFs. Bsed on these four prmeters, which depend only on the ccumultion length M, the CF nd CCF of the SK or TDK estimtors my be computed by using either the closed forms expressions provided by equtions (3.9 nd (3.7, respectively, or by direct numericl estimtion of the integrls given by eqution (3.6. Compred to the symmetricl thresholds originlly suggested in Pper I, the procedure described in this study properly tkes into ccount the intrinsic skewness of the probbility density functions of the SK nd TDK estimtors, which provides better overll RFI detection performnce for either smll or lrge ccumultion lengths. These theoreticlly estblished results re shown in Gry, Liu & Nit [2, 3] to be exctly obeyed by dt tken in the KSRBL spectrometer hrdwre implementtion of the lgorithm, where the improvement in RFI excision by use of these modified thresholds is confirmed. The modified thresholds become ever more importnt when smller number M of ccumultions is used. A simple procedure hs been written in IDL (Interctive Dt Lnguge for numericl clcultion of the thresholds for ny M, nd is vilble upon request from the uthors. References [1] Y. Dou, D. E. Gry, Z. Liu, G. M. Nit, S. -C. Bong, K. -S. Cho, Y. -D. Prk, & Y. -J. Moon, The Koren Solr Rdio Burst Loctor (KSRBL, PASP, 121, 512 [2] D. E. Gry, Z. Liu, & G. M. Nit, A Widebnd Spectrometer with RFI Detection, PASP, 122, 560 [3] D. E. Gry, Z. Liu, & G. M. Nit, Hrdwre Implementtion of n SK Spectrometer, in this proceedings, PoS(RFI [4] J. Heinrich, A Guide to the Person Type IV Distribution, Collider Detector t Fermilb internl note 6820, [5] M. G. Kendll & A. Sturt, The Advnced Theory of Sttistics Vol.1, Griffin, London 1958 [6] Y. Nghr, The PDF nd CF of Person type IV distributions nd the ML estimtion of the prmeters, Sttistics & Probbility Letters, 43, 251 [7] G. M. Nit, D. E. Gry, Z. Liu, G. J. Hurford, & S. M. White, Rdio Frequency Interference Excision Using Spectrl Domin Sttistics, PASP, 119, 805 [8] G. M. Nit & D. E. Gry, Sttistics of the Spectrl Kurtosis Estimtor, PASP, 122, 595 [9] K. Person, Contributions to the Mthemticl Theory of Evolution.-II. Skew Vrition in Homogeneous Mteril, Philos. Trns. R. Soc. London A, 186, 343 [10] E.W. Stcy, A Generliztion of the Gmm Distribution, Ann. Mth. Sttist. Assoc., 33, 1187 [11] R. Willink, A Closed-Form Expression For The Person Type IV Distribution Function, Austrl. NZ J. Stt., 50(2, 199 8
Tests for the Ratio of Two Poisson Rates
Chpter 437 Tests for the Rtio of Two Poisson Rtes Introduction The Poisson probbility lw gives the probbility distribution of the number of events occurring in specified intervl of time or spce. The Poisson
More informationCS667 Lecture 6: Monte Carlo Integration 02/10/05
CS667 Lecture 6: Monte Crlo Integrtion 02/10/05 Venkt Krishnrj Lecturer: Steve Mrschner 1 Ide The min ide of Monte Crlo Integrtion is tht we cn estimte the vlue of n integrl by looking t lrge number of
More informationContinuous Random Variables
STAT/MATH 395 A - PROBABILITY II UW Winter Qurter 217 Néhémy Lim Continuous Rndom Vribles Nottion. The indictor function of set S is rel-vlued function defined by : { 1 if x S 1 S (x) if x S Suppose tht
More informationChapter 9: Inferences based on Two samples: Confidence intervals and tests of hypotheses
Chpter 9: Inferences bsed on Two smples: Confidence intervls nd tests of hypotheses 9.1 The trget prmeter : difference between two popultion mens : difference between two popultion proportions : rtio of
More informationChapter 5 : Continuous Random Variables
STAT/MATH 395 A - PROBABILITY II UW Winter Qurter 216 Néhémy Lim Chpter 5 : Continuous Rndom Vribles Nottions. N {, 1, 2,...}, set of nturl numbers (i.e. ll nonnegtive integers); N {1, 2,...}, set of ll
More information1 Probability Density Functions
Lis Yn CS 9 Continuous Distributions Lecture Notes #9 July 6, 28 Bsed on chpter by Chris Piech So fr, ll rndom vribles we hve seen hve been discrete. In ll the cses we hve seen in CS 9, this ment tht our
More informationMonte Carlo method in solving numerical integration and differential equation
Monte Crlo method in solving numericl integrtion nd differentil eqution Ye Jin Chemistry Deprtment Duke University yj66@duke.edu Abstrct: Monte Crlo method is commonly used in rel physics problem. The
More informationTHE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.
THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS RADON ROSBOROUGH https://intuitiveexplntionscom/picrd-lindelof-theorem/ This document is proof of the existence-uniqueness theorem
More informationLecture 3 Gaussian Probability Distribution
Introduction Lecture 3 Gussin Probbility Distribution Gussin probbility distribution is perhps the most used distribution in ll of science. lso clled bell shped curve or norml distribution Unlike the binomil
More informationThe steps of the hypothesis test
ttisticl Methods I (EXT 7005) Pge 78 Mosquito species Time of dy A B C Mid morning 0.0088 5.4900 5.5000 Mid Afternoon.3400 0.0300 0.8700 Dusk 0.600 5.400 3.000 The Chi squre test sttistic is the sum of
More information3.4 Numerical integration
3.4. Numericl integrtion 63 3.4 Numericl integrtion In mny economic pplictions it is necessry to compute the definite integrl of relvlued function f with respect to "weight" function w over n intervl [,
More informationMath& 152 Section Integration by Parts
Mth& 5 Section 7. - Integrtion by Prts Integrtion by prts is rule tht trnsforms the integrl of the product of two functions into other (idelly simpler) integrls. Recll from Clculus I tht given two differentible
More informationNon-Linear & Logistic Regression
Non-Liner & Logistic Regression If the sttistics re boring, then you've got the wrong numbers. Edwrd R. Tufte (Sttistics Professor, Yle University) Regression Anlyses When do we use these? PART 1: find
More informationNUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.
NUMERICAL INTEGRATION 1 Introduction The inverse process to differentition in clculus is integrtion. Mthemticlly, integrtion is represented by f(x) dx which stnds for the integrl of the function f(x) with
More informationGoals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite
Unit #8 : The Integrl Gols: Determine how to clculte the re described by function. Define the definite integrl. Eplore the reltionship between the definite integrl nd re. Eplore wys to estimte the definite
More informationConstruction of Gauss Quadrature Rules
Jim Lmbers MAT 772 Fll Semester 2010-11 Lecture 15 Notes These notes correspond to Sections 10.2 nd 10.3 in the text. Construction of Guss Qudrture Rules Previously, we lerned tht Newton-Cotes qudrture
More informationSection 11.5 Estimation of difference of two proportions
ection.5 Estimtion of difference of two proportions As seen in estimtion of difference of two mens for nonnorml popultion bsed on lrge smple sizes, one cn use CLT in the pproximtion of the distribution
More informationReview of Calculus, cont d
Jim Lmbers MAT 460 Fll Semester 2009-10 Lecture 3 Notes These notes correspond to Section 1.1 in the text. Review of Clculus, cont d Riemnn Sums nd the Definite Integrl There re mny cses in which some
More informationMath 1B, lecture 4: Error bounds for numerical methods
Mth B, lecture 4: Error bounds for numericl methods Nthn Pflueger 4 September 0 Introduction The five numericl methods descried in the previous lecture ll operte by the sme principle: they pproximte the
More information7.2 The Definite Integral
7.2 The Definite Integrl the definite integrl In the previous section, it ws found tht if function f is continuous nd nonnegtive, then the re under the grph of f on [, b] is given by F (b) F (), where
More informationLecture 21: Order statistics
Lecture : Order sttistics Suppose we hve N mesurements of sclr, x i =, N Tke ll mesurements nd sort them into scending order x x x 3 x N Define the mesured running integrl S N (x) = 0 for x < x = i/n for
More informationStudent Activity 3: Single Factor ANOVA
MATH 40 Student Activity 3: Single Fctor ANOVA Some Bsic Concepts In designed experiment, two or more tretments, or combintions of tretments, is pplied to experimentl units The number of tretments, whether
More informationTHE EXISTENCE OF NEGATIVE MCMENTS OF CONTINUOUS DISTRIBUTIONS WALTER W. PIEGORSCH AND GEORGE CASELLA. Biometrics Unit, Cornell University, Ithaca, NY
THE EXISTENCE OF NEGATIVE MCMENTS OF CONTINUOUS DISTRIBUTIONS WALTER W. PIEGORSCH AND GEORGE CASELLA. Biometrics Unit, Cornell University, Ithc, NY 14853 BU-771-M * December 1982 Abstrct The question of
More informationNumerical integration
2 Numericl integrtion This is pge i Printer: Opque this 2. Introduction Numericl integrtion is problem tht is prt of mny problems in the economics nd econometrics literture. The orgniztion of this chpter
More information63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1
3 9. SEQUENCES AND SERIES 63. Representtion of functions s power series Consider power series x 2 + x 4 x 6 + x 8 + = ( ) n x 2n It is geometric series with q = x 2 nd therefore it converges for ll q =
More informationPi evaluation. Monte Carlo integration
Pi evlution y 1 1 x Computtionl Physics 2018-19 (Phys Dep IST, Lisbon) Fernndo Bro (311) Monte Crlo integrtion we wnt to evlute the following integrl: F = f (x) dx remember tht the expecttion vlue of the
More informationarxiv: v1 [astro-ph.im] 24 May 2010
Mon. Not. R. Astron. Soc., 5 (22 Printed November 28 (MN LATEX style file v2.2 The Generalized Spectral Kurtosis Estimator arxiv:5.437v [astro-ph.im] 24 May 2 G. M. Nita and D. E. Gary Center for Solar-Terrestrial
More informationProperties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives
Block #6: Properties of Integrls, Indefinite Integrls Gols: Definition of the Definite Integrl Integrl Clcultions using Antiderivtives Properties of Integrls The Indefinite Integrl 1 Riemnn Sums - 1 Riemnn
More informationEstimation of Binomial Distribution in the Light of Future Data
British Journl of Mthemtics & Computer Science 102: 1-7, 2015, Article no.bjmcs.19191 ISSN: 2231-0851 SCIENCEDOMAIN interntionl www.sciencedomin.org Estimtion of Binomil Distribution in the Light of Future
More informationMIXED MODELS (Sections ) I) In the unrestricted model, interactions are treated as in the random effects model:
1 2 MIXED MODELS (Sections 17.7 17.8) Exmple: Suppose tht in the fiber breking strength exmple, the four mchines used were the only ones of interest, but the interest ws over wide rnge of opertors, nd
More informationTangent Line and Tangent Plane Approximations of Definite Integral
Rose-Hulmn Undergrdute Mthemtics Journl Volume 16 Issue 2 Article 8 Tngent Line nd Tngent Plne Approximtions of Definite Integrl Meghn Peer Sginw Vlley Stte University Follow this nd dditionl works t:
More informationChapter 4 Contravariance, Covariance, and Spacetime Diagrams
Chpter 4 Contrvrince, Covrince, nd Spcetime Digrms 4. The Components of Vector in Skewed Coordintes We hve seen in Chpter 3; figure 3.9, tht in order to show inertil motion tht is consistent with the Lorentz
More informationAn approximation to the arithmetic-geometric mean. G.J.O. Jameson, Math. Gazette 98 (2014), 85 95
An pproximtion to the rithmetic-geometric men G.J.O. Jmeson, Mth. Gzette 98 (4), 85 95 Given positive numbers > b, consider the itertion given by =, b = b nd n+ = ( n + b n ), b n+ = ( n b n ) /. At ech
More informationP 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0)
1 Tylor polynomils In Section 3.5, we discussed how to pproximte function f(x) round point in terms of its first derivtive f (x) evluted t, tht is using the liner pproximtion f() + f ()(x ). We clled this
More informationAcceptance Sampling by Attributes
Introduction Acceptnce Smpling by Attributes Acceptnce smpling is concerned with inspection nd decision mking regrding products. Three spects of smpling re importnt: o Involves rndom smpling of n entire
More informationLecture 19: Continuous Least Squares Approximation
Lecture 19: Continuous Lest Squres Approximtion 33 Continuous lest squres pproximtion We begn 31 with the problem of pproximting some f C[, b] with polynomil p P n t the discrete points x, x 1,, x m for
More informationNew data structures to reduce data size and search time
New dt structures to reduce dt size nd serch time Tsuneo Kuwbr Deprtment of Informtion Sciences, Fculty of Science, Kngw University, Hirtsuk-shi, Jpn FIT2018 1D-1, No2, pp1-4 Copyright (c)2018 by The Institute
More informationUniversity of Texas MD Anderson Cancer Center Department of Biostatistics. Inequality Calculator, Version 3.0 November 25, 2013 User s Guide
University of Texs MD Anderson Cncer Center Deprtment of Biosttistics Inequlity Clcultor, Version 3.0 November 5, 013 User s Guide 0. Overview The purpose of the softwre is to clculte the probbility tht
More informationLecture 6: Singular Integrals, Open Quadrature rules, and Gauss Quadrature
Lecture notes on Vritionl nd Approximte Methods in Applied Mthemtics - A Peirce UBC Lecture 6: Singulr Integrls, Open Qudrture rules, nd Guss Qudrture (Compiled 6 August 7) In this lecture we discuss the
More informationA New Statistic Feature of the Short-Time Amplitude Spectrum Values for Human s Unvoiced Pronunciation
Xiodong Zhung A ew Sttistic Feture of the Short-Time Amplitude Spectrum Vlues for Humn s Unvoiced Pronuncition IAODOG ZHUAG 1 1. Qingdo University, Electronics & Informtion College, Qingdo, 6671 CHIA Abstrct:
More informationTHERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION
XX IMEKO World Congress Metrology for Green Growth September 9,, Busn, Republic of Kore THERMAL EXPANSION COEFFICIENT OF WATER FOR OLUMETRIC CALIBRATION Nieves Medin Hed of Mss Division, CEM, Spin, mnmedin@mityc.es
More informationJournal of Inequalities in Pure and Applied Mathematics
Journl of Inequlities in Pure nd Applied Mthemtics MOMENTS INEQUALITIES OF A RANDOM VARIABLE DEFINED OVER A FINITE INTERVAL PRANESH KUMAR Deprtment of Mthemtics & Computer Science University of Northern
More information8 Laplace s Method and Local Limit Theorems
8 Lplce s Method nd Locl Limit Theorems 8. Fourier Anlysis in Higher DImensions Most of the theorems of Fourier nlysis tht we hve proved hve nturl generliztions to higher dimensions, nd these cn be proved
More informationSection 4.8. D v(t j 1 ) t. (4.8.1) j=1
Difference Equtions to Differentil Equtions Section.8 Distnce, Position, nd the Length of Curves Although we motivted the definition of the definite integrl with the notion of re, there re mny pplictions
More informationNumerical Analysis: Trapezoidal and Simpson s Rule
nd Simpson s Mthemticl question we re interested in numericlly nswering How to we evlute I = f (x) dx? Clculus tells us tht if F(x) is the ntiderivtive of function f (x) on the intervl [, b], then I =
More informationPurpose of the experiment
Newton s Lws II PES 6 Advnced Physics Lb I Purpose of the experiment Exmine two cses using Newton s Lws. Sttic ( = 0) Dynmic ( 0) fyi fyi Did you know tht the longest recorded flight of chicken is thirteen
More informationThe Shortest Confidence Interval for the Mean of a Normal Distribution
Interntionl Journl of Sttistics nd Proility; Vol. 7, No. 2; Mrch 208 ISSN 927-7032 E-ISSN 927-7040 Pulished y Cndin Center of Science nd Eduction The Shortest Confidence Intervl for the Men of Norml Distriution
More informationUndergraduate Research
Undergrdute Reserch A Trigonometric Simpson s Rule By Ctherine Cusimno Kirby nd Sony Stnley Biogrphicl Sketch Ctherine Cusimno Kirby is the dughter of Donn nd Sm Cusimno. Originlly from Vestvi Hills, Albm,
More informationS. S. Dragomir. 2, we have the inequality. b a
Bull Koren Mth Soc 005 No pp 3 30 SOME COMPANIONS OF OSTROWSKI S INEQUALITY FOR ABSOLUTELY CONTINUOUS FUNCTIONS AND APPLICATIONS S S Drgomir Abstrct Compnions of Ostrowski s integrl ineulity for bsolutely
More informationExpectation and Variance
Expecttion nd Vrince : sum of two die rolls P(= P(= = 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12 P(=2) = 1/36 P(=3) = 1/18 P(=4) = 1/12 P(=5) = 1/9 P(=7) = 1/6 P(=13) =? 2 1/36 3 1/18 4 1/12 5 1/9 6 5/36 7 1/6
More informationNumerical Integration
Chpter 1 Numericl Integrtion Numericl differentition methods compute pproximtions to the derivtive of function from known vlues of the function. Numericl integrtion uses the sme informtion to compute numericl
More informationSession 13
780.20 Session 3 (lst revised: Februry 25, 202) 3 3. 780.20 Session 3. Follow-ups to Session 2 Histogrms of Uniform Rndom Number Distributions. Here is typicl figure you might get when histogrmming uniform
More informationDOING PHYSICS WITH MATLAB MATHEMATICAL ROUTINES
DOIG PHYSICS WITH MATLAB MATHEMATICAL ROUTIES COMPUTATIO OF OE-DIMESIOAL ITEGRALS In Cooper School of Physics, University of Sydney in.cooper@sydney.edu.u DOWLOAD DIRECTORY FOR MATLAB SCRIPTS mth_integrtion_1d.m
More information2D1431 Machine Learning Lab 3: Reinforcement Learning
2D1431 Mchine Lerning Lb 3: Reinforcement Lerning Frnk Hoffmnn modified by Örjn Ekeberg December 7, 2004 1 Introduction In this lb you will lern bout dynmic progrmming nd reinforcement lerning. It is ssumed
More informationPredict Global Earth Temperature using Linier Regression
Predict Globl Erth Temperture using Linier Regression Edwin Swndi Sijbt (23516012) Progrm Studi Mgister Informtik Sekolh Teknik Elektro dn Informtik ITB Jl. Gnesh 10 Bndung 40132, Indonesi 23516012@std.stei.itb.c.id
More information20 MATHEMATICS POLYNOMIALS
0 MATHEMATICS POLYNOMIALS.1 Introduction In Clss IX, you hve studied polynomils in one vrible nd their degrees. Recll tht if p(x) is polynomil in x, the highest power of x in p(x) is clled the degree of
More informationAbstract inner product spaces
WEEK 4 Abstrct inner product spces Definition An inner product spce is vector spce V over the rel field R equipped with rule for multiplying vectors, such tht the product of two vectors is sclr, nd the
More information221B Lecture Notes WKB Method
Clssicl Limit B Lecture Notes WKB Method Hmilton Jcobi Eqution We strt from the Schrödinger eqution for single prticle in potentil i h t ψ x, t = [ ] h m + V x ψ x, t. We cn rewrite this eqution by using
More informationRiemann Sums and Riemann Integrals
Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 2013 Outline 1 Riemnn Sums 2 Riemnn Integrls 3 Properties
More informationReview of Gaussian Quadrature method
Review of Gussin Qudrture method Nsser M. Asi Spring 006 compiled on Sundy Decemer 1, 017 t 09:1 PM 1 The prolem To find numericl vlue for the integrl of rel vlued function of rel vrile over specific rnge
More informationLecture 20: Numerical Integration III
cs4: introduction to numericl nlysis /8/0 Lecture 0: Numericl Integrtion III Instructor: Professor Amos Ron Scribes: Mrk Cowlishw, Yunpeng Li, Nthnel Fillmore For the lst few lectures we hve discussed
More information13: Diffusion in 2 Energy Groups
3: Diffusion in Energy Groups B. Rouben McMster University Course EP 4D3/6D3 Nucler Rector Anlysis (Rector Physics) 5 Sept.-Dec. 5 September Contents We study the diffusion eqution in two energy groups
More informationSection 6.1 INTRO to LAPLACE TRANSFORMS
Section 6. INTRO to LAPLACE TRANSFORMS Key terms: Improper Integrl; diverge, converge A A f(t)dt lim f(t)dt Piecewise Continuous Function; jump discontinuity Function of Exponentil Order Lplce Trnsform
More information5.7 Improper Integrals
458 pplictions of definite integrls 5.7 Improper Integrls In Section 5.4, we computed the work required to lift pylod of mss m from the surfce of moon of mss nd rdius R to height H bove the surfce of the
More informationRiemann Sums and Riemann Integrals
Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 203 Outline Riemnn Sums Riemnn Integrls Properties Abstrct
More informationA NOTE ON ESTIMATION OF THE GLOBAL INTENSITY OF A CYCLIC POISSON PROCESS IN THE PRESENCE OF LINEAR TREND
A NOTE ON ESTIMATION OF THE GLOBAL INTENSITY OF A CYCLIC POISSON PROCESS IN THE PRESENCE OF LINEAR TREND I WAYAN MANGKU Deprtment of Mthemtics, Fculty of Mthemtics nd Nturl Sciences, Bogor Agriculturl
More information1 Online Learning and Regret Minimization
2.997 Decision-Mking in Lrge-Scle Systems My 10 MIT, Spring 2004 Hndout #29 Lecture Note 24 1 Online Lerning nd Regret Minimiztion In this lecture, we consider the problem of sequentil decision mking in
More informationCredibility Hypothesis Testing of Fuzzy Triangular Distributions
666663 Journl of Uncertin Systems Vol.9, No., pp.6-74, 5 Online t: www.jus.org.uk Credibility Hypothesis Testing of Fuzzy Tringulr Distributions S. Smpth, B. Rmy Received April 3; Revised 4 April 4 Abstrct
More informationData Assimilation. Alan O Neill Data Assimilation Research Centre University of Reading
Dt Assimiltion Aln O Neill Dt Assimiltion Reserch Centre University of Reding Contents Motivtion Univrite sclr dt ssimiltion Multivrite vector dt ssimiltion Optiml Interpoltion BLUE 3d-Vritionl Method
More informationThe Regulated and Riemann Integrals
Chpter 1 The Regulted nd Riemnn Integrls 1.1 Introduction We will consider severl different pproches to defining the definite integrl f(x) dx of function f(x). These definitions will ll ssign the sme vlue
More informationNUMERICAL INTEGRATION
NUMERICAL INTEGRATION How do we evlute I = f (x) dx By the fundmentl theorem of clculus, if F (x) is n ntiderivtive of f (x), then I = f (x) dx = F (x) b = F (b) F () However, in prctice most integrls
More information31.2. Numerical Integration. Introduction. Prerequisites. Learning Outcomes
Numericl Integrtion 3. Introduction In this Section we will present some methods tht cn be used to pproximte integrls. Attention will be pid to how we ensure tht such pproximtions cn be gurnteed to be
More informationTesting categorized bivariate normality with two-stage. polychoric correlation estimates
Testing ctegorized bivrite normlity with two-stge polychoric correltion estimtes Albert Mydeu-Olivres Dept. of Psychology University of Brcelon Address correspondence to: Albert Mydeu-Olivres. Fculty of
More information1 The Lagrange interpolation formula
Notes on Qudrture 1 The Lgrnge interpoltion formul We briefly recll the Lgrnge interpoltion formul. The strting point is collection of N + 1 rel points (x 0, y 0 ), (x 1, y 1 ),..., (x N, y N ), with x
More informationCzechoslovak Mathematical Journal, 55 (130) (2005), , Abbotsford. 1. Introduction
Czechoslovk Mthemticl Journl, 55 (130) (2005), 933 940 ESTIMATES OF THE REMAINDER IN TAYLOR S THEOREM USING THE HENSTOCK-KURZWEIL INTEGRAL, Abbotsford (Received Jnury 22, 2003) Abstrct. When rel-vlued
More informationA REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007
A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H Thoms Shores Deprtment of Mthemtics University of Nebrsk Spring 2007 Contents Rtes of Chnge nd Derivtives 1 Dierentils 4 Are nd Integrls 5 Multivrite Clculus
More informationSUMMER KNOWHOW STUDY AND LEARNING CENTRE
SUMMER KNOWHOW STUDY AND LEARNING CENTRE Indices & Logrithms 2 Contents Indices.2 Frctionl Indices.4 Logrithms 6 Exponentil equtions. Simplifying Surds 13 Opertions on Surds..16 Scientific Nottion..18
More informationDiscrete Mathematics and Probability Theory Spring 2013 Anant Sahai Lecture 17
EECS 70 Discrete Mthemtics nd Proility Theory Spring 2013 Annt Shi Lecture 17 I.I.D. Rndom Vriles Estimting the is of coin Question: We wnt to estimte the proportion p of Democrts in the US popultion,
More informationCHM Physical Chemistry I Chapter 1 - Supplementary Material
CHM 3410 - Physicl Chemistry I Chpter 1 - Supplementry Mteril For review of some bsic concepts in mth, see Atkins "Mthemticl Bckground 1 (pp 59-6), nd "Mthemticl Bckground " (pp 109-111). 1. Derivtion
More information1B40 Practical Skills
B40 Prcticl Skills Comining uncertinties from severl quntities error propgtion We usully encounter situtions where the result of n experiment is given in terms of two (or more) quntities. We then need
More informationOrdinary differential equations
Ordinry differentil equtions Introduction to Synthetic Biology E Nvrro A Montgud P Fernndez de Cordob JF Urchueguí Overview Introduction-Modelling Bsic concepts to understnd n ODE. Description nd properties
More informationReinforcement Learning
Reinforcement Lerning Tom Mitchell, Mchine Lerning, chpter 13 Outline Introduction Comprison with inductive lerning Mrkov Decision Processes: the model Optiml policy: The tsk Q Lerning: Q function Algorithm
More informationPhysics 201 Lab 3: Measurement of Earth s local gravitational field I Data Acquisition and Preliminary Analysis Dr. Timothy C. Black Summer I, 2018
Physics 201 Lb 3: Mesurement of Erth s locl grvittionl field I Dt Acquisition nd Preliminry Anlysis Dr. Timothy C. Blck Summer I, 2018 Theoreticl Discussion Grvity is one of the four known fundmentl forces.
More informationMarkscheme May 2016 Mathematics Standard level Paper 1
M6/5/MATME/SP/ENG/TZ/XX/M Mrkscheme My 06 Mthemtics Stndrd level Pper 7 pges M6/5/MATME/SP/ENG/TZ/XX/M This mrkscheme is the property of the Interntionl Bcclurete nd must not be reproduced or distributed
More informationDiscrete Mathematics and Probability Theory Summer 2014 James Cook Note 17
CS 70 Discrete Mthemtics nd Proility Theory Summer 2014 Jmes Cook Note 17 I.I.D. Rndom Vriles Estimting the is of coin Question: We wnt to estimte the proportion p of Democrts in the US popultion, y tking
More informationNumerical Integration
Chpter 5 Numericl Integrtion Numericl integrtion is the study of how the numericl vlue of n integrl cn be found. Methods of function pproximtion discussed in Chpter??, i.e., function pproximtion vi the
More informationThe First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a).
The Fundmentl Theorems of Clculus Mth 4, Section 0, Spring 009 We now know enough bout definite integrls to give precise formultions of the Fundmentl Theorems of Clculus. We will lso look t some bsic emples
More informationf(x) dx, If one of these two conditions is not met, we call the integral improper. Our usual definition for the value for the definite integral
Improper Integrls Every time tht we hve evluted definite integrl such s f(x) dx, we hve mde two implicit ssumptions bout the integrl:. The intervl [, b] is finite, nd. f(x) is continuous on [, b]. If one
More informationMath 426: Probability Final Exam Practice
Mth 46: Probbility Finl Exm Prctice. Computtionl problems 4. Let T k (n) denote the number of prtitions of the set {,..., n} into k nonempty subsets, where k n. Argue tht T k (n) kt k (n ) + T k (n ) by
More informationMain topics for the First Midterm
Min topics for the First Midterm The Midterm will cover Section 1.8, Chpters 2-3, Sections 4.1-4.8, nd Sections 5.1-5.3 (essentilly ll of the mteril covered in clss). Be sure to know the results of the
More informationPhysics 116C Solution of inhomogeneous ordinary differential equations using Green s functions
Physics 6C Solution of inhomogeneous ordinry differentil equtions using Green s functions Peter Young November 5, 29 Homogeneous Equtions We hve studied, especilly in long HW problem, second order liner
More informationFor the percentage of full time students at RCC the symbols would be:
Mth 17/171 Chpter 7- ypothesis Testing with One Smple This chpter is s simple s the previous one, except it is more interesting In this chpter we will test clims concerning the sme prmeters tht we worked
More informationComparison Procedures
Comprison Procedures Single Fctor, Between-Subects Cse /8/ Comprison Procedures, One-Fctor ANOVA, Between Subects Two Comprison Strtegies post hoc (fter-the-fct) pproch You re interested in discovering
More informationDriving Cycle Construction of City Road for Hybrid Bus Based on Markov Process Deng Pan1, a, Fengchun Sun1,b*, Hongwen He1, c, Jiankun Peng1, d
Interntionl Industril Informtics nd Computer Engineering Conference (IIICEC 15) Driving Cycle Construction of City Rod for Hybrid Bus Bsed on Mrkov Process Deng Pn1,, Fengchun Sun1,b*, Hongwen He1, c,
More informationAPPROXIMATE INTEGRATION
APPROXIMATE INTEGRATION. Introduction We hve seen tht there re functions whose nti-derivtives cnnot be expressed in closed form. For these resons ny definite integrl involving these integrnds cnnot be
More information221A Lecture Notes WKB Method
A Lecture Notes WKB Method Hmilton Jcobi Eqution We strt from the Schrödinger eqution for single prticle in potentil i h t ψ x, t = [ ] h m + V x ψ x, t. We cn rewrite this eqution by using ψ x, t = e
More informationHybrid Group Acceptance Sampling Plan Based on Size Biased Lomax Model
Mthemtics nd Sttistics 2(3): 137-141, 2014 DOI: 10.13189/ms.2014.020305 http://www.hrpub.org Hybrid Group Acceptnce Smpling Pln Bsed on Size Bised Lomx Model R. Subb Ro 1,*, A. Ng Durgmmb 2, R.R.L. Kntm
More informationLecture 14: Quadrature
Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be rel-vlues nd smooth The pproximtion of n integrl by numericl
More informationMATH 144: Business Calculus Final Review
MATH 144: Business Clculus Finl Review 1 Skills 1. Clculte severl limits. 2. Find verticl nd horizontl symptotes for given rtionl function. 3. Clculte derivtive by definition. 4. Clculte severl derivtives
More informationConstruction and Selection of Single Sampling Quick Switching Variables System for given Control Limits Involving Minimum Sum of Risks
Construction nd Selection of Single Smpling Quick Switching Vribles System for given Control Limits Involving Minimum Sum of Risks Dr. D. SENHILKUMAR *1 R. GANESAN B. ESHA RAFFIE 1 Associte Professor,
More information