Sequential LND sensitivity test for binary response data

Sensitivity tests are used to make inferences about a sensitivity, a characteristic property of some products that cannot be observed directly. For binary response sensitivity data (dead or alive, explode or unexplode), the Langlie and Neyer are two well-known sensitivity tests. The priorities of the Langlie and Neyer tests are investigated in this paper. It is shown that the Langlie test has an advantage in getting an overlap, while the Neyer test has better estimation precision. Aiming at improving both the speed of getting an overlap and the estimation precision, we propose a new sensitivity test which replaces the first part of the Neyer test with the Langlie test. Our simulation studies indicate that the proposed test outperforms the Langlie, Neyer and Dror and Steinberg tests from the viewpoints of estimation precision and probability of obtaining an overlap.     

Estimation procedures for grouped data – a comparative study

Interval-censored data are very common in the reliability and lifetime data analysis. This paper investigates the performance of different estimation procedures for a special type of interval-censored data, i.e. grouped data, from three widely used lifetime distributions. The approaches considered here include the maximum likelihood estimation, the minimum distance estimation based on chi-square criterion, the moment estimation based on imputation (IM) method and an ad hoc estimation procedure. Although IM-based techniques are extensively used recently, we show that this method is not always effective. It is found that the ad hoc estimation procedure is equivalent to the minimum distance estimation with another distance metric and more effective in the simulation. The procedures of different approaches are presented and their performances are investigated by Monte Carlo simulation for various combinations of sample sizes and parameter settings. The numerical results provide guidelines to analyse grouped data for practitioners when they need to choose a good estimation approach.     

空间滞后分位数回归模型的工具变量估计及参数检验

本文提出一种遵循空间数据分布特征的空间分位数回归模型,并着重探讨该模型的估计方法和参数检验问题。本文构建了上述模型的一个工具变量估计法,通过数理证明建立了估计量的大样本理论,并基于估计量的渐近分布构造了模型的参数检验方法。本文还通过数值模拟方法和应用实例考察估计方法和参数检验方法的实际应用效果,数值模拟结果显示,估计方法和参数检验方法在有限样本条件下均可以达到较高的精确度和稳定性。在应用实例中,本文利用所构建的理论方法重新检验我国“资源诅咒”效应的存在性,实证结果体现了理论方法的应用价值。     

Robust estimation for partially linear single-index model

In this article, we investigate a new estimation approach for the partially linear single-index model based on modal regression method, where the non parametric function is estimated by penalized spline method. Moreover, we develop an expection maximum (EM)-type algorithm and establish the large sample properties of the proposed estimation method. A distinguishing characteristic of the newly proposed estimation is robust against outliers through introducing an additional tuning parameter which can be automatically selected using the observed data. Simulation studies and real data example are used to evaluate the finite-sample performance, and the results show that the newly proposed method works very well.     

A unified empirical likelihood approach for testing MCAR and subsequent estimation

For an estimation with missing data, a crucial step is to determine if the data are missing completely at random (MCAR), in which case a complete‐case analysis would suffice. Most existing tests for MCAR do not provide a method for a subsequent estimation once the MCAR is rejected. In the setting of estimating means, we propose a unified approach for testing MCAR and the subsequent estimation. Upon rejecting MCAR, the same set of weights used for testing can then be used for estimation. The resulting estimators are consistent if the missingness of each response variable depends only on a set of fully observed auxiliary variables and the true outcome regression model is among the user‐specified functions for deriving the weights. The proposed method is based on the calibration idea from survey sampling literature and the empirical likelihood theory.     

Existence,uniqueness and consistency of estimation of life characteristics of three-parameter Weibull distribution based on Type-II right censored data

In this paper, we propose a new method of estimation for the parameters and quantiles of the three-parameter Weibull distribution based on Type-II right censored data. The method, based on a data transformation, overcomes the problem of unbounded likelihood. In the proposed method, under mild conditions, the estimates always exist uniquely, and the estimators are also consistent over the entire parameter space. Through Monte Carlo simulations, we further show that the proposed method of estimation performs well compared to some prominent methods in terms of bias and root mean squared error in small-sample situations. Finally, two real data sets are used to illustrate the proposed method of estimation.     

Estimation of k-Factor GIGARCH Process: A Monte Carlo Study

