Curve crossing for the reflected Levy process at zero and infinity

Let $R_{t}=\sup_{0\leq s\leq t}X_{s}-X_{t}$ be a Levy process reflected in its maximum. We give necessary and sufficient conditions for finiteness of passage times above power law boundaries at infinity. Information as to when the expected passage time for $R_{t}$ is finite, is given. We also discuss the almost sure finiteness of $\limsup_{t\to 0}R_{t}/t^{\kappa}$, for each $\kappa\geq 0$.

The unsteady flow of a weakly compressible fluid in a thin porous layer III: three-dimensional computations

We describe a novel method for determining the pressure and velocity fields for a weakly compressible fluid flowing in a thin three-dimensional layer composed of an inhomogeneous, anisotropic porous medium, with vertical side walls and variable upper and lower boundaries, in the presence of vertical wells injecting and/or extracting fluid. Our approach uses the method of matched asymptotic expansions to derive expressions for all significant process quantities, the computation of which requires only the solution of linear, elliptic, two-dimensional boundary value and eigenvalue problems. In this article, we provide full implementation details and present numerical results demonstrating the efficiency and accuracy of our scheme.

The Carnegie curve

The Earth’s fair weather atmospheric electric field shows, in clean air, an average daily variation which follows universal time, globally independent of the measurement position. This single diurnal cycle variation (maximum around 19UT and minimum around 03UT) is widely known as the Carnegie curve, after the geophysical survey vessel of the Carnegie Institution of Washington on which the original measurement campaigns demonstrating the universal time variation were undertaken. The Carnegie curve’s enduring importance is in providing a reference variation against which atmospheric electricity measurements are still compared; it is believed to originate from regular daily variations in atmospheric electrification associated with the different global disturbed weather regions. Details of the instrumentation, measurement principles and data obtained on the Carnegie’s seventh and final cruise are reviewed here, also deriving new harmonic coefficients allowing calculation of the Carnegie curve for different seasons. The additional harmonic analysis now identifies changes in the phasing of the maximum and minimum in the Carnegie curve, which shows a systematic seasonal variation, linked to the solstices and equinoxes, respectively.

The M curve and the performance of Spanish international new ventures

In this paper we address three challenges. First, we discuss how international new ventures (INVs) are probably not explained by the Uppsala model as there is no time for learning about foreign markets in newly born and small firms. Only in the longer term can INVs develop experiential learning to overcome the liability of foreignness as they expand abroad. Second, we advance theoretically on previous research demonstrating that the multinationality−performance relationship of INVs follows a traditional S-shaped relationship, but they first experience a ‘born global illusion’ which leads to a non-traditional M curve. Third, using a panel data analysis for the period 1994–2008 we find empirically that Spanish INVs follow an inverted U curve in the very short term, where no learning takes place, but that experience gained over time yields an M-curve relationship once learning takes place.

A comparative review of dimension reduction methods in approximate Bayesian computation

Approximate Bayesian computation (ABC) methods make use of comparisons between simulated and observed summary statistics to overcome the problem of computationally intractable likelihood functions. As the practical implementation of ABC requires computations based on vectors of summary statistics, rather than full data sets, a central question is how to derive low-dimensional summary statistics from the observed data with minimal loss of information. In this article we provide a comprehensive review and comparison of the performance of the principal methods of dimension reduction proposed in the ABC literature. The methods are split into three nonmutually exclusive classes consisting of best subset selection methods, projection techniques and regularization. In addition, we introduce two new methods of dimension reduction. The first is a best subset selection method based on Akaike and Bayesian information criteria, and the second uses ridge regression as a regularization procedure. We illustrate the performance of these dimension reduction techniques through the analysis of three challenging models and data sets.

Constructing summary statistics for approximate Bayesian computation: semi-automatic approximate Bayesian computation

Many modern statistical applications involve inference for complex stochastic models, where it is easy to simulate from the models, but impossible to calculate likelihoods. Approximate Bayesian computation (ABC) is a method of inference for such models. It replaces calculation of the likelihood by a step which involves simulating artificial data for different parameter values, and comparing summary statistics of the simulated data with summary statistics of the observed data. Here we show how to construct appropriate summary statistics for ABC in a semi-automatic manner. We aim for summary statistics which will enable inference about certain parameters of interest to be as accurate as possible. Theoretical results show that optimal summary statistics are the posterior means of the parameters. Although these cannot be calculated analytically, we use an extra stage of simulation to estimate how the posterior means vary as a function of the data; and we then use these estimates of our summary statistics within ABC. Empirical results show that our approach is a robust method for choosing summary statistics that can result in substantially more accurate ABC analyses than the ad hoc choices of summary statistics that have been proposed in the literature. We also demonstrate advantages over two alternative methods of simulation-based inference.

A finite element method for nonlinear elliptic problems

