A Generalization of the Skew-Normal Distribution: The Beta Skew-Normal

We consider a new generalization of the skew-normal distribution introduced by Azzalini (1985 Azzalini , A. ( 1985 ). A class of distributions which includes the normal ones . Scand. J. Statis. 12 ( 2 ): 171 – 178 .[Web of Science ®] , [Google Scholar]). We denote this distribution Beta skew-normal (BSN) since it is a special case of the Beta generated distribution (Jones, 2004 Jones , M. C. ( 2004 ). Families of distributions of order statistics . Test 13 ( 1 ): 1 – 43 .[Crossref], [Web of Science ®] , [Google Scholar]). Some properties of the BSN are studied. We pay attention to some generalizations of the skew-normal distribution (Bahrami et al., 2009 Bahrami , W. , Agahi , H. , Rangin , H. ( 2009 ). A two-parameter Balakrishnan skew-normal distribution . J. Statist. Res. Iran 6 : 231 – 242 . [Google Scholar]; Sharafi and Behboodian, 2008 Sharafi , M. , Behboodian , J. ( 2008 ). The Balakrishnan skew-normal density . Statist. Pap. 49 : 769 – 778 .[Crossref], [Web of Science ®] , [Google Scholar]; Yadegari et al., 2008 Yadegari , I. , Gerami , A. , Khaledi , M. J. ( 2008 ). A generalization of the Balakrishnan skew-normal distribution . Statist. Probab. Lett. 78 : 1165 – 1167 .[Crossref], [Web of Science ®] , [Google Scholar]) and to their relations with the BSN.     

Is Each NPMLE of a Continuous Bivariate Distribution Function with Singly Right-Censored Data Really Inconsistent?

We consider non-parametric estimation of a continuous cdf of a random vector (X ₁, X ₂). With bivariate RC data, it is stated in van der Laan (1996 Van der Laan , M. J. ( 1996 ) Efficient estimation in the bivariate censoring model and repairing NPMLE . Ann. Statist. 24 : 596 – 627 .[Crossref], [Web of Science ®] , [Google Scholar], p. 598¹⁰, Ann. Statist.), Quale et al. (2006 Quale , C. M. , van der Laan , M. J. , Robins , J. R. ( 2006 ). Locally efficient estimation with bivariate right-censored data . JASA. 101 : 1076 – 1084 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar], JASA) etc. that “it is well known that the NPMLE for continuous data is inconsistent (Tsai et al. (1986 Tsai , W. Y. , Leurgans , S. , Crowley , J. ( 1986 ). Nonparametric estimation of a bivariate survival function in the presence of censoring . Ann. Statist. 14 : 1351 – 1365 .[Crossref], [Web of Science ®] , [Google Scholar])).” The claim is based on a result in Tsai et al. (1986 Tsai , W. Y. , Leurgans , S. , Crowley , J. ( 1986 ). Nonparametric estimation of a bivariate survival function in the presence of censoring . Ann. Statist. 14 : 1351 – 1365 .[Crossref], [Web of Science ®] , [Google Scholar], p.1352, Ann. Statist.) that if X ₁ is right censored but not X ₂, then common ways for defining one NPMLE lead to inconsistency. If X ₁ is right censored and X ₂ is type I right-censored (which includes the case in Tsai et al.), we present a consistent NPMLE. The result corrects a common misinterpretation of Tsai's example (Tsai et al., 1986 Tsai , W. Y. , Leurgans , S. , Crowley , J. ( 1986 ). Nonparametric estimation of a bivariate survival function in the presence of censoring . Ann. Statist. 14 : 1351 – 1365 .[Crossref], [Web of Science ®] , [Google Scholar], Ann. Statist.).     

A Comparison of Procedures for Controlling the False Discovery Rate in the Presence of Small Variance Genes: A Simulation Study

