Doubly robust empirical likelihood inference in covariate-missing data problems

Missing covariate data occurs often in regression analysis. We study methods for estimating the regression coefficients in an assumed conditional mean function when some covariates are completely observed but other covariates are missing for some subjects. We adopt the semiparametric perspective of Robins et al. [Estimation of regression coefficients when some regressors are not always observed. J Amer Statist Assoc. 1994;89:846–866] on regression analyses with missing covariates, in which they pioneered the use of two working models, the working propensity score model and the working conditional score model. A recent approach to missing covariate data analysis is the empirical likelihood method of Qin et al. [Empirical likelihood in missing data problems. J Amer Statist Assoc. 2009;104:1492–1503], which effectively combines unbiased estimating equations. In this paper, we consider an alternative likelihood approach based on the full likelihood of the observed data. This full likelihood-based method enables us to generate estimators for the vector of the regression coefficients that are (a) asymptotically equivalent to those of Qin et al. [Empirical likelihood in missing data problems. J Amer Statist Assoc. 2009;104:1492–1503] when the working propensity score model is correctly specified, and (b) doubly robust, like the augmented inverse probability weighting (AIPW) estimators of Robins et al. [Estimation of regression coefficients when some regressors are not always observed. J Am Statist Assoc. 1994;89:846–866]. Thus, the proposed full likelihood-based estimators improve on the efficiency of the AIPW estimators when the working propensity score model is correct but the working conditional score model is possibly incorrect, and also improve on the empirical likelihood estimators of Qin, Zhang and Leung [Empirical likelihood in missing data problems. J Amer Statist Assoc. 2009;104:1492–1503] when the reverse is true, that is, the working conditional score model is correct but the working propensity score model is possibly incorrect. In addition, we consider a regression method for estimation of the regression coefficients when the working conditional score model is correctly specified; the asymptotic variance of the resulting estimator is no greater than the semiparametric variance bound characterized by the theory of Robins et al. [Estimation of regression coefficients when some regressors are not always observed. J Amer Statist Assoc. 1994;89:846–866]. Finally, we compare the finite-sample performance of various estimators in a simulation study.     

Randomized Response Techniques Using Maximum Likelihood Estimator

The randomized response (RR) technique with two decks of cards proposed by Odumade and Singh (2009) can always be made more efficient than the RR techniques proposed by Warner (1965), Mangat and Singh (1990), and Mangat (1994) by adjusting the proportion of cards in the decks. Arnab et al. (2012) generalized Odumade and Singh strategy (2009) for complex survey designs and wider class of estimators. In this paper improvement of Arnab et al. (2012) estimator has been made by using maximum likelihood method.     

Optimal Bayesian variable sampling plans for exponential distributions based on modified type-II hybrid censored samples

For the conventional type-II hybrid censoring scheme (HCS) in Childs et al., a Bayesian variable sampling plan among the class of the maximum likelihood estimators was derived by Lin et al. under the loss function, which does not include the cost of experimental time. Instead of taking the conventional type-II hybrid censoring scheme, a persuasive argument leads to taking the modified type-II hybrid censoring scheme (MHCS) if the cost of experimental time is included in the loss function. In this article, we apply the decision-theoretic approach for the concerned acceptance sampling. With the type-II MHCS, based on a sufficient statistics, the optimal Bayesian sampling plan is derived under a general loss function. Furthermore, for the conjugate prior distribution, the closed-form formula of the Bayes decision rule can be obtained under the quadratic decision loss. Numerical study is given to demonstrate the performance of the proposed Bayesian sampling plan.     

A stratified estimation of a sensitive character by using geometric distribution

Abstract

In this paper, we consider the estimation of a sensitive character when the population is consisted of several strata; this is undertaken by applying Niharika et al.’s model which is using geometric distribution as a randomization device. A sensitive parameter is estimated for the case in which stratum size is known, and proportional and optimum allocation methods are taken into account. We extended the Niharika et al.’s model to the case of an unknown stratum size; a sensitive parameter is estimated by applying stratified double sampling to the Niharika et al.’s model. Finally, the efficiency of the proposed model is compared with that of Niharika et al. in terms of the estimator variance.     

Large deviations for the stochastic present value of aggregate claims in the nonstandard compound renewal risk model with widely upper Orthant dependent claims

Abstract

Recently, Jiang et al. (Statist. Probab. Lett. 101, 83–91) obtained the asymptotic formulas for the large deviations for the stochastic present value of aggregate claims in the renewal risk model with Pareto-type claims and stochastic return on investments, where the price process of the investment portfolio is described as a geometric Lévy process. In the paper, we extend the above results to a nonstandard compound renewal risk model with widely upper orthant dependent and dominatedly-varying-tailed claims.     

Mean and sensitivity estimation of a sensitive variable through additive scrambling

Abstract

