On goodness-of-fit tests for Aalen's additive model based on left-truncated right-censored or doubly censored data

ABSTRACT

Gandy and Jensen (2005 Gandy, A., Jensen, U. (2005). On goodness-of-fit tests for Aalen's additive risk model. Scan. J. Stat. 32:425–445.[Crossref], [Web of Science ®] , [Google Scholar]) proposed goodness-of-fit tests for Aalen's additive risk model. In this article, we demonstrate that the approach of Gandy and Jensen (2005 Gandy, A., Jensen, U. (2005). On goodness-of-fit tests for Aalen's additive risk model. Scan. J. Stat. 32:425–445.[Crossref], [Web of Science ®] , [Google Scholar]) can be applied to left-truncated right-censored (LTRC) data and doubly censored data. A simulation study is conducted to investigate the performance of the proposed tests. The proposed tests are illustrated using heart transplant data.     

MULTI-LEVEL FACTORIAL DESIGNS WITH MINIMUM NUMBERS OF LEVEL CHANGES

The order of experimental runs in a fractional factorial experiment is essential when the cost of level changes in factors is considered. The generalized foldover scheme given by [1] Coster, D. C. and Cheng, C. S. 1988. Minimum cost trend free run orders of fractional factorial designs. The Annals of Statistics, 16: 1188–1205. [Crossref], [Web of Science ®] , [Google Scholar]gives an optimal order to experimental runs in an experiment with specified defining contrasts. An experiment can be specified by a design requirement such as resolution or estimation of some interactions. To meet such a requirement, we can find several sets of defining contrasts. Applying the generalized foldover scheme to these sets of defining contrasts, we obtain designs with different numbers of level changes and then the design with minimum number of level changes. The difficulty is to find all the sets of defining contrasts. An alternative approach is investigated by [2] Cheng, C. S., Martin, R. J. and Tang, B. 1998. Two-level factorial designs with extreme numbers of level changes. The Annals of Statistics, 26: 1522–1539. [Crossref], [Web of Science ®] , [Google Scholar]for two-level fractional factorial experiments. In this paper, we investigate experiments with all factors in slevels.     

On Augmented Franses Tests for Seasonal Unit Roots

This article extends the results reported in del Barrio Castro, Osborn and Taylor (2012 del Barrio Castro, T., Osborn, D.R., Taylor, A. M.R. (2012). On augmented HEGY tests for seasonal unit roots. Econometric Theor. 18:1121–1143.[Crossref], [Web of Science ®] , [Google Scholar]) to the approach followed by Franses (1991a Franses, P. H. (1991a). Model selection and seasonality in time series. Tibergen Institute Series, 18. [Google Scholar],b Franses, P.H. (1991b). Seasonality, non-stationarity and the forecasting of monthly time series. Int. J. Forecast. 7:199–208.[Crossref], [Web of Science ®] , [Google Scholar]) to test for seasonal unit roots, providing the asymptotic representation to the seasonal unit roots tests proposed by Franses for a general number of seasons S.     

Bayesian estimation of sensitivity level and population proportion of a sensitive characteristic in a binary optional unrelated question RRT model

Sihm et al. (2016 Sihm, J. S., A. Chhabra, and S. N. Gupta. 2016. An optional unrelated question RRT model. Involve: A Journal of Mathematics 9 (2):195–209.[Crossref] , [Google Scholar]) proposed an unrelated question binary optional randomized response technique (RRT) model for estimating the proportion of population that possess a sensitive characteristic and the sensitivity level of the question. In our work, decision theoretic approach has been followed to obtain Bayes estimates of the two parameters along with their corresponding minimal Bayes posterior expected losses (BPEL) using beta prior and squared error loss function (SELF). Relative losses are also examined to compare the performances of the Bayes estimates with those of the classical estimates obtained by Sihm et al. (2016 Sihm, J. S., A. Chhabra, and S. N. Gupta. 2016. An optional unrelated question RRT model. Involve: A Journal of Mathematics 9 (2):195–209.[Crossref] , [Google Scholar]). The results obtained are illustrated with the help of real survey data using non informative prior.     

Quantile-based Chernoff distance for truncated random variables

Several probability distributions such as power-Pareto distribution (see Gilchrist 2000 Gilchrist, W. 2000. Statistical modelling with quantile functions. Boca Raton, FL: Chapman and Hall/CRC.[Crossref] , [Google Scholar] and Hankin and Lee 2006 Hankin, R. K. S., and A. Lee. 2006. A new family of non-negative distributions. Australian and New Zealand Journal of Statistics 48:67–78.[Crossref], [Web of Science ®] , [Google Scholar]), various forms of lambda distributions (see Ramberg and Schmeiser 1974 Ramberg, J. S., and B. W. Schmeiser. 1974. An appropriate method for generating asymmetric random variables. Communications of the ACM 17:78–82.[Crossref], [Web of Science ®] , [Google Scholar] and Freimer et al. 1988 Freimer, M., S. Mudholkar, G. Kollia, and C. T. Lin. 1988. A study of the generalized lambda family. Communications in Statistics - Theory and Methods 17:3547–67.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), Govindarajulu distribution (see Nair, Sankaran, and Vineshkumar 2012 Nair, U. N., P. G. Sankaran, and B. Vineshkumar. 2012. The Govindarajulu distribution: some properties and applications. Communications in Statistics—Theory and Methods 41:4391–406.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), etc., do not have manageable distribution functions, though they have tractable quantile functions. Hence, analytical study of the properties of Chernoff distance of two random variables associated with these distributions via traditional distribution function-based tool becomes difficult. To make this simple, in this paper, we introduce quantile-based Chernoff distance for (left or right) truncated random variables and study its various properties. Some useful bounds as well as characterization results are obtained.     

Note on posterior inference for the Bingham distribution

The properties of high-dimensional Bingham distributions have been studied by Kume and Walker (2014 Kume, A., and S. G. Walker. 2014. On the Bingham distribution with large dimension. Journal of Multivariate Analysis 124:345–52.[Crossref], [Web of Science ®] , [Google Scholar]). Fallaize and Kypraios (2016 Fallaize, C. J., and T. Kypraios. 2016. Exact Bayesian inference for the Bingham distribution. Statistics and Computing 26:349–60.[Crossref], [Web of Science ®] , [Google Scholar]) propose the Bayesian inference for the Bingham distribution and they use developments in Bayesian computation for distributions with doubly intractable normalizing constants (Møller et al. 2006 Møller, J., A. N. Pettitt, R. Reeves, and K. K. Berthelsen. 2006. An efficient Markov chain Monte Carlo method for distributions with intractable normalising constants. Biometrika 93 (2):451–458.[Crossref], [Web of Science ®] , [Google Scholar]; Murray, Ghahramani, and MacKay 2006 Murray, I., Z. Ghahramani, and D. J. C. MacKay. 2006. MCMC for doubly intractable distributions. In Proceedings of the 22nd annual conference on uncertainty in artificial intelligence (UAI-06), 359–66. AUAI Press. [Google Scholar]). However, they rely heavily on two Metropolis updates that they need to tune. In this article, we propose instead a model selection with the marginal likelihood.     

Quantile Regression with Left-Truncated and Right-Censored Data in a Reproducing Kernel Hilbert Space

Li et al. (2007 Li, Y., Liu, Y., Zhu, J. (2007). Quantile regression in reproducing kernel Hilbert spaces. J. Amer. Statist. Assoc. 102:255–268.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) developed an estimation method for quantile functions in a reproducing kernel Hilbert space for complete data, and Park and Kim (2011 Park, J., Kim, J. (2011). Quantile regression with an epsilon-insensitive loss in a reproducing kernel Hilbert space. Statist. Probab. Lett. 81:62–70.[Crossref], [Web of Science ®] , [Google Scholar]) proposed an estimation method using the ε-insensitive loss. This article extends these estimation methods to left-truncated and right-censored data. As a measure of goodness of fit, the check loss and the ε-insensitive loss were used to estimate the quantile function. The ε-insensitive loss can shrink the estimated coefficients toward zero; hence, it can reduce the variability of the estimates. Simulation studies show that the estimated quantile functions based on the ε-insensitive loss perform slightly better when ε is adequately chosen.     

Restricted estimation and testing of hypothesis in linear measurement errors models

ABSTRACT

In this article, the linear models with measurement error both in the response and in the covariates are considered. Following Shalabh et al. (2007 Shalabh, Garg, G., Misra, N. (2007). Restricted regression estimation in measurement error models. Comput. Stat. Data Anal. 52:1149–1166.[Crossref], [Web of Science ®] , [Google Scholar], 2009 Shalabh, Garg, G., Misra, N. (2009). Use of prior information in the consistent estimation of regression coefficients in measurement error models. J. Multivariate Anal. 100:1498–1520.[Crossref], [Web of Science ®] , [Google Scholar]), we propose several restricted estimators for the regression coefficients. The consistency and asymptotic normality of the restricted estimators are established. Furthermore, we also discuss the superiority of the restricted estimators to unrestricted estimators under Pitman closeness criterion. We also develop several variance estimators and establish their asymptotic distributions. Wald-type statistics are constructed for testing the linear restrictions. Finally, Monte Carlo simulations are conducted to illustrate the finite-sample properties of the proposed estimators.     

Bayes minimax ridge regression estimators

The problem of estimating of the vector β of the linear regression model y = Aβ + ? with ? ～ N_p(0, σ²I_p) under quadratic loss function is considered when common variance σ² is unknown. We first find a class of minimax estimators for this problem which extends a class given by Maruyama and Strawderman (2005 Maruyama, Y., and W. E. Strawderman. 2005. A new class of generalized Bayes minimax ridge regression estimators. Annals of Statistics 33:1753–70.[Crossref], [Web of Science ®] , [Google Scholar]) and using these estimators, we obtain a large class of (proper and generalized) Bayes minimax estimators and show that the result of Maruyama and Strawderman (2005 Maruyama, Y., and W. E. Strawderman. 2005. A new class of generalized Bayes minimax ridge regression estimators. Annals of Statistics 33:1753–70.[Crossref], [Web of Science ®] , [Google Scholar]) is a special case of our result. We also show that under certain conditions, these generalized Bayes minimax estimators have greater numerical stability (i.e., smaller condition number) than the least-squares estimator.     

ESTIMATION OF PARAMETERS IN ARFIMA PROCESSES: A SIMULATION STUDY

It is known that, in the presence of short memory components, the estimation of the fractional parameter d in an Autoregressive Fractionally Integrated Moving Average, ARFIMA(p, d, q), process has some difficulties (see [1] Smith, J., Taylor, N. and Yadav, S. 1997. Comparing the bias and misspecification in ARFIMA models. Journal of Time Series Analysis, 18(5): 507–527. [Crossref] , [Google Scholar]). In this paper, we continue the efforts made by Smith et al. [1] Smith, J., Taylor, N. and Yadav, S. 1997. Comparing the bias and misspecification in ARFIMA models. Journal of Time Series Analysis, 18(5): 507–527. [Crossref] , [Google Scholar] and Beveridge and Oickle [2] Beveridge, S. and Oickle, C. 1993. Estimating fractionally integrated time series models. Economics Letters, 43: 137–142.  [Google Scholar] by conducting a simulation study to evaluate the convergence properties of the iterative estimation procedure suggested by Hosking [3] Hosking, J. 1981. Fractional differencing. Biometrika, 68(1): 165–176. [Crossref], [Web of Science ®] , [Google Scholar]. In this context we consider some semiparametric approaches and a parametric method proposed by Fox-Taqqu[4] Fox, R. and Taqqu, M. S. 1986. Large-sample properties of parameter estimates for strongly dependent stationary gaussian time series. The Annals of Statistics, 14(2): 517–532. [Crossref], [Web of Science ®] , [Google Scholar]. We also investigate the method proposed by Robinson [5] Robinson, P. M. 1995a. Log-periodogram regression of time series with long range dependence. The Annals of Statistics, 23(3): 1048–1072. [Crossref], [Web of Science ®] , [Google Scholar] and a modification using the smoothed periodogram function.     

A Simple Corrected Score for Logistic Regression with Errors-in-Covariates

We develop a simple corrected score for logistic regression with errors-in-covariates. The new method is an extension of the consistent functional methods proposed by Huang and Wang (2001) and is closely related to the corrected score method by Nakamura (1990 Nakamura, T. (1990). Corrected score function for errors-in-variables models: Methodology and application to generalized linear models. Biometrika. 77:127–137.[Crossref], [Web of Science ®] , [Google Scholar]) and Stefanski (1989 Stefanski, L.A. (1989). Unbiased estimation of a nonlinear function a normal mean with application to measurement error models. Commun. Stat. Ser. A - Theory Methods. 18:4335–4358.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). The new method requires that the measurement error distribution is known, but does not require normality. The new method yields a consistent and asymptotically normal estimator under regularity conditions. We examine the finite-sample performance of the new estimator through simulation studies. Finally, we illustrate the new method by applying it to an AIDS study.     

Least absolute relative error estimation for functional quadratic multiplicative model

ABSTRACT

As an alternative to the functional quadratic model due to Yao and Müller (2010 Yao, F., Müller, H.-G. (2010). Functional quadratic regression. Biometrika 97:49–64.[Crossref], [Web of Science ®] , [Google Scholar]), we consider a functional quadratic multiplicative model. This multiplicative model provides a useful alternative when the relative error is considered for analyzing data with positive responses. The existing work for functional models are mainly based on absolute errors. The commonly used least squares criterion is just such an example. In many practical applications, however, people concern on the size of relative error rather than that of error itself. Therefore, the estimation procedure based on least absolute relative errors, which is proposed by Chen et al. (2010 Chen, K., Guo, S., Lin, Y., Ying, Z. (2010). Least absolute relative error estimation. J. Am. Stat. Assoc. 105:1104–1112.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) for the linear multiplicative model, is developed for functional quadratic multiplicative model. The asymptotic behaviors of the proposed estimators are established. Some simulation studies show that the estimation procedure has good prediction performance. Moreover, a real data set is analyzed for illustrating the proposed methods.     

Generalized linear regression models incorporating original outcome distributions

ABSTRACT

In this article, we propose an approach for incorporating continuous and discrete original outcome distributions into the usual exponential family regression models. The new approach is an extension of the works of Suissa (1991 Suissa, S. (1991). Binary methods for continuous outcomes: A parametric alternative. J. Clin. Epidemiol. 44:241–248.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) and Suissa and Blais (1995 Suissa, S., Blais, L. (1995). Binary regression with continuous outcomes. Stat. Med. 14:247–255.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]), which present methods to estimate the risk of an event defined in a sample subspace of an original continuous outcome variable. Simulation studies are presented in order to illustrate the performance of the developed methodology. Real data sets are analyzed by using the proposed models.     

Group Variable Selection with Oracle Property by Weight-Fused Adaptive Elastic Net Model for Strongly Correlated Data

This paper is the generalization of weight-fused elastic net (Fu and Xu, 2012 Fu, G., Xu, Q. (2012). Grouping variable selection by weight fused elastic net for multi-collinear data. Communications in Statistics-Simulation and Computation 41(2):205–221.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), which performs group variable selection by combining weight-fused LASSO(wfLasso) and elastic net (Zou and Hastie, 2005 Zou, H., Hastie, T. (2005). Regularization and variable selection via the elastic net. Journal of the Royal Statistical Society: Series B (Statistical Methodology) 67(2):301–320.[Crossref], [Web of Science ®] , [Google Scholar]) penalties. In this study, the elastic net penalty is replaced by adaptive elastic net penalty (AdaEnet) (Zou and Zhang, 2009 Zou, H., Zhang, H. (2009). On the adaptive elastic-net with a diverging number of parameters. Annals of Statistics 37(4):1733–1751.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]), and a new group variable selection algorithm with oracle property (Fan and Li, 2001 Fan, J., Li, R. (2001). Variable selection via nonconcave penalized likelihood and its oracle properties. Journal of the American Statistical Association 96(456):1348–1360.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]; Zou, 2006 Zou, H. (2006). The adaptive lasso and its oracle properties. Journal of the American Statistical Association 101(476):1418–1429.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) is obtained.     

