Comparing k Cumulative Incidence Functions Through Resampling Methods

Tests for the equality of k cumulative incidence functions in a competing risks model are proposed. Test statistics are based on a vector of processes related to the cumulative incidence functions. Since their asymptotic distributions appear very complicated and depend on the underlying distribution of the data, two resampling techniques, namely the well-known bootstrap method and the so-called random symmetrization method, are used to approximate the critical values of the tests. Without making any assumptions on the nature of dependence between the risks, the tests allow one to compare k risks simultaneously for k 2 under the random censorship model. Tests against ordered alternatives are also considered. Simulation studies indicate that the proposed tests perform very well with moderate sample size. A real application to cancer mortality data is given.     

The Distribution of the Record Position and its Applications

This is mainly an expository article on the positions of records in sequences of ordered elements. Such sequences are obtained, for example, when observing and ordering continuous iid random variables. In practice, records are of interest, for example, in meteorology and sports. A k-record is obtained when a new element is placed at position k counted from the top. The sequences of time points, when new k-records occur, are studied by elementary random walk methods. In the last section, it is shown that the time scale can be changed so that the time points of the k-records follow, approximately, a Poisson process.     

Estimation of the Jump Size Density in a Mixed Compound Poisson Process

In this paper, we consider a mixed compound Poisson process, that is, a random sum of independent and identically distributed (i.i.d.) random variables where the number of terms is a Poisson process with random intensity. We study nonparametric estimators of the jump density by specific deconvolution methods. Firstly, assuming that the random intensity has exponential distribution with unknown expectation, we propose two types of estimators based on the observation of an i.i.d. sample. Risks bounds and adaptive procedures are provided. Then, with no assumption on the distribution of the random intensity, we propose two non‐parametric estimators of the jump density based on the joint observation of the number of jumps and the random sum of jumps. Risks bounds are provided, leading to unusual rates for one of the two estimators. The methods are implemented and compared via simulations.     

Non-ergodic martingale estimating functions and related asymptotics

Well-known estimation methods such as conditional least squares, quasilikelihood and maximum likelihood (ML) can be unified via a single framework of martingale estimating functions (MEFs). Asymptotic distributions of estimates for ergodic processes use constant norm (e.g. square root of the sample size) for asymptotic normality. For certain non-ergodic-type applications, however, such as explosive autoregression and super-critical branching processes, one needs a random norm in order to get normal limit distributions. In this paper, we are concerned with non-ergodic processes and investigate limit distributions for a broad class of MEFs. Asymptotic optimality (within a certain class of non-ergodic MEFs) of the ML estimate is deduced via establishing a convolution theorem using a random norm. Applications to non-ergodic autoregressive processes, generalized autoregressive conditional heteroscedastic-type processes, and super-critical branching processes are discussed. Asymptotic optimality in terms of the maximum random limiting power regarding large sample tests is briefly discussed.     

Some novel results on pairwise quasi-asymptotical independence with applications to risk theory

In this article we obtain some novel results on pairwise quasi-asymptotically independent (pQAI) random variables. Concretely speaking, let X₁, …, X_n be n real-valued pQAI random variables, and W₁, …, W_n be another n non negative and arbitrarily dependent random variables, but independent of X₁, …, X_n. Under some mild conditions, we prove that W₁X₁, …, W_nX_n are still pQAI as well. Our result is in a general setting whether the primary random variables X₁, …, X_n are heavy-tailed or not. Finally, a special case of above result is applied to risk theory for investigating the finite-time ruin probability for a discrete-time risk model with a wide type of dependence structure.     

Consistency of the Subsample Bootstrap empirical process

In the classical Bootstrap approach the number of distinct observation in the resample is random. To overcome this hitch Rao et al. [Bootstrap by sequential resampling, J. Statist. Plan. Inference 64 (1997), pp. 257–281] have proposed a modified resampling procedure – the so-called Sequential Bootstrap or 0.632-Bootstrap – in which each resample has exactly the same number m ～eq ?0.632 n? of distinct observations. Motivated by this idea we introduce an akin procedure, the Subsample Bootstrap, where additionally even the size of each resample is equal. It will turn out that the Subsample Bootstrap empirical process is consistent for a wide class of Donsker classes.     

Mixing and moments properties of a non-stationary copula-based Markov process

Abstract

We provide conditions under which a non-stationary copula-based Markov process is geometric β-mixing and geometric ρ-mixing. Our results generalize some results of Beare who considers the stationary case. As a particular case we introduce a stochastic process, that we call convolution-based Markov process, whose construction is obtained by using the C-convolution operator which allows the increments to be dependent. Within this subclass of processes we characterize a modified version of the standard random walk where copulas and marginal distributions involved are in the same elliptical family. We study mixing and moments properties to identify the differences compared to the standard case.     

