On- and offline detection of structural breaks in thermal spraying processes

We investigate and develop methods for structural break detection, considering time series from thermal spraying process monitoring. Since engineers induce technical malfunctions during the processes, the time series exhibit structural breaks at known time points, giving us valuable information to conduct the investigations. First, we consider a recently developed robust online (also real-time) filtering (i.e. smoothing) procedure that comprises a test for local linearity. This test rejects when jumps and trend changes are present, so that it can also be useful to detect such structural breaks online. Second, based on the filtering procedure we develop a robust method for the online detection of ongoing trends. We investigate these two methods as to the online detection of structural breaks by simulations and applications to the time series from the manipulated spraying processes. Third, we consider a recently developed fluctuation test for constant variances that can be applied offline, i.e. after the whole time series has been observed, to control the spraying results. Since this test is not reliable when jumps are present in the time series, we suggest data transformation based on filtering and demonstrate that this transformation makes the test applicable.     

Selection procedures for sparse data

The problem of selecting the best of k populations is studied for data which are incomplete as some of the values have been deleted randomly. This situation is met in extreme value analysis where only data exceeding a threshold are observable. For increasing sample size we study the case where the probability that a value is observed tends to zero, but the sparse condition is satisfied, so that the mean number of observable values in each population is bounded away from zero and infinity as the sample size tends to infinity. The incomplete data are described by thinned point processes which are approximated by Poisson point processes. Under weak assumptions and after suitable transformations these processes converge to a Poisson point process. Optimal selection rules for the limit model are used to construct asymptotically optimal selection rules for the original sequence of models. The results are applied to extreme value data for high thresholds data.     

Efficient and flexible model-based clustering of jumps in diffusion processes

Jump–diffusion processes involving diffusion processes with discontinuous movements, called jumps, are widely used to model time-series data that commonly exhibit discontinuity in their sample paths. The existing jump–diffusion models have been recently extended to multivariate time-series data. The models are, however, still limited by a single parametric jump-size distribution that is common across different subjects. Such strong parametric assumptions for the shape and structure of a jump-size distribution may be too restrictive and unrealistic for multiple subjects with different characteristics. This paper thus proposes an efficient Bayesian nonparametric method to flexibly model a jump-size distribution while borrowing information across subjects in a clustering procedure using a nested Dirichlet process. For efficient posterior computation, a partially collapsed Gibbs sampler is devised to fit the proposed model. The proposed methodology is illustrated through a simulation study and an application to daily stock price data for companies in the S&P 100 index from June 2007 to June 2017.     

Analyses of the Markov modulated fluid flow with one-sided ph-type jumps using coupled queues and the completed graphs

In this paper, we analyze Markov modulated fluid flow processes with one-sided ph-type jumps using the completed graph and also through the limits of coupled queueing processes to be constructed. For the models, we derive various results on time-dependent distributions and distributions of first passage times, and present the Riccati equations that transform matrices of the first return times to 0 satisfy. The Riccati equations enable us to compute the transform matrices using Newton’s method which is known very fast and stable. Finally, we present some duality results between the model with ph-type downward jumps and the model with ph-type upward jumps. This paper contains extended results of Ahn (2009) and probabilistic interpretations given by the completed graphs.     

Contrast function estimation for the drift parameter of ergodic jump diffusion process

In this paper, we consider an ergodic diffusion process with jumps whose drift coefficient depends on an unknown parameter. We suppose that the process is discretely observed. We introduce an estimator based on a contrast function, which is efficient without requiring any conditions on the rate at which the step discretization goes to zero, and where we allow the observed process to have nonsummable jumps. This extends earlier results where the condition on the step discretization was needed and where the process was supposed to have summable jumps. In general situations, our contrast function is not explicit and one has to resort to some approximation. In the case of a finite jump activity, we propose explicit approximations of the contrast function such that the efficient estimation of the drift parameter is feasible. This extends the results obtained by Kessler in the case of continuous processes.     

Central Limit Theorems of Local Polynomial Threshold Estimator for Diffusion Processes with Jumps

Central limit theorems play an important role in the study of statistical inference for stochastic processes. However, when the non‐parametric local polynomial threshold estimator, especially local linear case, is employed to estimate the diffusion coefficients of diffusion processes, the adaptive and predictable structure of the estimator conditionally on the σ ‐field generated by diffusion processes is destroyed, so the classical central limit theorem for martingale difference sequences cannot work. In high‐frequency data, we proved the central limit theorems of local polynomial threshold estimators for the volatility function in diffusion processes with jumps by Jacod's stable convergence theorem. We believe that our proof procedure for local polynomial threshold estimators provides a new method in this field, especially in the local linear case.     

Integer autoregressive models with structural breaks

