Bivariate Continuous Distributions with Specified Conditional Hazard Functions

Analysis of covariance (ANCOVA) is the standard procedure for comparing several treatments when the response variable depends on one or more covariates. We consider the problem of testing the equality of treatment effects when the variances are not assumed to be equal. It is well known that classical F test is not robust with respect to the assumption of equal variances and may lead to misleading conclusions if the variances are not equal. Ananda (1998 Ananda , M. M. A. ( 1998 ). Bayesian and non-Bayesian solutions to analysis of covariance models under heteroscedasticity . J. Econometrics 86 : 177 – 192 .[Crossref], [Web of Science ®] , [Google Scholar]) developed a generalized F test for testing the equality of treatment effects. However, simulation studies show that the actual size of this test can be much higher than the nominal level when the sample sizes are small, particularly when the number of treatments is large. In this article, we develop a test using the parametric bootstrap approach of Krishnamoorthy et al. (2007 Krishnamoorthy , K. , Lu , F. , Mathew , T. ( 2007 ). A parametric bootstrap approach for ANOVA with unequal variances: Fixed and random models . Computat. Statist. Data Anal. 51 : 5731 – 5742 .[Crossref], [Web of Science ®] , [Google Scholar]). Our simulations show that the actual size of our proposed test is close to the nominal level, irrespective of the number of treatments and sample sizes. Our simulations also indicate that our proposed PB test is more robust, with respect to the assumption of normality, than the generalized F test. Therefore, our proposed PB test provides a satisfactory alternative to the generalized F test.     

A stratified unrelated question randomized response model using Neyman allocation

This paper applies stratified random sampling using Neyman allocation to Mangat et al. (1992 Mangat, N.S., Singh, R., Singh, S. (1992). An improved unrelated question randomized response strategy. Cal. Stat. Assoc. Bull. 42:277–281.[Crossref] , [Google Scholar]) unrelated question randomized response (RR) strategy for both completely truthful reporting and less than completely truthful reporting. It is shown that, for the prior information given, our new model is more efficient in terms of variance (in the case of completely truthful reporting) and mean square error (in terms of less than completely truthful reporting) than Kim and Elam's (2007 Kim, J.M., Elam, M.E. (2007). A stratified unrelated question randomized response model. Stat. Papers 48:215–233.[Crossref], [Web of Science ®] , [Google Scholar]) model. Numerical illustrations and graphs are also given in support of the present study.     

Bootstrap methods for multivariate hypothesis testing

The nonparametric and parametric bootstrap methods for multivariate hypothesis testing are developed. They are used to approximate the null distribution of the test statistics proposed by Duchesne and Francq (2015 Duchesne, P., Francq, C. (2015). Multivariate hypothesis testing using generalized and {2}-inverses—with applications. Statistics 49:475–496.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]), resulting in bootstrap testing procedures. In the problem of testing for the mean vector of a multivariate distribution, the asymptotic validity of the bootstrap methods is proved. The finite sample performance of the new solutions is demonstrated by means of Monte Carlo simulation studies. They indicate that for small-sample size, the bootstrap tests provide a better finite sample properties than the asymptotic tests considered by Duchesne and Francq (2015 Duchesne, P., Francq, C. (2015). Multivariate hypothesis testing using generalized and {2}-inverses—with applications. Statistics 49:475–496.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]).     

A Modification of a Percentile Estimation Procedure Based on Generalized Polya Urns

Adaptive designs find an important application in the estimation of unknown percentiles for an underlying dose-response curve. A nonparametric adaptive design was suggested by Mugno et al. (2004 Mugno, R.A., Zhus, W., Rosenberger, W.F. (2004). Adaptive urn designs for estimating several percentiles of a dose-response curve. Statist. Med. 23(13):2137–2150.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) to simultaneously estimate multiple percentiles of an unknown dose-response curve via generalized Polya urns. In this article, we examine the properties of the design proposed by Mugno et al. (2004 Mugno, R.A., Zhus, W., Rosenberger, W.F. (2004). Adaptive urn designs for estimating several percentiles of a dose-response curve. Statist. Med. 23(13):2137–2150.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) when delays in observing responses are encountered. Using simulations, we evaluate a modification of the design under varying group sizes. Our results demonstrate unbiased estimation with minimal loss in efficiency when compared to the original compound urn design.     

On generalized exponential chain ratio estimators under stratified two-phase random sampling