A MONTE CARLO COMPARISON OF VARIOUS ASYMPTOTIC APPROXIMATIONS TO THE DISTRIBUTION OF INSTRUMENTAL VARIABLES ESTIMATORS

ABSTRACT

We examine empirical relevance of three alternative asymptotic approximations to the distribution of instrumental variables estimators by Monte Carlo experiments. We find that conventional asymptotics provides a reasonable approximation to the actual distribution of instrumental variables estimators when the sample size is reasonably large. For most sample sizes, we find Bekker[11] Bekker, P. A. 1994. Alternative Approximations to the Distributions of Instrumental Variable Estimators. Econometrica, 62: 657–681. [Crossref], [Web of Science ®] , [Google Scholar] asymptotics provides reasonably good approximation even when the first stage R ² is very small. We conclude that reporting Bekker[11] Bekker, P. A. 1994. Alternative Approximations to the Distributions of Instrumental Variable Estimators. Econometrica, 62: 657–681. [Crossref], [Web of Science ®] , [Google Scholar] confidence interval would suffice for most microeconometric (cross-sectional) applications, and the comparative advantage of Staiger and Stock[5] Staiger, D. and Stock, J. H. 1997. Instrumental Variables Regression with Weak Instruments. Econometrica, 65: 556–586. [Crossref], [Web of Science ®] , [Google Scholar] asymptotic approximation is in applications with sample sizes typical in macroeconometric (time series) applications.     