The Significance Analysis of Microarrays (SAM; Tusher et al., 2001 Tusher , V. G. , Tibshirani , R. , Chu , G. ( 2001 ). Significance analysis of microarrys applied to the ionizing radiation response . Proceedings of the National Academy of Sciences 98 : 5116 – 5121 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) method is widely used in analyzing gene expression data while controlling the FDR by using resampling-based procedure in the microarray setting. One of the main components of the SAM procedure is the adjustment of the test statistic. The introduction of the fudge factor to the test statistic aims at deflating the large value of test statistics due to the small standard error of gene-expression. Lin et al. (2008 Lin , D. , Shkedy , Z. , Burzykowski , T. , Göhlmann , H. W. H. , De Bondt , A. , Perera , T. , Geerts , T. , Bijnens , L. ( 2008 ). Significance analysis of microarray (SAM) for comparisons of several treatments with one control . Biometric Journal, MCP 50 ( 5 ): 801 – 823 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) pointed out that the fudge factor does not effectively improve the power and the control of the FDR as compared to the SAM procedure without the fudge factor in the presence of small variance genes. Motivated by the simulation results presented in Lin et al. (2008 Lin , D. , Shkedy , Z. , Burzykowski , T. , Göhlmann , H. W. H. , De Bondt , A. , Perera , T. , Geerts , T. , Bijnens , L. ( 2008 ). Significance analysis of microarray (SAM) for comparisons of several treatments with one control . Biometric Journal, MCP 50 ( 5 ): 801 – 823 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]), in this article, we extend our study to compare several methods for choosing the fudge factor in the modified t-type test statistics and use simulation studies to investigate the power and the control of the FDR of the considered methods.     

Aspects of Dependence in Generalized Farlie-Gumbel-Morgenstern Distributions

In this article, we obtained a dependence measure for generalized Farlie-Gumbel-Morgenstern (FGM) family in view of Kochar and Gupta (1987 Kochar , S. G. , Gupta , R. P. ( 1987 ). Competitors of Kendall-tau test for testing independence against PQD . Biometrika 74 ( 3 ): 664 – 669 .[Crossref], [Web of Science ®] , [Google Scholar]) and then compared this measure with Spearman's rho and Kendall's tau in FGM family. Moreover, we evaluated the empirical power of 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 PQD . Biometrika 74 ( 3 ): 664 – 669 .[Crossref], [Web of Science ®] , [Google Scholar], 1990 Kochar , S. G. , Gupta , R. P. ( 1990 ). Distribution-free tests based on sub-sample extrema for testing against positive dependence . Australian Journal of Statistics 32 : 45 – 51 .[Crossref] , [Google Scholar]) based on exact distribution of a U-statistics. This is derived via a simulation study for sample of sizes n = 6, 8, 10, 12, 16, and 20. Also, we compared our simulation results with those achieved by Amini et al. (2010 Amini , M. , Jabbari , H. , Mohtashami Borzadaran , G. R. , Azadbakhsh , M. ( 2010 ). Power comparison of independence test for the Farlie-Gumbel-Moregenstern family . Communications of the Korean Statistical Society 17 ( 4 ): 493 – 505 .[Crossref] , [Google Scholar]) and Güven and Kotz (2008 Güven , B. , Kotz , S. ( 2008 ). Test of independence for generalized Farlie-Gumbel-Morgenstern distributions . Journal of Computational and Applied Mathematics 212 : 102 – 111 .[Crossref], [Web of Science ®] , [Google Scholar]).     

Improved Bar-Lev,Bobovitch, and Boukai Randomized Response Models

In this article, two new improved randomized response models have been proposed. The proposed models are found to be more efficient than the recent randomized response model studied by Bar-Lev et al. (2004 Bar-Lev , S. K. , Bobovitch , E. , Boukai , B. ( 2004 ). A note on randomized response models for quantitative data . Metrika 60 : 255 – 260 .[Crossref], [Web of Science ®] , [Google Scholar]). The relative efficiency of the proposed models has been studied with respect to the Bar-Lev et al. (2004 Bar-Lev , S. K. , Bobovitch , E. , Boukai , B. ( 2004 ). A note on randomized response models for quantitative data . Metrika 60 : 255 – 260 .[Crossref], [Web of Science ®] , [Google Scholar]) model under different situations.     