Sanaullah et al. (2014 Sanaullah, A., Ali, H.M., Noor ul Amin, M., Hanif, M. (2014). Generalized exponential chain ratio estimators under stratified two-phase random sampling. Appl. Math. Comput. 226:541–547.[Crossref], [Web of Science ®] , [Google Scholar]) have suggested generalized exponential chain ratio estimators under stratified two-phase sampling scheme for estimating the finite population mean. However, the bias and mean square error (MSE) expressions presented in that work need some corrections, and consequently the study based on efficiency comparison also requires corrections. In this article, we revisit Sanaullah et al. (2014 Sanaullah, A., Ali, H.M., Noor ul Amin, M., Hanif, M. (2014). Generalized exponential chain ratio estimators under stratified two-phase random sampling. Appl. Math. Comput. 226:541–547.[Crossref], [Web of Science ®] , [Google Scholar]) estimator and provide the correct bias and MSE expressions of their estimator. We also propose an estimator which is more efficient than several competing estimators including the classes of estimators in Sanaullah et al. (2014 Sanaullah, A., Ali, H.M., Noor ul Amin, M., Hanif, M. (2014). Generalized exponential chain ratio estimators under stratified two-phase random sampling. Appl. Math. Comput. 226:541–547.[Crossref], [Web of Science ®] , [Google Scholar]). Three real datasets are used for efficiency comparisons.     

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.     

Optimal allocations of coverage limits for two independent random losses of insurance policy

This article recasts the optimal allocations of coverage limits for two independent random losses. Under some regularity conditions on the two concerned probability density functions, we build the sufficient and necessary condition for the existence of the optimal allocation of coverage limits, and derive the optimal allocation whenever they do exist. The results supplement Lu and Meng (2011 Lu, Z.Y., Meng, L.L. (2011). Stochastic comparisons for allocations of upper limits and deductibles with applications. Insur.: Math. Econ. 48:338–343.[Crossref], [Web of Science ®] , [Google Scholar], Proposition 5.2) and Hu and Wang (2014 Hu, S., Wang, R. (2014). Stochastic comparisons and optimal allocation for policy limits and deductibles. Commun. Stat. – Theory Methods 43:151–164.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar], Theorem 5.1).     

A two-stage Land et al.'s randomized response model for estimating a rare sensitive attribute using Poisson distribution

The crux of this paper is to estimate the mean of the number of persons possessing a rare sensitive attribute based on the Mangat (1992 Mangat, N.S. (1992). Two stage reandomized response sampling procedure using unrelated question. J. Ind. Soc. Agric. Stat. 44(1):82–87. [Google Scholar]) randomization device by utilizing the Poisson distribution in survey sampling. It is shown that the proposed model is more efficient than Land et al. (2011 Land, M., Singh, S., Sedory, S.A. (2011). Estimation of a rare attribute using Poisson distribution. Statistics doi:10.1080/02331888.2010.524300[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) when the proportion of persons possessing a rare unrelated attribute is known. Properties of the proposed randomized response model have been studied along with recommendations. We have also extended the proposed model to stratified random sampling on the lines of Lee et al. (2013 Lee, G.S., Uhm, D., Kim, J.M. (2013). Estimation of a rare sensitive attribute in stratified sampling using Poisson distribution. Statistics 47(3):575–589.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]). It has been also shown that the proposed estimator is better than Lee et al.'s (2013 Lee, G.S., Uhm, D., Kim, J.M. (2013). Estimation of a rare sensitive attribute in stratified sampling using Poisson distribution. Statistics 47(3):575–589.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) estimator. Numerical illustrations are also given in support of the present study.     

Hill's estimator under weak dependence

In this article, we investigate the asymptotic normality of the Hill's estimator of the tail index parameter, when the observations are weakly dependent in the sense of Doukhan and Louhichi (1999 Doukhan, P., Louhichi, S. (1999). A new weak dependence condition and applications to moment inequalities. Stochastic Process. Appl. 84:313–342.[Crossref], [Web of Science ®] , [Google Scholar]) and are drawn from a strictly linear process. We show that the previous result on Hill estimator obtained by Rootzen et al. (1990 Rootzen, H., Leadbetter, M., De Haan, L. (1990). Tail and quantile estimation for strongly mixing stationary sequences. Technical report. No. 292, Center for Stochastic Processes, Department of Statistics, University of North Carolina, Chapel Hill. [Google Scholar]) and Resnick and Starica (1997 Resnick, S., Starica, C. (1997). Asymptotic behavior of Hill's estimator for autoregressive data. Commun. Statistics-stochastic Models 13:703–723.[Taylor &; Francis Online] , [Google Scholar]) for strong mixing can be extended to weak dependence.     