Even though integer-valued time series are common in practice, the methods for their analysis have been developed only in recent past. Several models for stationary processes with discrete marginal distributions have been proposed in the literature. Such processes assume the parameters of the model to remain constant throughout the time period. However, this need not be true in practice. In this paper, we introduce non-stationary integer-valued autoregressive (INAR) models with structural breaks to model a situation, where the parameters of the INAR process do not remain constant over time. Such models are useful while modelling count data time series with structural breaks. The Bayesian and Markov Chain Monte Carlo (MCMC) procedures for the estimation of the parameters and break points of such models are discussed. We illustrate the model and estimation procedure with the help of a simulation study. The proposed model is applied to the two real biometrical data sets.     

A blocked Gibbs sampler for NGG-mixture models via a priori truncation

We define a new class of random probability measures, approximating the well-known normalized generalized gamma (NGG) process. Our new process is defined from the representation of NGG processes as discrete measures where the weights are obtained by normalization of the jumps of Poisson processes and the support consists of independent identically distributed location points, however considering only jumps larger than a threshold \(\varepsilon \). Therefore, the number of jumps of the new process, called \(\varepsilon \)-NGG process, is a.s. finite. A prior distribution for \(\varepsilon \) can be elicited. We assume such a process as the mixing measure in a mixture model for density and cluster estimation, and build an efficient Gibbs sampler scheme to simulate from the posterior. Finally, we discuss applications and performance of the model to two popular datasets, as well as comparison with competitor algorithms, the slice sampler and a posteriori truncation.     

Basket CDS pricing with default intensities using a regime-switching shot-noise model

In this paper, we employ an intensity-based credit risk model with regime-switching to study the valuation of basket CDS in a homogeneous portfolio. We assume that the default intensities are described by some dependent regime-switching shot-noise processes and the individual jumps of the intensity are driven by a common factor. By using the conditional Laplace transform of the regime-switching shot-noise process, we obtain the closed form results for pricing the fair spreads of the basket CDS. We present some numerical examples to illustrate the effect of the model parameters on the fair spreads.     

Discontinuities in robust nonparametric regression with α-mixing dependence

The main idea of the paper is to introduce a robust regression estimation method under an α-mixing dependence assumption, staying free of any parametric model restrictions while also allowing for some sudden changes in the unknown regression function. The sudden changes in the model may correspond to discontinuity points (jumps) or higher order breaks (jumps in corresponding derivatives) as well. We firstly derive some important statistical properties for local polynomial M-smoother estimates and we will propose a statistical test to decide whether some given point of interest is significantly important for a change to occur or not. As the asymptotic distribution of the test statistic depends on quantities which are left unknown we also introduce a bootstrap algorithm which can be used to mimic the target distribution of interest. All necessary proofs are provided together with some experimental results from a simulation study and a real data example.     

Testing the Martingale Hypothesis

We propose new tests of the martingale hypothesis based on generalized versions of the Kolmogorov–Smirnov and Cramér–von Mises tests. The tests are distribution-free and allow for a weak drift in the null model. The methods do not require either smoothing parameters or bootstrap resampling for their implementation and so are well suited to practical work. The article develops limit theory for the tests under the null and shows that the tests are consistent against a wide class of nonlinear, nonmartingale processes. Simulations show that the tests have good finite sample properties in comparison with other tests particularly under conditional heteroscedasticity and mildly explosive alternatives. An empirical application to major exchange rate data finds strong evidence in favor of the martingale hypothesis, confirming much earlier research.     

A Stochastic Volatility Model With Conditional Skewness

We develop a discrete-time affine stochastic volatility model with time-varying conditional skewness (SVS). Importantly, we disentangle the dynamics of conditional volatility and conditional skewness in a coherent way. Our approach allows current asset returns to be asymmetric conditional on current factors and past information, which we term contemporaneous asymmetry. Conditional skewness is an explicit combination of the conditional leverage effect and contemporaneous asymmetry. We derive analytical formulas for various return moments that are used for generalized method of moments (GMM) estimation. Applying our approach to S&P500 index daily returns and option data, we show that one- and two-factor SVS models provide a better fit for both the historical and the risk-neutral distribution of returns, compared to existing affine generalized autoregressive conditional heteroscedasticity (GARCH), and stochastic volatility with jumps (SVJ) models. Our results are not due to an overparameterization of the model: the one-factor SVS models have the same number of parameters as their one-factor GARCH competitors and less than the SVJ benchmark.     

Semiparametric jump-preserving estimation for single-index models

Estimation of the single-index model with a discontinuous unknown link function is considered in this paper. Existed refined minimum average variance estimation (rMAVE) method can estimate the single-index parameter and unknown link function simultaneously by minimising the average pointwise conditional variance, where the conditional variance can be estimated using the local linear fit method with centred kernel function. When there are jumps in the link function, big biases around jumps can appear. For this reason, we embed the jump-preserving technique in the rMAVE method, then propose an adaptive jump-preserving estimation procedure for the single-index model. Concretely speaking, the conditional variance is obtained by the one among local linear fits with centred, left-sided and right-sided kernel functions who has minimum weighted residual mean squares. The resulting estimators can preserve the jumps well and also give smooth estimates of the continuity parts. Asymptotic properties are established under some mild conditions. Simulations and real data analysis show the proposed method works well.     