In this article, we discuss the parameter estimation for a k-factor generalized long-memory process with conditionally heteroskedastic noise. Two estimation methods are proposed. The first method is based on the conditional distribution of the process and the second is obtained as an extension of Whittle's estimation approach. For comparison purposes, Monte Carlo simulations are used to evaluate the finite sample performance of these estimation techniques, using four different conditional distribution functions.     

Relative Risk Regression for Binary Outcomes: Methods and Recommendations

Relative risks are often considered preferable to odds ratios for quantifying the association between a predictor and a binary outcome. Relative risk regression is an alternative to logistic regression where the parameters are relative risks rather than odds ratios. It uses a log link binomial generalised linear model, or log‐binomial model, which requires parameter constraints to prevent probabilities from exceeding 1. This leads to numerical problems with standard approaches for finding the maximum likelihood estimate (MLE), such as Fisher scoring, and has motivated various non‐MLE approaches. In this paper we discuss the roles of the MLE and its main competitors for relative risk regression. It is argued that reliable alternatives to Fisher scoring mean that numerical issues are no longer a motivation for non‐MLE methods. Nonetheless, non‐MLE methods may be worthwhile for other reasons and we evaluate this possibility for alternatives within a class of quasi‐likelihood methods. The MLE obtained using a reliable computational method is recommended, but this approach requires bootstrapping when estimates are on the parameter space boundary. If convenience is paramount, then quasi‐likelihood estimation can be a good alternative, although parameter constraints may be violated. Sensitivity to model misspecification and outliers is also discussed along with recommendations and priorities for future research.     

测量误差模型稳健方法及应用研究

回归模型在经济学、生物医学、流行病学、工农业生产等众多领域有着广泛的应用,而在实际数据收集时常常出现无法获得变量的精确数据或全部数据的情况,即常碰到测量误差数据、缺失数据等复杂数据情形。对于回归模型中存在测量误差的情况,如在参数估计时不加以修正,则易产生估计偏差,使得估计精度下降。对于数据缺失情形,如不采取合理的处理方法也会导致模型分析结果不佳。故此,本文研究含有测量误差数据时,解释变量具有随机缺失时的线性测量误差模型和部分线性测量误差模型的稳健参数估计问题。本文提出了一种在测量误差服从拉普拉斯分布时参数的损失修正估计,通过蒙特卡洛模拟和医学研究中的实证分析,显示本文所提的估计方法具有偏差小、精度高、稳健性强的优势。     

Smooth Semi‐nonparametric Analysis for Mixture Cure Models and Its Application to Breast Cancer

Mixture cure models are widely used when a proportion of patients are cured. The proportional hazards mixture cure model and the accelerated failure time mixture cure model are the most popular models in practice. Usually the expectation–maximisation (EM) algorithm is applied to both models for parameter estimation. Bootstrap methods are used for variance estimation. In this paper we propose a smooth semi‐nonparametric (SNP) approach in which maximum likelihood is applied directly to mixture cure models for parameter estimation. The variance can be estimated by the inverse of the second derivative of the SNP likelihood. A comprehensive simulation study indicates good performance of the proposed method. We investigate stage effects in breast cancer by applying the proposed method to breast cancer data from the South Carolina Cancer Registry.     

An improved method of estimation for the parameters of the Birnbaum–Saunders distribution

The Birnbaum–Saunders (BS) distribution is a positively skewed distribution, frequently used for analysing lifetime data. In this paper, we propose a simple method of estimation for the parameters of the two-parameter BS distribution by making use of some key properties of the distribution. Compared with the maximum likelihood estimators and the modified moment estimators, the proposed method has smaller bias, but having the same mean square errors as these two estimators. We also discuss some methods of construction of confidence intervals. The performance of the estimators is then assessed by means of Monte Carlo simulations. Finally, an example is used to illustrate the method of estimation developed here.     

小域估计方法在小品种农作物播种面积统计中的应用

 我国以省级目标总体开展的现行农产量抽样调查,着眼于农作物主要品种的省级推算,而小品种农作物的总体分布比较偏态,往往有效样本量相对不足,不能解决小品种农作物播种面积的推算问题,同时对分县的主要品种农作物播种面积进行的直接推算也不能满足精度要求。现阶段对小品种农作物播种面积的统计方法研究成为农村统计方法制度改革迫切需要研究的课题之一。本文选择了河北省张家口的蔚县,利用小域估计方法对小品种农作物播种面积进行了统计推断,从推断结果看得到了比较好的估计精度。利用蔚县为总体的实际数据进行的抽样仿真分析,从实证的角度阐述了小域估计方法对这一问题的有效性,而且分析结果也表明该方法可以显著提高估计效果。     

