Matrix mean squared error comparisons of some biased estimators with two biasing parameters

To deal with multicollinearity problem, the biased estimators with two biasing parameters have recently attracted much research interest. The aim of this article is to compare one of the last proposals given by Yang and Chang (2010 Yang, H., and X. Chang. 2010. A new two-parameter estimator in linear regression. Communications in Statistics: Theory and Methods 39 (6):923–34.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) with Liu-type estimator (Liu 2003 Liu, K. 2003. Using Liu-type estimator to combat collinearity. Communications in Statistics: Theory and Methods 32 (5):1009–20.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) and k ? d class estimator (Sakallioglu and Kaciranlar 2008 Sakallioglu, S., and S. Kaciranlar. 2008. A new biased estimator based on ridge estimation. Statistical Papers 49:669–89.[Crossref], [Web of Science ®] , [Google Scholar]) under the matrix mean squared error criterion. As well as giving these comparisons theoretically, we support the results with the extended simulation studies and real data example, which show the advantages of the proposal given by Yang and Chang (2010 Yang, H., and X. Chang. 2010. A new two-parameter estimator in linear regression. Communications in Statistics: Theory and Methods 39 (6):923–34.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) over the other proposals with increasing multicollinearity level.     

On weighted cumulative residual entropy

In analogy with the weighted Shannon entropy proposed by Belis and Guiasu (1968 Belis, M., Guiasu, S. (1968). A quantitative-qualitative measure of information in cybernetic systems. IEEE Trans. Inf. Th. IT-4:593–594.[Crossref], [Web of Science ®] , [Google Scholar]) and Guiasu (1986 Guiasu, S. (1986). Grouping data by using the weighted entropy. J. Stat. Plann. Inference 15:63–69.[Crossref], [Web of Science ®] , [Google Scholar]), we introduce a new information measure called weighted cumulative residual entropy (WCRE). This is based on the cumulative residual entropy (CRE), which is introduced by Rao et al. (2004 Rao, M., Chen, Y., Vemuri, B.C., Wang, F. (2004). Cumulative residual entropy: a new measure of information. IEEE Trans. Info. Theory 50(6):1220–1228.[Crossref], [Web of Science ®] , [Google Scholar]). This new information measure is “length-biased” shift dependent that assigns larger weights to larger values of random variable. The properties of WCRE and a formula relating WCRE and weighted Shannon entropy are given. Related studies of reliability theory is covered. Our results include inequalities and various bounds to the WCRE. Conditional WCRE and some of its properties are discussed. The empirical WCRE is proposed to estimate this new information measure. Finally, strong consistency and central limit theorem are provided.     

An objective Bayesian estimation in a two-period crossover design

Two-period crossover design is one of the commonly used designs in clinical trials. But, the estimation of treatment effect is complicated by the possible presence of carryover effect. It is known that ignoring the carryover effect when it exists can lead to poor estimates of the treatment effect. The classical approach by Grizzle (1965 Grizzle, J.E. (1965). The two-period change-over design and its use in clinical trials. Biometrics 21:467–480. See Grizzle (1974) for corrections.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) consists of two stages. First, a preliminary test is conducted on carryover effect. If the carryover effect is significant, analysis is based only on data from period one; otherwise, analysis is based on data from both periods. A Bayesian approach with improper priors was proposed by Grieve (1985 Grieve, A.P. (1985). A Bayesian analysis of the two-period crossover design for clinical trials. Biometrics 41:979–990.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) which uses a mixture of two models: a model with carryover effect and another without. The indeterminacy of the Bayes factor due to the arbitrary constant in the improper prior was addressed by assigning a minimally discriminatory value to the constant. In this article, we present an objective Bayesian estimation approach to the two-period crossover design which is also based on a mixture model, but using the commonly recommended Zellner–Siow g-prior. We provide simulation studies and a real data example and compare the numerical results with Grizzle (1965 Grizzle, J.E. (1965). The two-period change-over design and its use in clinical trials. Biometrics 21:467–480. See Grizzle (1974) for corrections.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar])’s and Grieve (1985 Grieve, A.P. (1985). A Bayesian analysis of the two-period crossover design for clinical trials. Biometrics 41:979–990.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar])’s approaches.     

Test of Independence for Baker’s Bivariate Distributions