On a Generalization of Bivariate Cauchy Distribution

This paper addresses a generalization of the bivariate Cauchy distribution discussed by Fang et al. (1990 Fang , K. T. , Kotz , S. , Ng , K. W. ( 1990 ). Symmetric Multivariate and Related Distributions . London : Chapman and Hall .[Crossref] , [Google Scholar]), derived from a trivariate normal distribution with a general correlation matrix. We obtain explicit expressions for the joint distribution function and joint density function, and show that they reduce in a special case to the corresponding expressions of Fang et al. (1990 Fang , K. T. , Kotz , S. , Ng , K. W. ( 1990 ). Symmetric Multivariate and Related Distributions . London : Chapman and Hall .[Crossref] , [Google Scholar]). Finally, we show that this generalized distribution is useful in determining the orthant probability of a bivariate skew-normal distribution of Azzalini and Dalla Valle (1996 Azzalini , A. , Dalla Valle , A. ( 1996 ). The multivariate skew-normal distribution . Biometrika 83 : 715 – 726 .[Crossref], [Web of Science ®] , [Google Scholar]).     

Bayesian Inference in Dynamic Disequilibrium Models: An Application to the Polish Credit Market

We propose a Bayesian approach for inference in a dynamic disequilibrium model. To circumvent the difficulties raised by the Maddala and Nelson (1974 Maddala , G. , Nelson , F. ( 1974 ). Maximum likelihood methods for models of markets in disequilibrium . Econometrica 42 ( 6 ): 1013 – 1030 .[Crossref], [Web of Science ®] , [Google Scholar]) specification in the dynamic case, we analyze a dynamic extended version of the disequilibrium model of Ginsburgh et al. (1980 Ginsburgh , V. , Tishler , A. , Zang , I. ( 1980 ). Alternative estimation methods for two regime models . European Economic Review 13 : 207 – 228 .[Crossref], [Web of Science ®] , [Google Scholar]). We develop a Gibbs sampler based on the simulation of the missing observations. The feasibility of the approach is illustrated by an empirical analysis of the Polish credit market, for which we conduct a specification search using the posterior deviance criterion of Spiegelhalter et al. (2002 Spiegelhalter , D. J. , Best , N. G. , Carlin , B. P. , Van Der Linde , A. ( 2002 ). Bayesian measures of complexity and fit . Journal of the Royal Statistical Society, Ser. B 64 : 583 – 639 .[Crossref] , [Google Scholar]).     

A Note on Importance Resampling for Multi-Dimensional Statistics

Johns (1988 Johns , M. V. (1988). Importance sampling for bootstrap confidence intervals. Journal of the American Statistical Association 83:709–714.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), Davison (1988 Davison , A. C. ( 1988 ). Discussion of paper by D. V. Hinkley . Journal of the Royal Statistical Society Series B 50 : 356 – 357 . [Google Scholar]), and Do and Hall (1991 Do , K. A. , Hall , P. ( 1991 ). On importance sampling for the bootstrap . Biometrika 78 : 161 – 167 .[Crossref], [Web of Science ®] , [Google Scholar]) used importance sampling for calculating bootstrap distributions of one-dimensional statistics. Realizing that their methods can not be extended easily to multi-dimensional statistics, Fuh and Hu (2004 Fuh , C. D. , Hu , I. ( 2004 ). Efficient importance sampling for events of moderate deviations with applications . Biometrika 91 : 471 – 490 .[Crossref], [Web of Science ®] , [Google Scholar]) proposed an exponential tilting formula for statistics of multi-dimension, which is optimal in the sense that the asymptotic variance is minimized for estimating tail probabilities of asymptotically normal statistics. For one-dimensional statistics, Hu and Su (2008 Hu , J. , Su , Z. ( 2008 ). Adaptive resampling algorithms for estimating bootstrap distributions . Journal of Statistical Planning and Inference 138 ( 6 ): 1763 – 1777 .[Crossref], [Web of Science ®] , [Google Scholar]) proposed a multi-step variance minimization approach that can be viewed as a generalization of the two-step variance minimization approach proposed by Do and Hall (1991 Do , K. A. , Hall , P. ( 1991 ). On importance sampling for the bootstrap . Biometrika 78 : 161 – 167 .[Crossref], [Web of Science ®] , [Google Scholar]). In this article, we generalize the approach of Hu and Su (2008 Hu , J. , Su , Z. ( 2008 ). Adaptive resampling algorithms for estimating bootstrap distributions . Journal of Statistical Planning and Inference 138 ( 6 ): 1763 – 1777 .[Crossref], [Web of Science ®] , [Google Scholar]) to multi-dimensional statistics, which applies to general statistics and does not resort to asymptotics. Empirical results on a real survival data set show that the proposed algorithm provides significant computational efficiency gains.     