This article focuses on reducing the additional variance due to randomization of the responses. The idea of additive scrambling and its inverse has been used along with (i) split sample approach and (ii) double response approach. Specifically, our proposal is based on Gupta et al. (2006) randomized response model. We selected this model for improvement because it provides estimator of mean and sensitivity level of a sensitive variable and is better than all of its competitors proposed earlier to it and even Gupta et al. (2006) sensitivity estimator is better than that of Gupta et al. (2010). Our suggested estimators are unbiased estimators and perform better than Gupta et al. (2006) estimator. The issue of privacy protection is also discussed.     

An FFT approach for option pricing under a regime-switching stochastic interest rate model

In this article, we investigate the pricing of European-style options under a Markovian regime-switching Hull–White interest rate model. The parameters of this model, including the mean-reversion level, the volatility of the stochastic interest rate, and the volatility of an asset’s value, are modulated by an observable, continuous-time, finite-state Markov chain. A closed-form expression for the characteristic function of the logarithmic terminal asset price is derived. Then, using the fast Fourier transform, a price of a European-style option is computed. In a two-state Markov chain case, numerical examples and empirical studies are presented to illustrate the practical implementation of the model.     

Estimation of population proportion of a qualitative character using randomized response technique in stratified random sampling

This paper proposes an efficient stratified randomized response model based on Chang et al.'s (2004) model. We have obtained the variance of the proposed estimator of π_s, the proportion of the respondents in the population belonging to a sensitive group, under proportional and Neyman allocations. It is shown that the estimator based on the proposed model is more efficient than the Chang et al.'s (2004) estimator under both proportional as well as Neyman allocations, Hong et al.'s (1994) estimator and Kim and Warde's (2004) estimator. Numerical illustration and pictorial representation are given in support of the present study.     

Modified ratio estimators using robust regression methods

When there is an outlier in the data set, the efficiency of traditional methods decreases. In order to solve this problem, Kadilar et al. (2007) adapted Huber-M method which is only one of robust regression methods to ratio-type estimators and decreased the effect of outlier problem. In this study, new ratio-type estimators are proposed by considering Tukey-M, Hampel M, Huber MM, LTS, LMS and LAD robust methods based on the Kadilar et al. (2007). Theoretically, we obtain the mean square error (MSE) for these estimators. We compared with MSE values of proposed estimators and MSE values of estimators based on Huber-M and OLS methods. As a result of these comparisons, we observed that our proposed estimators give more efficient results than both Huber M approach which was proposed by Kadilar et al. (2007) and OLS approach. Also, under all conditions, all of the other proposed estimators except Lad method are more efficient than robust estimators proposed by Kadilar et al. (2007). And, these theoretical results are supported with the aid of a numerical example and simulation by basing on data that includes an outlier.     

The generalized family of estimators of population mean using auxiliary information in double sampling

Abstract

We suggested the class of estimators of the population mean with its bias and mean square error. It has been shown that the suggested class is more efficient than the usual unbiased, ratio, product and regression estimators and estimators due to Bahl and Tuteja (1991), Singh et al. (2009), and Upadhyaya et al. (2011). In addition an empirical study also carried out to and founded that the members of suggested family also have improvement over Grover and Kaur (2011) and Shabbir and Gupta (2011) classes. Two-phase (double) sampling version of the proposed class was also given.     

A Study on EWMA charts with runs rules—the Markov chain approach

ABSTRACT

Runs rules are usually used with Shewhart-type charts to enhance the charts' sensitivities toward small and moderate shifts. Abbas et al. in 2011 took it a step further by proposing two runs rules schemes, applied to the exponentially weighted moving average (EWMA) chart and evaluated their average run length (ARL) performances using simulation. They showed that the proposed schemes are superior to the classical EWMA chart and other schemes being investigated. Besides pointing out some erroneous ARL and standard deviation of the run length (SDRL) computations in Abbas et al., this paper presents a Markov chain approach for computing the ARL, percentiles of the run length (RL) distribution and SDRL, for the two runs rules schemes of Abbas et al. Using Markov chain, we also propose two combined runs rules EWMA schemes to quicken the two schemes of Abbas et al. in responding to large shifts. The runs rules (basic and combined rules) EWMA schemes will be compared with some existing control charting methods, where the former charts are shown to prevail.     

On the Computation of Gauss Hypergeometric Functions

The pioneering study undertaken by Liang et al. in 2008 (Journal of the American Statistical Association, 103, 410–423) and the hundreds of papers citing that work make use of certain hypergeometric functions. Liang et al. and many others claim that the computation of the hypergeometric functions is difficult. Here, we show that the hypergeometric functions can in fact be reduced to simpler functions that can often be computed using a pocket calculator.     

Bayesian Approach in Nonparametric Count Regression with Binomial Kernel

