We investigate some aspects of the problem of the estimation of birth distributions (BDs) in multi-type Galton–Watson trees (MGWs) with unobserved types. More precisely, we consider two-type MGWs called spinal-structured trees. This kind of tree is characterized by a spine of special individuals whose BD $\nu$ is different from the other individuals in the tree (called normal, and whose BD is denoted by $\mu$). In this work, we show that even in such a very structured two-type population, our ability to distinguish the two types and estimate $\mu$ and $\nu$ is constrained by a trade-off between the growth-rate of the population and the similarity of $\mu$ and $\nu$. Indeed, if the growth-rate is too large, large deviation events are likely to be observed in the sampling of the normal individuals, preventing us from distinguishing them from special ones. Roughly speaking, our approach succeeds if $r\lt \mathfrak{D}(\mu,\nu)$, where r is the exponential growth-rate of the population and $\mathfrak{D}$ is a divergence measuring the dissimilarity between $\mu$ and $\nu$.

We discuss a variant, named ‘Rattle’, of the product replacement algorithm. Rattle is a Markov chain, that returns a random element of a black box group. The limiting distribution of the element returned is the uniform distribution. We prove that, if the generating sequence is long enough, the probability distribution of the element returned converges unexpectedly quickly to the uniform distribution.

We extend the construction principle of phase-type (PH) distributions to allow for inhomogeneous transition rates and show that this naturally leads to direct probabilistic descriptions of certain transformations of PH distributions. In particular, the resulting matrix distributions enable the carrying over of fitting properties of PH distributions to distributions with heavy tails, providing a general modelling framework for heavy-tail phenomena. We also illustrate the versatility and parsimony of the proposed approach in modelling a real-world heavy-tailed fire insurance dataset.

We consider a stochastic evolutionary model for a phenotype developing amongst n related species with unknown phylogeny. The unknown tree is modelled by a Yule process conditioned on n contemporary nodes. The trait value is assumed to evolve along lineages as an Ornstein-Uhlenbeck process. As a result, the trait values of the n species form a sample with dependent observations. We establish three limit theorems for the sample mean corresponding to three domains for the adaptation rate. In the case of fast adaptation, we show that for large n the normalized sample mean is approximately normally distributed. Using these limit theorems, we develop novel confidence interval formulae for the optimal trait value.

Financial data are as a rule asymmetric, although most econometric models are symmetric. This applies also to continuous-time models for high-frequency and irregularly spaced data. We discuss some asymmetric versions of the continuous-time GARCH model, concentrating then on the GJR-COGARCH model. We calculate higher-order moments and extend the first-jump approximation. These results are prerequisites for moment estimation and pseudo maximum likelihood estimation of the GJR-COGARCH model parameters, respectively, which we derive in detail.

We study nonparametric estimation problems for discrete-time stochastic networks of GeomX/G/∞ queues. We assume that we are only able to observe the external arrival and external departure processes at the nodes over a stretch of time. Based on such incomplete information of the system, we aim to construct estimators for the unknown general service time distributions at the nodes without imposing any parametric condition. We propose two different estimation approaches. The first approach is based on the construction of a so-called sequence of differences, and a crucial relation between the expected number of external departures at a node and specific sojourn time distributions in the network. The second approach directly utilizes the structure of the cross-covariance functions between external arrival and departure processes at the nodes. Both methods lead to deconvolution problems which we solve explicitly. A detailed simulation study illustrates the numerical performances of our estimators and shows their advantages and disadvantages.

Many complex systems can be modeled via Markov jump processes. Applications include chemical reactions, population dynamics, and telecommunication networks. Rare-event estimation for such models can be difficult and is often computationally expensive, because typically many (or very long) paths of the Markov jump process need to be simulated in order to observe the rare event. We present a state-dependent importance sampling approach to this problem that is adaptive and uses Markov chain Monte Carlo to sample from the zero-variance importance sampling distribution. The method is applicable to a wide range of Markov jump processes and achieves high accuracy, while requiring only a small sample to obtain the importance parameters. We demonstrate its efficiency through benchmark examples in queueing theory and stochastic chemical kinetics.

In Markov Chain Monte Carlo (MCMC) simulations, thermal equilibria quantities are estimated by ensemble average over a sample set containing a large number of correlated samples. These samples are selected in accordance with the probability distribution function, known from the partition function of equilibrium state. As the stochastic error of the simulation results is significant, it is desirable to understand the variance of the estimation by ensemble average, which depends on the sample size (i.e., the total number of samples in the set) and the sampling interval (i.e., cycle number between two consecutive samples). Although large sample sizes reduce the variance, they increase the computational cost of the simulation. For a given CPU time, the sample size can be reduced greatly by increasing the sampling interval, while having the corresponding increase in variance be negligible if the original sampling interval is very small. In this work, we report a few general rules that relate the variance with the sample size and the sampling interval. These results are observed and confirmed numerically. These variance rules are derived for the MCMC method but are also valid for the correlated samples obtained using other Monte Carlo methods. The main contribution of this work includes the theoretical proof of these numerical observations and the set of assumptions that lead to them.

We consider the problem of estimation in a partially observed discrete-time Galton-Watson branching process, focusing on the first two moments of the offspring distribution. Our study is motivated by modelling the counts of new cases at the onset of a stochastic epidemic, allowing for the facts that only a part of the cases is detected, and that the detection mechanism may affect the evolution of the epidemic. In this setting, the offspring mean is closely related to the spreading potential of the disease, while the second moment is connected to the variability of the mean estimators. Inference for branching processes is known for its nonstandard characteristics, as compared with classical inference. When, in addition, the true process cannot be directly observed, the problem of inference suffers significant further perturbations. We propose nonparametric estimators related to those used when the underlying process is fully observed, but suitably modified to take into account the intricate dependence structure induced by the partial observation and the interaction scheme. We show consistency, derive the limiting laws of the estimators, and construct asymptotic confidence intervals, all valid conditionally on the explosion set.