Baker (2008 Baker, R. (2008). An order-statistics-based method for constructing multivariate distributions with fixed marginals. Journal of Multivariate Analysis 99: 2312–2327.[Crossref], [Web of Science ®] , [Google Scholar]) introduced a new method for constructing multivariate distributions with given marginals based on order statistics. In this paper, we provide a test of independence for a pair of absolutely continuous random variables (X, Y) jointly distributed according to Baker’s bivariate distributions. Our purpose is to test the hypothesis that X and Y are independent versus the alternative that X and Y are positively (negatively) quadrant dependent. The asymptotic distribution of the proposed test statistic is investigated. Also, the powers of the proposed test and the class of distribution-free tests proposed by Kochar and Gupta (1987 Kochar, S. G., Gupta, R. P. (1987). Competitors of Kendall-tau test for testing independence against positive quadrant dependence. Biometrika 74(3): 664–666.[Crossref], [Web of Science ®] , [Google Scholar]) are compared empirically via a simulation study.     

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.     

An efficient use of moment's ratios of scrambling variables in a randomized response technique

This paper aimed at providing an efficient new unbiased estimator for estimating the proportion of a potentially sensitive attribute in survey sampling. The suggested randomization device makes use of the means, variances of scrambling variables, and the two scalars lie between “zero” and “one.” Thus, the same amount of information has been used at the estimation stage. The variance formula of the suggested estimator has been obtained. We have compared the proposed unbiased estimator with that of Kuk (1990 Kuk, A.Y.C. (1990). Asking sensitive questions inderectely. Biometrika 77:436–438.[Crossref], [Web of Science ®] , [Google Scholar]) and Franklin (1989 Franklin, L.A. (1989). A comparision of estimators for randomized response sampling with continuous distribution s from a dichotomous population. Commun. Stat. Theor. Methods 18:489–505.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]), and Singh and Chen (2009 Singh, S., Chen, C.C. (2009). Utilization of higher order moments of scrambling variables in randomized response sampling. J. Stat. Plann. Inference. 139:3377–3380.[Crossref], [Web of Science ®] , [Google Scholar]) estimators. Relevant conditions are obtained in which the proposed estimator is more efficient than Kuk (1990 Kuk, A.Y.C. (1990). Asking sensitive questions inderectely. Biometrika 77:436–438.[Crossref], [Web of Science ®] , [Google Scholar]) and Franklin (1989 Franklin, L.A. (1989). A comparision of estimators for randomized response sampling with continuous distribution s from a dichotomous population. Commun. Stat. Theor. Methods 18:489–505.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) and Singh and Chen (2009 Singh, S., Chen, C.C. (2009). Utilization of higher order moments of scrambling variables in randomized response sampling. J. Stat. Plann. Inference. 139:3377–3380.[Crossref], [Web of Science ®] , [Google Scholar]) estimators. The optimum estimator (OE) in the proposed class of estimators has been identified which finally depends on moments ratios of the scrambling variables. The variance of the optimum estimator has been obtained and compared with that of the Kuk (1990 Kuk, A.Y.C. (1990). Asking sensitive questions inderectely. Biometrika 77:436–438.[Crossref], [Web of Science ®] , [Google Scholar]) and Franklin (1989 Franklin, L.A. (1989). A comparision of estimators for randomized response sampling with continuous distribution s from a dichotomous population. Commun. Stat. Theor. Methods 18:489–505.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) estimator and Singh and Chen (2009 Singh, S., Chen, C.C. (2009). Utilization of higher order moments of scrambling variables in randomized response sampling. J. Stat. Plann. Inference. 139:3377–3380.[Crossref], [Web of Science ®] , [Google Scholar]) estimator. It is interesting to mention that the “optimum estimator” of the class of estimators due to Singh and Chen (2009 Singh, S., Chen, C.C. (2009). Utilization of higher order moments of scrambling variables in randomized response sampling. J. Stat. Plann. Inference. 139:3377–3380.[Crossref], [Web of Science ®] , [Google Scholar]) depends on the parameter π under investigation which limits the use of Singh and Chen (2009 Singh, S., Chen, C.C. (2009). Utilization of higher order moments of scrambling variables in randomized response sampling. J. Stat. Plann. Inference. 139:3377–3380.[Crossref], [Web of Science ®] , [Google Scholar]) OE in practice while the proposed OE in this paper is free from such a constraint. The proposed OE depends only on the moments ratios of scrambling variables. This is an advantage over the Singh and Chen (2009 Singh, S., Chen, C.C. (2009). Utilization of higher order moments of scrambling variables in randomized response sampling. J. Stat. Plann. Inference. 139:3377–3380.[Crossref], [Web of Science ®] , [Google Scholar]) estimator. Numerical illustrations are given in the support of the present study when the scrambling variables follow normal distribution. Theoretical and empirical results are very sound and quite illuminating in the favor of the present study.     