An optional randomized response model for estimating a rare sensitive attribute using Poisson distribution

The crux of this article is to estimate the mean of the number of persons possessing a rare sensitive attribute based on the Mangat (1991 Mangat, N.S. (1991). An optional randomized response sampling technique using non–stigmatized attribute. Statistica. 51(4):595–602. [Google Scholar]) randomization device by utilizing the Poisson distribution in simple random sampling and stratified sampling. Properties of the proposed randomized response (RR) model have been studied along with recommendations. It is also shown that the proposed model is more efficient than that of Land et al. (2011 Land, M., Singh, S., Sedory, S.A. (2011). Estimation of a rare attribute using Poisson distribution. Statistics. 46(3):351–360.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) in simple random sampling and that of Lee et al. (2013 Lee, G.S., Uhm, D., Kim, J.M. (2013). Estimation of a rare sensitive attribute in stratified sampling using Poisson distribution. Statistics. 47(3):575–589.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) in stratified random sampling when the proportion of persons possessing a rare unrelated attribute is known. Numerical illustrations are also given in support of the present study.     

Dynamic version of weighted cumulative residual entropy

Some extensions of Shannon entropy to the survival function have been recently proposed. Misagh et al. (2011 Misagh, F., Panahi, Y., Yari, G.H., Shahi, R. (2011, September). Weighted cumulative entropy and its estimation. In: Quality and Reliability (ICQR), 2011, IEEE International conference (pp. 477–480), IEEE.[Crossref] , [Google Scholar]) introduced weighted cumulative residual entropy (WCRE) that was studied more by Mirali et al. (2015 Mirali, M., Baratpour, S., Fakoor, V. (2015). On weighted cumulative residual entropy. Commun. Stat. Theory Methods. doi:10.1080103610926.2015.1053932.[Web of Science ®] , [Google Scholar]). In this article, the dynamic version of WCRE is proposed. Some relationships of this measure with well-known reliability measures and ageing classes are studied and some characterization results for exponential and Rayleigh distributions are provided. Also, a non parametric estimation of dynamic version of WCRE is introduced and its asymptotic behavior is investigated.     

Improved Exponential Type Estimators of the Mean of a Sensitive Variable in the Presence of Nonsensitive Auxiliary Information

Recently, Koyuncu et al. (2013 Koyuncu, N., Gupta, S., Sousa, R. (2014). Exponential type estimators of the mean of a sensitive variable in the presence of non-sensitive auxiliary information. Communications in Statistics- Simulation and Computation[PubMed], [Web of Science ®] , [Google Scholar]) proposed an exponential type estimator to improve the efficiency of mean estimator based on randomized response technique. In this article, we propose an improved exponential type estimator which is more efficient than the Koyuncu et al. (2013 Koyuncu, N., Gupta, S., Sousa, R. (2014). Exponential type estimators of the mean of a sensitive variable in the presence of non-sensitive auxiliary information. Communications in Statistics- Simulation and Computation[PubMed], [Web of Science ®] , [Google Scholar]) estimator, which in turn was shown to be more efficient than the usual mean estimator, ratio estimator, regression estimator, and the Gupta et al. (2012 Gupta, S., Shabbir, J., Sousa, R., Corte-Real, P. (2012). Regression estimation of the mean of a sensitive variable in the presence of auxiliary information. Communications in Statistics – Theory and Methods 41:2394–2404.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) estimator. Under simple random sampling without replacement (SRSWOR) scheme, bias and mean square error expressions for the proposed estimator are obtained up to first order of approximation and comparisons are made with the Koyuncu et al. (2013 Koyuncu, N., Gupta, S., Sousa, R. (2014). Exponential type estimators of the mean of a sensitive variable in the presence of non-sensitive auxiliary information. Communications in Statistics- Simulation and Computation[PubMed], [Web of Science ®] , [Google Scholar]) estimator. A simulation study is used to observe the performances of these two estimators. Theoretical findings are also supported by a numerical example with real data. We also show how to, extend the proposed estimator to the case when more than one auxiliary variable is available.     

Bayesian Estimation for Linear Functions of Poisson Rates