Computing the Exact Fisher Information Matrix of Periodic State-Space Models

The purpose of this article is to develop algorithms for computing the exact Fisher information matrix of periodic time-varying state-space models. We first present a relatively simple recursive algorithm which computes the elements of the exact information matrix without involving numerical differentiation, since all required derivatives are analytically evaluated. The proposed algorithm extends the procedure due to Cavanaugh and Shumway (1996 Cavanaugh , J. E. , Shumway , R. H. ( 1996 ). On computing the expected Fisher information matrix for state-space model parameters . Statist. Probab. Lett. 26 : 347 – 355 .[Crossref], [Web of Science ®] , [Google Scholar]) to the periodic state-space framework. Exploiting the approach used in Klein et al. (2000 Klein , A. , Mélard , G. , Zahaf , T. ( 2000 ). Construction of the exact Fisher information matrix of Gaussian time series models by means of matrix differential rules . Linear Alg. Applic. 321 : 209 – 232 .[Crossref], [Web of Science ®] , [Google Scholar]), a second algorithm is proposed in order to obtain the exact information matrix as a whole instead of element by element. The algorithms are first developed in a general framework and then specialized to the case of a periodic Gaussian vector autoregressive moving-average (PVARMA) model.     

Additive Mixed Effect Model for Clustered Doubly Censored Data

Cai and Zeng (2011 Cai, J. and Zeng, D. 2011. Additive mixed effect model for clustered failure time data. Biometrics, 67(4): 1340–1351. [Crossref], [PubMed] , [Google Scholar]) proposed an additive mixed effect model to analyze clustered right-censored data. In this article, we demonstrate that the approach of Cai and Zeng (2011 Cai, J. and Zeng, D. 2011. Additive mixed effect model for clustered failure time data. Biometrics, 67(4): 1340–1351. [Crossref], [PubMed] , [Google Scholar]) can be extended to clustered doubly censored data. Furthermore, when both left- and right-censoring variables are always observed, we propose alternative estimators using the approach of Cai and Cheng (2004 Cai, T. and Cheng, S. C. 2004. Semiparametric regression analysis for doubly censored data. Biometrika, 91: 277–290. [Crossref], [Web of Science ®] , [Google Scholar]). A simulation study is conducted to investigate the performance of the proposed estimators.     

The performance of the truncated mixed control chart

Abstract

When the mixed chart proposed by Aslam et al. (2015 Aslam, M., M. Azam, N. Khan, and C.-H. Jun. 2015. A mixed control chart to monitor the process. International Journal of Production Research 53 (15):4684–93. doi:10.1080/00207543.2015.1031354.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) is in use, the sample items are classified as defective or not defective and, depending on the number of defectives, the quality characteristic X of the sample items are also measured. In this case, an Xbar chart decides the state of the process. The previous conforming/non-conforming classification truncates the X distribution and, because of that, the mathematical development to obtain the ARLs is complex. Aslam et al. (2015 Aslam, M., M. Azam, N. Khan, and C.-H. Jun. 2015. A mixed control chart to monitor the process. International Journal of Production Research 53 (15):4684–93. doi:10.1080/00207543.2015.1031354.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) didn’t pay attention to the fact that the X distribution is truncated and, due to that, they obtained incorrect ARLs.     

