In this paper we introduce multitype branching processes with inhomogeneous Poisson immigration, and consider in detail the critical Markov case when the local intensity r(t) of the Poisson random measure is a regularly varying function. Various multitype limit distributions (conditional and unconditional) are obtained depending on the rate at which r(t) changes with time. The asymptotic behaviour of the first and second moments, and the probability of nonextinction are investigated.

In this paper we provide an introduction to statistical inference for the classical linear birth‒death process, focusing on computational aspects of the problem in the setting of discretely observed processes. The basic probabilistic properties are given in Section 2, focusing on computation of the transition functions. This is followed by a brief discussion of simulation methods in Section 3, and of frequentist methods in Section 4. Section 5 is devoted to Bayesian methods, from rejection sampling to Markov chain Monte Carlo and approximate Bayesian computation. In Section 6 we consider the time-inhomogeneous case. The paper ends with a brief discussion in Section 7.

The aging of population is perhaps the most important problem that developed countries must face in the near future. Dependency can be seen as a consequence of the process of gradual aging. In a health context, this contingency is defined as a lack of autonomy in performing basic activities of daily living that requires the care of another person or significant help. In Europe in general and in Spain in particular, this phenomena represents a problem with economic, political, social and demographic implications. The prevalence of dependency in the population, as well as its intensity and evolution over the course of a person’s life are issues of greatest importance that should be addressed. The aim of this work is the estimation of life expectancy free of dependency (LEFD) based on functional trajectories to enhance the regular estimation of health expectancy. Using information from the Spanish survey EDAD 2008, we estimate the number of years spent free of dependency for disabled people according to gender, dependency degree (moderate, severe, major) and the earlier or later onset of dependency compared to a central trend. The main findings are as follows: first, we show evidence that to estimate LEFD ignoring the information provided by the functional trajectories may lead to non-representative LEFD estimates; second, in general, dependency-free life expectancy is higher for women than for men. However, its intensity is higher in women with later onset on dependency; Third, the loss of autonomy is higher (and more abrupt) in men than in women. Finally, the diversity of patterns observed at later onset of dependency tends to a dependency extreme-pattern in both genders.

This work investigates analytically, the use of piezoelectric tiles placed on stairways for vibrational energy harvesting – harnessing electrical power from natural vibrational phenomena – from pedestrian footfalls. While energy harvesting from pedestrian traffic along flat pathways has been studied in the linear regime and realised in practical applications, the greater amounts of energy naturally expended in traversing stairways suggest better prospects for harvesting. Considering the characteristics of two types of commercially available piezoelectric tiles – Navy Type III and Navy Type V – analytical models for the coupled electromechanical system are formulated. The harvesting potential of the tiles is then studied under conditions of both deterministic and carefully developed random excitation profiles for three distinct cases: linear, monostable nonlinear and an array of monostable nonlinear tiles on adjacent steps with linear coupling between them. The results indicate enhanced power output when the tiles are: (1) placed on stairways, (2) uncoupled and (3) subjected to excitation profiles with stochastic frequency. In addition, the Navy Type V tiles are seen to outperform the Navy Type III tiles. Finally, the strongly nonlinear regime outperforms the linear one suggesting that the realisation of commercially available piezoelectric tiles with appropriately tailored nonlinear characteristics will likely have a significant impact on energy harvesting from pedestrian traffic.

We consider the pricing of European options under a modified Black–Scholes equation having fractional derivatives in the “spatial” (price) variable. To be specific, the underlying price is assumed to follow a geometric Koponen–Boyarchenko–Levendorski process. This pure jump Lévy process could better capture the real behaviour of market data. Despite many difficulties caused by the “globalness” of the fractional derivatives, we derive an explicit closed-form analytical solution by solving the fractional partial differential equation analytically, using the Fourier transform technique. Based on the newly derived formula, we also examine, in theory, many basic properties of the option price under the current model. On the other hand, for practical purposes, we impose a reliable implementation method for the current formula so that it can be easily used in the trading market. With the numerical results, the impact of different parameters on the option price are also investigated.

Taylor's law (TL) originated as an empirical pattern in ecology. In many sets of samples of population density, the variance of each sample was approximately proportional to a power of the mean of that sample. In a family of nonnegative random variables, TL asserts that the population variance is proportional to a power of the population mean. TL, sometimes called fluctuation scaling, holds widely in physics, ecology, finance, demography, epidemiology, and other sciences, and characterizes many classical probability distributions and stochastic processes such as branching processes and birth-and-death processes. We demonstrate analytically for the first time that a version of TL holds for a class of distributions with infinite mean. These distributions, a subset of stable laws, and the associated TL differ qualitatively from those of light-tailed distributions. Our results employ and contribute to the methodology of Albrecher and Teugels (2006) and Albrecher et al. (2010). This work opens a new domain of investigation for generalizations of TL.

We consider a class of Sevastyanov branching processes with nonhomogeneous Poisson immigration. These processes relax the assumption required by the Bellman–Harris process which imposes the lifespan and offspring of each individual to be independent. They find applications in studies of the dynamics of cell populations. In this paper we focus on the subcritical case and examine asymptotic properties of the process. We establish limit theorems, which generalize classical results due to Sevastyanov and others. Our key findings include a novel law of large numbers and a central limit theorem which emerge from the nonhomogeneity of the immigration process.