Analysis of Order Statistics from Distributions with Regularly Varying Tails

In this article, necessary conditions for comparing order statistics from distributions with regularly varying tails are discussed in terms of various stochastic orders. A necessary and sufficient condition for stochastically comparing tail behaviors of order statistics is derived. The main results generalize and recover some results in Kleiber (2002 Kleiber, C. (2002). Variability ordering of heavy-tailed distributions with applications to order statistics. Statist. Probab. Lett. 58:381388.[Crossref], [Web of Science ®] , [Google Scholar], 2004 Kleiber, C. (2004). Lorenz ordering of order statistics from log-logistic and related distributions. J. Statist. Plann. Infer. 120:2004.[Crossref], [Web of Science ®] , [Google Scholar]). Extensions to coherent systems are mentioned as well.     

Model selection using union-intersection principle for non nested models

In this article, we have extended the Vuong’s (1989 Vuong, Q.H. (1989). Likelihood ratio tests for model selection and non-nested hypothesis. Econometrica. 57:307–333.[Crossref], [Web of Science ®] , [Google Scholar]) model selection test to three models in accordance to union-intersection principle. Using the Kullback–Leibler criterion to measure the closeness of a model to the truth, we propose a simple likelihood ratio-based statistics for testing the null hypothesis that the competing models are equally close to the true data-generating process against the alternative hypothesis that at least one model is closer. We show that the distribution of the test statistic is asymptotically equal to the distribution of the maximum of dependent random variables with bivariate folded standard normal distribution. The density function of the maximum of dependent random variables with elliptically contoured distributions has been obtained by other researchers, but, not for distributions which do not belong to the elliptically contoured distributions family. In this article, the exact distribution of the maximum of dependent random variables with bivariate folded standard normal distribution is calculated as an asymptotic distribution of the proposed test statistic. The test is directional and is derived successively for the cases where the competing models are non nested and whether three, two, one, or none of them are misspecified.     

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.     

Bimodal symmetric-asymmetric power-normal families

This article proposes new symmetric and asymmetric distributions applying methods analogous as the ones in Kim (2005 Kim, H.J. (2005). On a class of two-piece skew-normal distributions. Statist.: J. Theoret. Appl. Statist. 39:537–553.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) and Arnold et al. (2009 Arnold, B.C., H.W. Gómez, and H.S. Salinas. (2009). On multiple constraint skewed models. Statist. J. Theoret. Appl. Statist. 43: 279–293.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) to the exponentiated normal distribution studied in Durrans (1992 Durrans, S.R. (1992). Distributions of fractional order statistics in hydrology. Water Resour. Res. 28:1649–1655.[Crossref], [Web of Science ®] , [Google Scholar]), that we call the power-normal (PN) distribution. The proposed bimodal extension, the main focus of the paper, is called the bimodal power-normal model and is denoted by BPN(α) model, where α is the asymmetry parameter. The authors give some properties including moments and maximum likelihood estimation. Two important features of the model proposed is that its normalizing constant has closed and simple form and that the Fisher information matrix is nonsingular, guaranteeing large sample properties of the maximum likelihood estimators. Finally, simulation studies and real applications reveal that the proposed model can perform well in both situations.     

Recurrence relations for single and product moments of generalized order statistics from Marshall–Olkin extended Pareto distributions

Marshall and Olkin (1997 Marshall, A.W., Olkin, I. (1997). A new method for adding a parameter to a family of distributions with application to the exponential and Weibull families. Biometrika 84(3):641–652.[Crossref], [Web of Science ®] , [Google Scholar]) introduced a new method of adding parameter to expand a family of distributions. Using this concept, in this article, the Marshall–Olkin extended Pareto distribution is introduced and some recurrence relations for single and product moments of generalized order statistics are studied. Also the results are deduced for record values and order statistics.     