We present a Galerkin method with piecewise polynomial continuous elements for fully nonlinear elliptic equations. A key tool is the discretization proposed in Lakkis and Pryer, 2011, allowing us to work directly on the strong form of a linear PDE. An added benefit to making use of this discretization method is that a recovered (finite element) Hessian is a byproduct of the solution process. We build on the linear method and ultimately construct two different methodologies for the solution of second order fully nonlinear PDEs. Benchmark numerical results illustrate the convergence properties of the scheme for some test problems as well as the Monge–Amp`ere equation and the Pucci equation.

A finite element method for second order nonvariational elliptic problems

We propose a numerical method to approximate the solution of second order elliptic problems in nonvariational form. The method is of Galerkin type using conforming finite elements and applied directly to the nonvariational (nondivergence) form of a second order linear elliptic problem. The key tools are an appropriate concept of “finite element Hessian” and a Schur complement approach to solving the resulting linear algebra problem. The method is illustrated with computational experiments on three linear and one quasi-linear PDE, all in nonvariational form.

Continuum sea ice rheology determined from subcontinuum mechanics

[1] A method is presented to calculate the continuum-scale sea ice stress as an imposed, continuum-scale strain-rate is varied. The continuum-scale stress is calculated as the area-average of the stresses within the floes and leads in a region (the continuum element). The continuum-scale stress depends upon: the imposed strain rate; the subcontinuum scale, material rheology of sea ice; the chosen configuration of sea ice floes and leads; and a prescribed rule for determining the motion of the floes in response to the continuum-scale strain-rate. We calculated plastic yield curves and flow rules associated with subcontinuum scale, material sea ice rheologies with elliptic, linear and modified Coulombic elliptic plastic yield curves, and with square, diamond and irregular, convex polygon-shaped floes. For the case of a tiling of square floes, only for particular orientations of the leads have the principal axes of strain rate and calculated continuum-scale sea ice stress aligned, and these have been investigated analytically. The ensemble average of calculated sea ice stress for square floes with uniform orientation with respect to the principal axes of strain rate yielded alignment of average stress and strain-rate principal axes and an isotropic, continuum-scale sea ice rheology. We present a lemon-shaped yield curve with normal flow rule, derived from ensemble averages of sea ice stress, suitable for direct inclusion into the current generation of sea ice models. This continuum-scale sea ice rheology directly relates the size (strength) of the continuum-scale yield curve to the material compressive strength.

Asymmetric output gap effects in Phillips Curve and mark-up pricing models: evidence for the US and the UK

A number of studies have found an asymmetric response of consumer price index inflation to the output gap in the US in simple Phillips curve models. We consider whether there are similar asymmetries in mark-up pricing models, that is, whether the mark-up over producers' costs also depends upon the sign of the (adjusted) output gap. The robustness of our findings to the price series is assessed, and also whether price-output responses in the UK are asymmetric.

Dependence of sea ice yield-curve shape on ice thickness

In this note, the authors discuss the contribution that frictional sliding of ice floes (or floe aggregates) past each other and pressure ridging make to the plastic yield curve of sea ice. Using results from a previous study that explicitly modeled the amount of sliding and ridging that occurs for a given global strain rate, it is noted that the relative contribution of sliding and ridging to ice stress depends upon ice thickness. The implication is that the shape and size of the plastic yield curve is dependent upon ice thickness. The yield-curve shape dependence is in addition to plastic hardening/weakening that relates the size of the yield curve to ice thickness. In most sea ice dynamics models the yield-curve shape is taken to be independent of ice thickness. The authors show that the change of the yield curve due to a change in the ice thickness can be taken into account by a weighted sum of two thickness-independent rheologies describing ridging and sliding effects separately. It would be straightforward to implement the thickness-dependent yield-curve shape described here into sea ice models used for global or regional ice prediction.

Abstract platforms of computation

Computational formalisms have been pushing the boundaries of the field of computing for the last 80 years and much debate has surrounded what computing entails; what it is, and what it is not. This paper seeks to explore the boundaries of the ideas of computation and provide a framework for enabling a constructive discussion of computational ideas. First, a review of computing is given, ranging from Turing Machines to interactive computing. Then, a variety of natural physical systems are considered for their computational qualities. From this exploration, a framework is presented under which all dynamical systems can be considered as instances of the class of abstract computational platforms. An abstract computational platform is defined by both its intrinsic dynamics and how it allows computation that is meaningful to an external agent through the configuration of constraints upon those dynamics. It is asserted that a platform’s computational expressiveness is directly related to the freedom with which constraints can be placed. Finally, the requirements for a formal constraint description language are considered and it is proposed that Abstract State Machines may provide a reasonable basis for such a language.

What will be the risk-free rate and benchmark yield curve following European monetary union?

Using a linear factor model, we study the behaviour of French, Germany, Italian and British sovereign yield curves in the run up to EMU. This allows us to determine which of these yield curves might best approximate a benchmark yield curve post EMU. We find that the best approximation for the risk free yield is the UK three month T-bill yield, followed by the German three month T-bill yield. As no one sovereign yield curve dominates all others, we find that a composite yield curve, consisting of French, Italian and UK bonds at different maturity points along the yield curve should be the benchmark post EMU.

The elliptic sine-Gordon equation: a nonlinear elliptic integrable PDE

In this paper, we summarise this recent progress to underline the features specific to this nonlinear elliptic case, and we give a new classification of boundary conditions on the semistrip that satisfy a necessary condition for yielding a boundary value problem can be effectively linearised. This classification is based on formulation the equation in terms of an alternative Lax pair.