Generalized Post-Widder inversion formula with application to statistics
|
|
- Gabriel Griffin
- 5 years ago
- Views:
Transcription
1 arxiv:5.9298v [math.st] 3 ov 25 Generalized Post-Widder inversion formula with application to statistics Denis Belomestny, Hilmar Mai 2, John Schoenmakers 3 August 24, 28 Abstract In this work we derive an inversion formula for the Laplace transform of a density observed on a curve in the comple domain, which generalizes the well known Post-Widder formula. We establish convergence of our inversion method and derive the corresponding convergence rates for the case of a Laplace transform of a smooth density. As an application we consider the problem of statistical inference for variance-mean miture models. We construct a nonparametric estimator for the miing density based on the generalized Post-Widder formula, derive bounds for its root mean square error and give a brief numerical eample. Keywords: Laplace transform, inversion formula, Post-Widder formula, variancemean mitures, density estimation. Introduction Let p be a probability density on R +, then the integral L(z) := e z p() d, Re z >, () eists and is called the Laplace transform of p. The Laplace transform is a popular tool for solving differential equations and convolution integral equations. Its inversion is of importance in many problems from e.g. physics, engineering and finance (c.f. [2] and [6] for various eamples). In general, the compleity of the inversion problem for L depends on the information available about the Laplace transform. If the Laplace transform is eplicitly given on its Duisburg-Essen University, Essen and IITP RAS, Moscow, Russia denis.belomestny@uni-due.de 2 CREST ESAE ParisTech, Paris, hilmar.mai@ensae.fr, (corresponding author) 3 Weierstrass Institute, Berlin, schoenma@wias-berlin.de
2 half-plane of convergence, the density p can be reconstructed using the so-called Bromvich contour integral (see, e.g. [7]) p() = 2πi c+i c i e z L(z) dz, >. In the real case, i.e. in the situation where the Laplace transform of p is known on the real ais only, the inversion of L is a well-known ill-posed problem (see for eample [4], [5] and references therein). One popular solution for this case is given by the well known Post-Widder formula which reads as follows (cf. [7]): ( ) p() = lim! ( ) + L () ( In some situations, the Laplace transform L can only be computed on some curve l in C, which is different from R + or {Re(z) = c} for some c >. In this paper we generalize the Post-Widder formula to the case of rather general curves l and derive the convergence rates of the resulting estimator. As an application of our results we consider the problem of estimating the miing density in a variance mean miture model (see e.g. [] and [3]). After constructing the estimator we derive bounds for its root mean square error (RMSE) and demonstrate its performance in a short numerical eample. An advantage of using the generalized Post-Widder formula here is that the resulting estimator can be evaluated without any numerical integration. The paper is organized as follows. In Section 2 we introduce the generalized Post- Widder inversion formula and discuss its convergence behavior. Section 3 is devoted to the statistical inference for variance mean mitures together with some numerical results. Finally, the proofs of our results are given in Section 5 to 7. 2 Generalized Post-Widder Laplace inversion In this section we will introduce a generalized Post-Widder inversion formula that etends the classical result by Post and Widder [7] to the situation when the Laplace transform of a continuous density on [, ] is given on a curve in the comple plane. Subsequently, we prove a convergence result and derive the rates of convergence for the resulting inverse Laplace transform. 2. Inversion formula and its kernel representation Let p be a continuous probability density on [, ) and let its Laplace transform L(z) be given on a curve: l := {z = y + ic(y) : y R + }, (2) such that c is piecewise smooth with c(y) = o(y) as y. In this setting the generalized Post-Widder formula can be described as follows. 2 ).
3 Definition 2. (Generalized Post-Widder formula). For any fied >, we introduce the generalized Post-Widder formula by p () := ( )! ( g ( )) + ( L () g ( )), (3) where L () denotes the th-derivative of the Laplace transform L and g(y) := y + ic(y). For fied >, we define the generalized Post-Widder kernel via K (t, ) := ( ( + ic )) + t e (+ic( ))t, t >. (4)! Our first result deals with the convergence of p to p as. Such a convergence follows from the properties of the generalized Post-Widder kernel K and a representation formula for p in terms of K and p. The latter representation is given by the following proposition. Proposition 2.2. It holds p () = p(t)k (t, ) dt. The following result states that K (t, ) converges to the delta function δ(t ) on (, ) for any fied >. Proposition 2.3. The following statements hold. (i) For r =,, 2, it holds t r ( + r/) ( + /) K (t, )dt = ( + i(/)c(/)) r (5) = + r ( c(/) r ir + O / + r c(/) ) 2,, >. / (6) Hence, in particular we have K (t, )dt = for all. (7) (ii) Let > be fied. For any δ (, ), there eists a natural number δ t r K (t, ) dt Ce (ln(+δ) δ)/8, > δ. { t δ, t } such that for r =,, 2, and some constant C not depending on and δ. 3
4 2.2 Convergence analysis By combining Proposition 2.2 and Proposition 2.3, the point-wise convergence of p to p follows for as stated in the following corollary. Corollary 2.4. For any fied and any continuous density p on [, ), we have lim p () = p(). (8) We may now sharpen the statement (8) under additional smoothness assumptions on the density p. In the following propositions we give eplicit convergence rates for p as. It turns out that the rates crucially depend on the growth behavior of the function c(y) as y. We henceforth assume that [ ] c 2 (y) γ := lim sup <. (9) y y The notation f(, ) = O (r(, )) for fied R and means in the sequel the usual O-notation where the actual order coefficient may depend on. We start with a local Lipschitz condition on p. Proposition 2.5. Let p be a locally Lipschitz continuous density on [, ) with Laplace transform () given on the curve (2). It then holds p () = p() + R (), where for and each >. R () = O ( /2 ) When the density p is differentiable, the rates of Proposition 2.5 can be improved as the following result shows. Proposition 2.6. Let p be a differentiable density on [, ) with Laplace transform () given on the curve (2). We then have for γ <, where Re [p ()] = p() + R (), R () = o( /2 ) for and each >. Further we have that p () = p() + O ( /2 ) for < γ <, and p () = p() + o( /2 ) for γ =. 4
5 We conclude this section by considering the Laplace inversion problem for a differentiable density p with locally Lipschitz derivative. It turns out that we can achieve the error term R of the order in this case. Proposition 2.7. Let p be a smooth density on [, ) such that its derivative p is locally Lipschitz, and its Laplace transform () is given on the curve (2). Then for γ <, where Re [p ()] = p() + R (), R () = O( ) for and each >. Moreover we have that p () = p() + O ( /2 ) for < γ <, and p () = p() + O ( ) for γ =. In the net section we discuss some applications of the generalized Post-Widder formula (3). 3 Application to statistical inference for variancemean mitures The problem of inverting a Laplace transform that is given on a curve l in the comple domain appears naturally in the contet of statistical inference for variance-mean miture models. In this section we apply our generalized Post-Widder Laplace inversion formula to estimate the miture density in a variance mean miture model. We start the construction of the estimator from the empirical characteristic function that can be written as the Laplace transform of the miture density evaluated on a certain curve in the comple plain. By inverting this Laplace transform we obtain a nonparametric estimator for the miing density p. Then we derive bounds for the RMSE and conclude by a numerical eample. 3. Variance-mean miture models A normal variance-mean miture model is defined as ( ) sµ q() := σ s ν σ p(s) ds, s where µ R, σ R +, ν is the density of a standard normal distribution and p is a miing density on R +. Variance-mean miture models play an important role in both theory and practice of statistics. In particular, such mitures appear as limit distributions in asymptotic theory for dependent random variables and they are useful for modeling data stemming from heavy-tailed and skewed distributions, see, e.g. [] and [3]. 5
6 As can be easily seen, the variance-mean miture distribution q coincides with the distribution of the random variable σ ξ X + µ ξ, where X is standard normal and ξ is a nonnegative random variable with the density p, which is independent of X. The class of variance-mean miture models is rather large. For eample, the class of the normal variance miture distributions (µ = ) can be described as follows: q is the density of a normal variance miture if and only if F[q]( u) is a completely monotone function in u. 3.2 Estimating the miing density Here we consider the problem of statistical inference for the miing density p based on a sample X,..., X n from the distribution q. The Fourier transform of the density q is given by Φ(u) := F[q](u) = e sψ(u) p(s) ds = L[p](ψ(u)) () with ψ(u) := iuµ + u 2 σ 2 /2 and from our data we can directly estimate the Fourier transform of q, e.g. by means of the so-called empirical Fourier transform: Φ n (u) := n n e iux k. () Then we end up with the problem of reconstructing the density p from its empirical Laplace transform observed on the curve k= l := {z = Re ψ(u) + i Im ψ(u) : u R + }, where we have Re[ψ(u)] = u 2 σ 2 /2 and Im[ψ(u)] = uµ. ote that l = {z = y + c(y) : y R + } with c(y) = µ 2y/σ. If σ than the function c is smooth and satisfies c(y) = o(y) as y. Moreover it holds [ ] c 2 (y) γ = lim sup = 2µ 2 /σ 2. (2) y y Hence if p is a differentiable density on [, ) such that p is locally Lipschitz, we can apply Proposition 2.7 to get the following asymptotic bound { O (/), µ =, p () p() = O (/ (3) ), µ, for p defined in (3) with g(y) = y iµ 2y/σ. Due to (), we have L(z) = Φ(ξ(z)), where ξ is the inverse of ψ on l. Without loss of generality we may assume that σ =, then ξ(z) = 2z µ 2 + iµ using the principal branch of the square root. So, for l we have that ξ (l) (z) = ( ) l (2 (l ))! ( ) 2z µ 2 2 l, 2 l (l )! 6
7 and by Faa di Bruno s formula it follows that for z l, L () (z) = = k,...,k, k +2k k = k +k k =k! k!... k!(!) k... (!) k Φ (k) (ξ(z)) ( ξ (l) (z) ) k l Φ (k) (ξ(z)) ( ) k ( 2z µ 2) 2 k F,k. (4) k= The coefficients F,k can be epressed as follows F,k := = k,...,k, k +2k k = k +k k =k! k!... k! k,...,k, k +2k ( k+)k k+ = k +k k k+ =k, = B,k (,..., (2 ( k))! 2 k ( k)! l= ( ) kl (2 (l ))! 2 l (l )!l!! k! k k+! ), k+ l= l= ( (2(l ))! ) kl 2 l (l )! l! where B,k stand for the partial Bell polynomials. In view of () and (4), we now introduce L () n (z) := Φ (k) n (ξ(z)) ( ) k ( 2z µ 2) 2 k F,k k= as an unbiased estimator for L () (z) at every z l. We so arrive at an empirical estimate for the miing density p: p n, () := ( )! = ( )! (g(/)) + L () n (g(/)) n (g(/)) + n j= e iξ(g(/))x j (ix j ) k ( ) k ( 2g(/) µ 2) 2 k F,k, (5) k= which obviously satisfies E [p n, ()] = p (). The coefficients F,k can be computed by evaluating the partial Bell polynomials B,k that are available in most computational algebra packages. Hence, we obtain an eplicit estimator for p that circumvents the use of numerical integration procedures as needed in other Laplace inversion techniques. 7
8 3.3 Convergence of the estimator Let us now analyze the variance of p n,. Theorem 3.. For some constant C >, depending on >, it holds that [ where β 2k := E X 2k]. Var [p n, ()] C n k β 2k, (6) k= Based on the estimate (6), we can derive upper bounds of the root mean square error (RMSE) for the density estimator p n,. Theorem 3.2. Fi some > and suppose that R () = p () p() = O ( ρ ). We then have the following bounds for the RMSE of p n,. (i) If β 2k A k k k for some A > and all natural k >, then ( ) RMSE(p n, ) = O ln ρ, n provided = where the constant C comes from the bound (6). ln n ln(ac) 2ρ ln ln n, (7) ln(ac) (ii) In the case β 2k A k k bk, k for some b >, it holds that which is achieved by choosing = RMSE(p n, ) lnρ ln n ln ρ n, (2ρ + ln n) (b ) ln ((2ρ + ln n) / (b )) 2ρ b. (8) Remark 3.3. Because of the inequality Var [p n, ()] Var [Re[p n, ()]], the results of Theorem 3.2 remain valid with p n, replaced by Re[p n, ]. 3.4 umerical eample Let the miing density p be the eponential density, i.e. p() = ep( ), then E [ X 2k] = 2k Γ( + k)γ(k + /2) A k k 2k π 8
9 p 4,2 p 4,3 ep(-) p 5 4,2 p 5 4,3 ep(-) Figure : Approimations p n, (real parts) for sample sizes n = (above), n = 5 (below) and different values of of the true eponential density p() = ep( ) in the normal mean variance miture model with parameters µ =. and σ =. 9
10 for some A >. Combining Proposition 2.7 with Theorem 3.2, we obtain RMSE(Re[p n, ]) ln ln(n) ln(n), n. In Figure one can see the result of numerical estimation of the underlying eponential density p() = ep( ) based on different numbers of terms in (5) and different sample sizes n. As can be observed, the estimation error increases as +. This effect can be eplained by noting that g(/) as + and so the variance increases for small (see the proof of Theorem 3.). 4 Acknowledgment The first author was supported by the Deutsche Forschungsgemeinschaft through the SFB 823 Statistical modelling of nonlinear dynamic processes and by a Russian Scientific Foundation grant (project 4-5-5). The second author received founding by the LABEX Louis Bachelier Finance et Croissance Durable in Paris. 5 Proofs In this section we gather the proofs of our results from Section Proof of Proposition 2.2 We start by proving the integral representation of p in terms of p and the generalized Post-Widder kernel K. We have by definition L(z) = differentiating -times results in L () (z) = e uz p(u)du, Re z > ( u) e uz p(u)du, yielding finally ( ( )) + p () = ( )! + ic ( u) e u( +ic( )) p(u) du = ( ( )) +! + ic (t) e t(+ic( )) p(t) dt ( ( )) + = + ic t e t(+ic( )) p(t) dt! = p(t)k (t, ) dt.
11 5.2 Proof of Proposition 2.3 (i): For r =,, 2,... we have t r K (t, )dt = =! ( + ic! ( + ic ( ( )) + t +r e t(+ic( )) dt )) r (+ic( )) z +r e z dz. ote that on the set { z : arc(z) = R e iθ, π/2 < θ < π/2 } it holds that z +r e z = R +r e R cos θ for R. Thus, by the Cauchy integral theorem, (+ic( )) z +r e z dz = t +r e t dt = ( + r)! from which (5) follows. Thus (5) holds for any integer integer r. The asymptotic epression (6) for r = and r = 2 can be seen from taking the logarithm of (5): ( + r/) ( + /) r ln ( + i(/)c(/)) r = ln ( + l/) r ln ( + i(/)c(/)) l= = r irc(/) / e ( + O + c(/) ) 2. / (ii): By Stirling s formula it follows that ( ( + i ) ( c )) + + K (t, ) = t e t e i(t)c( 2π ) ( + O(/)) and hence ( ( ) ) + = + ic 2π t e t e e i(t)c( ) ( + O(/)) ln K (t, ) = i(t)c ( ) ( + ( + ) ln + i ( )) c + 2 ln + ln t t + + O(/) 2π ( ( ) ) ( ) j j = ( + ) ic + j 2 ln 2π j= ( ) + ln t t + i(t)c + O(/) = 2 ln ( ) ) 2 (c 2π + O + O(/) ( ) ( ) + ln t t + + ic ( t) + ic. (9)
12 In particular, for t we have K (t, ) = ep 2 ln 2π + O(/) + ln } t {{ t + } ( ( ) ) 2 + O c < }{{} as. (2) Let us fi > and δ > arbitrarily. W.l.o.g. we may assume that δ <. Because c( ) for, there eist a number δ such that for any > δ and any t > with t δ, i.e. and so for r =,, 2, > δ, ( ( ) ) 2 ln t t + + O c < (ln t t + ), 2 [ K (t, ) ep 2 ln 2π + O(/) + ] (ln t t + ), 2 { t δ, t } t r K (t, ) dt C 2π { t δ, t } t r ( te t+) /2 dt. (2) ote that for t δ we have te t+ = e ln t t+ ln t t+ + e ln(+δ) δ 2 2, hence for r =,, 2, > δ, { t δ, t } t ( r te t+) /2 e (ln(+δ) δ)/4 e (ln(+δ) δ)/4 { t δ, t } { t δ, t } t r e (ln t t+)/4 dt t r e (ln t t+)/4 dt e (ln(+δ) δ)/4 t r+ e ( t+)/4 dt C 2 e (ln(+δ) δ)/4, where ln ( + δ) δ < and C 2 >. It net follows from (2) that for some constant C > t r K (t, ) dt Ce (ln(+δ) δ)/8 (22) for, > δ. { t δ, t } 5.3 Proof of Proposition 2.5 In order to derive the convergence rates in (8), we proceed with the following lemma. 2
13 Lemma 5.. For r =, 2,.. [ ( ( t r (r+)/2 Γ ((r + ) /2) K (t, ) dt C6 ep O c 2 r/2 Proof. Let us fi >. For δ > with < δ < we have by (2), +δ [ ( ( ) ) ] 2 t r K (t, ) dt C 2π ep O c where +δ (t ) r ep [ (ln t t + )] dt = δ Hence by (23) and Theorem 2.3-(ii) +δ t r K (t, ) dt C 6 (r+)/2 ) 2 )]. (t ) r ep [ (ln t t + )] dt, (23) u r ep [ (ln ( + u) u)] du 2 Γ ((r + ) /2) r/2 = 6(r+)/2 2 [ ( ( ep O c 2 Γ ((r + ) /2) (r+)/2. ) 2 )]. Similarly, δ [ ( ( ) )] t r K (t, ) dt = C 2π ep O c 2 2 δ ( t) r ep [ (ln t t + )] dt, where δ δ = ( t) r dt ep [ (ln t t + )] u r du ep [ (ln ( u) + u)] δ u r dt ep [ u 2 /6 ]. By applying Theorem 2.3-(ii) once again the statement of the lemma is proved. 3
14 Let us now fi >. By the Lipschitz assumption on p we thus have for fied < δ <, p () p() p(t) p() K (t, ) dt K t K (t, ) dt + K e (ln(+δ) δ)/8 due to assumption (9). t <δ [ )] K /2 ep O (c 2 ( )2 K 2 /2 5.4 Proof of Proposition 2.6 Let us fi >. By differentiability of p we may find for any ε > a δ ε with < δ ε < such that p(t) =: p() + (t ) (p () + E t ), with Et < ε for all t > with t < δ ε.due to Proposition 2.2 we then have p () p() = (p(t) p()) K (t, )dt = (p(t) p()) K (t, )dt t δ +p () (t )K (t, )dt + (t )Et K (t, )dt =: ( ) + ( ) 2 + ( ) 3. t <δ t <δ Since p is bounded we have by Theorem (2.3)-(ii) that ( ) D e (ln(+δε) δε)/8 (24) for some D > and > δ ε (cf. the proof of Theorem 2.3)-(ii). By Theorem (2.3) we have (t )K (t, )dt = (t )K (t, )dt (t )K (t, )dt t <δ ε t δ ε, t = ( ic(/) / + O + c(/) ) 2 + f,δε (25) / with f,δε D 2 e (ln(+δε) δε)/8 for > δ ε. et, by Lemma 5. and assumption (9) we have [ )] ε ( ) 3 ε (t K (t, ) dt D 3 ep O (c 2 ( )2 (26) t <δ /2 ε D 4 /2 4
15 From (24) (26) we gather that ( ) p () = p() + p () ic(/) / ow, since + O ( + c(/) / )2 + o( /2 ). (27) c(/) / = O ( /2), if < γ <, and c(/) / = o( /2 ) if γ =, (28) the statements follow by taking the real part of (27). 5.5 Proof of Proposition 2.7 Let us fi >. By the Lipschitz assumption on p we may find a δ with < δ < such that p(t) =: p() + (t )p () + 2 (t )2 R t, with Rt < R for some constant R > and for all t > with t < δ. Due to Proposition 2.2 we thus have p () p() = (p(t) p()) K (t, )dt = (p(t) p()) K (t, )dt + t <δ ( (t )p () + 2 (t )2 R t t δ ) K (t, )dt =: ( ) + ( ) 2. (29) Since p is bounded we have by Theorem 2.3-(ii) again that ( ) D e (ln(+δ) δ)/8 for some D > and > δ (cf. the proof of Theorem (2.3)-(ii)). ow let us consider ( ) 2. From Theorem 2.3 it follows similar to the proof of Proposition 2.6 that t <δ (t )K (t, )dt = ic(/) / + O ( + c(/) ) 2 + f,δ (3) / with f,δ D 2 e (ln(+δ) δ)/8, and by Lemma 5. we have that (t ) 2 Rt K (t, )dt R t 2 K (t, ) dt t <δ D [ )] 3 ep O (c 2 ( )2 (3) We thus get by (29), (3), (3), and assumption (9), ( ) ( p () = p() + p () ic(/) + O / + c(/) ) 2 + O( ) / from which the statements follow by taking the real part and taking (28) into account. 5
16 6 Proof of Proposition 3. For a generic constant C >, depending on and changing in this proof from line to line, we may write Var [p n, ()] C 2+2 n E (ix 2 ) k ( ) k ( 2 2g(/) µ 2) 2 k F,k k= C n E (ix ) k ( ) k ( 2 2g(/) µ 2) 2 k F,k C n C n k= k= k= [ E X 2k ] 2g(/) µ 2 k 2 F,k 2 [ k 2 F,kE 2 X 2k]. It is not difficult to see that for l =,..., k +, (2 (l ))! 2 l (l )! Cl l! and so from the definition of the Bell polynomials it follows that F,k C k B,k (!,..., ( k + )!) ( ) ( ) = C k ( k)! k k C k. 7 Proof of Proposition 3.2 (i): Without loss of generality we may assume that A >. Since for k, k k k, we get from (6), Var [p n, ()] C A k A (AC). (32) n A n k= By substituting according to (7) into (32) we obtain Var [p n, ()] A (AC) A n = A A ln 2ρ n, while for the squared bias we have ( ) R 2 () = O ( 2ρ ) = O ln 2ρ, n 6
17 hence (i) follows. (ii): In this case (6) yields Var [p n, ()] C n k= A k k (b )k C n (b ) (33) for another constant C >. From (8) it follows straightforwardly that ( ln (b )+2ρ 2ρ + ln n = ln ((2ρ + ln n) / (b )) ln (2ρ + ln n) (b ) ln ((2ρ + ln n) / (b )) 2ρ ) b 2ρ + ln n = = ln ln ((2ρ + ln n) / (b )) [ (2ρ + ln n) / (b ) ln ((2ρ + ln n) / (b )) 2ρ + ln n ln ((2ρ + ln n) / (b )) ln (2ρ + ln n) / (b ) ( )] ln ((2ρ + ln n) / (b )) 2ρ 2ρ + ln n ( ) ln ln n ln ((2ρ + ln n) / (b )) 2ρ + O ln n ( ) ln ln n + O. ln n ln ln ((2ρ + ln n) / (b )) = ln n (2ρ + ln n) ln ((2ρ + ln n) / (b )) On the other hand, from (8), So, for n, hence We thus obtain by (33), while which gives (ii). ln C (2ρ + ln n) ln C (b ) ln ((2ρ + ln n) / (b )) ln ( C (b )+2ρ) ( ) ln ln n ln n + O ln n ln ln ((2ρ + ln n) / (b )) (2ρ + ln n) ln ((2ρ + ln n) / (b )) }{{} + ln C (b ) ln ln ((2ρ + ln n) / (b )), }{{} C (b )+2ρ 2n Var(p n, ()) 2 2ρ ln2ρ ln n ln 2ρ n ( ln 2ρ ) R 2 () = O ( 2ρ ln n ) = O ln 2ρ n 7
18 References [] O Barndorff-ielsen, John Kent, and Michael Sørensen. ormal variance-mean mitures and z distributions. International Statistical Review, 5(2):45 59, 982., 3. [2] R.E. Bellman, R.E. Kalaba, and J.A. Lockett. umerical inversion of the Laplace transform: applications to biology, economics, engineering, and physics. Modern analytic and computational methods in science and mathematics. American Elsevier Pub. Co., 966. [3] icholas H Bingham and Rüdiger Kiesel. Semi-parametric modelling in finance: theoretical foundations. Quantitative Finance, 2(4):24 25, 22., 3. [4] V V Kryzhniy. Direct regularization of the inversion of real-valued Laplace transforms. Inverse Problems, 9(3):573, 23. [5] V V Kryzhniy. umerical inversion of the laplace transform: analysis via regularized analytic continuation. Inverse Problems, 22(2):579, 26. [6] M.M. Lavrent ev, V.G. Romanov, and S.P. Shishatski. Ill-posed Problems of Mathematical Physics and Analysis. Translations of Mathematical Monographs. American Mathematical Society, 986. [7] D.V. Widder. The Laplace transform. Princeton mathematical series. Princeton university press, 946., 2 8
Generalized Post-Widder inversion formula with application to statistics
Geeralized Post-Widder iversio formula with applicatio to statistics Deis Belomesty, Hilmar Mai, Joh Schoemakers 3 May 6, 7 Abstract I this work we derive a iversio formula for the Laplace trasform of
More information= 1 2 x (x 1) + 1 {x} (1 {x}). [t] dt = 1 x (x 1) + O (1), [t] dt = 1 2 x2 + O (x), (where the error is not now zero when x is an integer.
Problem Sheet,. i) Draw the graphs for [] and {}. ii) Show that for α R, α+ α [t] dt = α and α+ α {t} dt =. Hint Split these integrals at the integer which must lie in any interval of length, such as [α,
More informationTaylor Series and Asymptotic Expansions
Taylor Series and Asymptotic Epansions The importance of power series as a convenient representation, as an approimation tool, as a tool for solving differential equations and so on, is pretty obvious.
More informationProperties of the Scattering Transform on the Real Line
Journal of Mathematical Analysis and Applications 58, 3 43 (001 doi:10.1006/jmaa.000.7375, available online at http://www.idealibrary.com on Properties of the Scattering Transform on the Real Line Michael
More informationTail Approximation of the Skew-Normal by the Skew-Normal Laplace: Application to Owen s T Function and the Bivariate Normal Distribution
Journal of Statistical and Econometric ethods vol. no. 3 - ISS: 5-557 print version 5-565online Scienpress Ltd 3 Tail Approimation of the Skew-ormal by the Skew-ormal Laplace: Application to Owen s T Function
More informationMTH3101 Spring 2017 HW Assignment 4: Sec. 26: #6,7; Sec. 33: #5,7; Sec. 38: #8; Sec. 40: #2 The due date for this assignment is 2/23/17.
MTH0 Spring 07 HW Assignment : Sec. 6: #6,7; Sec. : #5,7; Sec. 8: #8; Sec. 0: # The due date for this assignment is //7. Sec. 6: #6. Use results in Sec. to verify that the function g z = ln r + iθ r >
More informationOn rate of convergence in distribution of asymptotically normal statistics based on samples of random size
Annales Mathematicae et Informaticae 39 212 pp. 17 28 Proceedings of the Conference on Stochastic Models and their Applications Faculty of Informatics, University of Debrecen, Debrecen, Hungary, August
More informationHarmonic Analysis Homework 5
Harmonic Analysis Homework 5 Bruno Poggi Department of Mathematics, University of Minnesota November 4, 6 Notation Throughout, B, r is the ball of radius r with center in the understood metric space usually
More information8-1 Exploring Exponential Models
8- Eploring Eponential Models Eponential Function A function with the general form, where is a real number, a 0, b > 0 and b. Eample: y = 4() Growth Factor When b >, b is the growth factor Eample: y =
More informationarxiv: v1 [math.cv] 4 Dec 2018
ON THE DE BRANGES LEMMA ALEKSEI KULIKOV arxiv:1812.1728v1 [math.cv] 4 Dec 218 Abstract. It is known that famous de Branges lemma about minimum of two entire functions of minimal type does not etend to
More informationSome New Inequalities for a Sum of Exponential Functions
Applied Mathematical Sciences, Vol. 9, 2015, no. 109, 5429-5439 HIKARI Ltd, www.m-hikari.com http://d.doi.org/10.12988/ams.2015.57471 Some New Inequalities for a Sum of Eponential Functions Steven G. From
More informationLOCAL TIMES OF RANKED CONTINUOUS SEMIMARTINGALES
LOCAL TIMES OF RANKED CONTINUOUS SEMIMARTINGALES ADRIAN D. BANNER INTECH One Palmer Square Princeton, NJ 8542, USA adrian@enhanced.com RAOUF GHOMRASNI Fakultät II, Institut für Mathematik Sekr. MA 7-5,
More informationDIFFERENTIAL EQUATIONS WITH A DIFFERENCE QUOTIENT
Electronic Journal of Differential Equations, Vol. 217 (217, No. 227, pp. 1 42. ISSN: 172-6691. URL: http://ejde.math.tstate.edu or http://ejde.math.unt.edu DIFFERENTIAL EQUATIONS WITH A DIFFERENCE QUOTIENT
More informationPiecewise Smooth Solutions to the Burgers-Hilbert Equation
Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang
More informationStatistical Skorohod embedding problem: optimality and asymptotic normality
Statistical Skorohod embedding problem: optimality and asymptotic normality Denis Belomestny and John Schoenmakers March 4, 5 Abstract Given a Brownian motion B, we consider the so-called statistical Skorohod
More informationAbout the Gamma Function
About the Gamma Function Notes for Honors Calculus II, Originally Prepared in Spring 995 Basic Facts about the Gamma Function The Gamma function is defined by the improper integral Γ) = The integral is
More informationSpatial smoothing using Gaussian processes
Spatial smoothing using Gaussian processes Chris Paciorek paciorek@hsph.harvard.edu August 5, 2004 1 OUTLINE Spatial smoothing and Gaussian processes Covariance modelling Nonstationary covariance modelling
More informationHomework 2 Solutions, MAP 6505, Fall 18. = lim
Homework 2 Solutions, MAP 655, Fall 18 1. Let {a k } R be a sequence. Investigate the convergence of the series a kδ (k) ( k) in the space of distributions D (R) and, as an etra credit, in the space of
More informationStatistics: Learning models from data
DS-GA 1002 Lecture notes 5 October 19, 2015 Statistics: Learning models from data Learning models from data that are assumed to be generated probabilistically from a certain unknown distribution is a crucial
More informationOn Asymptotic Expansions of the Scale Mixtures of Stable Distribution
On Asymptotic Epansions of the Scale Mitures of Stable Distribution Dhaifalla K. Al-Mutairi and Vu Kim Tuan Department of Statistics and Operations Research and Department of Mathematics and Computer Science,
More informationRates of Convergence to Self-Similar Solutions of Burgers Equation
Rates of Convergence to Self-Similar Solutions of Burgers Equation by Joel Miller Andrew Bernoff, Advisor Advisor: Committee Member: May 2 Department of Mathematics Abstract Rates of Convergence to Self-Similar
More informationNew asymptotic expansion for the Γ (x) function
New asymptotic epansion for the Γ function Gergő Nemes December 7, 2008 http://d.doi.org/0.3247/sl2math08.005 Abstract Using a series transformation, Stirling-De Moivre asymptotic series approimation to
More informationAnswer Key 1973 BC 1969 BC 24. A 14. A 24. C 25. A 26. C 27. C 28. D 29. C 30. D 31. C 13. C 12. D 12. E 3. A 32. B 27. E 34. C 14. D 25. B 26.
Answer Key 969 BC 97 BC. C. E. B. D 5. E 6. B 7. D 8. C 9. D. A. B. E. C. D 5. B 6. B 7. B 8. E 9. C. A. B. E. D. C 5. A 6. C 7. C 8. D 9. C. D. C. B. A. D 5. A 6. B 7. D 8. A 9. D. E. D. B. E. E 5. E.
More informationEvaluation of integrals by differentiation with respect to a parameter
December 8 Evaluation of integrals by differentiation with respect to a parameter Khristo N Boyadzhiev Department of Mathematics and Statistics, Ohio Northern University, Ada, OH 458, USA k-boyadzhiev@onuedu
More informationPhenomena in high dimensions in geometric analysis, random matrices, and computational geometry Roscoff, France, June 25-29, 2012
Phenomena in high dimensions in geometric analysis, random matrices, and computational geometry Roscoff, France, June 25-29, 202 BOUNDS AND ASYMPTOTICS FOR FISHER INFORMATION IN THE CENTRAL LIMIT THEOREM
More informationExact Calculation of Normalized Maximum Likelihood Code Length Using Fourier Analysis
Eact Calculation of ormalized Maimum Likelihood Code Length Using Fourier Analysis Atsushi Suzuki and Kenji Yamanishi The University of Tokyo Graduate School of Information Science and Technology Bunkyo,
More informationAdditive functionals of infinite-variance moving averages. Wei Biao Wu The University of Chicago TECHNICAL REPORT NO. 535
Additive functionals of infinite-variance moving averages Wei Biao Wu The University of Chicago TECHNICAL REPORT NO. 535 Departments of Statistics The University of Chicago Chicago, Illinois 60637 June
More informationLecture 2: CDF and EDF
STAT 425: Introduction to Nonparametric Statistics Winter 2018 Instructor: Yen-Chi Chen Lecture 2: CDF and EDF 2.1 CDF: Cumulative Distribution Function For a random variable X, its CDF F () contains all
More informationSobolev Spaces. Chapter 10
Chapter 1 Sobolev Spaces We now define spaces H 1,p (R n ), known as Sobolev spaces. For u to belong to H 1,p (R n ), we require that u L p (R n ) and that u have weak derivatives of first order in L p
More informationThe incomplete gamma functions. Notes by G.J.O. Jameson. These notes incorporate the Math. Gazette article [Jam1], with some extra material.
The incomplete gamma functions Notes by G.J.O. Jameson These notes incorporate the Math. Gazette article [Jam], with some etra material. Definitions and elementary properties functions: Recall the integral
More informationStatistical inference on Lévy processes
Alberto Coca Cabrero University of Cambridge - CCA Supervisors: Dr. Richard Nickl and Professor L.C.G.Rogers Funded by Fundación Mutua Madrileña and EPSRC MASDOC/CCA student workshop 2013 26th March Outline
More informationarxiv: v2 [math.st] 20 Nov 2016
A Lynden-Bell integral estimator for the tail inde of right-truncated data with a random threshold Nawel Haouas Abdelhakim Necir Djamel Meraghni Brahim Brahimi Laboratory of Applied Mathematics Mohamed
More informationNATIONAL UNIVERSITY OF SINGAPORE Department of Mathematics MA4247 Complex Analysis II Lecture Notes Part I
NATIONAL UNIVERSITY OF SINGAPORE Department of Mathematics MA4247 Comple Analsis II Lecture Notes Part I Chapter 1 Preliminar results/review of Comple Analsis I These are more detailed notes for the results
More informationM151B Practice Problems for Exam 1
M151B Practice Problems for Eam 1 Calculators will not be allowed on the eam. Unjustified answers will not receive credit. 1. Compute each of the following its: 1a. 1b. 1c. 1d. 1e. 1 3 4. 3. sin 7 0. +
More informationAPPROXIMATIONS FOR SOME FUNCTIONS OF PRIMES
APPROXIMATIONS FOR SOME FUNCTIONS OF PRIMES Fataneh Esteki B. Sc., University of Isfahan, 2002 A Thesis Submitted to the School of Graduate Studies of the University of Lethbridge in Partial Fulfillment
More informationON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM
ON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM OLEG ZUBELEVICH DEPARTMENT OF MATHEMATICS THE BUDGET AND TREASURY ACADEMY OF THE MINISTRY OF FINANCE OF THE RUSSIAN FEDERATION 7, ZLATOUSTINSKY MALIY PER.,
More informationHarmonic sets and the harmonic prime number theorem
Harmonic sets and the harmonic prime number theorem Version: 9th September 2004 Kevin A. Broughan and Rory J. Casey University of Waikato, Hamilton, New Zealand E-mail: kab@waikato.ac.nz We restrict primes
More informationMath 180A. Lecture 16 Friday May 7 th. Expectation. Recall the three main probability density functions so far (1) Uniform (2) Exponential.
Math 8A Lecture 6 Friday May 7 th Epectation Recall the three main probability density functions so far () Uniform () Eponential (3) Power Law e, ( ), Math 8A Lecture 6 Friday May 7 th Epectation Eample
More informationA GENERALIZATION OF POST-WIDDER OPERATORS
ANALELE ŞTIINŢIFICE ALE UNIVERSITĂŢII AL.I. CUZA DIN IAŞI S.N.) MATEMATICĂ Tomul LXII 16 f.1 A GENERALIZATION OF POST-WIDDER OPERATORS BASED ON -INTEGERS BY DIDEM AYDIN ALI ARAL and GÜLEN BAŞCANBAZ-TUNCA
More informationLecture 5: Rules of Differentiation. First Order Derivatives
Lecture 5: Rules of Differentiation First order derivatives Higher order derivatives Partial differentiation Higher order partials Differentials Derivatives of implicit functions Generalized implicit function
More informationOn the Bennett-Hoeffding inequality
On the Bennett-Hoeffding inequality of Iosif 1,2,3 1 Department of Mathematical Sciences Michigan Technological University 2 Supported by NSF grant DMS-0805946 3 Paper available at http://arxiv.org/abs/0902.4058
More informationOn lower limits and equivalences for distribution tails of randomly stopped sums 1
On lower limits and equivalences for distribution tails of randomly stopped sums 1 D. Denisov, 2 S. Foss, 3 and D. Korshunov 4 Eurandom, Heriot-Watt University and Sobolev Institute of Mathematics Abstract
More informationMath 411, Complex Analysis Definitions, Formulas and Theorems Winter y = sinα
Math 411, Complex Analysis Definitions, Formulas and Theorems Winter 014 Trigonometric Functions of Special Angles α, degrees α, radians sin α cos α tan α 0 0 0 1 0 30 π 6 45 π 4 1 3 1 3 1 y = sinα π 90,
More informationINTRODUCTION TO REAL ANALYTIC GEOMETRY
INTRODUCTION TO REAL ANALYTIC GEOMETRY KRZYSZTOF KURDYKA 1. Analytic functions in several variables 1.1. Summable families. Let (E, ) be a normed space over the field R or C, dim E
More informationSome commonly encountered sets and their notations
NATIONAL UNIVERSITY OF SINGAPORE DEPARTMENT OF MATHEMATICS (This notes are based on the book Introductory Mathematics by Ng Wee Seng ) LECTURE SETS & FUNCTIONS Some commonly encountered sets and their
More informationBLOWUP THEORY FOR THE CRITICAL NONLINEAR SCHRÖDINGER EQUATIONS REVISITED
BLOWUP THEORY FOR THE CRITICAL NONLINEAR SCHRÖDINGER EQUATIONS REVISITED TAOUFIK HMIDI AND SAHBI KERAANI Abstract. In this note we prove a refined version of compactness lemma adapted to the blowup analysis
More informationGeneralized Euler constants
Math. Proc. Camb. Phil. Soc. 2008, 45, Printed in the United Kingdom c 2008 Cambridge Philosophical Society Generalized Euler constants BY Harold G. Diamond AND Kevin Ford Department of Mathematics, University
More informationEstimation of the characteristics of a Lévy process observed at arbitrary frequency
Estimation of the characteristics of a Lévy process observed at arbitrary frequency Johanna Kappus Institute of Mathematics Humboldt-Universität zu Berlin kappus@math.hu-berlin.de Markus Reiß Institute
More informationOn Laplace Finite Marchi Fasulo Transform Of Generalized Functions. A.M.Mahajan
On Laplace Finite Marchi Fasulo Transform Of Generalized Functions A.M.Mahajan Department of Mathematics, Walchand College of Arts and Science, Solapur-43 006 email:-ammahajan9@gmail.com Abstract : In
More informationPart IB Complex Analysis
Part IB Complex Analysis Theorems Based on lectures by I. Smith Notes taken by Dexter Chua Lent 2016 These notes are not endorsed by the lecturers, and I have modified them (often significantly) after
More informationHow large is the class of operator equations solvable by a DSM Newton-type method?
This is the author s final, peer-reviewed manuscript as accepted for publication. The publisher-formatted version may be available through the publisher s web site or your institution s library. How large
More informationRemarks on Bronštein s root theorem
Remarks on Bronštein s root theorem Guy Métivier January 23, 2017 1 Introduction In [Br1], M.D.Bronštein proved that the roots of hyperbolic polynomials (1.1) p(t, τ) = τ m + m p k (t)τ m k. which depend
More informationMath 212a (Fall 2013) Yum-Tong Siu 1. Tauberian Theorems
Math 22a (Fall 203) Yum-Tong Siu Tauberian Theorems For a series a n of comple numbers, the convergence statement a n = A means that the n-th partial sum s n = n k=0 a k as a sequence of comple numbers
More informationTWO PROOFS OF THE PRIME NUMBER THEOREM. 1. Introduction
TWO PROOFS OF THE PRIME NUMBER THEOREM PO-LAM YUNG Introduction Let π() be the number of primes The famous prime number theorem asserts the following: Theorem (Prime number theorem) () π() log as + (This
More informationProblem Set 9 Solutions
8.4 Problem Set 9 Solutions Total: 4 points Problem : Integrate (a) (b) d. ( 4 + 4)( 4 + 5) d 4. Solution (4 points) (a) We use the method of partial fractions to write A B (C + D) = + +. ( ) ( 4 + 5)
More informationVolterra Integral Equations of the First Kind with Jump Discontinuous Kernels
Volterra Integral Equations of the First Kind with Jump Discontinuous Kernels Denis Sidorov Energy Systems Institute, Russian Academy of Sciences e-mail: contact.dns@gmail.com INV Follow-up Meeting Isaac
More informationGrouped Network Vector Autoregression
Statistica Sinica: Supplement Grouped Networ Vector Autoregression Xuening Zhu 1 and Rui Pan 2 1 Fudan University, 2 Central University of Finance and Economics Supplementary Material We present here the
More informationarxiv: v2 [math.pr] 15 Nov 2016
STEIN S METHOD, MANY INTERACTING WORLDS AND QUANTUM MECHANICS Ian W. McKeague, Erol Peköz, and Yvik Swan Columbia University, Boston University, and Université de Liège arxiv:606.0668v2 [math.pr] 5 Nov
More informationLECTURE-15 : LOGARITHMS AND COMPLEX POWERS
LECTURE-5 : LOGARITHMS AND COMPLEX POWERS VED V. DATAR The purpose of this lecture is twofold - first, to characterize domains on which a holomorphic logarithm can be defined, and second, to show that
More informationUNCERTAINTY FUNCTIONAL DIFFERENTIAL EQUATIONS FOR FINANCE
Surveys in Mathematics and its Applications ISSN 1842-6298 (electronic), 1843-7265 (print) Volume 5 (2010), 275 284 UNCERTAINTY FUNCTIONAL DIFFERENTIAL EQUATIONS FOR FINANCE Iuliana Carmen Bărbăcioru Abstract.
More informationGiven the vectors u, v, w and real numbers α, β, γ. Calculate vector a, which is equal to the linear combination α u + β v + γ w.
Selected problems from the tetbook J. Neustupa, S. Kračmar: Sbírka příkladů z Matematiky I Problems in Mathematics I I. LINEAR ALGEBRA I.. Vectors, vector spaces Given the vectors u, v, w and real numbers
More informationConvergence rates in weighted L 1 spaces of kernel density estimators for linear processes
Alea 4, 117 129 (2008) Convergence rates in weighted L 1 spaces of kernel density estimators for linear processes Anton Schick and Wolfgang Wefelmeyer Anton Schick, Department of Mathematical Sciences,
More informationON THE MOMENT DETERMINACY OF PRODUCTS OF NON-IDENTICALLY DISTRIBUTED RANDOM VARIABLES JORDAN S T OYA N OV (NEWCASTLE AND LJUBLJANA)
PROBABILITY AND MATHEMATICAL STATISTICS Vol. 36, Fasc. 1 2016), pp. 21 33 ON THE MOMENT DETERMINACY OF PRODUCTS OF NON-IDENTICALLY DISTRIBUTED RANDOM VARIABLES BY GWO DONG L I N TAIPEI) AND JORDAN S T
More information2. Higher-order Linear ODE s
2. Higher-order Linear ODE s 2A. Second-order Linear ODE s: General Properties 2A-1. On the right below is an abbreviated form of the ODE on the left: (*) y + p()y + q()y = r() Ly = r() ; where L is the
More informationMath 565: Introduction to Harmonic Analysis - Spring 2008 Homework # 2, Daewon Chung
Math 565: Introduction to Harmonic Analysis - Spring 8 Homework #, Daewon Chung. Show that, and that Hχ [a, b] : lim H χ [a, b] a b. H χ [a, b] : sup > H χ [a, b] a b. [Solution] Let us pick < min a, b
More informationAbelian and Tauberian Theorems for Solutions of the One-Dimensional Heat Equation with Respect to Proximate Order
International Journal of Mathematical Analysis Vol. 11, 217, no. 4, 173-188 HIKARI Ltd, www.m-hikari.com https://doi.org/1.12988/ijma.217.714 Abelian and Tauberian Theorems for Solutions of the One-Dimensional
More informationMATH39001 Generating functions. 1 Ordinary power series generating functions
MATH3900 Generating functions The reference for this part of the course is generatingfunctionology by Herbert Wilf. The 2nd edition is downloadable free from http://www.math.upenn. edu/~wilf/downldgf.html,
More informationOn the Convergence of the Summation Formulas Constructed by Using a Symbolic Operator Approach
On the Convergence of the Summation Formulas Constructed by Using a Symbolic Operator Approach Tian-Xiao He 1, Leetsch C. Hsu 2, and Peter J.-S. Shiue 3 1 Department of Mathematics and Computer Science
More informationCharacterizations on Heavy-tailed Distributions by Means of Hazard Rate
Acta Mathematicae Applicatae Sinica, English Series Vol. 19, No. 1 (23) 135 142 Characterizations on Heavy-tailed Distributions by Means of Hazard Rate Chun Su 1, Qi-he Tang 2 1 Department of Statistics
More informationComplex Variables. Cathal Ormond
Complex Variables Cathal Ormond Contents 1 Introduction 3 1.1 Definition: Polar Form.............................. 3 1.2 Definition: Length................................ 3 1.3 Definitions.....................................
More informationChapter 3 Single Random Variables and Probability Distributions (Part 1)
Chapter 3 Single Random Variables and Probability Distributions (Part 1) Contents What is a Random Variable? Probability Distribution Functions Cumulative Distribution Function Probability Density Function
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationScattered Data Interpolation with Polynomial Precision and Conditionally Positive Definite Functions
Chapter 3 Scattered Data Interpolation with Polynomial Precision and Conditionally Positive Definite Functions 3.1 Scattered Data Interpolation with Polynomial Precision Sometimes the assumption on the
More informationMathematical Tripos Part III Michaelmas 2017 Distribution Theory & Applications, Example sheet 1 (answers) Dr A.C.L. Ashton
Mathematical Tripos Part III Michaelmas 7 Distribution Theory & Applications, Eample sheet (answers) Dr A.C.L. Ashton Comments and corrections to acla@damtp.cam.ac.uk.. Construct a non-zero element of
More informationOn probabilities of large and moderate deviations for L-statistics: a survey of some recent developments
UDC 519.2 On probabilities of large and moderate deviations for L-statistics: a survey of some recent developments N. V. Gribkova Department of Probability Theory and Mathematical Statistics, St.-Petersburg
More informationOn the Bennett-Hoeffding inequality
On the Bennett-Hoeffding inequality Iosif 1,2,3 1 Department of Mathematical Sciences Michigan Technological University 2 Supported by NSF grant DMS-0805946 3 Paper available at http://arxiv.org/abs/0902.4058
More informationSYLLABUS FOR ENTRANCE EXAMINATION NANYANG TECHNOLOGICAL UNIVERSITY FOR INTERNATIONAL STUDENTS A-LEVEL MATHEMATICS
SYLLABUS FOR ENTRANCE EXAMINATION NANYANG TECHNOLOGICAL UNIVERSITY FOR INTERNATIONAL STUDENTS A-LEVEL MATHEMATICS STRUCTURE OF EXAMINATION PAPER. There will be one -hour paper consisting of 4 questions..
More informationCOMPLETE MONOTONICITIES OF FUNCTIONS INVOLVING THE GAMMA AND DIGAMMA FUNCTIONS. 1. Introduction
COMPLETE MONOTONICITIES OF FUNCTIONS INVOLVING THE GAMMA AND DIGAMMA FUNCTIONS FENG QI AND BAI-NI GUO Abstract. In the article, the completely monotonic results of the functions [Γ( + 1)] 1/, [Γ(+α+1)]1/(+α),
More informationPointwise convergence rates and central limit theorems for kernel density estimators in linear processes
Pointwise convergence rates and central limit theorems for kernel density estimators in linear processes Anton Schick Binghamton University Wolfgang Wefelmeyer Universität zu Köln Abstract Convergence
More informationExplicit Bounds for the Distribution Function of the Sum of Dependent Normally Distributed Random Variables
Explicit Bounds for the Distribution Function of the Sum of Dependent Normally Distributed Random Variables Walter Schneider July 26, 20 Abstract In this paper an analytic expression is given for the bounds
More information1969 AP Calculus BC: Section I
969 AP Calculus BC: Section I 9 Minutes No Calculator Note: In this eamination, ln denotes the natural logarithm of (that is, logarithm to the base e).. t The asymptotes of the graph of the parametric
More informationReal Analysis Math 125A, Fall 2012 Sample Final Questions. x 1+x y. x y x. 2 (1+x 2 )(1+y 2 ) x y
. Define f : R R by Real Analysis Math 25A, Fall 202 Sample Final Questions f() = 3 + 2. Show that f is continuous on R. Is f uniformly continuous on R? To simplify the inequalities a bit, we write 3 +
More informationDIFFERENTIATING UNDER THE INTEGRAL SIGN
DIFFEENTIATING UNDE THE INTEGAL SIGN KEITH CONAD I had learned to do integrals by various methods shown in a book that my high school physics teacher Mr. Bader had given me. [It] showed how to differentiate
More informationUnbiased simulation of distributions with explicitly known integral transforms
Unbiased simulation of distributions with explicitly known integral transforms Denis Belomestny, Nan Chen, and Yiwei Wang Abstract In this paper, we propose an importance-sampling based method to obtain
More informationPre-Calculus and Trigonometry Capacity Matrix
Pre-Calculus and Capacity Matri Review Polynomials A1.1.4 A1.2.5 Add, subtract, multiply and simplify polynomials and rational epressions Solve polynomial equations and equations involving rational epressions
More informationSmall BGK waves and nonlinear Landau damping (higher dimensions)
Small BGK waves and nonlinear Landau damping higher dimensions Zhiwu Lin and Chongchun Zeng School of Mathematics Georgia Institute of Technology Atlanta, GA, USA Abstract Consider Vlasov-Poisson system
More informationConvergence of greedy approximation I. General systems
STUDIA MATHEMATICA 159 (1) (2003) Convergence of greedy approximation I. General systems by S. V. Konyagin (Moscow) and V. N. Temlyakov (Columbia, SC) Abstract. We consider convergence of thresholding
More informationENEE 621 SPRING 2016 DETECTION AND ESTIMATION THEORY THE PARAMETER ESTIMATION PROBLEM
c 2007-2016 by Armand M. Makowski 1 ENEE 621 SPRING 2016 DETECTION AND ESTIMATION THEORY THE PARAMETER ESTIMATION PROBLEM 1 The basic setting Throughout, p, q and k are positive integers. The setup With
More information10.4: WORKING WITH TAYLOR SERIES
04: WORKING WITH TAYLOR SERIES Contributed by OpenSta Mathematics at OpenSta CNX In the preceding section we defined Taylor series and showed how to find the Taylor series for several common functions
More informationUpper Bounds for the Ruin Probability in a Risk Model With Interest WhosePremiumsDependonthe Backward Recurrence Time Process
Λ4flΛ4» ν ff ff χ Vol.4, No.4 211 8fl ADVANCES IN MATHEMATICS Aug., 211 Upper Bounds for the Ruin Probability in a Risk Model With Interest WhosePremiumsDependonthe Backward Recurrence Time Process HE
More informationFerromagnetic Ising models
Stat 36 Stochastic Processes on Graphs Ferromagnetic Ising models Amir Dembo Lecture 3-4 - 0/5/2007 Introduction By ferromagnetic Ising models we mean the study of the large-n behavior of the measures
More informationIRRATIONAL FACTORS SATISFYING THE LITTLE FERMAT THEOREM
International Journal of Number Theory Vol., No. 4 (005) 499 5 c World Scientific Publishing Company IRRATIONAL FACTORS SATISFYING THE LITTLE FERMAT THEOREM Int. J. Number Theory 005.0:499-5. Downloaded
More informationFourier transform of tempered distributions
Fourier transform of tempered distributions 1 Test functions and distributions As we have seen before, many functions are not classical in the sense that they cannot be evaluated at any point. For eample,
More information4 Inverse function theorem
Tel Aviv University, 2014/15 Analysis-III,IV 71 4 Inverse function theorem 4a What is the problem................ 71 4b Simple observations before the theorem..... 72 4c The theorem.....................
More informationTools from Lebesgue integration
Tools from Lebesgue integration E.P. van den Ban Fall 2005 Introduction In these notes we describe some of the basic tools from the theory of Lebesgue integration. Definitions and results will be given
More informationMA3111S COMPLEX ANALYSIS I
MA3111S COMPLEX ANALYSIS I 1. The Algebra of Complex Numbers A complex number is an expression of the form a + ib, where a and b are real numbers. a is called the real part of a + ib and b the imaginary
More informationStochastic Volatility and Correction to the Heat Equation
Stochastic Volatility and Correction to the Heat Equation Jean-Pierre Fouque, George Papanicolaou and Ronnie Sircar Abstract. From a probabilist s point of view the Twentieth Century has been a century
More informationB.Tech. Theory Examination (Semester IV) Engineering Mathematics III
Solved Question Paper 5-6 B.Tech. Theory Eamination (Semester IV) 5-6 Engineering Mathematics III Time : hours] [Maimum Marks : Section-A. Attempt all questions of this section. Each question carry equal
More informationStatistical Geometry Processing Winter Semester 2011/2012
Statistical Geometry Processing Winter Semester 2011/2012 Linear Algebra, Function Spaces & Inverse Problems Vector and Function Spaces 3 Vectors vectors are arrows in space classically: 2 or 3 dim. Euclidian
More informationOn the Stability of a Differential-Difference Analogue of a Two-Dimensional Problem of Integral Geometry
Filomat 3:3 18, 933 938 https://doi.org/1.98/fil183933b Published b Facult of Sciences and Mathematics, Universit of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat On the Stabilit of a Differential-Difference
More information