Bivariate Kumaraswamy distribution: properties and a new method to generate bivariate classes

In this paper, we introduce a bivariate Kumaraswamy (BVK) distribution whose marginals are Kumaraswamy distributions. The cumulative distribution function of this bivariate model has absolutely continuous and singular parts. Representations for the cumulative and density functions are presented and properties such as marginal and conditional distributions, product moments and conditional moments are obtained. We show that the BVK model can be obtained from the Marshall and Olkin survival copula and obtain a tail dependence measure. The estimation of the parameters by maximum likelihood is discussed and the Fisher information matrix is determined. We propose an EM algorithm to estimate the parameters. Some simulations are presented to verify the performance of the direct maximum-likelihood estimation and the proposed EM algorithm. We also present a method to generate bivariate distributions from our proposed BVK distribution. Furthermore, we introduce a BVK distribution which has only an absolutely continuous part and discuss some of its properties. Finally, a real data set is analysed for illustrative purposes.     

Penalized optimal scoring for the classification of multi-dimensional functional data

Many fields of research need to classify individual systems based on one or more data series, which are obtained by sampling an unknown continuous curve with noise. In other words, the underlying process is an unknown function which the observed variables represent only imperfectly. Although functional logistic regression has many attractive features for this classification problem, this method is applicable only when the number of individuals to be classified (or available to estimate the model) is large compared to the number of curves sampled per individual.To overcome this limitation, we use penalized optimal scoring to construct a new method for the classification of multi-dimensional functional data. The proposed method consists of two stages. First, the series of observed discrete values available for each individual are expressed as a set of continuous curves. Next, the penalized optimal scoring model is estimated on the basis of these curves. A similar penalized optimal scoring method was described in my previous work, but this model is not suitable for the analysis of continuous functions. In this paper we adopt a Gaussian kernel approach to extend the previous model. The high accuracy of the new method is demonstrated on Monte Carlo simulations, and used to predict defaulting firms on the Japanese Stock Exchange.     

Asymptotic normality of numbers of observations in random regions determined by order statistics

In this paper, we study the joint limiting behaviour of numbers of observations that fall into regions determined by order statistics and Borel sets. We show that suitably centred and normed versions of these numbers are asymptotically multivariate normal under some conditions. We consider two cases: one where the population distribution function is discontinuous and the other where it is continuous and the order statistics are extreme. Finally, we compare results obtained for the two cases with their analogues for absolutely continuous distribution function and central-order statistics.     

A new absolute continuous bivariate generalized exponential distribution

The generalized exponential is the most commonly used distribution for analyzing lifetime data. This distribution has several desirable properties and it can be used quite effectively to analyse several skewed life time data. The main aim of this paper is to introduce absolutely continuous bivariate generalized exponential distribution using the method of Block and Basu (1974). In fact, the Block and Basu exponential distribution will be extended to the generalized exponential distribution. We call the new proposed model as the Block and Basu bivariate generalized exponential distribution, then, discuss its different properties. In this case the joint probability distribution function and the joint cumulative distribution function can be expressed in compact forms. The model has four unknown parameters and the maximum likelihood estimators cannot be obtained in explicit form. To compute the maximum likelihood estimators directly, one needs to solve a four dimensional optimization problem. The EM algorithm has been proposed to compute the maximum likelihood estimations of the unknown parameters. One data analysis is provided for illustrative purposes. Finally, we propose some generalizations of the proposed model and compare their models with each other.     

Effect of censoring on adjusting for covariates in comparison of survival times

The increase in the variance of the estimate of treatment effect which results from omitting a dichotomous or continuous covariate is quantified as a function of censoring. The efficiency of not adjusting for a covariate is measured by the ratio of the variance obtained with and without adjustment for the covariate. The variance is derived using the Weibull proportional hazards model. Under random censoring, the efficiency of not adjusting for a continuous covariate is an increasing function of the percentage of censored observations.     

On the Estimation of Integrated Volatility With Jumps and Microstructure Noise

In this article, we propose a nonparametric procedure to estimate the integrated volatility of an Itô semimartingale in the presence of jumps and microstructure noise. The estimator is based on a combination of the preaveraging method and threshold technique, which serves to remove microstructure noise and jumps, respectively. The estimator is shown to work for both finite and infinite activity jumps. Furthermore, asymptotic properties of the proposed estimator, such as consistency and a central limit theorem, are established. Simulations results are given to evaluate the performance of the proposed method in comparison with other alternative methods.     