Inverse Sampling with Unequal Selection Probabilities

In some real situations the population of interest is divided into two groups, of which one contains only a few units. In other cases, the population may be considered as subdivided into two group', for example, if only a few units display a value of the variable of interest which is highly different from zero, while all the other units show a value equal to or near zero. In both cases, inverse sampling is more efficient than classical fixed sample-size designs to obtain the parameter estimators for the whole population as well as for its groups (e.g., Salehi and Seber, 2004 Salehi , M. M. , Seber , G. A. F. ( 2004 ). A general inverse sampling scheme and its application to adaptive cluster sampling . Austral. NZ J. Statist. 46 : 483 – 494 .[Crossref], [Web of Science ®] , [Google Scholar]). In fact, in this design the procedure selection continues until a prefixed number of units with the characteristic of interest is sampled. Since it is not known a priori to which group the population units belong, the sample size is a random variable. Christman and Lan (2001 Christman , M. C. , Lan , F. ( 2001 ). Inverse adaptive cluster sampling . Biometrics 57 : 1096 – 1105 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) and Salehi and Seber (2001 Salehi , M. M. , Seber , G. A. F. ( 2001 ). A new proof of Murthy's estimator which applies to sequential sampling . Austral. NZ J. Statist. 43 : 281 – 286 . [Google Scholar] 2004 Salehi , M. M. , Seber , G. A. F. ( 2004 ). A general inverse sampling scheme and its application to adaptive cluster sampling . Austral. NZ J. Statist. 46 : 483 – 494 .[Crossref], [Web of Science ®] , [Google Scholar]) considered inverse sampling designs when all the population units have equal selection probabilities. In this article, we consider the general case in which the units may have unequal probabilities of being included in the sample. In fact, in many real situations different units may have different selection probabilities because of some inherent features of the sampling procedure, or in order to obtain better estimates. We derive unbiased estimators of the totals of the two groups, their variance and the corresponding unbiased variance estimators in inverse sampling with replacement. Finally, we derive similar results for more complex designs, where the selection procedure stops before observing the prefixed number of units from the rare group.     

Relations for Moments of Progressively Type-II Censored Order Statistics from Log-Logistic Distribution with Applications to Inference

In this article, we establish several recurrence relations for the single and product moments of progressively Type-II right censored order statistics from a log-logistic distribution. The use of these relations in a systematic recursive manner would enable the computation of all the means, variances and covariances of progressively Type-II right censored order statistics from the log-logistic distribution for all sample sizes n, effective sample sizes m, and all progressive censoring schemes (R ₁,…, R _m ). The results established here generalize the corresponding results for the usual order statistics due to Balakrishnan and Malik (1987 Balakrishnan , N. , Malik , H. J. ( 1987 ). Moments of order statistics from truncated log-logistic distribution . J. Statist. Plann. Infer. 17 : 251 – 267 .[Crossref], [Web of Science ®] , [Google Scholar]) and Balakrishnan et al. (1987 Balakrishnan , N. , Malik , H. J. , Puthenpura , S. ( 1987 ). Best linear unbiased estimation of location and scale parameters of the log-logistic distribution . Commun. Statist. Theor. Meth. 16 : 3477 – 3495 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). The moments so determined are then utilized to derive best linear unbiased estimators for the scale- and location-scale log-logistic distributions. A comparison of these estimates with the maximum likelihood estimates is made through Monte Carlo simulation. The best linear unbiased predictors of progressively censored failure times is then discussed briefly. Finally, a numerical example is presented to illustrate all the methods of inference developed here.     

Efficient Ratio and Product Estimators in Stratified Random Sampling