Complete moment convergence of moving-average process generated by a class of random variables

In this article, we establish the complete moment convergence of a moving-average process generated by a class of random variables satisfying the Rosenthal-type maximal inequality and the week mean dominating condition. On the one hand, we give the correct proof for the case p = 1 in Ko (2015 Ko, M.H. (2015). Complete moment convergence of moving average process generated by a class of random variables. J. Inequalities Appl. 2015(1):1–9. Article ID 225.[Crossref], [Web of Science ®] , [Google Scholar]); on the other hand, we also consider the case αp = 1 which was not considered in Ko (2015 Ko, M.H. (2015). Complete moment convergence of moving average process generated by a class of random variables. J. Inequalities Appl. 2015(1):1–9. Article ID 225.[Crossref], [Web of Science ®] , [Google Scholar]). The results obtained in this article generalize some corresponding ones for some dependent sequences.     

On the optimal allocation of active redundancies in series system

This paper studies the allocations of two non identical active redundancies in series systems in terms of the reversed hazard rate order and hazard rate order, which generalizes some results built in Valdés and Zequeira (2003 Valdés, J. E., and R. I. Zequeira 2003. On the optimal allocation of an active redundancy in a two-component series system. Stat. Probab. Lett. 63:325–32.[Crossref], [Web of Science ®] , [Google Scholar], 2006 Valdés, J. E., and R. I. Zequeira 2006. On the optimal allocation of two active redundancies in a two-component series system. Oper. Res. Lett. 34:49–52.[Crossref], [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 fiducial generators

Fiducial inference has been gaining presence recently and it is the intention of the present article to look at the notion of fiducial generators; meaning procedures to simulate parameter values that in some sense correspond to simulations from some implicit fiducial distribution. It is well known that when the distribution has group structure, stemming from the natural pivotal associated, a fiducial may be obtained. It is in the non group distributions that there appears to be still room for finding a fiducial distribution. Recently some general procedures have been proposed for dealing with generalized fiducials, but these depend on certain choices for a structural equation or a fiducial equation, as in Hannig (2009 Hannig, J. (2009). On generalized fiducial inference. Stat. Sin. 19:491–544.[Web of Science ®] , [Google Scholar]) or Taraldsen and Lindqvist (2013 Taraldsen, G., Lindqvist, B.H. (2013). Fiducial theory and optimal inference. Ann. Stat. 41(1):323–341.[Crossref], [Web of Science ®] , [Google Scholar]), respectively. A brief presentation is made of an earlier approach to fiducial inference for multivariate parameters, as in Brillinger (1962 Brillinger, D.R. (1962). Examples bearing on the definition of fiducial probability with a bibliography. Ann. Math. Stat. 33(4):1349–1355.[Crossref] , [Google Scholar]), and the implied fiducial generator introduced in Engen and Lillegård (1997 Engen, S., Lillegård, M. (1997). Stochastic simulation conditioned on sufficient statistics. Biometrika 84(1):235–240.[Crossref], [Web of Science ®] , [Google Scholar]), trying to connect them. Three interesting non group distributions are seen; two of them, the truncated exponential and the two-parameter gamma, already reported in literature. A third non group distribution is analyzed; the inverse Gaussian, connecting the fiducial that results following Brillinger (1962 Brillinger, D.R. (1962). Examples bearing on the definition of fiducial probability with a bibliography. Ann. Math. Stat. 33(4):1349–1355.[Crossref] , [Google Scholar]), with a result pertaining confidence limits for the shape parameter in Hsieh (1990 Hsieh, H.K. (1990). Inferences on the coefficient of variation of an inverse-Gaussian distribution. Commun. Stat. - Theory Methods 19(5):1589–1605.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]). In the three cases, comparisons are made with the Bayesian posteriors that have been known to be close numerically. Some discussion is made on the issue of singularities of the fiducial density and its connection with densities that do not integrate to unity. As to the case of discrete observables, some comments are made for the Bernoulli distribution, only.     

On Complete Convergence for Weighted Sums of NA Random Fields

In this paper, we prove the complete convergence for the weighted sums of negatively associated random variables with multidimensional indices. The main result generalizes Theorem 2.1 in Kuczmaszewska and Lagodowski (2011 Kuczmaszewska, A., Lagodowski, Z.A. (2011). Convergence rates in the SLLN for some classes of dependent random field. J. Math. Anal. Appl. 380:571–584.[Crossref], [Web of Science ®] , [Google Scholar]) to the case of weighted sums.     

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.     

A one-sample location test based on weighted averaging of two test statistics when the dimension and the sample size are large

We discuss a one-sample location test that can be used when the dimension and the sample size are large. It is well-known that the power of Hotelling’s test decreases when the dimension is close to the sample size. To address this loss of power, some non exact approaches were proposed, e.g., Dempster (1958 Dempster, A.P. (1958). A high dimensional two sample significance test. Ann. Math. Stat. 29:995–1010.[Crossref] , [Google Scholar], 1960 Dempster, A.P. (1960). A significance test for the separation of two highly multivariate small samples. Biometrics 16:41–50.[Crossref], [Web of Science ®] , [Google Scholar]), Bai and Saranadasa (1996 Bai, Z.D., Saranadasa, H. (1996). Effect of high dimension: by an example of a two sample problem. Stat. Sin. 6:311–329.[Web of Science ®] , [Google Scholar]), and Srivastava and Du (2008 Srivastava, M.S., Du, M. (2008). A test for the mean vector with fewer observations than the dimension. J. Multivariate Anal. 99:386–402.[Crossref], [Web of Science ®] , [Google Scholar]). In this article, we focus on Hotelling’s test and Dempster’s test. The comparative merits and demerits of these two tests vary according to the local parameters. In particular, we consider the situation where it is difficult to determine which test should be used, that is, where the two tests are asymptotically equivalent in terms of local power. We propose a new statistic based on the weighted averaging of Hotelling’s T²-statistic and Dempster’s statistic that can be applied in such a situation. Our weight is determined on the basis of the maximum local asymptotic power on a restricted parameter space that induces local asymptotic equivalence between Hotelling’s test and Dempster’s test. Numerical results show that our test is more stable than Hotelling’s T²-statistic and Dempster’s statistic in most parameter settings.     

Sample size estimation for comparing rates of change in K-group repeated count outcomes

Sample size estimation for comparing the rates of change in two-arm repeated measurements has been investigated by many investigators. In contrast, the literature has paid relatively less attention to sample size estimation for studies with multi-arm repeated measurements where the design and data analysis can be more complex than two-arm trials. For continuous outcomes, Jung and Ahn (2004 Jung, S., Ahn, C. (2004). K-sample test and sample size calculation for comparing slopes in data with repeated measurements. Biometrical J. 46(5):554–564.[Crossref], [Web of Science ®] , [Google Scholar]) and Zhang and Ahn (2013 Zhang, S., Ahn, C. (2013). Sample size calculation for comparing time-averaged responses in k-group repeated measurement studies. Comput. Stat. Data Anal. 58:283–291.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) have presented sample size formulas to compare the rates of change and time-averaged responses in multi-arm trials, using the generalized estimating equation (GEE) approach. To our knowledge, there has been no corresponding development for multi-arm trials with count outcomes. We present a sample size formula for comparing the rates of change in multi-arm repeated count outcomes using the GEE approach that accommodates various correlation structures, missing data patterns, and unbalanced designs. We conduct simulation studies to assess the performance of the proposed sample size formula under a wide range of designing configurations. Simulation results suggest that empirical type I error and power are maintained close to their nominal levels. The proposed method is illustrated using an epileptic clinical trial example.     

Stress–strength models with more than two states under exponential distribution

In the stress–strength models, analysis is based on the reliability of the system where the system is either in operational state or in failure state. Ery?lmaz (2011 Ery?lmaz, S. (2011). A new perspective to stress–strength models. Ann. Inst. Stat. Math. 63(1):101–115.[Crossref], [Web of Science ®] , [Google Scholar]) introduced the stress–strength reliability in a different framework assigning more than two states to the system depending on the difference between strength and stress values. Unlike Ery?lmaz (2011 Ery?lmaz, S. (2011). A new perspective to stress–strength models. Ann. Inst. Stat. Math. 63(1):101–115.[Crossref], [Web of Science ®] , [Google Scholar]), the present article deals with the ratio of the strength and stress values when the stress and strength follow independent exponential distributions. This article presents in detail the estimation aspect of the multistate stress–strength reliability function.