Quantile function expansion using regularly varying functions
|
|
- Nora Campbell
- 6 years ago
- Views:
Transcription
1 Quantile function expansion using regularly varying functions arxiv: v [math.st] 9 Aug 07 Thomas Fung a, an Eugene Seneta b a Department of Statistics, Macquarie University, NSW 09, Australia b School of Mathematics an Statistics, University of Syney, NSW 006, Australia August 0, 07 Abstract We present a simple result that allows us to evaluate the asymptotic orer of the remainer of a partial asymptotic expansion of the quantile function hu as u 0 + or. This is focusse on important univariate istributions when h has no simple close form, with a view to assessing asymptotic rate of ecay to zero of tail epenence in the context of bivariate copulas. The Introuction motivates the stuy in terms of the stanar Normal. The Normal, Skew-Normal an Gamma are use as initial examples. Finally, we iscuss approximation to the lower quantile of the Variance-Gamma an Skew-Slash istributions. Keywors: Asymptotic expansion; asymptotic tail epenence; Quantile function; Regularly varying functions; Skew-Slash istribution; Variance-Gamma istribution. Corresponing Author. Honorary Associate, University of Syney. aress: thomas.fung@mq.eu.au Thomas Fung.
2 Introuction This paper is motivate by the nee for a generally applicable proceure to stuy the asymptotic behaviour as u 0 + of F i u, i =, an Cu, u = P X F u, X F u where the F u s are the inverse of continuous an strictly increasing cf s i F i u, i =,, an a bivariate copula function respectively. A ranom vector X = X, X with marginal inverse istribution function F i u, i =, has coefficient of lower tail epenence λ L if the limit λ L = lim u 0 + λ L u exists, where λ L u = P X F u X F u = P X F u, X F u P Cu, u X F =. u u If λ L = 0 then X is sai to be asymptotically inepenent in the lower tail. In this situation in particular the asymptotic rate of approach to the limit 0 of λ L u is an inication of the strength of asymptotic inepenence. The classical case is the bivariate Normal with correlation coefficient ρ as iscusse in Embrechts, McNeil an Straumann 0. It was shown in Fung an Seneta 0 that if ρ λu = Φ Φ u u θ Lu + ρ where Φx, < x < is the cf of the stanar Normal, an Lu is a slowly varying function as u 0 +, then λ L u uθ Lu θ +. Fung an Seneta 0 also showe in their Theorem 3 that λu u ρ +ρ Lu where + ρ ρ Lu 4π log u +ρ. ρ
3 The proof within their Theorem epene heavily on the very specific asymptotic relation as x between the cf Φ an the corresponing stanar Normal pf f. We were unable to use in that paper, for this purpose, the expression for the quantile function Φ u of the stanar Normal istribution yu = logu 4π log u, 3 given for example by Lefor an Tawn 997 from a truncate expansion of Φ u, since we i not know the asymptotic orer of the reminer Φ u yu. We were able to prove that Φ u yu as u 0 +, but inasmuch as this asserts only that Φ u = yu + o = log u + o, we can only be sure of the ominant term of a truncate asymptotic expansion, an this was inaequate to procee from. Since However our general result below, when applie to the stanar Normal, gives Φ u = yu + O Φy y π e y log u. as y, from λu kue yu ρ +ρ +O log u where ku = log u π +ρ ρ ; = kue ρ +ρy u +O log u 3
4 = kue ρ +ρy u +O log u = kue ρ +ρy u+o as u 0 +, since y uo = O log u Thus, since y u log u as u 0 +. log u λu kue +ρ ρy u an the right han sie simplifies to the right han sie of. In this note we present a simple result that allows us to evaluate the asymptotic orer of the ifference between yu an hu as u 0 + or, where hu is the quantile function corresponing to a cumulative istribution function g for important istributions where hu has no close form, an yu is an asymptotic close form expression. This is the general result which we iscuss in Section. In Section 3, we will illustrate our results by consiering the quantile function for the Generalise Gamma-type tail which has various commonly use istributions such as Normal, Skew-Normal, Gamma, Variance-Gamma an a Skew-Slash as special cases. Detaile extreme value structure of such istributions is important in a financial mathematics context. These iniviual examples will be iscusse in Section 4. Main Result Our main result is summarise into the following Theorem. Theorem. Suppose that gx is a strictly positive continuous an strictly increasing cumulative istribution function cf on, A], A < 0. Suppose further that some function yu as u 0 + satisfies ygx = x + O ζx 4 4
5 as x, such that ζx 0, x. Suppose finally that ζx = ψ x with ψw = w ρ Lw for some constant ρ 0 an function Lw, w > 0, slowly varying at 0. If hu, u 0, ga] is the inverse function of g, then hu = yu + Oζyu. Proof. We begin with the fact that ghu = u yghu = yu hu + Oζhu = yu, using 4 Thus hu yu as u 0+ an hence ζhu ζyu = ψ /hu ψ /yu 5 by the Uniform Convergence Theorem of slowly varying function of Seneta 976, Theorem.. From 5, ζhu ζyu O ζhu = O ζyu as u 0+. Finally, hu + Oζhu = yu hu = yu + Oζhu hu = yu + Oζhu hu = yu + Oζyu. 5
6 Note that the formulation of the theorem requires, in the case ρ = 0 only, that the function Lw 0, w 0. The conition require in Theorem is simply that the correction term in 4 is relate to a regular varying function which is quite general an shoul apply to a wie range of istributions. The result in Theorem not only ensures hu an yu will be asymptotically equivalent, it will also stipulate how accurate yu will be for hu. We express as a corollary to Theorem the corresponing result for the upper tail quantiles. Corollary. Suppose that gx is a strictly positive continuous an strictly increasing cumulative istribution function cf on [A,, A > 0. Suppose further that some function yu as u satisfies ygx = x + O ζx as x, such that ζx 0, x. If ζx = ψ x with ψw = w ρ Lw for some constant ρ 0 an function Lw, w > 0, slowly varying at 0. If hu = g u, u [ga, is the inverse function of g i.e. the upper tail quantile function then hu = yu + Oζyu. In the next section, we will illustrate our results by consiering the quantile function corresponing to a cf that has a Generalise Gamma-type tail behaviour. 6
7 3 Quantile for Generalise Gamma-type tail behaviour 3. Lower Tail We are intereste in approximation to the quantile function for the Generalise Gammatype tail. We consier a cf g to have a Generalise Gamma-type lower tail behaviour if g can be expresse as gx = a x b e c x + O, x < 0, 6 x e for some constants a, c,, e > 0 an b R. Several istributions of current interest have such tails, an we iscuss them as special cases in the next section. Suppose that an approximation to the quantile function hu = g u to be b yu = c c u b b a log log c b u b, for small u > 0. 7 This can be obtaine via the recursive metho for an inverse function see Chapter.4 of De Bruijn 96 for instance on g. Then ygx = = = { b c b c log c log b log c b [ b c c b x + log a x b e c x +O x e a a x b e c x + O b b x e c b x e c b x + O x e log ca b x b e c b x + O x e c b x + log + O x e log c b x + log ca b b x + log + O x e ]} 7
8 = { [ b c c b x + log log log c b x + O x e ca b b x + O c b x x e ]} { [ b = c ]} log x c b x + O + O x x e log x = x + O + O x e+ x x + O log x, if e x = x + O, otherwise. x +e c log b x As yu c log u, the behaviour of the quantile function h is yu + O hu = yu + O log c log u log u log u e + = yu + O log log u, if e log u, otherwise, 8 by Theorem. We next expan yu to obtain an expansion for hu with an error term of appropriate lower orer. After some algebra, yu = c [ a b log b log u b c log u log u + b 3 log a b ] 3 +O. log u 3 c log u + b 3 + b 3 log u log [ log a b c log u c b ] 8
9 As a result, when e, h becomes log log u hu =yu + O log u = c a b log b log u b c c log u c log c c u b c 3 log + O ; 9 c u log u On the other han, when > e, h becomes hu =yu + O log u e + = c log u b + O log u e + c log u c. a b log b c c log u c 3. Upper Tail Similarly, a cf g is sai to have a Generalise Gamma-type upper tail behaviour if g can be expresse as gx = ax b e cx + O, x, 0 x e 9
10 for some constants a, c,, e > 0 an b R, which woul suggest the following approximation to the upper tail quantile function: b yu = c c u b b a log log c u b b, as u. Then by Corollary, the upper tail quantile function h can be expresse as hu yu + O = yu + O log log u log u c log u + = log u + c +O log u e +, if e;, if < e; b log log u c c log u + b log log u + O c 3 c log u log u e + b log log u c c log u + log log u log u b log a b c c c log u, if e; b log + c c log u, if > e. a b c as u. 4 Applications In this section, we will illustrate our results with application to the Normal, Skew-Normal, Gamma, Variance-Gamma an Skew-Slash istributions. 0
11 4. Stanar Normal From Feller 968 Chapter VII Lemma, the form of g is gx = e x + O, for x <. π x x If we compare this with 6, we have a = π, b =, c = that the approximation y becomes: an = e =. This means yu = { log πu logu } = log u 4π log u 3 an is the same as 3. As in this case we have = e, we can obtain the behaviour of the quantile function by using using 8 an 9 an get hu =yu + O log u = log log u log u + log 4π + log u log u + + O. 8 log u 3 log u 3 If we use only the ominating term of yu in 3 an so put y u = log u then after some algebra we have ygx = x + Oζx where ζx = log x x. As a result, hu = y u + O log u
12 by using Theorem. We can see that y u is a less accurate approximation to hu than yu, an the expression for hu that it gives only improves on y u by giving the correct orer of the ifference hu y u. We are aware that there exists a whole fiel of literature on fining an efficient an accurate approximation in the numerical sense for the Normal quantile functions, such as Abramowitz an Stegun 964, Beasley an Springer 977 to the more recent Voutier 00. Soranzo an Epure 04 provie a substantial bibliography on this subject. But that is not our focus an so our methoology will most likely be outperforme by the more sophisticate approximation in the literature. Take the approximation to the Normal quantile in Section.. in Voutier 00 for example which gives y V u = c 3 log u + c c + log u + c 0 log u + log u + 0 for e 37 / < u < where c 3 = , c 0 = , c = , c = , 0 = , = To see how the expansions perform against the approximations, we plot the stanar Normal quantile function, enote as Φ u, in R against y V, y an y an the results are shown in Fig.. From Fig., we can see that h V gives almost ientical result to the built-in Normal quantile function in R an is better than our expansion yu. We also see that y is an orer worse than the other methos.
13 Figure : Comparison in R of Φ u, y V, y an y 4. Skew-Normal The quantile function of the Skew-Normal istribution is another example that is covere by our result. The istribution was first introuce by Azzalini 985 an has evelope into an extensive theory presente in a recent monograph by Azzalini an Capitanio 04. A ranom variable is sai to have a stanar Skew-Normal istribution if its ensity function is fx = π e x Φλx, x R where λ R controls the skewness of the istribution. When λ = 0, it reuces to the stanar Normal as special case so we shall exclue the case of λ = 0 from our subsequent iscussion. Using Lemma of Capitanio 00 or Azzalini an Capitanio 04, pp. 5 3
14 53, we have e πλ+λ gx = x +λ x + O x, λ > 0; π x e x + O x, λ < 0, x. 4 This means if we compare 4 with 6, we have a = πλ+λ, b =, c = + λ / an = e = when λ > 0; a = / π, b =, c = /, = e = when λ < 0. This means that the approximation y becomes: { yu = = +λ { log [ log u π/ logu +λ uπλ+λ log +λ u ]}, λ > 0; }, λ < 0; log πλu logu/ + λ +λ, if λ > 0; log 5 u 4π log u, if λ < 0. As in this case we have = e, we can obtain the behaviour of the quantile function by using 8 an 9 to get hu =yu + O log u log u + +λ = + +λ +λ log u 3 + O log u + + O + 8 log u 3 log log u + logπλ +λ +λ log u +λ +λ log u log u 3 log log u + log π log u log u log u 3, λ > 0;, λ < 0. 6 The expression 6 justifies the use of expression 5 above in Theorem of Fung an Seneta 06. 4
15 4.3 Gamma Using, one can also etermine the accuracy of the gamma-like upper tail. Suppose that X Γα, β an let its pf be Then fx = βα Γα xα e βx, x > 0. gx = βα Γα xα e βx + O, 7 x as x. This means that if we compare 7 with 0, we have a = βα Γα, b = α, c = β an = e =. By using an approximation to the upper quantile function when u woul be α yu = β β α log log u α β α /Γα β α = β log uγα uβ α α log α α u α As in this case we have = e again, we have for the upper quantile function hu = yu + O log log u log u = β log uγα uβ α α log = β as u by using. α α log u + α β + O log log u log u log log u log log u β log Γα + O log u In Section 4 of Fung an Seneta 0, via their Theorem, the authors obtaine hu = yu + o, u with yu = log u. β 5
16 4.4 Variance-Gamma an Skew-Slash Finally, we propose an approximation to the lower tail quantile function of the Variance- Gamma VG an a Skew-Slash istribution, as they can have similar tail structure. A ranom variable X is sai to have a skew VG istribution introuce in Maan, Carr an Chang 998 an further stuie in Schoutens 003, Seneta 004 an Tjetjep an Seneta 006, if X is efine by a Normal variance-mean mixture as X Y Nµ + θy, σ Y, where µ, θ R an σ > 0. The istribution of Y is Γ ν, ν with ν > 0 such that EY =. When µ = 0 an σ =, from Schoutens 003, we have gx = ν ν Γ ν x ν e ν +θ +θ x ν + θ ν = a x ν e ν +θ +θ x + ν θ + θ + O x + O x 8 as x, where a > 0 is efine by the above.. By comparing with 6, a is efine as above, b = ν, c = ν + θ + θ an = e = which in turn woul suggest the following approximation to the lower quantile function: yu = u log + ν θ + θ a [ log u log u, ν θ + θ + ν θ + θ ν +θ +θ ν ν ] ν as u 0 + by using 7. We can then obtain the behaviour of the quantile function as = hu u log + ν θ + θ a [ log u + ν θ + θ ν +θ +θ ν ν ] ν + O log u 6
17 = log u ν θ + θ ν ν + θ + θ log ν + ν θ + θ a ν ν ν +θ +θ + O log u as u 0 +, by using 8 an 9. A Skew-Slash istribution on the other han was first propose in multivariate form in Arslan 008 an further stuie in Ling an Peng 05. A ranom variable X is sai to have a Skew-Slash istribution if X is efine by a Normal variance-mean mixture as X Y Nµ + θ/y, σ Y, where µ, θ R an σ > 0. The istribution of Y is Betaλ, ; that is, its pf is: fy = λy λ, 0 < y < ; = 0 otherwise. Here we only consier the case of θ > 0 so that the tail behaviour is similar to that of the Variance-Gamma in 8. When µ = 0, σ = an θ > 0, from Ling an Peng 05, Lemma. an equation, we have gx = λθλ x λ+ e θ x + O, x as x. By comparing this with 6, we have a = λθλ, b = λ +, c = θ, an = e =, which in turn suggests the following approximation to the lower quantile function: λ + yu = θ log θ λ+ u λ+ λθ λ / log θ λ+ u λ+ log u θ, as u 0 + by using 7. Ling an Peng 05, Lemma 3., equation 3 use the fact that hu = yu + o, u 0 + with yu = log u/θ. We can obtain the behaviour of the 7
18 quantile function as λ + hu = θ = log u θ log θ λ+ u λ+ λθ λ / log θ + λ + θ λ+ u λ+ as u 0 + by using 8 an 9 once again. log λθ λ θ + O log u + O log u. Acknowlegements The authors thank the referee for a careful reaing of the original version, an a list of suggestions which have resulte in a much improve paper. References Abramowitz, M. an Stegun, I.E. e., Hanbook of Mathematical Functions With Formulas, Graphs an Mathematical Tables. National Bureau of Stanars, Washington. Arslan, O. 008 An alternative multivariate skew-slash istribution. Statistics an Probability Letters 78: Azzalini, A. 985 A class of istributions which inclues the normal ones. Scan. J. Statist. :7 78. Azzalini, A. an Capitanio, A. 04 The Skew-Normal an Relate Families. IMS Monographs. Cambrige University Press, Cambrige. Beasley, J. D. an Springer S. G The percentage points of the Normal Distribution, Applie Statistics 6, 8. Capitanio, A. 00 On the approximation of the tail probability of the scalar skewnormal istribution. Metron 68:
19 De Bruijn, N.G. 96 Asymptotic Methos in Analysis. n E. North-Hollan Publishing Co., Amsteram. Embrechts, P., McNeil, A. an Straumann, D. 00. Correlation an epenency in risk management: properties an pitfalls. In, Dempster, M. an Moffatt, H., es. Risk Management: Value at Risk an Beyon. Cambrige University Press, Cambrige, pp Maan, D.B., Carr, P an Chang, E.C. 998 The variance gamma process an option pricing. European Finance : Feller, W. 968 An Introuction to Probability Theory an Its Applications. Wiley, New York. Fung T. an Seneta, E. 0 The bivariate normal copula function is regularly varying. Statistics an Probability Letters 8: Fung, T. an Seneta, E. 06 Tail asymptotics for the bivariate skew normal. Journal of Multivariate analysis 44:9 38. Lefor, A.W. an Tawn J.A. 997 Moelling epenence with joint tail regions. Journal of the Royal Statistical Society: Series B 59: Ling, C. an Peng, Z. 05 Tail epenence for two skew slash istributions. Statistics an Its Interface 8: Schoutens, W.003 Lévy Processes in Finance: Pricing Financial Derivatives. Wiley, New York. Seneta, E. 004 Fitting the variance-gamma moel to financial ata. Journal of Applie Probability 4A Heye Festschrift: Seneta, E. 976 Regularly Varying Functions. Lecture Notes in Mathematics 508. Springer, Berlin. 9
20 Soranzo, A. an Epure, E. 04. Very simply explicitly invertible approximations of normal cumulative an normal quantile function. Applie Mathematical Sciences 8, Tjetjep, A. an Seneta E. 006 Skewe normal variance-mean moels for asset pricing an the metho of moments. International Statistical Review 74:09 6. Voutier, P.M. 00. A new approximation to the normal istribution quantile function. arxiv: [stat.co] 0
Modeling of Dependence Structures in Risk Management and Solvency
Moeling of Depenence Structures in Risk Management an Solvency University of California, Santa Barbara 0. August 007 Doreen Straßburger Structure. Risk Measurement uner Solvency II. Copulas 3. Depenent
More informationarxiv: v4 [math.pr] 27 Jul 2016
The Asymptotic Distribution of the Determinant of a Ranom Correlation Matrix arxiv:309768v4 mathpr] 7 Jul 06 AM Hanea a, & GF Nane b a Centre of xcellence for Biosecurity Risk Analysis, University of Melbourne,
More informationA simple tranformation of copulas
A simple tranformation of copulas V. Durrleman, A. Nikeghbali & T. Roncalli Groupe e Recherche Opérationnelle Créit Lyonnais France July 31, 2000 Abstract We stuy how copulas properties are moifie after
More informationTopic 7: Convergence of Random Variables
Topic 7: Convergence of Ranom Variables Course 003, 2016 Page 0 The Inference Problem So far, our starting point has been a given probability space (S, F, P). We now look at how to generate information
More informationComputing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions
Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5
More informationA Modification of the Jarque-Bera Test. for Normality
Int. J. Contemp. Math. Sciences, Vol. 8, 01, no. 17, 84-85 HIKARI Lt, www.m-hikari.com http://x.oi.org/10.1988/ijcms.01.9106 A Moification of the Jarque-Bera Test for Normality Moawa El-Fallah Ab El-Salam
More informationMODELLING DEPENDENCE IN INSURANCE CLAIMS PROCESSES WITH LÉVY COPULAS ABSTRACT KEYWORDS
MODELLING DEPENDENCE IN INSURANCE CLAIMS PROCESSES WITH LÉVY COPULAS BY BENJAMIN AVANZI, LUKE C. CASSAR AND BERNARD WONG ABSTRACT In this paper we investigate the potential of Lévy copulas as a tool for
More informationLeast-Squares Regression on Sparse Spaces
Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction
More informationSome Examples. Uniform motion. Poisson processes on the real line
Some Examples Our immeiate goal is to see some examples of Lévy processes, an/or infinitely-ivisible laws on. Uniform motion Choose an fix a nonranom an efine X := for all (1) Then, {X } is a [nonranom]
More information1. Aufgabenblatt zur Vorlesung Probability Theory
24.10.17 1. Aufgabenblatt zur Vorlesung By (Ω, A, P ) we always enote the unerlying probability space, unless state otherwise. 1. Let r > 0, an efine f(x) = 1 [0, [ (x) exp( r x), x R. a) Show that p f
More informationTail dependence for two skew slash distributions
1 Tail dependence for two skew slash distributions Chengxiu Ling 1, Zuoxiang Peng 1 Faculty of Business and Economics, University of Lausanne, Extranef, UNIL-Dorigny, 115 Lausanne, Switzerland School of
More informationTime-of-Arrival Estimation in Non-Line-Of-Sight Environments
2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor
More informationMark J. Machina CARDINAL PROPERTIES OF "LOCAL UTILITY FUNCTIONS"
Mark J. Machina CARDINAL PROPERTIES OF "LOCAL UTILITY FUNCTIONS" This paper outlines the carinal properties of "local utility functions" of the type use by Allen [1985], Chew [1983], Chew an MacCrimmon
More informationBivariate distributions characterized by one family of conditionals and conditional percentile or mode functions
Journal of Multivariate Analysis 99 2008) 1383 1392 www.elsevier.com/locate/jmva Bivariate istributions characterize by one family of conitionals an conitional ercentile or moe functions Barry C. Arnol
More informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
More informationSeparation of Variables
Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical
More informationALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS
ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS MARK SCHACHNER Abstract. When consiere as an algebraic space, the set of arithmetic functions equippe with the operations of pointwise aition an
More informationConvergence of Random Walks
Chapter 16 Convergence of Ranom Walks This lecture examines the convergence of ranom walks to the Wiener process. This is very important both physically an statistically, an illustrates the utility of
More informationA simple model for the small-strain behaviour of soils
A simple moel for the small-strain behaviour of soils José Jorge Naer Department of Structural an Geotechnical ngineering, Polytechnic School, University of São Paulo 05508-900, São Paulo, Brazil, e-mail:
More informationConstruction of the Electronic Radial Wave Functions and Probability Distributions of Hydrogen-like Systems
Construction of the Electronic Raial Wave Functions an Probability Distributions of Hyrogen-like Systems Thomas S. Kuntzleman, Department of Chemistry Spring Arbor University, Spring Arbor MI 498 tkuntzle@arbor.eu
More informationChapter 6: Energy-Momentum Tensors
49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.
More informationarxiv:math/ v1 [math.pr] 19 Apr 2001
Conitional Expectation as Quantile Derivative arxiv:math/00490v math.pr 9 Apr 200 Dirk Tasche November 3, 2000 Abstract For a linear combination u j X j of ranom variables, we are intereste in the partial
More informationLogarithmic spurious regressions
Logarithmic spurious regressions Robert M. e Jong Michigan State University February 5, 22 Abstract Spurious regressions, i.e. regressions in which an integrate process is regresse on another integrate
More informationPROBABILITY THEORY EXPONENTIAL FUNCTIONAL OF LÉVY PROCESSES: GENERALIZED WEIERSTRASS PRODUCTS AND WIENER-HOPF FACTORIZATION P. PATIE AND M.
PROBABILITY THEORY EXPONENTIAL FUNCTIONAL OF LÉVY PROCESSES: GENERALIZED WEIERSTRASS PRODUCTS AND WIENER-HOPF FACTORIZATION P. PATIE AND M. SAVOV Abstract. In this note, we state a representation of the
More informationThe derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)
Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)
More informationNOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,
NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which
More informationA LIMIT THEOREM FOR RANDOM FIELDS WITH A SINGULARITY IN THE SPECTRUM
Teor Imov r. ta Matem. Statist. Theor. Probability an Math. Statist. Vip. 81, 1 No. 81, 1, Pages 147 158 S 94-911)816- Article electronically publishe on January, 11 UDC 519.1 A LIMIT THEOREM FOR RANDOM
More informationLie symmetry and Mei conservation law of continuum system
Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive
More informationSurvey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013
Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing
More informationThis module is part of the. Memobust Handbook. on Methodology of Modern Business Statistics
This moule is part of the Memobust Hanbook on Methoology of Moern Business Statistics 26 March 2014 Metho: Balance Sampling for Multi-Way Stratification Contents General section... 3 1. Summary... 3 2.
More informationAgmon Kolmogorov Inequalities on l 2 (Z d )
Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,
More informationarxiv: v1 [math.pr] 4 Feb 2016
Mittag-Leffler Lévy Processes Arun Kumar an N. S. Upahye *Inian Statistical Institute, Chennai Center, Taramani, Chennai-636, Inia an **Department of Mathematics, Inian Institute of Technology Maras, Chennai
More informationResearch Article U-Statistic for Multivariate Stable Distributions
Hinawi Probability an Statistics Volume 2017, Article ID 3483827, 12 pages https://oiorg/101155/2017/3483827 Research Article U-Statistic for Multivariate Stable Distributions Mahi Teimouri, Saei Rezakhah,
More informationSturm-Liouville Theory
LECTURE 5 Sturm-Liouville Theory In the three preceing lectures I emonstrate the utility of Fourier series in solving PDE/BVPs. As we ll now see, Fourier series are just the tip of the iceberg of the theory
More informationLecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012
CS-6 Theory Gems November 8, 0 Lecture Lecturer: Alesaner Mąry Scribes: Alhussein Fawzi, Dorina Thanou Introuction Toay, we will briefly iscuss an important technique in probability theory measure concentration
More informationHyperbolic Moment Equations Using Quadrature-Based Projection Methods
Hyperbolic Moment Equations Using Quarature-Base Projection Methos J. Koellermeier an M. Torrilhon Department of Mathematics, RWTH Aachen University, Aachen, Germany Abstract. Kinetic equations like the
More informationarxiv:hep-th/ v1 3 Feb 1993
NBI-HE-9-89 PAR LPTHE 9-49 FTUAM 9-44 November 99 Matrix moel calculations beyon the spherical limit arxiv:hep-th/93004v 3 Feb 993 J. Ambjørn The Niels Bohr Institute Blegamsvej 7, DK-00 Copenhagen Ø,
More informationTHE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE
Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek
More informationLectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs
Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent
More informationMath Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors
Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+
More informationunder the null hypothesis, the sign test (with continuity correction) rejects H 0 when α n + n 2 2.
Assignment 13 Exercise 8.4 For the hypotheses consiere in Examples 8.12 an 8.13, the sign test is base on the statistic N + = #{i : Z i > 0}. Since 2 n(n + /n 1) N(0, 1) 2 uner the null hypothesis, the
More informationensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y
Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay
More informationA note on asymptotic formulae for one-dimensional network flow problems Carlos F. Daganzo and Karen R. Smilowitz
A note on asymptotic formulae for one-imensional network flow problems Carlos F. Daganzo an Karen R. Smilowitz (to appear in Annals of Operations Research) Abstract This note evelops asymptotic formulae
More informationImproving Estimation Accuracy in Nonrandomized Response Questioning Methods by Multiple Answers
International Journal of Statistics an Probability; Vol 6, No 5; September 207 ISSN 927-7032 E-ISSN 927-7040 Publishe by Canaian Center of Science an Eucation Improving Estimation Accuracy in Nonranomize
More informationIntroduction. A Dirichlet Form approach to MCMC Optimal Scaling. MCMC idea
Introuction A Dirichlet Form approach to MCMC Optimal Scaling Markov chain Monte Carlo (MCMC quotes: Metropolis et al. (1953, running coe on the Los Alamos MANIAC: a feasible approach to statistical mechanics
More informationSpring 2016 Network Science
Spring 206 Network Science Sample Problems for Quiz I Problem [The Application of An one-imensional Poisson Process] Suppose that the number of typographical errors in a new text is Poisson istribute with
More informationarxiv: v1 [math.st] 20 Jun 2017
On the joint asymptotic istribution of the restricte arxiv:1706.06632v1 [math.st] 20 Jun 2017 estimators in multivariate regression moel Sévérien Nkurunziza an Youzhi Yu Abstract The main Theorem of Jain
More informationDEGREE DISTRIBUTION OF SHORTEST PATH TREES AND BIAS OF NETWORK SAMPLING ALGORITHMS
DEGREE DISTRIBUTION OF SHORTEST PATH TREES AND BIAS OF NETWORK SAMPLING ALGORITHMS SHANKAR BHAMIDI 1, JESSE GOODMAN 2, REMCO VAN DER HOFSTAD 3, AND JÚLIA KOMJÁTHY3 Abstract. In this article, we explicitly
More informationThe Entropy of Random Finite Sets
The Entropy of Ranom Finite Sets Mohamma Rezaeian an Ba-Ngu Vo Department of Electrical an Electronic Engineering, University of Melbourne, Victoria, 300, Australia rezaeian, b.vo@ee.unimelb.eu.au Abstract
More informationRole of parameters in the stochastic dynamics of a stick-slip oscillator
Proceeing Series of the Brazilian Society of Applie an Computational Mathematics, v. 6, n. 1, 218. Trabalho apresentao no XXXVII CNMAC, S.J. os Campos - SP, 217. Proceeing Series of the Brazilian Society
More informationModelling and simulation of dependence structures in nonlife insurance with Bernstein copulas
Moelling an simulation of epenence structures in nonlife insurance with Bernstein copulas Prof. Dr. Dietmar Pfeifer Dept. of Mathematics, University of Olenburg an AON Benfiel, Hamburg Dr. Doreen Straßburger
More informationA New Discrete Family of Reduced Modified Weibull Distribution
International Journal of Statistical Distributions an Applications 2017; 3(3): 25-31 http://www.sciencepublishinggroup.com/j/ijsa oi: 10.11648/j.ijs.20170303.11 ISSN: 2472-3487 (Print); ISSN: 2472-3509
More information6 General properties of an autonomous system of two first order ODE
6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationJointly continuous distributions and the multivariate Normal
Jointly continuous istributions an the multivariate Normal Márton alázs an álint Tóth October 3, 04 This little write-up is part of important founations of probability that were left out of the unit Probability
More informationEnergy behaviour of the Boris method for charged-particle dynamics
Version of 25 April 218 Energy behaviour of the Boris metho for charge-particle ynamics Ernst Hairer 1, Christian Lubich 2 Abstract The Boris algorithm is a wiely use numerical integrator for the motion
More informationLower bounds on Locality Sensitive Hashing
Lower bouns on Locality Sensitive Hashing Rajeev Motwani Assaf Naor Rina Panigrahy Abstract Given a metric space (X, X ), c 1, r > 0, an p, q [0, 1], a istribution over mappings H : X N is calle a (r,
More informationEquilibrium in Queues Under Unknown Service Times and Service Value
University of Pennsylvania ScholarlyCommons Finance Papers Wharton Faculty Research 1-2014 Equilibrium in Queues Uner Unknown Service Times an Service Value Laurens Debo Senthil K. Veeraraghavan University
More informationThe Uniqueness of the Overall Assurance Interval for Epsilon in DEA Models by the Direction Method
Available online at http://nrm.srbiau.ac.ir Vol., No., Summer 5 Journal of New Researches in Mathematics Science an Research Branch (IAU) he Uniqueness of the Overall Assurance Interval for Epsilon in
More informationTail dependence in bivariate skew-normal and skew-t distributions
Tail dependence in bivariate skew-normal and skew-t distributions Paola Bortot Department of Statistical Sciences - University of Bologna paola.bortot@unibo.it Abstract: Quantifying dependence between
More informationMonotonicity for excited random walk in high dimensions
Monotonicity for excite ranom walk in high imensions Remco van er Hofsta Mark Holmes March, 2009 Abstract We prove that the rift θ, β) for excite ranom walk in imension is monotone in the excitement parameter
More informationSpurious Significance of Treatment Effects in Overfitted Fixed Effect Models Albrecht Ritschl 1 LSE and CEPR. March 2009
Spurious Significance of reatment Effects in Overfitte Fixe Effect Moels Albrecht Ritschl LSE an CEPR March 2009 Introuction Evaluating subsample means across groups an time perios is common in panel stuies
More informationELEC3114 Control Systems 1
ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.
More informationarxiv: v1 [physics.class-ph] 20 Dec 2017
arxiv:1712.07328v1 [physics.class-ph] 20 Dec 2017 Demystifying the constancy of the Ermakov-Lewis invariant for a time epenent oscillator T. Pamanabhan IUCAA, Post Bag 4, Ganeshkhin, Pune - 411 007, Inia.
More informationLEGENDRE TYPE FORMULA FOR PRIMES GENERATED BY QUADRATIC POLYNOMIALS
Ann. Sci. Math. Québec 33 (2009), no 2, 115 123 LEGENDRE TYPE FORMULA FOR PRIMES GENERATED BY QUADRATIC POLYNOMIALS TAKASHI AGOH Deicate to Paulo Ribenboim on the occasion of his 80th birthay. RÉSUMÉ.
More informationRandom products and product auto-regression
Filomat 7:7 (013), 1197 103 DOI 10.98/FIL1307197B Publishe by Faculty of Sciences an Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Ranom proucts an prouct auto-regression
More informationThermal conductivity of graded composites: Numerical simulations and an effective medium approximation
JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University
More informationNotes on Lie Groups, Lie algebras, and the Exponentiation Map Mitchell Faulk
Notes on Lie Groups, Lie algebras, an the Exponentiation Map Mitchell Faulk 1. Preliminaries. In these notes, we concern ourselves with special objects calle matrix Lie groups an their corresponing Lie
More informationd dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1
Lecture 5 Some ifferentiation rules Trigonometric functions (Relevant section from Stewart, Seventh Eition: Section 3.3) You all know that sin = cos cos = sin. () But have you ever seen a erivation of
More informationSYNCHRONOUS SEQUENTIAL CIRCUITS
CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents
More informationHIGH QUANTILE ESTIMATION AND THE PORT METHODOLOGY
REVSTAT Statistical Journal Volume 7, Number 3, November 2009, 245 264 HIGH QUANTILE ESTIMATION AND THE PORT METHODOLOGY Authors: Lígia Henriques-Rorigues Instituto Politécnico e Tomar an CEAUL, Portugal
More informationOptimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations
Optimize Schwarz Methos with the Yin-Yang Gri for Shallow Water Equations Abessama Qaouri Recherche en prévision numérique, Atmospheric Science an Technology Directorate, Environment Canaa, Dorval, Québec,
More informationTHE EFFICIENCIES OF THE SPATIAL MEDIAN AND SPATIAL SIGN COVARIANCE MATRIX FOR ELLIPTICALLY SYMMETRIC DISTRIBUTIONS
THE EFFICIENCIES OF THE SPATIAL MEDIAN AND SPATIAL SIGN COVARIANCE MATRIX FOR ELLIPTICALLY SYMMETRIC DISTRIBUTIONS BY ANDREW F. MAGYAR A issertation submitte to the Grauate School New Brunswick Rutgers,
More informationLower Bounds for the Smoothed Number of Pareto optimal Solutions
Lower Bouns for the Smoothe Number of Pareto optimal Solutions Tobias Brunsch an Heiko Röglin Department of Computer Science, University of Bonn, Germany brunsch@cs.uni-bonn.e, heiko@roeglin.org Abstract.
More informationFLUCTUATIONS IN THE NUMBER OF POINTS ON SMOOTH PLANE CURVES OVER FINITE FIELDS. 1. Introduction
FLUCTUATIONS IN THE NUMBER OF POINTS ON SMOOTH PLANE CURVES OVER FINITE FIELDS ALINA BUCUR, CHANTAL DAVID, BROOKE FEIGON, MATILDE LALÍN 1 Introuction In this note, we stuy the fluctuations in the number
More informationTransforms. Convergence of probability generating functions. Convergence of characteristic functions functions
Transforms For non-negative integer value ranom variables, let the probability generating function g X : [0, 1] [0, 1] be efine by g X (t) = E(t X ). The moment generating function ψ X (t) = E(e tx ) is
More informationCalculus of Variations
16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t
More information. ISSN (print), (online) International Journal of Nonlinear Science Vol.6(2008) No.3,pp
. ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.6(8) No.3,pp.195-1 A Bouneness Criterion for Fourth Orer Nonlinear Orinary Differential Equations with Delay
More informationFinancial Time Series Analysis Week 5
Financial Time Series Analysis Week 5 25 Estimation in AR moels Central Limit Theorem for µ in AR() Moel Recall : If X N(µ, σ 2 ), normal istribute ranom variable with mean µ an variance σ 2, then X µ
More informationBinary Discrimination Methods for High Dimensional Data with a. Geometric Representation
Binary Discrimination Methos for High Dimensional Data with a Geometric Representation Ay Bolivar-Cime, Luis Miguel Corova-Roriguez Universia Juárez Autónoma e Tabasco, División Acaémica e Ciencias Básicas
More informationDiscrete Mathematics
Discrete Mathematics 309 (009) 86 869 Contents lists available at ScienceDirect Discrete Mathematics journal homepage: wwwelseviercom/locate/isc Profile vectors in the lattice of subspaces Dániel Gerbner
More informationStable Polynomials over Finite Fields
Rev. Mat. Iberoam., 1 14 c European Mathematical Society Stable Polynomials over Finite Fiels Domingo Gómez-Pérez, Alejanro P. Nicolás, Alina Ostafe an Daniel Saornil Abstract. We use the theory of resultants
More informationA Review of Multiple Try MCMC algorithms for Signal Processing
A Review of Multiple Try MCMC algorithms for Signal Processing Luca Martino Image Processing Lab., Universitat e València (Spain) Universia Carlos III e Mari, Leganes (Spain) Abstract Many applications
More informationThe Role of Models in Model-Assisted and Model- Dependent Estimation for Domains and Small Areas
The Role of Moels in Moel-Assiste an Moel- Depenent Estimation for Domains an Small Areas Risto Lehtonen University of Helsini Mio Myrsylä University of Pennsylvania Carl-Eri Särnal University of Montreal
More informationThe Press-Schechter mass function
The Press-Schechter mass function To state the obvious: It is important to relate our theories to what we can observe. We have looke at linear perturbation theory, an we have consiere a simple moel for
More informationOn the First Passage Time and Leapover Properties of Lévy Motions
On the First Passage Time an Leapover Properties of Lévy Motions T. Koren a, A.V. Chechkin b, an J. Klafter a a School of Chemistry, Tel Aviv University, Tel Aviv 69978 Israel b Institute for Theoretical
More informationSurvey-weighted Unit-Level Small Area Estimation
Survey-weighte Unit-Level Small Area Estimation Jan Pablo Burgar an Patricia Dörr Abstract For evience-base regional policy making, geographically ifferentiate estimates of socio-economic inicators are
More informationTMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments
Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary
More informationChaos, Solitons and Fractals Nonlinear Science, and Nonequilibrium and Complex Phenomena
Chaos, Solitons an Fractals (7 64 73 Contents lists available at ScienceDirect Chaos, Solitons an Fractals onlinear Science, an onequilibrium an Complex Phenomena journal homepage: www.elsevier.com/locate/chaos
More informationAdvanced Partial Differential Equations with Applications
MIT OpenCourseWare http://ocw.mit.eu 18.306 Avance Partial Differential Equations with Applications Fall 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.eu/terms.
More informationTail dependence coefficient of generalized hyperbolic distribution
Tail dependence coefficient of generalized hyperbolic distribution Mohalilou Aleiyouka Laboratoire de mathématiques appliquées du Havre Université du Havre Normandie Le Havre France mouhaliloune@gmail.com
More informationEstimation of the Maximum Domination Value in Multi-Dimensional Data Sets
Proceeings of the 4th East-European Conference on Avances in Databases an Information Systems ADBIS) 200 Estimation of the Maximum Domination Value in Multi-Dimensional Data Sets Eleftherios Tiakas, Apostolos.
More informationMonotonicity of facet numbers of random convex hulls
Monotonicity of facet numbers of ranom convex hulls Gilles Bonnet, Julian Grote, Daniel Temesvari, Christoph Thäle, Nicola Turchi an Florian Wespi arxiv:173.31v1 [math.mg] 7 Mar 17 Abstract Let X 1,...,
More informationUniqueness of limit cycles of the predator prey system with Beddington DeAngelis functional response
J. Math. Anal. Appl. 290 2004 113 122 www.elsevier.com/locate/jmaa Uniqueness of limit cycles of the preator prey system with Beington DeAngelis functional response Tzy-Wei Hwang 1 Department of Mathematics,
More information19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control
19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior
More informationTopic Modeling: Beyond Bag-of-Words
Hanna M. Wallach Cavenish Laboratory, University of Cambrige, Cambrige CB3 0HE, UK hmw26@cam.ac.u Abstract Some moels of textual corpora employ text generation methos involving n-gram statistics, while
More informationTest of Hypotheses in a Time Trend Panel Data Model with Serially Correlated Error Component Disturbances
HE UNIVERSIY OF EXAS A SAN ANONIO, COLLEGE OF BUSINESS Working Paper SERIES Date September 25, 205 WP # 000ECO-66-205 est of Hypotheses in a ime ren Panel Data Moel with Serially Correlate Error Component
More informationarxiv: v1 [stat.co] 23 Jun 2012
Noname manuscript No. (will be inserte by the eitor) Moments Calculation For the Doubly Truncate Multivariate Normal Density Manjunath B G Stefan Wilhelm arxiv:1206.5387v1 [stat.co] 23 Jun 2012 This version:
More informationFluid Models with Jumps
Flui Moels with Jumps Elena I. Tzenova, Ivo J.B.F. Aan, Viyahar G. Kulkarni August 2, 2004 Abstract In this paper we stuy a general stochastic flui moel with a single infinite capacity buffer, where the
More informationChapter 4. Electrostatics of Macroscopic Media
Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1
More information