Multivariate Bayesian U-type asymmetric designs for non parametric response surface prediction under correlated errors

The objective of this paper is to study U-type designs for Bayesian non parametric response surface prediction under correlated errors. The asymptotic Bayes criterion is developed in terms of the asymptotic approach of Mitchell et al. (1994 Mitchell, T., Sacks, J., Ylvisaker, D. (1994). Asymptotic Bayes criteria for nonparametric response surface design. Ann. Stat. 22:634–651.[Crossref], [Web of Science ®] , [Google Scholar]) for a more general covariance kernel proposed by Chatterjee and Qin (2011 Chatterjee, K., Qin, H. (2011). Generalized discrete discrepancy and its applications in experimental designs. J. Stat. Plann. Inference 141:951–960.[Crossref], [Web of Science ®] , [Google Scholar]). A relationship between the asymptotic Bayes criterion and other criteria, such as orthogonality and aberration, is then developed. A lower bound for the criterion is also obtained, and numerical results show that this lower bound is tight. The established results generalize those of Yue et al. (2011 Yue, R.X., Qin, H., Chatterjee, K. (2011). Optimal U-type design for Bayesian nonparametric multiresponse prediction. J. Stat. Plann. Inference 141:2472–2479.[Crossref], [Web of Science ®] , [Google Scholar]) from symmetrical case to asymmetrical U-type designs.     