Option pricing under the Merton model of the short rate in subdiffusive Brownian motion regime

In this paper, we propose an extension of the Merton short rate model, which reflects the subdiffusive nature of the short rate dynamics. The subdiffusive property is manifested by the random periods of time, during which the asset price does not change. We derive explicit formulas for European call and put options and present some simulation results for the case of α stable. Moreover, we discuss the implied volatility of this model.     

The covariance of the backward and forward recurrence times in a renewal process: the stationary case and asymptotics for the ordinary case

In a Poisson process, it is well-known that the forward and backward recurrence times at a given time point t are independent random variables. In a renewal process, although the joint distribution of these quantities is known (asymptotically), it seems that very few results regarding their covariance function exist. In the present paper, we study this covariance and, in particular, we state both necessary and sufficient conditions for it to be positive, zero or negative in terms of reliability classifications and the coefficient of variation of the underlying inter-renewal and the associated equilibrium distribution. Our results apply either for an ordinary renewal process in the steady state or for a stationary process.     

Modified Normal-based Approximation to the Percentiles of Linear Combination of Independent Random Variables with Applications

A modified normal-based approximation for calculating the percentiles of a linear combination of independent random variables is proposed. This approximation is applicable in situations where expectations and percentiles of the individual random variables can be readily obtained. The merits of the approximation are evaluated for the chi-square and beta distributions using Monte Carlo simulation. An approximation to the percentiles of the ratio of two independent random variables is also given. Solutions based on the approximations are given for some classical problems such as interval estimation of the normal coefficient of variation, survival probability, the difference between or the ratio of two binomial proportions, and for some other problems. Furthermore, approximation to the percentiles of a doubly noncentral F distribution is also given. For all the problems considered, the approximation provides simple satisfactory solutions. Two examples are given to show applications of the approximation.     

A Bayesian approach to zero-inflated data in extremes

Abstract

The generalized extreme value (GEV) distribution is known as the limiting result for the modeling of maxima blocks of size n, which is used in the modeling of extreme events. However, it is possible for the data to present an excessive number of zeros when dealing with extreme data, making it difficult to analyze and estimate these events by using the usual GEV distribution. The Zero-Inflated Distribution (ZID) is widely known in literature for modeling data with inflated zeros, where the inflator parameter w is inserted. The present work aims to create a new approach to analyze zero-inflated extreme values, that will be applied in data of monthly maximum precipitation, that can occur during months where there was no precipitation, being these computed as zero. An inference was made on the Bayesian paradigm, and the parameter estimation was made by numerical approximations of the posterior distribution using Markov Chain Monte Carlo (MCMC) methods. Time series of some cities in the northeastern region of Brazil were analyzed, some of them with predominance of non-rainy months. The results of these applications showed the need to use this approach to obtain more accurate and with better adjustment measures results when compared to the standard distribution of extreme value analysis.     

Modelling columnarity of pyramidal cells in the human cerebral cortex

For modelling the location of pyramidal cells in the human cerebral cortex, we suggest a hierarchical point process in that exhibits anisotropy in the form of cylinders extending along the z-axis. The model consists first of a generalised shot noise Cox process for the xy-coordinates, providing cylindrical clusters, and next of a Markov random field model for the z-coordinates conditioned on the xy-coordinates, providing either repulsion, aggregation or both within specified areas of interaction. Several cases of these hierarchical point processes are fitted to two pyramidal cell data sets, and of these a final model allowing for both repulsion and attraction between the points seem adequate. We discuss how the final model relates to the so-called minicolumn hypothesis in neuroscience.     

Bayesian inference and prediction of the Pareto distribution based on ordered ranked set sampling

In this paper, order statistics from independent and non identically distributed random variables is used to obtain ordered ranked set sampling (ORSS). Bayesian inference of unknown parameters under a squared error loss function of the Pareto distribution is determined. We compute the minimum posterior expected loss (the posterior risk) of the derived estimates and compare them with those based on the corresponding simple random sample (SRS) to assess the efficiency of the obtained estimates. Two-sample Bayesian prediction for future observations is introduced by using SRS and ORSS for one- and m-cycle. A simulation study and real data are applied to show the proposed results.     

Regression Diagnostics of the Semiparametric Proportional Rate Model for Irregularly Spaced Repeated Measurements

