Researcher profile

E. Di Nardo

E. Di Nardo contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

11 published item(s)

preprint2022arXiv

kStatistics: Unbiased Estimates of Joint Cumulant Products from the Multivariate Faà Di Bruno's Formula

kStatistics is a package in R that serves as a unified framework for estimating univariate and multivariate cumulants as well as products of univariate and multivariate cumulants of a random sample, using unbiased estimators with minimum variance. The main computational machinery of kStatistics is an algorithm for computing multi-index partitions. The same algorithm underlies the general-purpose multivariate Faà di Bruno's formula, which has been therefore included in the last release of the package. This formula gives the coefficients of formal power series compositions as well as the partial derivatives of multivariable function compositions. One of the most significant applications of this formula is the possibility to generate many well-known polynomial families as special cases. So, in the package, there are special functions for generating very popular polynomial families, such as the Bell polynomials. However further families can be obtained, for suitable choices of the formal power series involved in the composition or when suitable symbolic strategies are employed. In both cases, we give examples on how to modify the R codes of the package to accomplish this task. Future developments are addressed at the end of the paper.

preprint2021arXiv

Input-output behaviour of a model neuron with alternating drift

The input-output behaviour of the Wiener neuronal model subject to alternating input is studied under the assumption that the effect of such an input is to make the drift itself of an alternating type. Firing densities and related statistics are obtained via simulations of the sample-paths of the process in the following three cases: the drift changes occur during random periods characterized by (i) exponential distribution, (ii) Erlang distribution with a preassigned shape parameter, and (iii) deterministic distribution. The obtained results are compared with those holding for the Wiener neuronal model subject to sinusoidal input

preprint2021arXiv

On the connection between orthant probabilities and the first passage time problem

This article describes a new Monte Carlo method for the evaluation of the orthant probabilities by sampling first passage times of a non-singular Gaussian discrete time-series across an absorbing boundary. This procedure makes use of a simulation of several time-series sample paths, aiming to record their first crossing instants. Thus, the computation of the orthant probabilities is traced back to the accurate simulation of a non-singular Gaussian discrete-time series. Moreover, if the simulation is also efficient, this method is shown to be more speedy than the others proposed in the literature. As example, we make use of the Davies-Harte algorithm in the evaluation of the orthant probabilities associated to the ARFIMA$(0,d,0)$ model. Test results are presented that compare this method with currently available software.

preprint2021arXiv

Simulation of first-passage times for alternating Brownian motions

The first-passage-time problem for a Brownian motion with alternating infinitesimal moments through a constant boundary is considered under the assumption that the time intervals between consecutive changes of these moments are described by an alternating renewal process. Bounds to the first-passage-time density and distribution function are obtained, and a simulation procedure to estimate first-passage-time densities is constructed. Examples of applications to problems in environmental sciences and mathematical finance are also provided.

preprint2021arXiv

Symbolic solutions of some linear recurrences

A symbolic method for solving linear recurrences of combinatorial and statistical interest is introduced. This method essentially relies on a representation of polynomial sequences as moments of a symbol that looks as the framework of a random variable with no reference to any probability space. We give several examples of applications and state an explicit form for the class of linear recurrences involving Sheffer sequences satisfying a special initial condition. The results here presented can be easily implemented in a symbolic software.

preprint2020arXiv

Polynomial traces and elementary symmetric functions in the latent roots of a non-central Wishart matrix

Hypergeometric functions and zonal polynomials are the tools usually addressed in the literature to deal with the expected value of the elementary symmetric functions in non-central Wishart latent roots. The method here proposed recovers the expected value of these symmetric functions by using the umbral operator applied to the trace of suitable polynomial matrices and their cumulants. The employment of a suitable linear operator in place of hypergeometric functions and zonal polynomials was conjectured by de Waal in 1972. Here we show how the umbral operator accomplishes this task and consequently represents an alternative tool to deal with these symmetric functions. When special formal variables are plugged in the variables, the evaluation through the umbral operator deletes all the monomials in the latent roots except those contributing in the elementary symmetric functions. Cumulants further simplify the computations taking advantage of the convolution structure of the polynomial trace. Open problems are addressed at the end of the paper.

preprint2013arXiv

On a representation of time space-harmonic polynomials via symbolic Lévy processes

In this paper, we review the theory of time space-harmonic polynomials developed by using a symbolic device known in the literature as the classical umbral calculus. The advantage of this symbolic tool is twofold. First a moment representation is allowed for a wide class of polynomial stochastic involving the Lévy processes in respect to which they are martingales. This representation includes some well-known examples such as Hermite polynomials in connection with Brownian motion. As a consequence, characterizations of many other families of polynomials having the time space-harmonic property can be recovered via the symbolic moment representation. New relations with Kailath-Segall polynomials are stated. Secondly the generalization to the multivariable framework is straightforward. Connections with cumulants and Bell polynomials are highlighted both in the univariate case and in the multivariate one. Open problems are addressed at the end of the paper.

preprint2012arXiv

A new family of time-space harmonic polynomials with respect to Lévy processes

By means of a symbolic method, a new family of time-space harmonic polynomials with respect to Lévy processes is given. The coefficients of these polynomials involve a formal expression of Lévy processes by which many identities are stated. We show that this family includes classical families of polynomials such as Hermite polynomials. Poisson-Charlier polynomials result to be a linear combinations of these new polynomials, when they have the property to be time-space harmonic with respect to the compensated Poisson process. The more general class of Lévy-Sheffer polynomials is recovered as a linear combination of these new polynomials, when they are time-space harmonic with respect to Lévy processes of very general form. We show the role played by cumulants of Lévy processes so that connections with boolean and free cumulants are also stated.

preprint2012arXiv

Multivariate Bernoulli and Euler polynomials via Lévy processes

By a symbolic method, we introduce multivariate Bernoulli and Euler polynomials as powers of polynomials whose coefficients involve multivariate Lévy processes. Many properties of these polynomials are stated straightforwardly thanks to this representation, which could be easily implemented in any symbolic manipulation system. A very simple relation between these two families of multivariate polynomials is provided.

preprint2010arXiv

A new algorithm for computing the multivariate Faà di Bruno's formula

A new algorithm for computing the multivariate Faà di Bruno's formula is provided. We use a symbolic approach based on the classical umbral calculus that turns the computation of the multivariate Faà di Bruno's formula into a suitable multinomial expansion. We propose a MAPLE procedure whose computational times are faster compared with the ones existing in the literature. Some illustrative applications are also provided.

preprint2010arXiv

Cumulants and convolutions via Abel polynomials

We provide an unifying polynomial expression giving moments in terms of cumulants, and viceversa, holding in the classical, boolean and free setting. This is done by using a symbolic treatment of Abel polynomials. As a by-product, we show that in the free cumulant theory the volume polynomial of Pitman and Stanley plays the role of the complete Bell exponential polynomial in the classical theory. Moreover via generalized Abel polynomials we construct a new class of cumulants, including the classical, boolean and free ones, and the convolutions linearized by them. Finally, via an umbral Fourier transform, we state a explicit connection between boolean and free convolution.