This paper suggests an efficient class of ratio and product estimators for estimating the population mean in stratified random sampling using auxiliary information. It is interesting to mention that, in addition to many, Koyuncu and Kadilar (2009 Koyuncu , N. , Kadilar , C. ( 2009 ). Ratio and product estimators in stratified random sampling . J. Statist. Plann. Infer. 139 : 2552 – 2558 .[Crossref], [Web of Science ®] , [Google Scholar]), Kadilar and Cingi (2003 Kadilar , C. , Cingi , H. ( 2003 ). Ratio estimator in stratified sampling . Biometr. J. 45 : 218 – 225 .[Crossref], [Web of Science ®] , [Google Scholar], 2005 Kadilar , C. , Cingi , H. ( 2005 ). A new estimator in stratified random sampling . Commun. Statist. Theor. Meth. 34 : 597 – 602 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), and Singh and Vishwakarma (2007 Singh , H. P. , Vishwakarma , G. K. ( 2007 ). Modified exponential ratio and product estimators for finite population mean in double sampling . Austr. J. Statist. 36 ( 3 ): 217 – 225 . [Google Scholar]) estimators are identified as members of the proposed class of estimators. The expressions of bias and mean square error (MSE) of the proposed estimators are derived under large sample approximation in general form. Asymptotically optimum estimator (AOE) in the class is identified alongwith its MSE formula. It has been shown that the proposed class of estimators is more efficient than combined regression estimator and Koyuncu and Kadilar (2009 Koyuncu , N. , Kadilar , C. ( 2009 ). Ratio and product estimators in stratified random sampling . J. Statist. Plann. Infer. 139 : 2552 – 2558 .[Crossref], [Web of Science ®] , [Google Scholar]) estimator. Moreover, theoretical findings are supported through a numerical example.     

The Identification of Potential Longitudinal Biomarkers and Measurements of Effectiveness for Biomarkers as Surrogates in Multivariate Survival Data

A proposed method based on frailty models is used to identify longitudinal biomarkers or surrogates for a multivariate survival. This method is an extention of earlier models by Wulfsohn and Tsiatis (1997 Wulfsohn , M. S. , Tsiatis , A. A. ( 1997 ). A joint model for survival and longitudinal data measured with error . Biometrics 53 ( 1 ): 330 – 339 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) and Song et al. (2002 Song , X. , Davidian , M. , Tsiatis , A. A. ( 2002 ). A Semiparametric likelihood approach to joint modeling of longitudinal and time-to-event data . Biometrics 58 ( 4 ): 742 – 753 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]). In this article, similar to Henderson et al. (2002 Henderson , R. , Diggle , P. J. , Dobson , A. ( 2002 ). Identification and efficacy of longitudinal markers for survival . Biostatistics 3 ( 1 ): 33 – 50 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]), a joint likelihood function combines the likelihood functions of the longitudinal biomarkers and the multivariate survival times. We use simulations to explore how the number of individuals, the number of time points per individual and the functional form of the random effects from the longitudianl biomarkers influence the power to detect the association of a longitudinal biomarker and the multivariate survival time. The proposed method is illustrate by using the gastric cancer data.     

Upper Bounds for the Ruin Probabilities of the Entrance-Based Risk Model

Li et al. (2005 Li , Z. , Zhu , J. , Chen , F. ( 2005 ). Study of a risk model based on the entrance process . Statist. Probab. Lett. 72 ( 1 ): 1 – 10 .[Crossref], [Web of Science ®] , [Google Scholar]) proposed a risk model based on the entrance process and studied the asymptotic behavior of the surplus when time goes to infinity. This article considers the ruin problem in that model. Some simple characteristics (stochastic intensity, compensator, mean process, etc.) of the risk process and other related processes are also considered. Under small claim condition, exponential upper bounds for the ruin probability are obtained.     

Semiparametric analysis of transformation models with doubly censored data