Fubini theorem for non additive measures in the framework of g-expectation

Since the seminal paper of Ghirardato (1997 Ghirardato, P. 1997. On the independence for non-additive measures, with a Fubini theorem. Journal of Economic Theory 73:261–91.[Crossref], [Web of Science ®] , [Google Scholar]), it is known that Fubini theorem for non additive measures can be available only for functions as “slice-comonotonic” in the framework of product algebra. Later, inspired by Ghirardato (1997 Ghirardato, P. 1997. On the independence for non-additive measures, with a Fubini theorem. Journal of Economic Theory 73:261–91.[Crossref], [Web of Science ®] , [Google Scholar]), Chateauneuf and Lefort (2008 Chateauneuf, A., and J. P. Lefort. 2008. Some Fubini theorems on product σ-algebras for non-additive measures. International Journal of Approximate Reasoning 48:686–96.[Crossref], [Web of Science ®] , [Google Scholar]) obtained some Fubini theorems for non additive measures in the framework of product σ-algebra. In this article, we study Fubini theorem for non additive measures in the framework of g-expectation. We give some different assumptions that provide Fubini theorem in the framework of g-expectation.     

A two-stage design for multivariate estimation of proportions

ABSTRACT

In this article, the problem of estimation of the several proportions of two inter-dependent sensitive characteristics prevailing in a given population is considered. A new two-stage randomized response model has been proposed by extending Lee et al.’s (2013 Lee, C., Sedory, S., Singh, S (2013). Estimation at least seven measures of qualitative variables from a single sample using randomized response technique. Stat. Probab. Lett. 83(1):399–409.[Crossref], [Web of Science ®] , [Google Scholar]) model. The expressions for biases, variances, and mean square errors of different estimators have been provided. The relative efficiency comparisons of the proposed estimators has been made, numerically by considering the different practicable choices of the design parameters. It was observed that the proposed extended version performs efficiently than simple and crossed models of Lee et al. (2013 Lee, C., Sedory, S., Singh, S (2013). Estimation at least seven measures of qualitative variables from a single sample using randomized response technique. Stat. Probab. Lett. 83(1):399–409.[Crossref], [Web of Science ®] , [Google Scholar]).     

