BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Nico M. Temme

Nico M. Temme appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.CAmath.NAmath.PRmath.STNumerical AnalysisStatistics TheoryMethodologyPopulations and Evolutionq-fin.CP

Catalog footprint

What is connected

19works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2022arXiv

Solution of integrals with fractional Brownian motion for different Hurst indices

In this paper, we will evaluate integrals that define the conditional expectation, variance and characteristic function of stochastic processes with respect to fractional Brownian motion (fBm) for all relevant Hurst indices, i.e. $H \in (0,1)$. The fractional Ornstein-Uhlenbeck (fOU) process, for example, gives rise to highly nontrivial integration formulas that need careful analysis when considering the whole range of Hurst indices. We will show that the classical technique of analytic continuation, from complex analysis, provides a way of extending the domain of validity of an integral, from $H\in(1/2,1)$, to the larger domain, $H\in(0,1)$. Numerical experiments for different Hurst indices confirm the robustness and efficiency of the integral formulations presented here. Moreover, we provide accurate and highly efficient financial option pricing results for processes that are related to the fOU process, with the help of Fourier cosine expansions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A faster and more accurate algorithm for calculating population genetics statistics requiring sums of Stirling numbers of the first kind

Stirling numbers of the first kind are used in the derivation of several population genetics statistics, which in turn are useful for testing evolutionary hypotheses directly from DNA sequences. Here, we explore the cumulative distribution function of these Stirling numbers, which enables a single direct estimate of the sum, using representations in terms of the incomplete beta function. This estimator enables an improved method for calculating an asymptotic estimate for one useful statistic, Fu's $F_s$. By reducing the calculation from a sum of terms involving Stirling numbers to a single estimate, we simultaneously improve accuracy and dramatically increase speed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Asymptotic expansions of Jacobi polynomials and of the nodes and weights of Gauss-Jacobi quadrature for large degree and parameters in terms of elementary functions

Asymptotic approximations of Jacobi polynomials are given in terms of elementary functions for large degree $n$ and parameters $α$ and $β$. From these new results, asymptotic expansions of the zeros are derived and methods are given to obtain the coefficients in the expansions. These approximations can be used as initial values in iterative methods for computing the nodes of Gauss--Jacobi quadrature for large degree and parameters. The performance of the asymptotic approximations for computing the nodes and weights of these Gaussian quadratures is illustrated with numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Numerical evaluation of Airy-type integrals arising in uniform asymptotic analysis

We describe a method to evaluate integrals that arise in the asymptotic analysis when two saddle points may be close together. These integrals, which appear in problems from optics, acoustics or quantum mechanics as well as in a wide class of special functions, can be transformed into Airy-type integrals and we use the trapezoidal rule to compute these integrals numerically. The quadrature method, which remains valid when two saddle points coalesce, is illustrated with numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

An improved algorithm and a Fortran 90 module for computing the conical function $P^m_{-1/2+iτ}(x)$

In this paper we describe an algorithm and a Fortran 90 module ({\bf Conical}) for the computation of the conical function $P^m_{-\tfrac12+iτ}(x)$ for $x>-1$, $m \ge 0$, $τ>0$. These functions appear in the solution of Dirichlet problems for domains bounded by cones; because of this, they are involved in a large number of applications in Engineering and Physics. In the Fortran 90 module, the admissible parameter ranges for computing the conical functions in standard IEEE double precision arithmetic are restricted to $(x,m,τ) \in (-1,1) \times [0,\,40] \times [0,\,100]$ and $(x,m,τ) \in (1,100) \times [0,\,100] \times [0,\,100]$. Based on tests of the three-term recurrence relation satisfied by these functions and direct comparison with Maple, we claim a relative accuracy close to $10^{-12}$ in the full parameter range, although a mild loss of accuracy can be found at some points of the oscillatory region of the conical functions. The relative accuracy increases to $10^{-13}\,-\,10^{-14}$ in the region of the monotonic regime of the functions where integral representations are computed ($-1<x<0$).

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Efficient and accurate algorithms for the computation and inversion of the incomplete gamma function ratios

Algorithms for the numerical evaluation of the incomplete gamma function ratios $P(a,x)=γ(a,x)/Γ(a)$ and $Q(a,x)=Γ(a,x)/Γ(a)$ are described for positive values of $a$ and $x$. Also, inversion methods are given for solving the equations $P(a,x)=p$, $Q(a,x)=q$, with $0<p,q<1$. Both the direct computation and the inversion of the incomplete gamma function ratios are used in many problems in statistics and applied probability. The analytical approach from earlier literature is summarized and new initial estimates are derived for starting the inversion algorithms. The performance of the associated software to our algorithms (the Fortran 90 module {\bf IncgamFI}) is analyzed and compared with earlier published algorithms.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

New Series Expansions of the Gauss Hypergeometric Function

The Gauss hypergeometric function ${}_2F_1(a,b,c;z)$ can be computed by using the power series in powers of $z, z/(z-1), 1-z, 1/z, 1/(1-z),(z-1)/z$. With these expansions ${}_2F_1(a,b,c;z)$ is not completely computable for all complex values of $z$. As pointed out in Gil, {\it et al.} [2007, \S2.3], the points $z=e^{\pm iπ/3}$ are always excluded from the domains of convergence of these expansions. Bühring [1987] has given a power series expansion that allows computation at and near these points. But, when $b-a$ is an integer, the coefficients of that expansion become indeterminate and its computation requires a nontrivial limiting process. Moreover, the convergence becomes slower and slower in that case. In this paper we obtain new expansions of the Gauss hypergeometric function in terms of rational functions of $z$ for which the points $z=e^{\pm iπ/3}$ are well inside their domains of convergence . In addition, these expansion are well defined when $b-a$ is an integer and no limits are needed in that case. Numerical computations show that these expansions converge faster than Bühring's expansion for $z$ in the neighborhood of the points $e^{\pm iπ/3}$, especially when $b-a$ is close to an integer number.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Uniform Asymptotic Methods for Integrals