A class of absolutely continuous bivariate distributions

Block and Basu bivariate exponential distribution is one of the most popular absolutely continuous bivariate distributions. Extensive work has been done on the Block and Basu bivariate exponential model over the past several decades. Interestingly it is observed that the Block and Basu bivariate exponential model can be extended to the Weibull model also. We call this new model as the Block and Basu bivariate Weibull model. We consider different properties of the Block and Basu bivariate Weibull model. The Block and Basu bivariate Weibull model has four unknown parameters and the maximum likelihood estimators cannot be obtained in closed form. To compute the maximum likelihood estimators directly, one needs to solve a four dimensional optimization problem. We propose to use the EM algorithm for computing the maximum likelihood estimators of the unknown parameters. The proposed EM algorithm can be carried out by solving one non-linear equation at each EM step. Our method can be also used to compute the maximum likelihood estimators for the Block and Basu bivariate exponential model. One data analysis has been preformed for illustrative purpose.     

Mixed Discrete and Continuous Cox Regression Model

The Cox (1972) regression model is extended to include discrete and mixed continuous/discrete failure time data by retaining the multiplicative hazard rate form of the absolutely continuous model. Application of martingale arguments to the regression parameter estimating function show the Breslow (1974) estimator to be consistent and asymptotically Gaussian under this model. A computationally convenient estimator of the variance of the score function can be developed, again using martingale arguments. This estimator reduces to the usual hypergeometric form in the special case of testing equality of several survival curves, and it leads more generally to a convenient consistent variance estimator for the regression parameter. A small simulation study is carried out to study the regression parameter estimator and its variance estimator under the discrete Cox model special case and an application to a bladder cancer recurrence dataset is provided.     

Modeling and simulation of a nonhomogeneous poisson process having cyclic behavior

In this paper we develop a unified approach to modeling and simulation of a nonhomogeneous Poisson process whose rate function exhibits cyclic behavior as well as a long-term evolutionary trend. The approach can be applied whether the oscillation frequency of the cyclic behavior is known or unknown. To model such a process, we use an exponential rate function whose exponent includes both a polynomial and a trigonometric component.Maximum likelihood estimates of the unknown continuous parameters of this function are obtained numerically, and the degree of the polynomial component is determined by a likelihood ratio test. If the oscillation frequency is unknown, then an initial estimate of this parameter is obtained via spectral analysis of the observed series of events; initial estimates of the remaining trigonometric (respectively, polynomial) parameters are computed from a standard maximum likelihood (respectively, moment-matching) procedure for an exponential-trigonometric (respectively, exponential-polynomial) rate function. To simulate the fitted process by the method of thinning, we present (a) a procedure for constructing an optimal piecewise linear majorizing rate function; and(b)a "piecewise thinning" simulation procedure based on the inverse transform method for generating events from a piecewise linear rate function. These procedures are applied to the storm-arrival process observed at an off-shore drilling site.     

Mixed-spectra analysis for stationary random fields

We consider a, discrete time, weakly stationary bidimensional process, for which the spectral measure is the sum of an absolutely continuous measure, a discrete measure of finite order and a finite number of absolutely continuous measures on several lines. In this paper we are interested in estimating the spectral density of the absolutely continuous measure and of the density on the lines. For this aim, by using the double kernel method, we construct consistent estimators of these densities and we study their asymptotic behaviors in term of the mean squared error with rate.     

Maximum likelihood estimation for reflected Ornstein-Uhlenbeck processes with jumps

In this paper, we investigate the maximum likelihood estimation for the reflected Ornstein-Uhlenbeck processes with jumps based on continuous observations. We derive likelihood functions by using semimartingale theory. From this we get explicit formulas for estimators. The strong consistence and asymptotic normality of estimators are proved by using the method of stochastic integration.     

Fast Censored Linear Regression

Weighted log‐rank estimating function has become a standard estimation method for the censored linear regression model, or the accelerated failure time model. Well established statistically, the estimator defined as a consistent root has, however, rather poor computational properties because the estimating function is neither continuous nor, in general, monotone. We propose a computationally efficient estimator through an asymptotics‐guided Newton algorithm, in which censored quantile regression methods are tailored to yield an initial consistent estimate and a consistent derivative estimate of the limiting estimating function. We also develop fast interval estimation with a new proposal for sandwich variance estimation. The proposed estimator is asymptotically equivalent to the consistent root estimator and barely distinguishable in samples of practical size. However, computation time is typically reduced by two to three orders of magnitude for point estimation alone. Illustrations with clinical applications are provided.     