Further Results for a General Family of Bivariate Copulas

A bivariate family of copulas has been initiated by Cuadras-Augé (1981 Cuadras, C.M., Augé, J. (1981). A continuous general multivariate distribution and its properties. Commun. Statist. (A) Theor. Meth. 10:339–353.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) and Marshall (1996 Marshall, A.W. (1996). Copulas, marginals, and joint distributions. In: Distributions with fixed marginals and related topics. IMS Lecture Notes Monogr. Ser. 28:213–222.[Crossref] , [Google Scholar]). Recently, Durante (2007 Durante, F. (2007). A new family of symmetric bivariate copulas. C. R. Math. Acad. Sci. Paris 344:195–198.[Crossref], [Web of Science ®] , [Google Scholar]) considered this family as a general family of symmetric bivariate copulas indexed by a generator function and studied some of its dependence properties. In this article, we obtain and describe further aspects of dependence for this family. For example, we have proved that the family has positive likelihood ratio dependence structure if and only if the family reduces to some well-known copulas. We also derive several proper forms for the generator function of this family. Considering a multivariate extension of the bivariate family of copulas provided by Durante et al. (2007 Durante, F., Quesada-Molina, J.J., Flores, M. (2007). On a family of multivariate copulas for aggregation processes. Inform. Sci. 177(24):5715–5724.[Crossref], [Web of Science ®] , [Google Scholar]), some dependence properties are studied. Finally, some positive dependence stochastic orderings for two random vectors having a copula from the proposed families, are discussed.     

A quantile-based study on ageing intensity function

Jiang, Ji, and Xiao (2003 Jiang, R., P. Ji, and X. Xiao. 2003. Aging property of unimodal failure rate models. Reliability Engineering and System Safety 79(1):113–6.[Crossref], [Web of Science ®] , [Google Scholar]) has introduced a quantitative measure known as the ageing intensity function for evaluating the ageing properties of a component/system. In recent years, there has been a great interest on the study of quantile function, an equivalent alternative to the distribution function approach. Unlike the distribution function approach, the quantile method possess some unique properties (see Gilchrist 2000 Gilchrist, W. 2000. Statistical modelling with quantile functions. Boca Raton, Florida: Chapman and Hall/CRC.[Crossref] , [Google Scholar], Nair, Sankaran, and Balakrishnan 2013 Nair N. U., P. G. Sankaran, N. Balakrishnan. 2013. Quantile-based reliability concepts. In: Quantile-Based Reliability Analysis. Statistics for Industry and Technology. New York, NY: Birkhäuser.[Crossref] , [Google Scholar]). Motivated with this, in the present paper we introduce a quantile-based ageing intensity function and study its various ageing properties. We also study some stochastic comparison of random variables based on the proposed measure.