Let {Zt}t≥0 be a Lévy process with Lévy measure ν, and let τ(t)=∫0tr(u) d u, where {r(t)}t≥0 is a positive ergodic diffusion independent from Z. Based upon discrete observations of the time-changed Lévy process Xt≔Zτt during a time interval [0,T], we study the asymptotic properties of certain estimators of the parameters β(φ)≔∫φ(x)ν(d x), which in turn are well known to be the building blocks of several nonparametric methods such as sieve-based estimation and kernel estimation. Under uniform boundedness of the second moments of r and conditions on φ necessary for the standard short-term ergodic property limt→ 0 E φ(Zt)/t = β(φ) to hold, consistency and asymptotic normality of the proposed estimators are ensured when the time horizon T increases in such a way that the sampling frequency is high enough relative to T.

Start with a necklace consisting of one white bead and one black bead, and add new beads one at a time by inserting each new bead between a randomly chosen adjacent pair of old beads, with the proviso that the new bead will be white if and only if both beads of the adjacent pair are black. Let Wn denote the number of white beads when the total number of beads is n. We show that EWn = n/3 and, with c2 = 2/45, that (Wn − n/3) / c√n is asymptotically standard normal. We find that, for all r ≥ 1 and n > 2r, the rth cumulant of the distribution of Wn is of the form nhr. We find the expected numbers of gaps of given length between white beads, and examine the asymptotics of the longest gaps.

Based on the concept of multipower variation we establish a class of easily computable and robust estimators for the integrated volatility, especially including the squared integrated volatility, in Lévy-type stochastic volatility models. We derive consistency and feasible distributional results for the estimators. Furthermore, we discuss the applications to time-changed CGMY, normal inverse Gaussian, and hyperbolic models with and without leverage, where the time-changes are based on integrated Cox-Ingersoll-Ross or Ornstein-Uhlenbeck-type processes. We deduce which type of market microstructure does not affect the estimates.

We investigate a sequence of Galton-Watson branching processes with immigration, where the offspring mean tends to its critical value 1 and the offspring variance tends to 0. It is shown that the fluctuation limit is an Ornstein-Uhlenbeck-type process. As a consequence, in contrast to the case in which the offspring variance tends to a positive limit, it transpires that the conditional least-squares estimator of the offspring mean is asymptotically normal. The norming factor is n3/2, in contrast to both the subcritical case, in which it is n1/2, and the nearly critical case with positive limiting offspring variance, in which it is n.

This paper introduces a new estimation technique for discretely observed diffusion processes. Transform functions are applied to transform the data to obtain good and easily calculated estimators of both the drift and diffusion coefficients. Consistency and asymptotic normality of the resulting estimators is investigated. Power transforms are used to estimate the parameters of affine diffusions, for which explicit estimators are obtained.

This paper investigates the probabilistic behaviour of the eigenvalue of the empirical transition matrix of a Markov chain which is of largest modulus other than 1, loosely called the second-largest eigenvalue. A central limit theorem is obtained for nonmultiple eigenvalues of the empirical transition matrix. When the Markov chain is actually a sequence of independent observations the distribution of the second-largest eigenvalue is determined and a test for independence is developed. The independence case is considered in more detail when the Markov chain has only two states, and some applications are given.

This paper has as its main theme the fitting in practice of the variance-gamma distribution, which allows for skewness, by moment methods. This fitting procedure allows for possible dependence of increments in log returns, while retaining their stationarity. It is intended as a step in a partial synthesis of some ideas of Madan, Carr and Chang (1998) and of Heyde (1999). Standard estimation and hypothesis-testing theory depends on a large sample of observations which are independently as well as identically distributed and consequently may give inappropriate conclusions in the presence of dependence.

The Bernstein-von Mises theorem, concerning the convergence of suitably normalized and centred posterior density to normal density, is proved for a certain class of linearly parametrized parabolic stochastic partial differential equations (SPDEs) as the number of Fourier coefficients in the expansion of the solution increases to infinity. As a consequence, the Bayes estimators of the drift parameter, for smooth loss functions and priors, are shown to be strongly consistent, asymptotically normal and locally asymptotically minimax (in the Hajek-Le Cam sense), and asymptotically equivalent to the maximum likelihood estimator as the number of Fourier coefficients become large. Unlike in the classical finite dimensional SDEs, here the total observation time and the intensity of noise remain fixed.

We consider the estimation of Markov transition matrices by Bayes’ methods. We obtain large and moderate deviation principles for the sequence of Bayesian posterior distributions.

We consider a parametrization of the Heath-Jarrow-Morton (HJM) family of term structure of interest rate models that allows a finite-dimensional Markovian representation of the stochastic dynamics. This parametrization results from letting the volatility function depend on time to maturity and on two factors: the instantaneous spot rate and one fixed-maturity forward rate. Our main purpose is an estimation methodology for which we have to model the observations under the historical probability measure. This leads us to consider as an additional third factor the market price of interest rate risk, that connects the historical and the HJM martingale measures. Assuming that the information comes from noisy observations of the fixed-maturity forward rate, the purpose is to estimate recursively, on the basis of this information, the three Markovian factors as well as the parameters in the model, in particular those in the volatility function. This leads to a nonlinear filtering problem, for the solution of which we describe an approximation methodology, based on time discretization and quantization. We prove the convergence of the approximate filters for each of the observed trajectories.