Counting process techniques have been successfully introduced to semiparametric inference of repeated measurements. Cheng and Wei (2000 Cheng , S. C. , Wei , L. J. ( 2000 ). Inference for a semiparametric model with panel data . Biometrika 87 : 89 – 97 .[Crossref], [Web of Science ®] , [Google Scholar]) proposed a simple inference procedure for the semiparametric proportional rate model, which reduces to relative risk regression models for binary data. While the baseline mean functions are completely unspecified, it still requires several assumptions for valid inference. In this article, a goodness-of-fit test for it is proposed based on cumulative residuals. Theoretical justification is provided and an illustration with a dataset from a clinical trial is given. Results of simulation studies to evaluate finite sample performance are also provided.     

An investigation of the spatial correlations for relative purchasing power in Baden–Württemberg

The relative purchasing power—i.e., the purchasing power per inhabitant—is one of the key characteristics for businesses deciding on site selection. Apart from that it also plays a major role in regional planning, pricing policy and market research. In this study we investigate the spatial correlations for relative purchasing power of the townships in Baden–Württemberg. In particular, changes in relative purchasing power are analysed for three different time intervals, 1987–1993, 1993–1998 and 1998–2004, by means of distance-dependent characteristics like the mark–correlation function, the Simpson indices α(r) and β(r) and by tests on random labelling. It is shown that there are positive correlations for small distances between different townships but that these positive correlations become weaker over the years until they are almost nonexistent (in the sense that hypotheses of random labelling are no longer rejected). A conclusion from this loss of spatial correlations with time is that the relative purchasing power might become more and more purely random. This means that the relative purchasing power in a township is less and less influenced by the relative purchasing power of townships nearby. We further analysed these changes in the Bodensee–Oberschwaben and Stuttgart regions to compare the development of the relative purchasing power in both urban and rural environments. This analysis was done in close cooperation with W. Brachat–Schwarz and W. Walla of the Statistical Office of Baden–Württemberg. S.E. is supported by a grant of the graduate college 1100 at Ulm University.     

Adaptive Deconvolution on the Non‐negative Real Line

In this paper, we consider the problem of adaptive density or survival function estimation in an additive model defined by Z=X+Y with X independent of Y, when both random variables are non‐negative. This model is relevant, for instance, in reliability fields where we are interested in the failure time of a certain material that cannot be isolated from the system it belongs. Our goal is to recover the distribution of X (density or survival function) through n observations of Z, assuming that the distribution of Y is known. This issue can be seen as the classical statistical problem of deconvolution that has been tackled in many cases using Fourier‐type approaches. Nonetheless, in the present case, the random variables have the particularity to be supported. Knowing that, we propose a new angle of attack by building a projection estimator with an appropriate Laguerre basis. We present upper bounds on the mean squared integrated risk of our density and survival function estimators. We then describe a non‐parametric data‐driven strategy for selecting a relevant projection space. The procedures are illustrated with simulated data and compared with the performances of a more classical deconvolution setting using a Fourier approach. Our procedure achieves faster convergence rates than Fourier methods for estimating these functions.     

Line transects of two-dimensional random fields: Estimation and design

In line-transect sampling data are collected continuously along the path travelled. These data are then used to estimate some characteristic of an entire region. We model the characteristic of the region as a linear functional of Z(), a realization of a random field, and our data as the values of Z(x) for x along the path. We study how to estimate the characteristic of interest and how to choose the path to accurately estimate it.     

Locking the Wiener Process to its Level-Crossing Time

We consider the specific transformation of a Wiener process {X(t), t ≥ 0} in the presence of an absorbing barrier a that results when this process is “time-locked” with respect to its first passage time T _a through a criterion level a, and the evolution of X(t) is considered backwards (retrospectively) from T _a . Formally, we study the random variables defined by Y(t) ≡ X(T _a  ? t) and derive explicit results for their density and mean, and also for their asymptotic forms. We discuss how our results can aid interpretations of time series “response-locked” to their times of crossing a criterion level.     

Non uniform exponential bounds on normal approximation by Stein’s method and monotone size bias couplings

It is known that the normal approximation is applicable for sums of non negative random variables, W, with the commonly employed couplings. In this work, we use the Stein’s method to obtain a general theorem of non uniform exponential bound on normal approximation base on monotone size bias couplings of W. Applications of the main result to give the bound on normal approximation for binomial random variable, the number of bulbs on at the terminal time in the lightbulb process, and the number of m runs are also provided.     

Correlation in a Bayesian framework

The authors consider the correlation between two arbitrary functions of the data and a parameter when the parameter is regarded as a random variable with given prior distribution. They show how to compute such a correlation and use closed form expressions to assess the dependence between parameters and various classical or robust estimators thereof, as well as between p‐values and posterior probabilities of the null hypothesis in the one‐sided testing problem. Other applications involve the Dirichlet process and stationary Gaussian processes. Using this approach, the authors also derive a general nonparametric upper bound on Bayes risks.