We give an overview of basic methods that can be used for obtaining asymptotic expansions of integrals: Watson's lemma, Laplace's method, the saddle point method, and the method of stationary phase. Certain developments in the field of asymptotic analysis will be compared with De Bruijn's book {\em Asymptotic Methods in Analysis}. The classical methods can be modified for obtaining expansions that hold uniformly with respect to additional parameters. We give an overview of examples in which special functions, such as the complementary error function, Airy functions, and Bessel functions, are used as approximations in uniform asymptotic expansions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Large Degree Asymptotics of Generalized Bessel Polynomials

Asymptotic expansions are given for large values of $n$ of the generalized Bessel polynomials $Y_n^μ(z)$. The analysis is based on integrals that follow from the generating functions of the polynomials. A new simple expansion is given that is valid outside a compact neighborhood of the origin in the $z-$plane. New forms of expansions in terms of elementary functions valid in sectors not containing the turning points $z=\pm i/n$ are derived, and a new expansion in terms of modified Bessel functions is given. Earlier asymptotic expansions of the generalized Bessel polynomials by Wong and Zhang (1997) and Dunster (2001) are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The tail of the maximum of Brownian motion minus a parabola

We analyze the tail behavior of the maximum N of Brownian motion minus a parabola and give an asymptotic expansion for P(N>x) as x tends to infinity. This extends a first order result on the tail behavior, which can be deduced from Huesler and Piterbarg (1999). We also point out the relation between certain results in Groeneboom (2010) and Janson, Louchard and Martin-Löf (2010).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A sampling theory for asymmetric communities

We introduce the first analytical model of asymmetric community dynamics to yield Hubbell's neutral theory in the limit of functional equivalence among all species. Our focus centers on an asymmetric extension of Hubbell's local community dynamics, while an analogous extension of Hubbell's metacommunity dynamics is deferred to an appendix. We find that mass-effects may facilitate coexistence in asymmetric local communities and generate unimodal species abundance distributions indistinguishable from those of symmetric communities. Multiple modes, however, only arise from asymmetric processes and provide a strong indication of non-neutral dynamics. Although the exact stationary distributions of fully asymmetric communities must be calculated numerically, we derive approximate sampling distributions for the general case and for nearly neutral communities where symmetry is broken by a single species distinct from all others in ecological fitness and dispersal ability. In the latter case, our approximate distributions are fully normalized, and novel asymptotic expansions of the required hypergeometric functions are provided to make evaluations tractable for large communities. Employing these results in a Bayesian analysis may provide a novel statistical test to assess the consistency of species abundance data with the neutral hypothesis.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Asymptotics of a ${}_3F_2$ polynomial associated with the Catalan-Larcombe-French sequence

The large $n$ behaviour of the hypergeometric polynomial $$\FFF{-n}{\sfrac12}{\sfrac12}{\sfrac12-n}{\sfrac12-n}{-1}$$ is considered by using integral representations of this polynomial. This ${}_3F_2$ polynomial is associated with the Catalan-Larcombe-French sequence. Several other representations are mentioned, with references to the literature, and another asymptotic method is described by using a generating function of the sequence. The results are similar to those obtained by Clark (2004) who used a binomial sum for obtaining an asymptotic expansion.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

The ABC of Hyper Recursions

Each family of Gauss hypergeometric functions $$ f_n={}_2F_1(a+ε_1n, b+ε_2n ;c+ε_3n; z), $$ for fixed $ε_j=0,\pm1$ (not all $ε_j$ equal to zero) satisfies a second order linear difference equation of the form $$ A_nf_{n-1}+B_nf_n+C_nf_{n+1}=0. $$ Because of symmetry relations and functional relations for the Gauss functions, many of the 26 cases (for different $ε_j$ values) can be transformed into each other. We give a set of basic equations from which all other equations can be obtained. For each basic equation, we study the existence of minimal solutions and the character of $f_n$ (minimal or dominant) as $n\to \pm\infty$. A second independent solution is given in each basic case which is dominant when $f_n$ is minimal and vice-versa. In this way, satisfactory pairs of linearly independent solutions for each of the 26 second order linear difference equations can be obtained.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Large Parameter Cases of the Gauss Hypergeometric Function

We consider the asymptotic behaviour of the Gauss hypergeometric function when several of the parameters a, b, c are large. We indicate which cases are of interest for orthogonal polynomials (Jacobi, but also Krawtchouk, Meixner, etc.), which results are already available and which cases need more attention. We also consider a few examples of 3F2-functions of unit argument, to explain which difficulties arise in these cases, when standard integrals or differential equations are not available.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

The Askey Scheme for Hypergeometric Orthogonal Polynomials Viewed from Asymptotic Analysis

Many limits are known for hypergeometric orthogonal polynomials that occur in the Askey scheme. We show how asymptotic representations can be derived by using the generating functions of the polynomials. For example, we discuss the asymptotic representation of the Meixner-Pollaczek, Jacobi, Meixner, and Krawtchouk polynomials in terms of Laguerre polynomials.

0 citations0 reviewsNo code linkedOpen access