ON THE BEHRENS–FISHER PROBLEM1

In this article the Behrens-Fisher problem is reformulated in terms of a structural model of inference. For this version of the problem a solution is obtained which is valid for arbitrary absolutely continuous error distributions. These results are further discussed for the standard normal distribution and for some other special cases with not normally distributed populations.     

Spectral analysis with replicated time series

A doubly stochastic process {x(b,t);b?B,t?Z} is considered, with (B,β,P_β) being a probability space so that for each b, {X(b,t);t ? Z} is a stationary process with an absolutely continuous spectral distribution. The population spectrum is defined as f(ω) = E_B[Q(b,ω)] with Q(b,ω) being the spectral density function of X(b,t). The aim of this paper is to estimate f(ω) by means of a random sample b₁,…,b_r from (B,β,P_β). For each b₁? B, the processes X(b₁,t) are observed at the same times t=1,…,N. Thus, r time series (x(b₁,t)} are available in order to estimate f(ω). A model for each individual periodogram, which involves f(ω), is formulated. It has been proven that a certain family of linear stationary processes follows the above model In this context, a kernel estimator is proposed in order to estimate f(ω). The bias, variance and asymptotic distribution of this estimator are investigated under certain conditions.     

Improvement of the Liu Estimator in Weighted Mixed Regression

ABSTRACT

The maximum likelihood approach to the proportional hazards model is considered. The purpose is to find a general approach to the analysis of the proportional hazards model, whether the baseline distribution is absolutely continuous, discrete, or a mixture. The advantage is that ties are treated without pain, while the performance for continuous data is almost the same as Cox's partial likelihood. The potential disadvantage with many nuisance parameters is taken care of by profiling them out for risk sets containing only one failure.     

Empirical Characteristic Function Estimation and Its Applications

This paper reviews the method of model-fitting via the empirical characteristic function. The advantage of using this procedure is that one can avoid difficulties inherent in calculating or maximizing the likelihood function. Thus it is a desirable estimation method when the maximum likelihood approach encounters difficulties but the characteristic function has a tractable expression. The basic idea of the empirical characteristic function method is to match the characteristic function derived from the model and the empirical characteristic function obtained from data. Ideas are illustrated by using the methodology to estimate a diffusion model that includes a self-exciting jump component. A Monte Carlo study shows that the finite sample performance of the proposed procedure offers an improvement over a GMM procedure. An application using over 72 years of DJIA daily returns reveals evidence of jump clustering.     

Markov bridges and enlarged filtrations

The semimartingale decomposition of a Markov chain, whose value at some future time is known, is obtained by considering an enlarged filtration.     

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.     

The modified Fréchet distribution and its properties

The Fréchet distribution is an absolutely continuous model which has wide applicability in extreme value theory. In this paper, we propose a new three-parameter model, so-called the modified Fréchet distribution, to extend the Fréchet distribution. By using the Lambert function, we obtain some properties of the new distribution. We provide a simulation study to illustrate the performance of the maximum likelihood estimates. The flexibility of the introduced distribution is illustrated by means of a real data set. We use some goodness-of-fit statistics to verify the adequacy of the proposed model. We prove empirically that it is appropriate for lifetime applications.     

Estimating the Spot Covariation of Asset Prices—Statistical Theory and Empirical Evidence

ABSTRACT

We propose a new estimator for the spot covariance matrix of a multi-dimensional continuous semimartingale log asset price process, which is subject to noise and nonsynchronous observations. The estimator is constructed based on a local average of block-wise parametric spectral covariance estimates. The latter originate from a local method of moments (LMM), which recently has been introduced by Bibinger et al.. We prove consistency and a point-wise stable central limit theorem for the proposed spot covariance estimator in a very general setup with stochastic volatility, leverage effects, and general noise distributions. Moreover, we extend the LMM estimator to be robust against autocorrelated noise and propose a method to adaptively infer the autocorrelations from the data. Based on simulations we provide empirical guidance on the effective implementation of the estimator and apply it to high-frequency data of a cross-section of Nasdaq blue chip stocks. Employing the estimator to estimate spot covariances, correlations, and volatilities in normal but also unusual periods yields novel insights into intraday covariance and correlation dynamics. We show that intraday (co-)variations (i) follow underlying periodicity patterns, (ii) reveal substantial intraday variability associated with (co-)variation risk, and (iii) can increase strongly and nearly instantaneously if new information arrives. Supplementary materials for this article are available online.