The probability matching prior for linear functions of Poisson parameters is derived. A comparison is made between the confidence intervals obtained by Stamey and Hamilton (2006 Stamey, J., Hamilton, C. (2006). A note on confidence intervals for a linear function of Poisson rates. Commun. Statist. Simul. &; Computat. 35(4):849–856.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]), and the intervals derived by us when using the Jeffreys’ and probability matching priors. The intervals obtained from the Jeffreys’ prior are in some cases fiducial intervals (Krishnamoorthy and Lee, 2010 Krishnamoorthy, K., Lee, M. (2010). Inference for functions of parameters in discrete distributions based on fiducial approach: Binomial and Poisson cases. J. Statist. Plann. Infere. 140(5):1182–1192.[Crossref], [Web of Science ®] , [Google Scholar]). A weighted Monte Carlo method is used for the probability matching prior. The power and size of the test, using Bayesian methods, is compared to tests used by Krishnamoorthy and Thomson (2004 Krishnamoorthy, K., Thomson, J. (2004). A more powerful test for comparing two Poisson means. J. Statist. Plann. Infere. 119(1):23–35.[Crossref], [Web of Science ®] , [Google Scholar]). The Jeffreys’, probability matching and two other priors are used.     

On extreme order statistics from heterogeneous beta distributions with applications

This paper treats the problem of stochastic comparisons for the extreme order statistics arising from heterogeneous beta distributions. Some sufficient conditions involved in majorization-type partial orders are provided for comparing the extreme order statistics in the sense of various magnitude orderings including the likelihood ratio order, the reversed hazard rate order, the usual stochastic order, and the usual multivariate stochastic order. The results established here strengthen and extend those including Kochar and Xu (2007 Kochar, S.C., Xu, M. (2007). Stochastic comparisons of parallel systems when components have proportional hazard rates. Probab. Eng. Inf. Sci. 21:597–609.[Crossref], [Web of Science ®] , [Google Scholar]), Mao and Hu (2010 Mao, T., Hu, T. (2010). Equivalent characterizations on orderings of order statistics and sample ranges. Probab. Eng. Inf. Sci. 24:245–262.[Crossref], [Web of Science ®] , [Google Scholar]), Balakrishnan et al. (2014 Balakrishnan, N., Barmalzan, G., Haidari, A. (2014). On usual multivariate stochastic ordering of order statistics from heterogeneous beta variables. J. Multivariate Anal. 127:147–150.[Crossref], [Web of Science ®] , [Google Scholar]), and Torrado (2015 Torrado, N. (2015). On magnitude orderings between smallest order statistics from heterogeneous beta distributions. J. Math. Anal. Appl. 426:824–838.[Crossref], [Web of Science ®] , [Google Scholar]). A real application in system assembly and some numerical examples are also presented to illustrate the theoretical results.     

Some extensions of circular balanced and circular strongly balanced repeated measurements designs

Repeated measurement designs are widely used in medicine, pharmacology, animal sciences, and psychology. In this paper the works of Iqbal and Tahir (2009 Iqbal, I., and M. H. Tahir. 2009. Circular strongly balanced repeated measurements designs. Communications in Statistics—Theory and Methods 38:3686–96.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) and Iqbal, Tahir, and Ghazali (2010 Iqbal, I., M. H. Tahir, and S. S. A. Ghazali. 2010. Circular first- and second-order balanced repeated measurements designs. Communications in Statistics—Theory and Methods 39:228–40.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) are generalized for the construction of circular-balanced and circular strongly balanced repeated measurements designs through the method of cyclic shifts for three periods.     

An Approximation Model of the Collective Risk Model with INAR(1) Claim Process

Cossette et al. (2010 Cossette, H., Marceau, E., Maume-Deschamps, V. (2010). Discerte-time risk models based on time series for count random variables. ASTIN Bull. 40:123–150.[Crossref], [Web of Science ®] , [Google Scholar], 2011 Cossette, H., Marceau, E., Toureille, F. (2011). risk models based on time series for count random variables. Insur. Math. Econ. 48:19–28.[Crossref], [Web of Science ®] , [Google Scholar]) gave a novel collective risk model where the total numbers of claims satisfy the first-order integer-valued autoregressive process. For a risk model, it is interesting to investigate the upper bound of ruin probability. However, the loss increments of the above model are dependent; it is difficult to derive the upper bound of ruin probability. In this article, we propose an approximation model with stationary independent increments. The upper bound of ruin probability and the adjustment coefficient are derived. The approximation model is illustrated via four simulated examples. Results show that the gap of the approximation model and dependent model can be ignored by adjusting values of parameters.     

Using the Inverse Maximum Ratio-Λ as a Technique to Quantify Surface Uniformity