Double censoring arises when T represents an outcome variable that can only be accurately measured within a certain range, [L, U], where L and U are the left- and right-censoring variables, respectively. In this note, using Martingale arguments of Chen et al. [3 Chen, K., Jin, Z. and Ying, Z. 2002. Semiparametric analysis of transformation models with censored data. Biometrika, 89: 659–668. [Crossref], [Web of Science ®] [Google Scholar]], we propose an estimator (denoted by ?β) for estimating regression coefficients of transformation model when L is always observed. Under Cox proportional hazards model, the proposed estimator is equivalent to the partial likelihood estimator for left-truncated and right-censored data if the left-censoring variables L were regarded as left-truncated variables. In this case, the estimator ?β can be obtained by the standard software. A simulation study is conducted to investigate the performance of ?β. For the purpose of comparison, the simulation study also includes the estimator proposed by Cai and Cheng [2 Cai, T. and Cheng, S. 2004. Semiparametric regression analysis for doubly censored data. Biometrika, 91: 277–290. [Crossref], [Web of Science ®] [Google Scholar]] for the case when L and U are always observed.     

On Universal Optimality of Circular Weakly Neighbor Balanced Designs Under an Interference Model

In this article, the optimal design problem in a fixed effects interference model with left-neighbor effects is studied. It is known (Druilhet, 1999 Druilhet , P. ( 1999 ). Optimality of circular neighbor balanced designs . J. Statist. Plann. Infer. 81 : 141 – 152 .[Crossref], [Web of Science ®] , [Google Scholar]) that circular neighbor balanced designs (CNBDs) are universally optimal in such a model. We prove the universal optimality of circular weakly neighbor balanced designs (CWNBDs), which require a smaller number of blocks than CNBDs. CWNBDs with the number of blocks smaller than the number of treatments belong to the class of partially neighbor balanced designs (PNBDs) defined by Wilkinson et al. (1983 Wilkinson , G. N. , Eckert , R. , Hancock , T. W. , Mayo , O. ( 1983 ). Nearest neighbour (NN) analysis of field experiments . J. Roy. Statist. Soc. Ser. B 45 : 151 – 210 . [Google Scholar]). We give a construction method for some CWNBDs, with examples.     

Step-Up Procedure for Multiple Comparison with a Control for Multivariate Normal Means

In this study, we consider the multiple comparison with a control for multivariate normal means. Specifically, we construct a step-up procedure by referring to Dunnett and Tamhane (1992 Dunnett , C. W. , Tamhane , A. C. ( 1992 ). A step-up multiple test procedure . Journal of the American Statistical Association 87 : 162 – 170 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). We derive recursive formulae for determining the critical values of the step-up procedure for a specified significance level. Then we formulate the power of the test. Finally, we compare the step-up procedure with the single-step procedure proposed by Nakamura and Imada (2005 Nakamura , T. , Imada , T. ( 2005 ). Multiple comparison procedure of Dunnett's type for multivariate normal means . Journal of the Japanese Society of Computational Statistics 18 : 21 – 32 .[Crossref] , [Google Scholar]) and the step-down procedure proposed by Imada and Douke (2007 Imada , T. , Douke , H. ( 2007 ). Step down procedure for comparing several treatments with a control based on multivariate normal response . Biometrical Journal 49 ( 1 ): 18 – 29 .[Crossref], [Web of Science ®] , [Google Scholar]) in terms of numerical examples regarding the power of the test.     

Choosing Ridge Parameter for Regression Problems

ABSTRACT

Hoerl and Kennard (1970a Hoerl , A. E. , Kennard , R. W. ( 1970a ). Ridge regression: biased estimation for non-orthogonal problems . Tech. 12 : 55 – 67 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) introduced the ridge regression estimator as an alternative to the ordinary least squares estimator in the presence of multicollinearity. In this article, a new approach for choosing the ridge parameter (K), when multicollinearity among the columns of the design matrix exists, is suggested and evaluated by simulation techniques, in terms of mean squared errors (MSE). A number of factors that may affect the properties of these methods have been varied. The MSE from this approach has shown to be smaller than using Hoerl and Kennard (1970a Hoerl , A. E. , Kennard , R. W. ( 1970a ). Ridge regression: biased estimation for non-orthogonal problems . Tech. 12 : 55 – 67 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) in almost all situations.