Asymptotic Normality in Monte Carlo Integration
|
|
- Clementine Warner
- 5 years ago
- Views:
Transcription
1 MATHEMATICS OF COMPUTATION, VOLUME 30, NUMBER 136 OCTOBER 1976, PAGES Asymptotic Nomality in Monte Calo Integation By Masashi Okamoto Abstact. To estimate a multiple integal of a function ove the unit cube, Habe poposed two Monte Calo estimatos /'j and J'2 based on 2N and 4N obsevations, espec- 2 2» tively, of the function. He also consideed estimatos Dy and D2 of the vaiances of/j and J'2, espectively. This pape shows that all these estimatos ae asymptotically nomally distibuted as N tends to infinity. 1. Intoduction. Monte Calo integation is a method to estimate the value of a definite integal of a given eal-valued function ove a finite egion (say, a cube) by obseving the value of the function only at a finite numbe of points in the egion which ae chosen suitably and stochastically. Kitagawa [4] poposed seveal estimating methods but he was concened mainly with the case when the function has a cetain pio distibution. Habe [1], [2] poposed a mesh estimato of the integal and then impoved it by means of the idea of "antithetic vaiâtes" due to Hammesley and Moton [3]. SpecificaUy, let /be a eal-valued function defined ove the unit cube Gs in the s- space and set / = jg f. Let A ( = 1,..., A) be a family of conguent subcubes aising by patitioning Gs so that the inteval [0, 1] on the x'-axis is divided into n equal subintevals fo each i = 1,..., s, N = nx ns. Let x be a andom point in A chosen independently fo each and let x' = 2c,. - x, c stands fo the cente of A. Then <1J> and '!-* /(*,) 1 f(x) + f(x') (1.2) 2 = A? Z-2- = 1 ae unbiased estimatos of/ based on N and 2A obsevations of/, espectively. To estimate the vaiances of/, and J2 we need eplications of obsevations. Let z be anothe andom point in A, the z being chosen independently othe and also of the x's. Define similaly z' 2c - z. Then (1.3),, 1 ± /(*,) + /(z) 1 *kx 2 and (1 4) /-' I V f(x) + f(x') + f(z) + f{z') 2 N x A of each Received August 26, 1975; evised Januay 27, AMS (MOS) subject classifications (1970). Pimay 65D30; Seconday 62E20, 60F05. Key wods and phases. Monte Calo integation, Habe's estimatos, asymptotic nomality. Lyapunov condition. 831 Copyight Ameican Mathematical Society
2 832 MASASHI OKAMOTO ae again unbiased estimatos of/ based on 2N and AN obsevations, espectively. using these andom points, the vaiances of J\ and f2 can be estimated unbiasedly by (L5) ^= 1 íf(x)-f(z)}2 and 4M '6) P' ' ^,?, j 2-1 j espectively. =\ Fo k = 1 and 2 let Ck denote the set of all eal-valued functions defined ove Gs and having continuous fcth ode patial deivatives. In the sequel we say just 'W o " to indicate that n fo evey i = I,..., s. Put n = min(jj,..., ns). Habe poved that if/e C1, then (1.7) va(/j) = T2N(Jt) + ottnn2)-1) as N -* -, By and also that if/ C2, then (1.8) va(/2) = t2n(j2) + o^a«4)"1) asjv^«2, L v ± f Z-^V TiV^l^ 12A ff-, î JGs\bxi) ' (1.9) T2.,_L_ j2 f i- f /-^-V "l 2' 1440A)2-4 Jc,^(xi)V + 5 ± _l_ ; pv_,,./=! (n.«.)2-^vaxw Futhemoe, since J'k (k = 1,2) is the aveage of two independent ealizations of Jk, its vaiance is half as lage as that of Jk. In estimating / by J'k, Habe used Dk as a measue of the eo, assuming that J'k is appoximately nomally distibuted. The pupose of this pape is to show, fist, that Jk (k = 1,2) is asymptotically nomally distibuted with mean / and vaiance T2N(Jk) as N >. This implies the asymptotic nomality of J'k. Next, it is shown that Dk(k= nomally distibuted with mean vax(j'k) and a pope vaiance. 1,2) is asymptotically 2. Asymptotic Nomality of /, and J2. Befoe stating the theoems we intoduce some notations. Let I=E{f(x)}=Nf /, K = K» fy - X, - C f,=f(c), t'-^ic) dx' and ff-jül-icu dxdx'
3 MONTE CARLO INTEGRATION 833 x = (x1,..., Xs), fo /,/ = 1,..., s, and = 1,..., N. The L2-nom will be denoted by. Theoem 1. Assume /EC1. As N,JX is asymptotically nomally distibuted with mean I and vaiance T2N(Jy) defined in (1.9). Poof. Since / = 2* j I/N, it holds that (2.1) A-' = Z {/(*, )-/,} A Taylo seies expansion of f(x) aound the point c gives N =i N (2.2) f(x)=f+ 5'//+/?l, í=i the emainde RX has the following popety in view of the unifom continuity of df/dx' in Gs: fo any e > 0, \Rl\ < e 6 fo evey, povided N is sufficiently lage. Again, by a Taylo expansion we have \R2\ < e/n fo evey, povided N is sufficiently lage. Substitution of (2.2) and (2.3) into (2.1) yields (2-4) J!-I = SN+RN, (2.5) ^=Â7ZtôX ^=Í (*l-*2)- / = 1 i= 1 = 1 Since X/,) = /, to pove the theoem it emains only to veify the following thee popositions: (OvaiS^-T^/^asA' «>, the symbol ~ means that the atio of the two sides tends to one, (ii) the sequence SN satisfies the Lyapunov condition of the cental limit theoem (see, e.g., Loève [5, p. 275]), and (iii) vai(rn) = o((nn2)'1) as N -* o. In fact, in the ight-hand side of the identity /, -/ SN - E(SN) RN - E(RN) TN TN TN the fist tem conveges to the standadized nomal distibution because of (i) and (ii), as the second tem conveges in pobability to zeo because of (iii) and the Chebyshev inequality. Now, pat (i) is essentially equivalent to Theoem 3 of Habe [1] o (1.7); but a poof is given hee fo the completeness of the poof. By definition, o' has a unifom distibution ove the inteval [- 1/(2«,), 1/(2«,.)] independently of each othe; and hence, va(s') = l/(12n2), which implies that
4 834 MASASH1 OKAMOTO vaís»,) = E i -T(^')2~^i) as TV-> oc. 12A2 ~\ ~\ nf Re (ii). Define XN = 2J=, «',//, then (2.6) - E(X2N)=N2va:(SN). Fom the inequality \ô'\ < l/(2n() it follows that (2.7) pu*,i» < (p/>/y~ f J& ax'" 0(Nn~3). By (i),(2.6) and (2.7) [Z^l^l3] /[Z^Ä)]372 = OÍTV-1/2) = 0(1), which is the Lyapunov condition fo 5^. Re (iii). Since R2 ae constant, va(/^) = va(/?l)//v2 < I>(/?2,.)//V2 < e2/(/v«2). This implies (iii), since e can be made abitaily small. Remak. Q.E.D. The asymptotic nomality may be poved by applying the cental limit theoem diectly to 7, in the fom (1.1), not indiectly to SN. This appoach, howeve, equies asymptotic expansions of E {f2(x)} and I up to the tems of ode n~2 so that a stonge assumption / G C2 is needed instead of /G C1. This emak is valid also with Theoems 2, 3 and 4 fo which a much stonge assumption / C* (k = 4, 4 and 6, espectively) is equied. Theoem 2. Assume /G C2.,4s A, /2 is asymptotically nomally distibuted with mean I and vaiance t2n(j2) defined in (1.9). Poof. Similaly as (2.1) (2.8) ^2-í = Z {f(x)+f(x')-2i}. =l Expanding f(x) and f(x') in Taylo seies aound the point c, we find (2.9) f(x) =f+í, Kf + \ i KKf + RU> i= 1 L ij= 1 f{x')=f-± Kfl + \ f S'M+K, i=l i,j=l \Rl\ and \R\\ ae bounded fom above by e 5 2 fo evey, povided N is sufficiently lage. Taking one moe tem in the expansion (2.3), we have (2.10) I=f+t 2fï+R2, í=i 24h? \R2\ < e/n2 fo evey. Substitution of (2.9) and (2.10) into (2.8) yields
5 MONTE CARLO INTEGRATION 835 ^2 ~~ I ~ $N + ^JV ' (2.11) ^=^E ^v =i \i,/=i /=i 12n? J 1 " =l (Äi + /i;-2ä2). Just as fo Theoem 1 we have only to pove the following thee popositions: (i) va^) ~ 2N(J2) as N -* -, (ii) SN satisfies the Lyapunov condition, and (iii) va(rn) = o((an4)-1) as N. Pat (i) is equivalent to the main theoem in Habe [2] o (1.8), while (ii) and (iii) can be veified by using easoning simila to that in the poof of Theoem 1. Q.E.D. Coollay 1. Fo k = 1 o 2, assume fgck. As N, Jk follows asymptotically a nomal distibution with mean I and vaiance (2.12) T%(j'k) = i 2N(Jk). 3. Asymptotic Nomality of Dx and D2. Though the asymptotic nomality of Dx o D2 may not be so impotant in pactice as that of any estimato of/ itself, it can be poved along a simila line of aguments fo the latte. Fist let us conside Theoem 3. Assume / G C1. As N, D\ is asymptotically nomally distibuted with mean va(/j) and vaiance 72( )2)= 1 Ulf /-V-Y + 10 T l f ibf bf\ N ' 2880A3) éik}jo'\9x*) tk^(nñ)2)g\^bxi) 2 W Poof. Similaly to (2.2) it is the case that i f(z)=f+'lsífí+r2> i=l = ttî,- > ) = *, -<V. \R2\ < «HM fo evey. Substitution of (2.2) and (3.1) into (1.5) yields (3-2) D2 = TN +RN, (3.3) 1 N / s tn = 2 ZÍZU). AN1 =\ \ =l and t = «', - V RN = -^ Z \(Rl-R2)2+^Rl-R2)t<f\ ANZ =\ ( 1=1 ]
6 836 MASASHI OKAMOTO Since E(D\) = va(/',), a poof of Theoem 3 can be educed to the following: (i) va(tn) ~ t2n(d2) as N-+, (ii) TN satisfies the Lyapunov condition, and (iii) va^) = o((n3n4yi) as N >. Re (i). Fist we have Since fo evey i, j and -, (3.4) (TN) = ^Z (ZvfJ- EW'^i- 6n( 15/14 36Wj.f and since the expectations of any othe monomial of the tj's of ode 2 o 4 vanish, we find i I I \ I <M-m ^m i fi \ 2 This implies Zvif-)2=^\iz(^ +ioz/ Substitution of the last fomula into (3.4) poves (i). Re (ii). Define Then a staightfowad and hence, calculation gives Y,E\XN\3 = 0(Nn~6); [pl^l3] /[Z^ÁV,)]3/2 = OÍA"1/2) = 0(1). Re (iii). Fom (3.3) it follows that l*jvi < ~h Z (H*ii + iifid2 + 2(iiô,ii + iifii) i Z l/íll; 4A2 * n i ) and hence, vaí/í^) < fi/?2,) = o((a V)"1). Q.E.D. Remak. Though t^/',) defined in (2.12) is indeed the leading tem in the asymptotic expansion of va(/j) in n~l, the phase "mean vai(j[)" in the statement
7 MONTE CARLO INTEGRATION 837 of Theoem 3 cannot be eplaced by "mean 2N(fxy\ in geneal. The eason is that the diffeence va(/',) - t2n(j\) is o((nn2y~l) and hence is not negligible in geneal as compaed with N(D2) which is of ode (A3'2«2)-1. Similaly, we can pove the following: Theoem 4. Assume f EC2. As N, D22 is asymptotically nomally distibuted with mean va(/2) and vaiance t2n(d\) which is a linea combination of seven tems of the fom y i f / a2/ a2/ \2 i,/,lb=i N3(ntnfttkn,)2 JGs\dx'dx> bxkdx' J some of i,j, k and I ae constained to be equal. The oot-squae tansfomation of the andom vaiable D\ (k = 1,2) induces the following: Coollay 2. Fo k = 1 o 2, assume /G Ck. As N, Dk is asymptotically nomally distibuted with mean o(fk) and vaiance lát2n(d2.)/yax(j'k), o(j'k) stands fo the standad deviation of J'k. Depatment of Applied Mathematics Faculty of Engineeing Science Osaka Univesity Toyonaka, Osaka 560, Japan 1. S. HABER, "A modified Monte-Calo quadatue," Math. Comp., v. 20, 1966, pp MR 35 # S. HABER, "A modified Monte-Calo quadatue. II," Math. Comp., v. 21, 1967, pp MR 38 # J. M. HAMMERSLEY & K. W. MORTON, "A new Monte Calo technique: Antithetic vaiâtes," Poc. Cambidge Philos. Soc, v. 52, 1956, pp MR 18, T. KITAGAWA, "Random integations," Bull. Math. Statist., v. 4, 1950, pp MR 14, M. LOÈVE, Pobability Theoy, 3d ed., Van Nostand, Pinceton, N.J., MR 34 # 3596.
A Multivariate Normal Law for Turing s Formulae
A Multivaiate Nomal Law fo Tuing s Fomulae Zhiyi Zhang Depatment of Mathematics and Statistics Univesity of Noth Caolina at Chalotte Chalotte, NC 28223 Abstact This pape establishes a sufficient condition
More information3.1 Random variables
3 Chapte III Random Vaiables 3 Random vaiables A sample space S may be difficult to descibe if the elements of S ae not numbes discuss how we can use a ule by which an element s of S may be associated
More informationJournal of Inequalities in Pure and Applied Mathematics
Jounal of Inequalities in Pue and Applied Mathematics COEFFICIENT INEQUALITY FOR A FUNCTION WHOSE DERIVATIVE HAS A POSITIVE REAL PART S. ABRAMOVICH, M. KLARIČIĆ BAKULA AND S. BANIĆ Depatment of Mathematics
More informationON LACUNARY INVARIANT SEQUENCE SPACES DEFINED BY A SEQUENCE OF MODULUS FUNCTIONS
STUDIA UNIV BABEŞ BOLYAI, MATHEMATICA, Volume XLVIII, Numbe 4, Decembe 2003 ON LACUNARY INVARIANT SEQUENCE SPACES DEFINED BY A SEQUENCE OF MODULUS FUNCTIONS VATAN KARAKAYA AND NECIP SIMSEK Abstact The
More informationOn the Poisson Approximation to the Negative Hypergeometric Distribution
BULLETIN of the Malaysian Mathematical Sciences Society http://mathusmmy/bulletin Bull Malays Math Sci Soc (2) 34(2) (2011), 331 336 On the Poisson Appoximation to the Negative Hypegeometic Distibution
More informationON INDEPENDENT SETS IN PURELY ATOMIC PROBABILITY SPACES WITH GEOMETRIC DISTRIBUTION. 1. Introduction. 1 r r. r k for every set E A, E \ {0},
ON INDEPENDENT SETS IN PURELY ATOMIC PROBABILITY SPACES WITH GEOMETRIC DISTRIBUTION E. J. IONASCU and A. A. STANCU Abstact. We ae inteested in constucting concete independent events in puely atomic pobability
More informationStanford University CS259Q: Quantum Computing Handout 8 Luca Trevisan October 18, 2012
Stanfod Univesity CS59Q: Quantum Computing Handout 8 Luca Tevisan Octobe 8, 0 Lectue 8 In which we use the quantum Fouie tansfom to solve the peiod-finding poblem. The Peiod Finding Poblem Let f : {0,...,
More informationNumerical approximation to ζ(2n+1)
Illinois Wesleyan Univesity Fom the SelectedWoks of Tian-Xiao He 6 Numeical appoximation to ζ(n+1) Tian-Xiao He, Illinois Wesleyan Univesity Michael J. Dancs Available at: https://woks.bepess.com/tian_xiao_he/6/
More informationq i i=1 p i ln p i Another measure, which proves a useful benchmark in our analysis, is the chi squared divergence of p, q, which is defined by
CSISZÁR f DIVERGENCE, OSTROWSKI S INEQUALITY AND MUTUAL INFORMATION S. S. DRAGOMIR, V. GLUŠČEVIĆ, AND C. E. M. PEARCE Abstact. The Ostowski integal inequality fo an absolutely continuous function is used
More informationHypothesis Test and Confidence Interval for the Negative Binomial Distribution via Coincidence: A Case for Rare Events
Intenational Jounal of Contempoay Mathematical Sciences Vol. 12, 2017, no. 5, 243-253 HIKARI Ltd, www.m-hikai.com https://doi.og/10.12988/ijcms.2017.7728 Hypothesis Test and Confidence Inteval fo the Negative
More information6 Matrix Concentration Bounds
6 Matix Concentation Bounds Concentation bounds ae inequalities that bound pobabilities of deviations by a andom vaiable fom some value, often its mean. Infomally, they show the pobability that a andom
More informationPearson s Chi-Square Test Modifications for Comparison of Unweighted and Weighted Histograms and Two Weighted Histograms
Peason s Chi-Squae Test Modifications fo Compaison of Unweighted and Weighted Histogams and Two Weighted Histogams Univesity of Akueyi, Bogi, v/noduslód, IS-6 Akueyi, Iceland E-mail: nikolai@unak.is Two
More informationA generalization of the Bernstein polynomials
A genealization of the Benstein polynomials Halil Ouç and Geoge M Phillips Mathematical Institute, Univesity of St Andews, Noth Haugh, St Andews, Fife KY16 9SS, Scotland Dedicated to Philip J Davis This
More informationLecture 28: Convergence of Random Variables and Related Theorems
EE50: Pobability Foundations fo Electical Enginees July-Novembe 205 Lectue 28: Convegence of Random Vaiables and Related Theoems Lectue:. Kishna Jagannathan Scibe: Gopal, Sudhasan, Ajay, Swamy, Kolla An
More informationMEASURES OF BLOCK DESIGN EFFICIENCY RECOVERING INTERBLOCK INFORMATION
MEASURES OF BLOCK DESIGN EFFICIENCY RECOVERING INTERBLOCK INFORMATION Walte T. Fedee 337 Waen Hall, Biometics Unit Conell Univesity Ithaca, NY 4853 and Tey P. Speed Division of Mathematics & Statistics,
More informationFunctions Defined on Fuzzy Real Numbers According to Zadeh s Extension
Intenational Mathematical Foum, 3, 2008, no. 16, 763-776 Functions Defined on Fuzzy Real Numbes Accoding to Zadeh s Extension Oma A. AbuAaqob, Nabil T. Shawagfeh and Oma A. AbuGhneim 1 Mathematics Depatment,
More informationarxiv: v1 [math.co] 4 May 2017
On The Numbe Of Unlabeled Bipatite Gaphs Abdullah Atmaca and A Yavuz Ouç axiv:7050800v [mathco] 4 May 207 Abstact This pape solves a poblem that was stated by M A Haison in 973 [] This poblem, that has
More informationMoment-free numerical approximation of highly oscillatory integrals with stationary points
Moment-fee numeical appoximation of highly oscillatoy integals with stationay points Sheehan Olve Abstact We pesent a method fo the numeical quadatue of highly oscillatoy integals with stationay points.
More informationarxiv: v1 [math.nt] 28 Oct 2017
ON th COEFFICIENT OF DIVISORS OF x n axiv:70049v [mathnt] 28 Oct 207 SAI TEJA SOMU Abstact Let,n be two natual numbes and let H(,n denote the maximal absolute value of th coefficient of divisos of x n
More informationSurveillance Points in High Dimensional Spaces
Société de Calcul Mathématique SA Tools fo decision help since 995 Suveillance Points in High Dimensional Spaces by Benad Beauzamy Januay 06 Abstact Let us conside any compute softwae, elying upon a lage
More informationIntroduction to Mathematical Statistics Robert V. Hogg Joeseph McKean Allen T. Craig Seventh Edition
Intoduction to Mathematical Statistics Robet V. Hogg Joeseph McKean Allen T. Caig Seventh Edition Peason Education Limited Edinbugh Gate Halow Essex CM2 2JE England and Associated Companies thoughout the
More informationChapter 3: Theory of Modular Arithmetic 38
Chapte 3: Theoy of Modula Aithmetic 38 Section D Chinese Remainde Theoem By the end of this section you will be able to pove the Chinese Remainde Theoem apply this theoem to solve simultaneous linea conguences
More informationOn the ratio of maximum and minimum degree in maximal intersecting families
On the atio of maximum and minimum degee in maximal intesecting families Zoltán Lóánt Nagy Lale Özkahya Balázs Patkós Máté Vize Septembe 5, 011 Abstact To study how balanced o unbalanced a maximal intesecting
More informationJournal of Number Theory
Jounal of umbe Theoy 3 2 2259 227 Contents lists available at ScienceDiect Jounal of umbe Theoy www.elsevie.com/locate/jnt Sums of poducts of hypegeometic Benoulli numbes Ken Kamano Depatment of Geneal
More informationKOEBE DOMAINS FOR THE CLASSES OF FUNCTIONS WITH RANGES INCLUDED IN GIVEN SETS
Jounal of Applied Analysis Vol. 14, No. 1 2008), pp. 43 52 KOEBE DOMAINS FOR THE CLASSES OF FUNCTIONS WITH RANGES INCLUDED IN GIVEN SETS L. KOCZAN and P. ZAPRAWA Received Mach 12, 2007 and, in evised fom,
More informationSOME GENERAL NUMERICAL RADIUS INEQUALITIES FOR THE OFF-DIAGONAL PARTS OF 2 2 OPERATOR MATRICES
italian jounal of pue and applied mathematics n. 35 015 (433 44) 433 SOME GENERAL NUMERICAL RADIUS INEQUALITIES FOR THE OFF-DIAGONAL PARTS OF OPERATOR MATRICES Watheq Bani-Domi Depatment of Mathematics
More informationAlternative Tests for the Poisson Distribution
Chiang Mai J Sci 015; 4() : 774-78 http://epgsciencecmuacth/ejounal/ Contibuted Pape Altenative Tests fo the Poisson Distibution Manad Khamkong*[a] and Pachitjianut Siipanich [b] [a] Depatment of Statistics,
More informationResults on the Commutative Neutrix Convolution Product Involving the Logarithmic Integral li(
Intenational Jounal of Scientific and Innovative Mathematical Reseach (IJSIMR) Volume 2, Issue 8, August 2014, PP 736-741 ISSN 2347-307X (Pint) & ISSN 2347-3142 (Online) www.acjounals.og Results on the
More informationNew problems in universal algebraic geometry illustrated by boolean equations
New poblems in univesal algebaic geomety illustated by boolean equations axiv:1611.00152v2 [math.ra] 25 Nov 2016 Atem N. Shevlyakov Novembe 28, 2016 Abstact We discuss new poblems in univesal algebaic
More information2 K. ENTACHER seies called Es classes, Koobov[4] developed the theoy of good lattice points. Recently, in a seies of papes, Lache et al. [6, 7, 8, 9]
BIT 37:(4) (997), 845{860. QUASI-MONTE CARLO METHODS FOR NUMERICAL INTEGRATION OF MULTIVARIATE HAAR SERIES KARL ENTACHER Depatment of Mathematics, Univesity of Salzbug, Hellbunnest. 34 A-5020 Salzbug,
More informationChapter 3 Optical Systems with Annular Pupils
Chapte 3 Optical Systems with Annula Pupils 3 INTRODUCTION In this chapte, we discuss the imaging popeties of a system with an annula pupil in a manne simila to those fo a system with a cicula pupil The
More informationNOTE. Some New Bounds for Cover-Free Families
Jounal of Combinatoial Theoy, Seies A 90, 224234 (2000) doi:10.1006jcta.1999.3036, available online at http:.idealibay.com on NOTE Some Ne Bounds fo Cove-Fee Families D. R. Stinson 1 and R. Wei Depatment
More informationCENTRAL INDEX BASED SOME COMPARATIVE GROWTH ANALYSIS OF COMPOSITE ENTIRE FUNCTIONS FROM THE VIEW POINT OF L -ORDER. Tanmay Biswas
J Koean Soc Math Educ Se B: Pue Appl Math ISSNPint 16-0657 https://doiog/107468/jksmeb01853193 ISSNOnline 87-6081 Volume 5, Numbe 3 August 018, Pages 193 01 CENTRAL INDEX BASED SOME COMPARATIVE GROWTH
More informationJENSEN S INEQUALITY FOR DISTRIBUTIONS POSSESSING HIGHER MOMENTS, WITH APPLICATION TO SHARP BOUNDS FOR LAPLACE-STIELTJES TRANSFORMS
J. Austal. Math. Soc. Se. B 40(1998), 80 85 JENSEN S INEQUALITY FO DISTIBUTIONS POSSESSING HIGHE MOMENTS, WITH APPLICATION TO SHAP BOUNDS FO LAPLACE-STIELTJES TANSFOMS B. GULJAŠ 1,C.E.M.PEACE 2 and J.
More informationCentral Coverage Bayes Prediction Intervals for the Generalized Pareto Distribution
Statistics Reseach Lettes Vol. Iss., Novembe Cental Coveage Bayes Pediction Intevals fo the Genealized Paeto Distibution Gyan Pakash Depatment of Community Medicine S. N. Medical College, Aga, U. P., India
More informationRegularity for Fully Nonlinear Elliptic Equations with Neumann Boundary Data
Communications in Patial Diffeential Equations, 31: 1227 1252, 2006 Copyight Taylo & Fancis Goup, LLC ISSN 0360-5302 pint/1532-4133 online DOI: 10.1080/03605300600634999 Regulaity fo Fully Nonlinea Elliptic
More information6 PROBABILITY GENERATING FUNCTIONS
6 PROBABILITY GENERATING FUNCTIONS Cetain deivations pesented in this couse have been somewhat heavy on algeba. Fo example, detemining the expectation of the Binomial distibution (page 5.1 tuned out to
More informationΣk=1. g r 3/2 z. 2 3-z. g 3 ( 3/2 ) g r 2. = 1 r = 0. () z = ( a ) + Σ. c n () a = ( a) 3-z -a. 3-z. z - + Σ. z 3, 5, 7, z ! = !
09 Maclauin Seies of Completed Riemann Zeta 9. Maclauin Seies of Lemma 9.. ( Maclauin seies of gamma function ) When is the gamma function, n is the polygamma function and B n,kf, f, ae Bell polynomials,
More informationOn a quantity that is analogous to potential and a theorem that relates to it
Su une quantité analogue au potential et su un théoème y elatif C R Acad Sci 7 (87) 34-39 On a quantity that is analogous to potential and a theoem that elates to it By R CLAUSIUS Tanslated by D H Delphenich
More informationarxiv: v1 [math.nt] 12 May 2017
SEQUENCES OF CONSECUTIVE HAPPY NUMBERS IN NEGATIVE BASES HELEN G. GRUNDMAN AND PAMELA E. HARRIS axiv:1705.04648v1 [math.nt] 12 May 2017 ABSTRACT. Fo b 2 and e 2, let S e,b : Z Z 0 be the function taking
More informationUsing Laplace Transform to Evaluate Improper Integrals Chii-Huei Yu
Available at https://edupediapublicationsog/jounals Volume 3 Issue 4 Febuay 216 Using Laplace Tansfom to Evaluate Impope Integals Chii-Huei Yu Depatment of Infomation Technology, Nan Jeon Univesity of
More informationOn the global uniform asymptotic stability of time-varying dynamical systems
Stud. Univ. Babeş-Bolyai Math. 59014), No. 1, 57 67 On the global unifom asymptotic stability of time-vaying dynamical systems Zaineb HajSalem, Mohamed Ali Hammami and Mohamed Mabouk Abstact. The objective
More informationQuasi-Randomness and the Distribution of Copies of a Fixed Graph
Quasi-Randomness and the Distibution of Copies of a Fixed Gaph Asaf Shapia Abstact We show that if a gaph G has the popety that all subsets of vetices of size n/4 contain the coect numbe of tiangles one
More informationSOME SOLVABILITY THEOREMS FOR NONLINEAR EQUATIONS
Fixed Point Theoy, Volume 5, No. 1, 2004, 71-80 http://www.math.ubbcluj.o/ nodeacj/sfptcj.htm SOME SOLVABILITY THEOREMS FOR NONLINEAR EQUATIONS G. ISAC 1 AND C. AVRAMESCU 2 1 Depatment of Mathematics Royal
More informationMultiple Criteria Secretary Problem: A New Approach
J. Stat. Appl. Po. 3, o., 9-38 (04 9 Jounal of Statistics Applications & Pobability An Intenational Jounal http://dx.doi.og/0.785/jsap/0303 Multiple Citeia Secetay Poblem: A ew Appoach Alaka Padhye, and
More informationBEST CONSTANTS FOR UNCENTERED MAXIMAL FUNCTIONS. Loukas Grafakos and Stephen Montgomery-Smith University of Missouri, Columbia
BEST CONSTANTS FOR UNCENTERED MAXIMAL FUNCTIONS Loukas Gafakos and Stehen Montgomey-Smith Univesity of Missoui, Columbia Abstact. We ecisely evaluate the oeato nom of the uncenteed Hady-Littlewood maximal
More informationA NOTE ON ROTATIONS AND INTERVAL EXCHANGE TRANSFORMATIONS ON 3-INTERVALS KARMA DAJANI
A NOTE ON ROTATIONS AND INTERVAL EXCHANGE TRANSFORMATIONS ON 3-INTERVALS KARMA DAJANI Abstact. Wepove the conjectue that an inteval exchange tansfomation on 3-intevals with coesponding pemutation (1; 2;
More informationPROBLEM SET #1 SOLUTIONS by Robert A. DiStasio Jr.
POBLM S # SOLUIONS by obet A. DiStasio J. Q. he Bon-Oppenheime appoximation is the standad way of appoximating the gound state of a molecula system. Wite down the conditions that detemine the tonic and
More informationOn the ratio of maximum and minimum degree in maximal intersecting families
On the atio of maximum and minimum degee in maximal intesecting families Zoltán Lóánt Nagy Lale Özkahya Balázs Patkós Máté Vize Mach 6, 013 Abstact To study how balanced o unbalanced a maximal intesecting
More informationA New Method of Estimation of Size-Biased Generalized Logarithmic Series Distribution
The Open Statistics and Pobability Jounal, 9,, - A New Method of Estimation of Size-Bied Genealized Logaithmic Seies Distibution Open Access Khushid Ahmad Mi * Depatment of Statistics, Govt Degee College
More informationand the initial value R 0 = 0, 0 = fall equivalence classes ae singletons fig; i = 1; : : : ; ng: (3) Since the tansition pobability p := P (R = j R?1
A CLASSIFICATION OF COALESCENT PROCESSES FOR HAPLOID ECHANGE- ABLE POPULATION MODELS Matin Mohle, Johannes Gutenbeg-Univesitat, Mainz and Seik Sagitov 1, Chalmes and Gotebogs Univesities, Gotebog Abstact
More information( ) [ ] [ ] [ ] δf φ = F φ+δφ F. xdx.
9. LAGRANGIAN OF THE ELECTROMAGNETIC FIELD In the pevious section the Lagangian and Hamiltonian of an ensemble of point paticles was developed. This appoach is based on a qt. This discete fomulation can
More informationApplication of Parseval s Theorem on Evaluating Some Definite Integrals
Tukish Jounal of Analysis and Numbe Theoy, 4, Vol., No., -5 Available online at http://pubs.sciepub.com/tjant/// Science and Education Publishing DOI:.69/tjant--- Application of Paseval s Theoem on Evaluating
More informationBrief summary of functional analysis APPM 5440 Fall 2014 Applied Analysis
Bief summay of functional analysis APPM 5440 Fall 014 Applied Analysis Stephen Becke, stephen.becke@coloado.edu Standad theoems. When necessay, I used Royden s and Keyzsig s books as a efeence. Vesion
More informationRandom Variables and Probability Distribution Random Variable
Random Vaiables and Pobability Distibution Random Vaiable Random vaiable: If S is the sample space P(S) is the powe set of the sample space, P is the pobability of the function then (S, P(S), P) is called
More informationk. s k=1 Part of the significance of the Riemann zeta-function stems from Theorem 9.2. If s > 1 then 1 p s
9 Pimes in aithmetic ogession Definition 9 The Riemann zeta-function ζs) is the function which assigns to a eal numbe s > the convegent seies k s k Pat of the significance of the Riemann zeta-function
More information556: MATHEMATICAL STATISTICS I
556: MATHEMATICAL STATISTICS I CHAPTER 5: STOCHASTIC CONVERGENCE The following efinitions ae state in tems of scala anom vaiables, but exten natually to vecto anom vaiables efine on the same obability
More informationEncapsulation theory: radial encapsulation. Edmund Kirwan *
Encapsulation theoy: adial encapsulation. Edmund Kiwan * www.edmundkiwan.com Abstact This pape intoduces the concept of adial encapsulation, wheeby dependencies ae constained to act fom subsets towads
More informationChapter 13 Gravitation
Chapte 13 Gavitation In this chapte we will exploe the following topics: -Newton s law of gavitation, which descibes the attactive foce between two point masses and its application to extended objects
More informationA Multivariate Normal Law for Turing s Formulae
A Multivaiate Nomal Law fo Tuing s Fomulae Zhiyi Zhang Depatment of Mathematics and Statistics Univesity of Noth Caolina at Chalotte Chalotte, NC 28223 Abstact This pape establishes a sufficient condition
More informationThe second law of thermodynamics - II.
Januay 21, 2013 The second law of themodynamics - II. Asaf Pe e 1 1. The Schottky defect At absolute zeo tempeatue, the atoms of a solid ae odeed completely egulaly on a cystal lattice. As the tempeatue
More informationOn Polynomials Construction
Intenational Jounal of Mathematical Analysis Vol., 08, no. 6, 5-57 HIKARI Ltd, www.m-hikai.com https://doi.og/0.988/ima.08.843 On Polynomials Constuction E. O. Adeyefa Depatment of Mathematics, Fedeal
More informationWhat Form of Gravitation Ensures Weakened Kepler s Third Law?
Bulletin of Aichi Univ. of Education, 6(Natual Sciences, pp. - 6, Mach, 03 What Fom of Gavitation Ensues Weakened Keple s Thid Law? Kenzi ODANI Depatment of Mathematics Education, Aichi Univesity of Education,
More informationTHE NUMBER OF TWO CONSECUTIVE SUCCESSES IN A HOPPE-PÓLYA URN
TH NUMBR OF TWO CONSCUTIV SUCCSSS IN A HOPP-PÓLYA URN LARS HOLST Depatment of Mathematics, Royal Institute of Technology S 100 44 Stocholm, Sweden -mail: lholst@math.th.se Novembe 27, 2007 Abstact In a
More informationKR- 21 FOR FORMULA SCORED TESTS WITH. Robert L. Linn, Robert F. Boldt, Ronald L. Flaugher, and Donald A. Rock
RB-66-4D ~ E S [ B A U R L C L Ii E TI KR- 21 FOR FORMULA SCORED TESTS WITH OMITS SCORED AS WRONG Robet L. Linn, Robet F. Boldt, Ronald L. Flaughe, and Donald A. Rock N This Bulletin is a daft fo inteoffice
More informationUnobserved Correlation in Ascending Auctions: Example And Extensions
Unobseved Coelation in Ascending Auctions: Example And Extensions Daniel Quint Univesity of Wisconsin Novembe 2009 Intoduction In pivate-value ascending auctions, the winning bidde s willingness to pay
More informationarxiv: v1 [math.ca] 31 Aug 2009
axiv:98.4578v [math.ca] 3 Aug 9 On L-convegence of tigonometic seies Bogdan Szal Univesity of Zielona Góa, Faculty of Mathematics, Compute Science and Econometics, 65-56 Zielona Góa, ul. Szafana 4a, Poland
More informationA pathway to matrix-variate gamma and normal densities
Linea Algeba and its Applications 396 005 317 38 www.elsevie.com/locate/laa A pathway to matix-vaiate gamma and nomal densities A.M. Mathai Depatment of Mathematics and Statistics, McGill Univesity, 805
More informationResearch Article On Alzer and Qiu s Conjecture for Complete Elliptic Integral and Inverse Hyperbolic Tangent Function
Abstact and Applied Analysis Volume 011, Aticle ID 697547, 7 pages doi:10.1155/011/697547 Reseach Aticle On Alze and Qiu s Conjectue fo Complete Elliptic Integal and Invese Hypebolic Tangent Function Yu-Ming
More informationA STABILITY RESULT FOR p-harmonic SYSTEMS WITH DISCONTINUOUS COEFFICIENTS. Bianca Stroffolini. 0. Introduction
Electonic Jounal of Diffeential Equations, Vol. 2001(2001), No. 02, pp. 1 7. ISSN: 1072-6691. URL: http://ejde.math.swt.edu o http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) A STABILITY RESULT
More informationConcomitants of Multivariate Order Statistics With Application to Judgment Poststratification
JASA jasa v.00/1/09 Pn:0/06/006; 16:6 F:jasatm05511.tex; (R p. 1 Concomitants of Multivaiate Ode Statistics With Application to Judgment Poststatification Xinlei WANG, Lynne STOKES, Johan LIM, and Min
More informationOn the Quasi-inverse of a Non-square Matrix: An Infinite Solution
Applied Mathematical Sciences, Vol 11, 2017, no 27, 1337-1351 HIKARI Ltd, wwwm-hikaicom https://doiog/1012988/ams20177273 On the Quasi-invese of a Non-squae Matix: An Infinite Solution Ruben D Codeo J
More informationSolving Some Definite Integrals Using Parseval s Theorem
Ameican Jounal of Numeical Analysis 4 Vol. No. 6-64 Available online at http://pubs.sciepub.com/ajna///5 Science and Education Publishing DOI:.69/ajna---5 Solving Some Definite Integals Using Paseval s
More informationSTUDY OF SOLUTIONS OF LOGARITHMIC ORDER TO HIGHER ORDER LINEAR DIFFERENTIAL-DIFFERENCE EQUATIONS WITH COEFFICIENTS HAVING THE SAME LOGARITHMIC ORDER
UNIVERSITATIS IAGELLONICAE ACTA MATHEMATICA doi: 104467/20843828AM170027078 542017, 15 32 STUDY OF SOLUTIONS OF LOGARITHMIC ORDER TO HIGHER ORDER LINEAR DIFFERENTIAL-DIFFERENCE EQUATIONS WITH COEFFICIENTS
More informationBerkeley Math Circle AIME Preparation March 5, 2013
Algeba Toolkit Rules of Thumb. Make sue that you can pove all fomulas you use. This is even bette than memoizing the fomulas. Although it is best to memoize, as well. Stive fo elegant, economical methods.
More informationSecret Exponent Attacks on RSA-type Schemes with Moduli N = p r q
Secet Exponent Attacks on RSA-type Schemes with Moduli N = p q Alexande May Faculty of Compute Science, Electical Engineeing and Mathematics Univesity of Padebon 33102 Padebon, Gemany alexx@uni-padebon.de
More informationBoundedness for Marcinkiewicz integrals associated with Schrödinger operators
Poc. Indian Acad. Sci. (Math. Sci. Vol. 24, No. 2, May 24, pp. 93 23. c Indian Academy of Sciences oundedness fo Macinkiewicz integals associated with Schödinge opeatos WENHUA GAO and LIN TANG 2 School
More informationGROWTH ESTIMATES THROUGH SCALING FOR QUASILINEAR PARTIAL DIFFERENTIAL EQUATIONS
Annales Academiæ Scientiaum Fennicæ Mathematica Volumen 32, 2007, 595 599 GROWTH ESTIMATES THROUGH SCALING FOR QUASILINEAR PARTIAL DIFFERENTIAL EQUATIONS Teo Kilpeläinen, Henik Shahgholian and Xiao Zhong
More informationarxiv: v1 [math.ca] 12 Mar 2015
axiv:503.0356v [math.ca] 2 Ma 205 AN APPLICATION OF FOURIER ANALYSIS TO RIEMANN SUMS TRISTRAM DE PIRO Abstact. We develop a method fo calculating Riemann sums using Fouie analysis.. Poisson Summation Fomula
More informationarxiv: v1 [math.co] 1 Apr 2011
Weight enumeation of codes fom finite spaces Relinde Juius Octobe 23, 2018 axiv:1104.0172v1 [math.co] 1 Ap 2011 Abstact We study the genealized and extended weight enumeato of the - ay Simplex code and
More informationMitscherlich s Law: Sum of two exponential Processes; Conclusions 2009, 1 st July
Mitschelich s Law: Sum of two exponential Pocesses; Conclusions 29, st July Hans Schneebege Institute of Statistics, Univesity of Elangen-Nünbeg, Gemany Summay It will be shown, that Mitschelich s fomula,
More informationSolution to HW 3, Ma 1a Fall 2016
Solution to HW 3, Ma a Fall 206 Section 2. Execise 2: Let C be a subset of the eal numbes consisting of those eal numbes x having the popety that evey digit in the decimal expansion of x is, 3, 5, o 7.
More informationEstimation of the Correlation Coefficient for a Bivariate Normal Distribution with Missing Data
Kasetsat J. (Nat. Sci. 45 : 736-74 ( Estimation of the Coelation Coefficient fo a Bivaiate Nomal Distibution with Missing Data Juthaphon Sinsomboonthong* ABSTRACT This study poposes an estimato of the
More informationLET a random variable x follows the two - parameter
INTERNATIONAL JOURNAL OF MATHEMATICS AND SCIENTIFIC COMPUTING ISSN: 2231-5330, VOL. 5, NO. 1, 2015 19 Shinkage Bayesian Appoach in Item - Failue Gamma Data In Pesence of Pio Point Guess Value Gyan Pakash
More informationA Power Method for Computing Square Roots of Complex Matrices
JOURNAL OF MATHEMATICAL ANALYSIS AND APPLICATIONS 13, 39345 1997 ARTICLE NO. AY975517 A Powe Method fo Computing Squae Roots of Complex Matices Mohammed A. Hasan Depatment of Electical Engineeing, Coloado
More informationTemporal-Difference Learning
.997 Decision-Making in Lage-Scale Systems Mach 17 MIT, Sping 004 Handout #17 Lectue Note 13 1 Tempoal-Diffeence Leaning We now conside the poblem of computing an appopiate paamete, so that, given an appoximation
More informationFractional Zero Forcing via Three-color Forcing Games
Factional Zeo Focing via Thee-colo Focing Games Leslie Hogben Kevin F. Palmowski David E. Robeson Michael Young May 13, 2015 Abstact An -fold analogue of the positive semidefinite zeo focing pocess that
More informationLocalization of Eigenvalues in Small Specified Regions of Complex Plane by State Feedback Matrix
Jounal of Sciences, Islamic Republic of Ian (): - () Univesity of Tehan, ISSN - http://sciencesutaci Localization of Eigenvalues in Small Specified Regions of Complex Plane by State Feedback Matix H Ahsani
More informationSUFFICIENT CONDITIONS FOR MAXIMALLY EDGE-CONNECTED AND SUPER-EDGE-CONNECTED GRAPHS DEPENDING ON THE CLIQUE NUMBER
Discussiones Mathematicae Gaph Theoy 39 (019) 567 573 doi:10.7151/dmgt.096 SUFFICIENT CONDITIONS FOR MAXIMALLY EDGE-CONNECTED AND SUPER-EDGE-CONNECTED GRAPHS DEPENDING ON THE CLIQUE NUMBER Lutz Volkmann
More informationRevision of Lecture Eight
Revision of Lectue Eight Baseband equivalent system and equiements of optimal tansmit and eceive filteing: (1) achieve zeo ISI, and () maximise the eceive SNR Thee detection schemes: Theshold detection
More informationExceptional regular singular points of second-order ODEs. 1. Solving second-order ODEs
(May 14, 2011 Exceptional egula singula points of second-ode ODEs Paul Gaett gaett@math.umn.edu http://www.math.umn.edu/ gaett/ 1. Solving second-ode ODEs 2. Examples 3. Convegence Fobenius method fo solving
More informationON THE INVERSE SIGNED TOTAL DOMINATION NUMBER IN GRAPHS. D.A. Mojdeh and B. Samadi
Opuscula Math. 37, no. 3 (017), 447 456 http://dx.doi.og/10.7494/opmath.017.37.3.447 Opuscula Mathematica ON THE INVERSE SIGNED TOTAL DOMINATION NUMBER IN GRAPHS D.A. Mojdeh and B. Samadi Communicated
More informationIntegral operator defined by q-analogue of Liu-Srivastava operator
Stud. Univ. Babeş-Bolyai Math. 582013, No. 4, 529 537 Integal opeato defined by q-analogue of Liu-Sivastava opeato Huda Aldweby and Maslina Daus Abstact. In this pape, we shall give an application of q-analogues
More informationOn the integration of the equations of hydrodynamics
Uebe die Integation de hydodynamischen Gleichungen J f eine u angew Math 56 (859) -0 On the integation of the equations of hydodynamics (By A Clebsch at Calsuhe) Tanslated by D H Delphenich In a pevious
More informationEnergy Levels Of Hydrogen Atom Using Ladder Operators. Ava Khamseh Supervisor: Dr. Brian Pendleton The University of Edinburgh August 2011
Enegy Levels Of Hydogen Atom Using Ladde Opeatos Ava Khamseh Supeviso: D. Bian Pendleton The Univesity of Edinbugh August 11 1 Abstact The aim of this pape is to fist use the Schödinge wavefunction methods
More informationMeasure Estimates of Nodal Sets of Polyharmonic Functions
Chin. Ann. Math. Se. B 39(5), 08, 97 93 DOI: 0.007/s40-08-004-6 Chinese Annals of Mathematics, Seies B c The Editoial Office of CAM and Spinge-Velag Belin Heidelbeg 08 Measue Estimates of Nodal Sets of
More informationSemicanonical basis generators of the cluster algebra of type A (1)
Semicanonical basis geneatos of the cluste algeba of type A (1 1 Andei Zelevinsky Depatment of Mathematics Notheasten Univesity, Boston, USA andei@neu.edu Submitted: Jul 7, 006; Accepted: Dec 3, 006; Published:
More informationarxiv: v2 [astro-ph] 16 May 2008
New Anomalies in Cosmic Micowave Backgound Anisotopy: Violation of the Isotopic Gaussian Hypothesis in Low-l Modes Shi Chun, Su and M.-C., Chu Depatment of Physics and Institute of Theoetical Physics,
More information9.1 The multiplicative group of a finite field. Theorem 9.1. The multiplicative group F of a finite field is cyclic.
Chapte 9 Pimitive Roots 9.1 The multiplicative goup of a finite fld Theoem 9.1. The multiplicative goup F of a finite fld is cyclic. Remak: In paticula, if p is a pime then (Z/p) is cyclic. In fact, this
More informationMathematisch-Naturwissenschaftliche Fakultät I Humboldt-Universität zu Berlin Institut für Physik Physikalisches Grundpraktikum.
Mathematisch-Natuwissenschaftliche Fakultät I Humboldt-Univesität zu Belin Institut fü Physik Physikalisches Gundpaktikum Vesuchspotokoll Polaisation duch Reflexion (O11) duchgefüht am 10.11.2009 mit Vesuchspatne
More information