Techniques used in variability assessment are subsequently used to draw conclusions regarding the “spread”/uniformity of data curves. Due to the limitations of these techniques, they are not adequate for circumstances where data manifest with multiple peaks. Examples of these manifestations (in three-dimensional space) include under-foot pressure distributions recorded for different types of footwear (Becerro-de-Bengoa-Vallejo et al., 2014 Biau, D.J. (2011). In brief: Standard deviation and standard error. Clinical Orthopaedics and Related Research 469(9):2661–2664.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]; Cibulka et al., 1994 Cibulka, M.T., Sinacore, D.R., Mueller, M.J. (1994). Shin splints and forefoot contact running: A case report. Journal of Orthopaedic &; Sports Physical Therapy 20(2):98–102.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]; Davies et al., 2003 Davies, M.B., Betts, R.P., Scott, I.R. (2003). Optical plantar pressure analysis following internal fixation for displaced intra-articular os calcis fractures. Foot &; Ankle International 24(11):851–856.[PubMed], [Web of Science ®] , [Google Scholar]), surface textures and interfaces designed to impact friction, and and and molecular surface structures such as viral epitopes (Torras and Garcia-Valls, 2004 Torras, C., Garcia-Valls, R. (2004). Quantification of membrane morphology by interpretation of scanning electron microscopy images. Journal of Membrane Science 233(1–2):119–127.[Crossref], [Web of Science ®] , [Google Scholar]; Pacejka, 1997; Fustaffson, 1997). This article proposes a technique for generating a single variable – Λ that will quantify the uniformity of such surfaces. We define and validate this technique using several mathematical and graphical models.     

Tuning Parameter Selection in Penalized Frailty Models

The penalized likelihood approach of Fan and Li (2001 Fan, J., Li, R. (2001). Variable selection via nonconcave penalized likelihood and its oracle properties. Journal of the American Association 96:1348–1360.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar], 2002 Fan, J., Li, R. (2002). Variable selection for Cox’s proportional hazards model and frailty model. The Annals of Statistics 30:74–99.[Crossref], [Web of Science ®] , [Google Scholar]) differs from the traditional variable selection procedures in that it deletes the non-significant variables by estimating their coefficients as zero. Nevertheless, the desirable performance of this shrinkage methodology relies heavily on an appropriate selection of the tuning parameter which is involved in the penalty functions. In this work, new estimates of the norm of the error are firstly proposed through the use of Kantorovich inequalities and, subsequently, applied to the frailty models framework. These estimates are used in order to derive a tuning parameter selection procedure for penalized frailty models and clustered data. In contrast with the standard methods, the proposed approach does not depend on resampling and therefore results in a considerable gain in computational time. Moreover, it produces improved results. Simulation studies are presented to support theoretical findings and two real medical data sets are analyzed.     

Multiple dependent state repetitive sampling plans based on one-sided process capability indices

In this research, multiple dependent state and repetitive group sampling are used to design a variable sampling plan based on one-sided process capability indices, which consider the quality of the current lot as well as the quality of the preceding lots. The sample size and critical values of the proposed plan are determined by minimizing the average sample number while satisfying the producer's risk and consumer's risk at corresponding quality levels. In addition, comparisons are made with the existing sampling plans [Pearn and Wu (2006a Pearn, W. L., and C. W. Wu. 2006a. Critical acceptance values and sample sizes of a variables sampling plan for very low fraction of defectives. Omega: International Journal of Management Science 34 (1):90–101.[Crossref], [Web of Science ®] , [Google Scholar]), Yen et al. (2015 Yen, C. H., C. H. Chang, and M. Aslam. 2015. Repetitive variable acceptance sampling plan for one-sided specification. Journal of Statistical Computation and Simulation 85 (6):1102–16.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar])] in terms of average sample number and operating characteristic curve. Finally, an example is provided to illustrate the proposed plan.     

Buffered Autoregressive Models With Conditional Heteroscedasticity: An Application to Exchange Rates

This article introduces a new model called the buffered autoregressive model with generalized autoregressive conditional heteroscedasticity (BAR-GARCH). The proposed model, as an extension of the BAR model in Li et al. (2015 Li, G.D., Guan, B., Li, W.K., and Yu, P. L.H. (2015), “Hysteretic Autoregressive Time Series Models,” Biometrika, 102, 717–723.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]), can capture the buffering phenomena of time series in both the conditional mean and variance. Thus, it provides us a new way to study the nonlinearity of time series. Compared with the existing AR-GARCH and threshold AR-GARCH models, an application to several exchange rates highlights the importance of the BAR-GARCH model.