Recently, Kokonendji et al. have adapted the well-known Nadaraya–Watson kernel estimator for estimating the count function m in the context of nonparametric discrete regression. The authors have also investigated the bandwidth selection using the cross-validation method. In this article, we propose a Bayesian approach in the context of nonparametric count regression for estimating the bandwidth and the variance of the model error, which has not been estimated in Kokonendji et al. The model error is considered as Gaussian with mean of zero and a variance of σ². The Bayes estimates cannot be obtained in closed form and then, we use the well-known Markov chain Monte Carlo (MCMC) technique to compute the Bayes estimates under the squared errors loss function. The performance of this proposed approach and the cross-validation method are compared through simulation and real count data.     

Goodness-of-fit testing for the marginal distribution of regime-switching models with an application to electricity spot prices

This paper complements a recently published study (Janczura and Weron in AStA-Adv Stat Anal 96(3):385–407, 2012) on efficient estimation of Markov regime-switching models. Here, we propose a new goodness-of-fit testing scheme for the marginal distribution of such models. We consider models with an observable (like threshold autoregressions) as well as a latent state process (like Markov regime-switching). The test is based on the Kolmogorov–Smirnov supremum-distance statistic and the concept of the weighted empirical distribution function. The motivation for this research comes from a recent stream of literature in energy economics concerning electricity spot price models. While the existence of distinct regimes in such data is generally unquestionable (due to the supply stack structure), the actual goodness-of-fit of the models requires statistical validation. We illustrate the proposed scheme by testing whether commonly used Markov regime-switching models fit deseasonalized electricity prices from the NEPOOL (US) day-ahead market.     

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.     

Calibration Estimator of a Rare Sensitive Attribute with Poisson Distribution

In this paper we consider the calibration procedure for a rare sensitive attribute with Poisson distribution which suggested by Land et al. (2012) using auxiliary information associated with the variable of interest. In the calibration procedure, we can use auxiliary information such as socio-demographical variables for the respondents of rare sensitive attribute questions from an external source, and then this estimator can be improved with respect to the problems of non coverage or non response. From the efficiency comparison study, we show that the calibrated Poisson RR estimators are more efficient than that of Land et al. (2012), when the known population cell and marginal counts of auxiliary information are used for the calibration procedure.     

A segmented regime-switching model with its application to stock market indices

This paper evaluates the ability of a Markov regime-switching log-normal (RSLN) model to capture the time-varying features of stock return and volatility. The model displays a better ability to depict a fat tail distribution as compared with using a log-normal model, which means that the RSLN model can describe observed market behavior better. Our major objective is to explore the capability of the model to capture stock market behavior over time. By analyzing the behavior of calibrated regime-switching parameters over different lengths of time intervals, the change-point concept is introduced and an algorithm is proposed for identifying the change-points in the series corresponding to the times when there are changes in parameter estimates. This algorithm for identifying change-points is tested on the Standard and Poor's 500 monthly index data from 1971 to 2008, and the Nikkei 225 monthly index data from 1984 to 2008. It is evident that the change-points we identify match the big events observed in the US stock market and the Japan stock market (e.g., the October 1987 stock market crash), and that the segmentations of stock index series, which are defined as the periods between change-points, match the observed bear–bull market phases.     

The IPW estimator for the joint distribution function of the gap times from recurrent event data

Lee et al. in 2016 proposed a nonparametric estimator of the joint distribution of the gap time between transplant and the first infection and the following gap times between consecutive infections. In this article, we propose an alternative estimator based on the inverse-probability weighted (IPW) approach. Asymptotic properties of the proposed estimator are established . Simulation results indicate that the IPW estimator performs as well as the estimator proposed by Lee et al. We also propose an IPW estimator for estimating the joint distribution function of the gap times between consecutive recurrent events beyond the first episode.     

The generalizations of strong law of large numbers for asymptotic even–odd Markov chains indexed by a homogeneous tree

Yang et al. (Yang et al., J. Math. Anal. Appl., 410 (2014), 179–189.) have obtained the strong law of large numbers and asymptotic equipartition property for the asymptotic even–odd Markov chains indexed by a homogeneous tree. In this article, we are going to study the strong law of large numbers and the asymptotic equipartition property for a class of non homogeneous Markov chains indexed by a homogeneous tree which are the generalizations of above results. We also provide an example showing that our generalizations are not trivial.     

Transformed symmetric generalized autoregressive moving average models

Cordeiro and Andrade [Transformed generalized linear models. J Stat Plan Inference. 2009;139:2970–2987] incorporated the idea of transforming the response variable to the generalized autoregressive moving average (GARMA) model, introduced by Benjamin et al. [Generalized autoregressive moving average models. J Am Stat Assoc. 2003;98:214–223], thus developing the transformed generalized autoregressive moving average (TGARMA) model. The goal of this article is to develop the TGARMA model for symmetric continuous conditional distributions with a possible nonlinear structure for the mean that enables the fitting of a wide range of models to several time series data types. We derive an iterative process for estimating the parameters of the new model by maximum likelihood and obtain a simple formula to estimate the parameter that defines the transformation of the response variable. Furthermore, we determine the moments of the original dependent variable which generalize previous published results. We illustrate the theory by means of real data sets and evaluate the results developed through simulation studies.