Lévy Processes,Phase-Type Distributions,and Martingales

Lévy processes are defined as processes with stationary independent increments and have become increasingly popular as models in queueing, finance, etc.; apart from Brownian motion and compound Poisson processes, some popular examples are stable processes, variance gamma processes, CGMY Lévy processes (tempered stable processes), NIG (normal inverse Gaussian) Lévy processes, and hyperbolic Lévy processes. We consider here a dense class of Lévy processes, compound Poisson processes with phase-type jumps in both directions and an added Brownian component. Within this class, we survey how to explicitly compute a number of quantities that are traditionally studied in the area of Lévy processes, in particular two-sided exit probabilities and associated Laplace transforms, the closely related scale function, one-sided exit probabilities and associated Laplace transforms coming up in queueing problems, and similar quantities for a Lévy process with reflection in 0. The solutions are in terms of roots to polynomials, and the basic equations are derived by purely probabilistic arguments using martingale optional stopping; a particularly useful martingale is the so-called Kella-Whitt martingale. Also, the relation to fluid models with a Brownian component is discussed.     

Forecasting energy futures volatility with threshold augmented heterogeneous autoregressive jump models

Abstract

This study forecasts the volatility of two energy futures markets (oil and gas), using high-frequency data. We, first, disentangle volatility into continuous volatility and jumps. Second, we apply wavelet analysis to study the relationship between volume and the volatility measures for different horizons. Third, we augment the heterogeneous autoregressive (HAR) model by nonlinearly including both jumps and volume. We then propose different empirical extensions of the HAR model. Our study shows that oil and gas volatilities nonlinearly depend on public information (jumps), private information (continuous volatility), and trading volume. Moreover, our threshold augmented HAR model with heterogeneous jumps and continuous volatility outperforms HAR model in forecasting volatility.     

Identification of average marginal effects under misspecification when covariates are normal

A previously known result in the econometrics literature is that when covariates of an underlying data generating process are jointly normally distributed, estimates from a nonlinear model that is misspecified as linear can be interpreted as average marginal effects. This has been shown for models with exogenous covariates and separability between covariates and errors. In this paper, we extend this identification result to a variety of more general cases, in particular for combinations of separable and nonseparable models under both exogeneity and endogeneity. So long as the underlying model belongs to one of these large classes of data generating processes, our results show that nothing else must be known about the true DGP—beyond normality of observable data, a testable assumption—in order for linear estimators to be interpretable as average marginal effects. We use simulation to explore the performance of these estimators using a misspecified linear model and show they perform well when the data are normal but can perform poorly when this is not the case.     

Asymptotic properties for multipower variation of semimartingales and Gaussian integral processes with jumps

This paper presents limit theorems of realized multipower variation for semimartingales and Gaussian integral processes with jumps observed in high frequency. In particular, we obtain a central limit theorem of realized multipower variation for semimartingale where some of the powers equal one and the others are less one. Convergence in probability and central limit theorems of realized threshold bipower variation for Gaussian integral processes with jumps are also obtained. These results provide new statistical tools to analyze and test the long memory effect in high frequency situation.     

Measuring serial dependence in categorical time series

The analysis of time-indexed categorical data is important in many fields, e.g., in telecommunication network monitoring, manufacturing process control, ecology, etc. Primary interest is in detecting and measuring serial associations and dependencies in such data. For cardinal time series analysis, autocorrelation is a convenient and informative measure of serial association. Yet, for categorical time series analysis an analogous convenient measure and corresponding concepts of weak stationarity have not been provided. For two categorical variables, several ways of measuring association have been suggested. This paper reviews such measures and investigates their properties in a serial context. We discuss concepts of weak stationarity of a categorical time series, in particular of stationarity in association measures. Serial association and weak stationarity are studied in the class of discrete ARMA processes introduced by Jacobs and Lewis (J. Time Ser. Anal. 4(1):19–36, 1983). An intrinsic feature of a time series is that, typically, adjacent observations are dependent. The nature of this dependence among observations of a time series is of considerable practical interest. Time series analysis is concerned with techniques for the analysis of this dependence. (Box et al. 1994p. 1)     

Inference for cluster point processes with over- or under-dispersed cluster sizes

Cluster point processes comprise a class of models that have been used for a wide range of applications. While several models have been studied for the probability density function of the offspring displacements and the parent point process, there are few examples of non-Poisson distributed cluster sizes. In this paper, we introduce a generalization of the Thomas process, which allows for the cluster sizes to have a variance that is greater or less than the expected value. We refer to this as the cluster sizes being over- and under-dispersed, respectively. To fit the model, we introduce minimum contrast methods and a Bayesian MCMC algorithm. These are evaluated in a simulation study. It is found that using the Bayesian MCMC method, we are in most cases able to detect over- and under-dispersion in the cluster sizes. We use the MCMC method to fit the model to nerve fiber data, and contrast the results to those of a fitted Thomas process.