Importance sampling has become an important tool for the computation of extreme quantiles and tail-based risk measures. For estimation of such nonlinear functionals of the underlying distribution, the standard efficiency analysis is not necessarily applicable. In this paper we therefore study importance sampling algorithms by considering moderate deviations of the associated weighted empirical processes. Using a delta method for large deviations, combined with classical large deviation techniques, the moderate deviation principle is obtained for importance sampling estimators of two of the most common risk measures: value at risk and expected shortfall.

Model and parameter uncertainties are common whenever some parametric model is selected to value a derivative instrument. Combining the Monte Carlo method with the Smolyak interpolation algorithm, we propose an accurate efficient numerical procedure to quantify the uncertainty embedded in complex derivatives. Except for the value function being sufficiently smooth with respect to the model parameters, there are no requirements on the payoff or candidate models. Numerical tests carried out quantify the uncertainty of Bermudan put options and down-and-out put options under the Heston model, with each model parameter specified in an interval.

Amethod for non-rigid image registration that is suitable for large deformations is presented. Conventional registration methods embed the image in a B-spline object, and the image is evolved by deforming the B-spline object. In this work, we represent the image using B-spline and deform the image using a composition approach. We also derive a computationally efficient algorithm for calculating the B-spline coefficients and gradients of the image by adopting ideas from signal processing using image filters. We demonstrate the application of our method on several different types of 2D and 3D images and compare it with existing methods.

We discuss modelling and simulation of volumetric rainfall in a catchment of the Murray–Darling Basin – an important food production region in Australia that was seriously affected by a recent prolonged drought. Consequently, there has been sustained interest in development of improved water management policies. In order to model accumulated volumetric catchment rainfall over a fixed time period, it is necessary to sum weighted rainfall depths at representative sites within each sub-catchment. Since sub-catchment rainfall may be highly correlated, the use of a Gamma distribution to model rainfall at each site means that catchment rainfall is expressed as a sum of correlated Gamma random variables. We compare four different models and conclude that a joint probability distribution for catchment rainfall constructed by using a copula of maximum entropy is the most effective.

In this paper we study a discrete-time optimal switching problem on a finite horizon. The underlying model has a running reward, terminal reward, and signed (positive and negative) switching costs. Using optimal stopping theory for discrete-parameter stochastic processes, we extend a well-known explicit dynamic programming method for computing the value function and the optimal strategy to the case of signed switching costs.

The advantages and limitations of frequency domain and time domain methods for estimating the interannual variability arising from day-to-day weather events are summarized. A modification of the time domain method is developed and its application in examining a precondition for the frequency domain method is demonstrated. A combined estimation procedure is proposed: it takes advantage of the strengths of both methods. The estimation procedures are tested with sets of synthetic data and are applied to long time series of three meteorological parameters. The impacts of the different methods on tests of potential long-range predictability for seasonal means are also discussed.

The dynamics of host-macroparasite infections pose considerable challenges for stochastic modelling because of the need to take into account a large number of relevant factors and many nonlinear interactions between them. This paper focuses attention on the infection transmission process and the effects of specific modelling assumptions about the mechanisms involved. Some dramatically simplified linear models are considered; they are based on multidimensional linear birth and death processes, and are designed to illuminate qualitative effects of interest. Both single and compound infections are allowed. It is shown that such simple models can generate and increase dispersion of parasite counts, even among homogeneous hosts.

The paper considers one of the standard processes for modeling returns in finance, the stochastic volatility process with regularly varying innovations. The aim of the paper is to show how point process techniques can be used to derive the asymptotic behavior of the sample autocorrelation function of this process with heavy-tailed marginal distributions. Unlike other non-linear models used in finance, such as GARCH and bilinear models, sample autocorrelations of a stochastic volatility process have attractive asymptotic properties. Specifically, in the infinite variance case, the sample autocorrelation function converges to zero in probability at a rate that is faster the heavier the tails of the marginal distribution. This behavior is analogous to the asymptotic behavior of the sample autocorrelations of independent identically distributed random variables.

This paper considers how students might be taught statistical consulting. Desirable characteristics of a good consultant are described. While skill in statistical techniques is considered, it is just one of several areas of expertise that are discussed. Suggestions are made as to what should be included in a consulting course to equip a student with these characteristics.

A general class of Markovian non-Gaussian bifurcating models for cell lineage data is presented. Examples include bifurcating autoregression, random coefficient autoregression, bivariate exponential, bivariate gamma, and bivariate Poisson models. Quasi-likelihood estimation for the model parameters and large-sample properties of the estimates are discussed.

Using the kernel representation of a continuous-time Lévy-driven ARMA (autoregressive moving average) process, we extend the class of nonnegative Lévy-driven Ornstein–Uhlenbeck processes employed by Barndorff-Nielsen and Shephard (2001) to allow for nonmonotone autocovariance functions. We also consider a class of fractionally integrated Lévy-driven continuous-time ARMA processes obtained by a simple modification of the kernel of the continuous-time ARMA process. Asymptotic properties of the kernel and of the autocovariance function are derived.