Consistent estimation of parameters and quantiles of the three-parameter gamma distribution based on Type-II right-censored data

In this paper, we propose a method of estimation of parameters and quantiles of the three-parameter gamma distribution based on Type-II right-censored data. In the proposed method, under mild conditions, the estimates always exist uniquely, and the estimators have consistency over the entire parameter space. Through Monte Carlo simulations, we further show that the proposed method performs well compared with another prominent method of estimation in terms of bias and root mean-squared error in small-sample situations. Finally, two real data sets are used for illustrating the proposed method.     

Bayes estimation of Moran–Downton bivariate exponential distribution based on censored samples

In this paper, we discuss a Bayesian estimation procedure for the parameters in a Moran–Downton bivariate exponential distribution based on complete and censored samples. A Markov-chain Monte Carlo method is used to obtain the Bayes estimates of the parameters. An intensive simulation experiment is conducted to study the performance of the proposed Bayesian estimation procedure. Discussions and suggestions are provided based on the simulation results. A numerical example is presented to illustrate the Bayesian estimation procedure developed here and some concluding remarks are provided.     

Mass imputation for two-phase sampling

Two-phase sampling is a cost-effective method of data collection using outcome-dependent sampling for the second-phase sample. In order to make efficient use of auxiliary information and to improve domain estimation, mass imputation can be used in two-phase sampling. Rao and Sitter (1995) introduce mass imputation for two-phase sampling and its variance estimation under simple random sampling in both phases. In this paper, we extend the Rao–Sitter method to general sampling design. The proposed method is further extended to mass imputation for categorical data. A limited simulation study is performed to examine the performance of the proposed methods.     

A Comparison of Utilized and Theoretical Covariance Weighting Matrices on the Estimation Performance of Quadratic Inference Functions

The quadratic inference function (QIF) method is increasingly popular for the marginal analysis of correlated data due to its advantages over generalized estimating equations. Asymptotic theory is used to derive analytical results from the QIF, and we, therefore, study three asymptotically equivalent weighting matrices in terms of finite-sample parameter estimation. Furthermore, to improve small-sample estimation, we study modifications to the estimation procedure. Examples are presented via simulations and application. Results show that although theoretical weighting matrices work best, the proposed estimation procedure, in which initial estimates are held constant within the matrix of estimated empirical covariances, is preferable in practice.     

A nonparametric plug-in rule for selecting optimal block lengths for block bootstrap methods

In this paper, we consider the problem of empirical choice of optimal block sizes for block bootstrap estimation of population parameters. We suggest a nonparametric plug-in principle that can be used for estimating ‘mean squared error’-optimal smoothing parameters in general curve estimation problems, and establish its validity for estimating optimal block sizes in various block bootstrap estimation problems. A key feature of the proposed plug-in rule is that it can be applied without explicit analytical expressions for the constants that appear in the leading terms of the optimal block lengths. Furthermore, we also discuss the computational efficacy of the method and explore its finite sample properties through a simulation study.     

Some Uniform Limit Results in Additive Regression Model

In the setting of additive regression model components estimation, we establish some uniform limit results in probability. Our results allow to give an asymptotic simultaneous 100% confidence band for these components. These results are stated in the framework of i.i.d random vectors when the marginal integration estimation method is used.     

A new orthogonality-based estimation for varying-coefficient partially linear models

Varying coefficient partially linear models are usually used for longitudinal data analysis, and an interest is mainly to improve efficiency of regression coefficients. By the orthogonality estimation technology and the quadratic inference function method, we propose a new orthogonality-based estimation method to estimate parameter and nonparametric components in varying coefficient partially linear models with longitudinal data. The proposed procedure can separately estimate the parametric and nonparametric components, and the resulting estimators do not affect each other. Under some mild conditions, we establish some asymptotic properties of the resulting estimators. Furthermore, the finite sample performance of the proposed procedure is assessed by some simulation experiments.