Robustness of the shrinkage estimator for the relative potency in the combination of multivariate bioassays

ABSTRACT

This article investigates the robustness of the shrinkage Bayesian estimator for the relative potency parameter in the combinations of multivariate bioassays proposed in Chen et al. (1999 Chen, D.G., Carter, E.M., Hubert, J.J., Kim, P.T. (1999). Empirical Bayesian estimation for combinations of multivariate bioassays. Biometrics 55(4):1035–1043. [Google Scholar]), which incorporated prior information on the model parameters based on Jeffreys’ rules. This investigation is carried out for the families of t-distribution and Cauchy-distribution based on the characteristics of bioassay theory since the t-distribution approaches the normal distribution which is the most commonly used distribution in the applications of bioassay as the degrees of freedom increases and the t-distribution approaches the Cauchy-distribution as the degrees of freedom approaches 1 which is also an important distribution in bioassay. A real data is used to illustrate the application of this investigation. This analysis further supports the application of the shrinkage Bayesian estimator to the theory of bioassay along with the empirical Bayesian estimator.     

On Complete Convergence for Arrays of Dependent Random Variables

A complete convergence theorem for an array of rowwise independent random variables was established by Sung et al. (2005 Sung , S. H. , Volodin , A. I. , Hu , T.-C. ( 2005 ). More on complete convergence for arrays . Statist. Probab. Lett. 71 : 303 – 311 .[Crossref], [Web of Science ®] , [Google Scholar]). This result has been generalized and extended by Kruglov et al. (2006 Kruglov , V. M. , Volodin , A. I. , Hu , T.-C. ( 2006 ). On complete convergence for arrays . Statist. Probab. Lett. 76 : 1631 – 1640 .[Crossref], [Web of Science ®] , [Google Scholar]) and Chen et al. (2007 Chen , P. , Hu , T.-C. , Liu , X. , Volodin , A. ( 2007 ). On complete convergence for arrays of rowwise negatively associated random variables . Theor. Probab. Appl. 52 : 393 – 397 . [Google Scholar]). In this article, we extend the results of Sung et al. (2005 Sung , S. H. , Volodin , A. I. , Hu , T.-C. ( 2005 ). More on complete convergence for arrays . Statist. Probab. Lett. 71 : 303 – 311 .[Crossref], [Web of Science ®] , [Google Scholar]), Kruglov et al. (2006 Kruglov , V. M. , Volodin , A. I. , Hu , T.-C. ( 2006 ). On complete convergence for arrays . Statist. Probab. Lett. 76 : 1631 – 1640 .[Crossref], [Web of Science ®] , [Google Scholar]), and Chen et al. (2007 Chen , P. , Hu , T.-C. , Liu , X. , Volodin , A. ( 2007 ). On complete convergence for arrays of rowwise negatively associated random variables . Theor. Probab. Appl. 52 : 393 – 397 . [Google Scholar]) to an array of dependent random variables satisfying Hoffmann-Jørgensen type inequalities.     

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.     

Some Fundamental Properties of a Multivariate von Mises Distribution

In application areas like bioinformatics, multivariate distributions on angles are encountered which show significant clustering. One approach to statistical modeling of such situations is to use mixtures of unimodal distributions. In the literature (Mardia et al., 2012 Mardia , K. V. , Kent , J. T. , Zhang , Z. , Taylor , C. , Hamelryck , T. ( 2012 ). Mixtures of concentrated multivariate sine distributions with applications to bioinformatics . J. Appl. Stat. 39 : 2475 – 2492 .[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]), the multivariate von Mises distribution, also known as the multivariate sine distribution, has been suggested for components of such models, but work in the area has been hampered by the fact that no good criteria for the von Mises distribution to be unimodal were available. In this article we study the question about when a multivariate von Mises distribution is unimodal. We give sufficient criteria for this to be the case and show examples of distributions with multiple modes when these criteria are violated. In addition, we propose a method to generate samples from the von Mises distribution in the case of high concentration.     

Marginal correlation from an extended random-effects model for repeated and overdispersed counts

Vangeneugden et al. [15 Vangeneugden, T., Molenberghs, G., Laenen, A., Geys, H., Beunckens, C. and Sotto, C. 2007. Marginal correlation in longitudinal binary data based on generalized linear mixed models, Tech. Rep., Hasselt University. submitted for publication [Google Scholar]] derived approximate correlation functions for longitudinal sequences of general data type, Gaussian and non-Gaussian, based on generalized linear mixed-effects models (GLMM). Their focus was on binary sequences, as well as on a combination of binary and Gaussian sequences. Here, we focus on the specific case of repeated count data, important in two respects. First, we employ the model proposed by Molenberghs et al. [13 Molenberghs, G., Verbeke, G. and Demétrio, C. G.B. 2007. An extended random-effects approach to modeling repeated, overdispersed count data. Lifetime Data Anal., 13: 513–531. [Crossref], [PubMed], [Web of Science ®] [Google Scholar]], which generalizes at the same time the Poisson-normal GLMM and the conventional overdispersion models, in particular the negative-binomial model. The model flexibly accommodates data hierarchies, intra-sequence correlation, and overdispersion. Second, means, variances, and joint probabilities can be expressed in closed form, allowing for exact intra-sequence correlation expressions. Next to the general situation, some important special cases such as exchangeable clustered outcomes are considered, producing insightful expressions. The closed-form expressions are contrasted with the generic approximate expressions of Vangeneugden et al. [15 Vangeneugden, T., Molenberghs, G., Laenen, A., Geys, H., Beunckens, C. and Sotto, C. 2007. Marginal correlation in longitudinal binary data based on generalized linear mixed models, Tech. Rep., Hasselt University. submitted for publication [Google Scholar]]. Data from an epileptic-seizures trial are analyzed and correlation functions derived. It is shown that the proposed extension strongly outperforms the classical GLMM.     

Estimating Incidence Rate on Current Status Data with Application to a Phase IV Cancer Trial

New drug discovery in the pediatrics has dramatically improved survival, but with long- term adverse events. This motivates the examination of adverse outcomes such as long-term toxicity in a phase IV trial. An ideal approach to monitor long-term toxicity is to systematically follow the survivors, which is generally not feasible. Instead, cross-sectional surveys are conducted in Hudson et al. (2007 Hudson , M. M. , Rai , S. N. , Nunez , C. , Merchant , T. E. , Marina , N. M. , Zalamea , N. , Cox , C. , Phipps , S. , Pompeu , R. , Rosenthal , D. ( 2007 ). Noninvasive evaluation of late anthracycline cardiac toxicity in childhood cancer survivors . J. Clin. Oncol. 25 : 3635 – 3643 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]), with one of the objectives to estimate the cumulative incidence rates along with specific interest in fixed-term (5 or 10 year) rates. We present inference procedures based on current status data to our motivating example with very interesting findings.     

Multivariate Generalized Distributions of Order k

Abstract

The study of multivariate distributions of order k, two of which are the multivariate negative binomial of order k and the multinomial of the same order, was introduced in Philippou et al. (Philippou, A. N., Antzoulakos, D. L., Tripsiannis, G. A. (1988 Philippou, A. N., Antzoulakos, D. L. and Tripsiannis, G. A. 1988. Multivariate distributions of order k. Statistics and Probability Letters, 7(3): 207–216.  [Google Scholar]). Multivariate distributions of order k. Statistics and Probability Letters 7(3):207–216.), and Philippou et al. (Philippou, A. N., Antzoulakos, D. L., Tripsiannis, G. A. (1990 Philippou, A. N., Antzoulakos, D. L. and Tripsiannis, G. A. 1990. Multivariate distributions of order k, part II. Statistics and Probability Letters, 10(1): 29–35.  [Google Scholar]). Multivariate distributions of order k, part II. Statistics and Probability Letters 10(1):29–35.). Recently, an order k (or cluster) generalized negative binomial distribution and a multivariate negative binomial distribution were derived in Sen and Jain (Sen, K., Jain, R. (1996 Sen, K. and Jain, R. 1996. “Cluster generalized negative binomial distribution”. In Probability Models and Statistics Medhi Festschrift, A. J., on the Occasion of his 70th Birthday Edited by: Borthakur, A. C. 227–241. New Delhi: New Age International Publishers.  [Google Scholar]). Cluster generalized negative binomial distribution. In: Borthakur et al. A. C., Eds.; Probability Models and Statistics Medhi Festschrift, A. J., on the Occasion of his 70th Birthday. New Age International Publishers: New Delhi, 227–241.) and Sen and Jain (Sen, K., Jain, R. (1997 Sen, K. and Jain, R. 1997. A multivariate generalized Polya-Eggenberger probability model-first passage approach. Communications in Statistics—Theory and Methods, 26: 871–884. [Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). A multivariate generalized Polya-Eggenberger probability model-first passage approach. Communications in Statistics-Theory and Methods 26:871–884.), respectively. In this paper, all four distributions are generalized to a multivariate generalized negative binomial distribution of order k by means of an appropriate sampling scheme and a first passage event. This new distribution includes as special cases several known and new multivariate distributions of order k, and gives rise in the limit to multivariate generalized logarithmic, Poisson and Borel-Tanner distributions of the same order. Applications are indicated.     

The Effect of Migration on a Similarity Index

Fewster and Buckland (2001 Fewster , R. M. , Buckland , S. T. ( 2001 ). Similarity indices for spatial ecological data . Biometrics 57 : 495 – 501 . [Google Scholar]) defined a similarity index between two communities by allowing changes between sites to reduce the influence of local discrepancies. The similarity index of Fewster and Buckland is calculated to attain the maximum similarity between two communities in the presence of migration. Instead of maximizing similarity, we propose random migration to measure the similarity of two communities with two types of stochastic migration. The similarity values based on the proposed methods can be treated as the expected value of similarity under migration. We use computer simulation and empirical examples to demonstrate our approach.     

Optimality of ratio type estimation methods for population mean in the presence of missing data

This article proposes various Searls-type ratio imputation methods (STRIM) on the lines of Ahmed et al. (2006 Ahmed, M. S., O. Al-Titi, Z. Al-Rawi, and W. Abu-Dayyeh. 2006. Estimation of a population mean using different imputation methods. Stat. Trans. 7 (6):1247–1264. [Google Scholar]). It is a well-known fact that the optimal ratio type estimator attains the MSE of regression estimator (or optimal difference estimator) but while using Searls-type transformation (STT) (Searls (1964 Searls, D. T. 1964. The utilization of a known coefficient of variation in the estimation procedure. J. Am. Stat. Assoc. 59:1225–1226.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar])) this may not always happen. These STRIM are shown to perform better than the imputation procedures of Ahmed et al. (2006 Ahmed, M. S., O. Al-Titi, Z. Al-Rawi, and W. Abu-Dayyeh. 2006. Estimation of a population mean using different imputation methods. Stat. Trans. 7 (6):1247–1264. [Google Scholar]). The STRIM may even outperform the Searls type difference imputation methods (STDIM) proposed by us in our earlier work, Bhushan and Pandey (2016 Bhushan, S., and A. P. Pandey. 2016. Optimal imputation of the missing data for estimation of population mean. Journal of Statistics and Management System 19 (6):755–69.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]). This study is concluded with the numerical study along with the theoretical comparison.     

Adaptive tangent distance classifier on recognition of handwritten digits

Simard et al. [16 Simard, P. Y., LeCun, Y., Denker, J. S. and Victorri, B. 2000. Transformation invariance in pattern recognition: Tangent distance and tangent propagation. J. Imaging Syst. Technol., 11: 181–197.  [Google Scholar] 17 Sona, D., Sperduti, A. and Starita, A. 1997. A constructive learning algorithm for discriminant tangent models. Advances in Neural Information Processing Systems. 1997, Cambridge, MA. Edited by: Mozer, M. C., Jordan, M. I. and Petsche, T. Vol. 9, pp.786–792. MIT Press.  [Google Scholar]] proposed a transformation distance called “tangent distance” (TD) which can make pattern recognition be efficient. The key idea is to construct a distance measure which is invariant with respect to some chosen transformations. In this research, we provide a method using adaptive TD based on an idea inspired by “discriminant adaptive nearest neighbor” [7 Hastie, T., Tibshirani, R. and Friedman, J. 2009. The Elements of Statistical Learning, Data Mining, Inference, and Prediction, 2, New York, Berlin, Heidelberg: Springer. Available at http://www-stat.stanford.edu/ElemStatLearn [Google Scholar]]. This method is relatively easy compared with many other complicated ones. A real handwritten recognition data set is used to illustrate our new method. Our results demonstrate that the proposed method gives lower classification error rates than those by standard implementation of neural networks and support vector machines and is as good as several other complicated approaches.     

Measures of predictor sensitivity for order-insensitive partitioning of multiple correlation

Lindeman et al. [12 Lindeman, R. H., Merenda, P. F. and Gold, R. Z. 1980. Introduction to Bivariate and Multivariate Analysis, Glenview, IL: Scott Foresman.  [Google Scholar]] provide a unique solution to the relative importance of correlated predictors in multiple regression by averaging squared semi-partial correlations obtained for each predictor across all p! orderings. In this paper, we propose a series of predictor sensitivity statistics that complement the variance decomposition procedure advanced by Lindeman et al. [12 Lindeman, R. H., Merenda, P. F. and Gold, R. Z. 1980. Introduction to Bivariate and Multivariate Analysis, Glenview, IL: Scott Foresman.  [Google Scholar]]. First, we detail the logic of averaging over orderings as a technique of variance partitioning. Second, we assess predictors by conditional dominance analysis, a qualitative procedure designed to overcome defects in the Lindeman et al. [12 Lindeman, R. H., Merenda, P. F. and Gold, R. Z. 1980. Introduction to Bivariate and Multivariate Analysis, Glenview, IL: Scott Foresman.  [Google Scholar]] variance decomposition solution. Third, we introduce a suite of indices to assess the sensitivity of a predictor to model specification, advancing a series of sensitivity-adjusted contribution statistics that allow for more definite quantification of predictor relevance. Fourth, we describe the analytic efficiency of our proposed technique against the Budescu conditional dominance solution to the uneven contribution of predictors across all p! orderings.     

Steady-state Gibbs sampler estimation for lung cancer data

This paper is based on the application of a Bayesian model to a clinical trial study to determine a more effective treatment to lower mortality rates and consequently to increase survival times among patients with lung cancer. In this study, Qian et al. [13 J. Qian, D.K. Stangl, and S. George, A Weibull model for survival data: Using prediction to decide when to stop a clinical trial, in Bayesian Biostatistics, D. Berry and D. Stangl, eds., Marcel Dekker, New York, 1996, pp. 187–205. [Google Scholar]] strived to determine if a Weibull survival model can be used to decide whether to stop a clinical trial. The traditional Gibbs sampler was used to estimate the model parameters. This paper proposes to use the independent steady-state Gibbs sampling (ISSGS) approach, introduced by Dunbar et al. [3 M. Dunbar, H.M. Samawi, R. Vogel, and L. Yu, A more efficient Gibbs sampler estimation using steady state simulation: Application to public health studies, J. Stat. Simul. Comput. 10.1080/00949655.2013.770857.[Taylor &; Francis Online] , [Google Scholar]], to improve the original Gibbs sampler in multidimensional problems. It is demonstrated that ISSGS provides accuracy with unbiased estimation and improves the performance and convergence of the Gibbs sampler in this application.     

Optimum Multivariate Stratified Sampling Designs with Travel Cost: A Multiobjective Integer Nonlinear Programming Approach

The problem of optimum allocation in stratified sampling and its solution is well known in sampling literature for univariate populations (see Cochran, 1977 Cochran , W. G. ( 1977 ). Sampling Techniques. , 3rd ed. New York : Wiley . [Google Scholar]; Sukhatme et al., 1984 Sukhatme , P. V. , Sukhatme , B. V. , Sukhatme , S. , Ashok , C. ( 1984 ). Sampling Theory of Surveys With Applications. , 3rd ed. Ames , and New Delhi : Iowa State University Press and Indian Society of Agricultural Statistics . [Google Scholar]). In multivariate populations where more than one characteristics are to be studied on every selected unit of the population the problem of finding an optimum allocation becomes more complex due to conflicting behaviour of characteristics. Various authors such as Dalenius (1953 Dalenius , T. ( 1953 ). The multivariate sampling problem . Skandinavisk Actuarietidskrift 36 : 92 – 102 . [Google Scholar], 1957 Dalenius , T. ( 1957 ). Sampling in Sweden. Contributions to the Methods and Theories of Sample Survey Practice . Stockholm : Almqvist and Wicksell . [Google Scholar]), Ghosh (1958 Ghosh , S. P. ( 1958 ). A note on stratified random sampling with multiple characters . Calcutta Statistical Association Bulletin 8 : 81 – 89 . [Google Scholar]), Yates (1960 Yates , F. ( 1960 ). Sampling Methods for Censuses and Surveys. , 3rd ed. London : Charles Griffin . [Google Scholar]), Aoyama (1963 Aoyama , H. ( 1963 ). Stratified random sampling with optimum allocation for multivariate populations . Annals of the Institute of Statistical Mathematics 14 : 251 – 258 .[Crossref], [Web of Science ®] , [Google Scholar]), Gren (1964 Gren , J. ( 1964 ). Some methods of sample allocation in multivariate stratified sampling . Przeglad Statystyczny 11 : 361 – 369 (in Polish) . [Google Scholar], 1966 Gren , J. ( 1966 ). Some application of non-linear programming in sampling methods . Przeglad Statystyczny 13 : 203 – 217 (in Polish) . [Google Scholar]), Folks and Antle (1965 Folks , J. L. , Antle , C. E. ( 1965 ). Optimum allocation of sampling units to the strata when there are r responses of interest . Journal of American Statistical Association 60 : 225 – 233 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), Hartley (1965 Hartley , H. O. (1965). Multiple purpose optimum allocation in stratified sampling. Proc. Amer. Statist. Assoc. Social Statist. Sec. 258–261. [Google Scholar]), Kokan and Khan (1967 Kokan , A. R. , Khan , S. U. ( 1967 ). Optimum allocation in multivariate surveys: An analytical solution . Journal of Royal Statistical Society, Ser. B 29 : 115 – 125 . [Google Scholar]), Chatterjee (1972 Chatterjee , S. ( 1972 ). A study of optimum allocation in multivariate stratified surveys . Skandinavisk Actuarietidskrift 55 : 73 – 80 . [Google Scholar]), Ahsan and Khan (1977 Ahsan , M. J. , Khan , S. U. ( 1977 ). Optimum allocation in multivariate stratified random sampling using prior information . Journal of Indian Statistical Association 15 : 57 – 67 . [Google Scholar], 1982 Ahsan , M. J. , Khan , S. U. ( 1982 ). Optimum allocation in multivariate stratified random sampling with overhead cost . Metrika 29 : 71 – 78 .[Crossref] , [Google Scholar]), Chromy (1987 Chromy , J. R. ( 1987 ). Design optimization with multiple objectives. Proceedings of the Survey Research Methods, 194–199 . [Google Scholar]), Wywial (1988 Wywial , J. ( 1988 ). Minimizing the spectral radius of means vector from sample variance-covariance matrix sample allocation between strata. Prace Naukowe Akademii Ekonomicznej we Wroclawiu 404:223–235 (in Polish) . [Google Scholar]), Bethel (1989 Bethel , J. ( 1989 ). Sample allocation in multivariate surveys . Survey Methodology 15 : 47 – 57 . [Google Scholar]), Kreienbrock (1993 Kreienbrock , L. ( 1993 ). Generalized measures of dispersion to solve the allocation problem in multivariate stratified random sampling . Communication in Statistics—Theory and Methds 22 : 219 – 239 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), Jahan et al. (1994 Jahan , N. , Khan , M. G. M. , Ahsan , M. J. ( 1994 ). A generalized compromise allocation . Journal of the Indian Statistical Association 32 : 95 – 101 . [Google Scholar]), Khan et al. (1997 Khan , M. G. M. , Ahsan , M. J. , Jahan , N. ( 1997 ). Compromise allocation in multivariate stratified sampling: An integer solution . Naval Research Logistics 44 : 69 – 79 .[Crossref], [Web of Science ®] , [Google Scholar]), Khan et al. (2003 Khan , M. G. M. , Khan , E. A. , Ahsan , M. J. ( 2003 ). An optimal multivariate stratified sampling design using dynamic programming . Australian & New Zealand J. Statist. 45 : 107 – 113 .[Crossref], [Web of Science ®] , [Google Scholar]), Ahsan et al. (2005 Ahsan , M. J. , Najmussehar, Khan , M. G. M. ( 2005 ). Mixed allocation in stratified sampling . Aligarh Journal of Statistics 25 : 87 – 97 . [Google Scholar]), Díaz-García and Ulloa (2006 Díaz-García , J. A. , Ulloa , C. L. ( 2006 ). Optimum allocation in multivariate stratified sampling: Multi-objective programming. Comunicación Técnica No. I-06-07/28-03-206 (PE/CIMAT), Guanajuato, México . [Google Scholar], 2008 Díaz-García , J. A. , Ulloa , C. L. ( 2008 ). Multi-objective optimization for optimum allocation in multivariate stratified sampling . Survey Methodology 34 : 215 – 222 .[Web of Science ®] , [Google Scholar]), Ahsan et al. (2009 Ansari , A. H. , Najmussehar, Ahsan , M. J. ( 2009 ). On multiple response stratified random sampling design . International Journal of Statistical Sciences , Kolkata, India, 1(1):1–11 . [Google Scholar]) etc. used different compromise criteria to work out a compromise allocation that is optimum for all characteristics in some sense.

Almost all the previous authors used some function of the sampling variances of the estimators of various characteristics to be measured as an objective that is to be minimized for a fixed cost given as a linear function of sample allocations. Because the variances are not unit free it is more logical to consider the minimization of some function of squared coefficient of variations as an objective. Previously this concept was used by Kozok (2006 Kozok , M. ( 2006 ). On sample allocation in multivariate surveys . Communication in Statistics—Simulation and Computation 35 : 901 – 910 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]).

Furthermore, investigators have to approach the sampled units in order to get the observations. This involves some travel cost. Usually this cost is neglected while constructing a cost function. This travel cost may be significant in some surveys. For example if the strata consist of some geographically difficult-to-approach areas.

The authors problem of optimum allocation in multivariate stratified sampling is discussed with an objective to minimize simultaneously the coefficients of variation of the estimators of various characteristics under a cost constraint that includes the measurement as well as travel cost. The formulated problem of obtaining an optimum compromise allocation turns out to be a multiobjective all-integer nonlinear programming problem. Three different approaches are considered: the value function approach, ∈ –constraint method, and Distance–based method, to obtain compromise allocations. The cost function considered also includes the travel cost within stratum to reach the selected units. Additional restrictions are placed on the sample sizes to avoid oversampling and ensure the availability of the estimates of the strata variances. Numerical examples are also presented to illustrate the computational details of the proposed methods.     

Marginal Correlation from Logit- and Probit-Beta-Normal Models for Hierarchical Binary Data

In hierarchical data settings, be it of a longitudinal, spatial, multi-level, clustered, or otherwise repeated nature, often the association between repeated measurements attracts at least part of the scientific interest. Quantifying the association frequently takes the form of a correlation function, including but not limited to intraclass correlation. Vangeneugden et al. (2010 Vangeneugden, T., Molenberghs, G., Laenen, A., Geys, H., Beunckens, C., Sotto, C. (2010). Marginal correlation in longitudinal binary data based on generalized linear mixed models. Communi. Stati. Theory &; Methods. 39:3542–3557. [Google Scholar]) derived approximate correlation functions for longitudinal sequences of general data type, Gaussian and non-Gaussian, based on generalized linear mixed-effects models. Here, we consider the extended model family proposed by Molenberghs et al. (2010 Molenberghs, G., Verbeke, G., Demétrio, C., Vieira, A. (2010). A family of generalized linear models for repeated measures with normal and conjugate random effects. Stat. Sci. 25:325–347.[Crossref], [Web of Science ®] , [Google Scholar]). This family flexibly accommodates data hierarchies, intra-sequence correlation, and overdispersion. The family allows for closed-form means, variance functions, and correlation function, for a variety of outcome types and link functions. Unfortunately, for binary data with logit link, closed forms cannot be obtained. This is in contrast with the probit link, for which such closed forms can be derived. It is therefore that we concentrate on the probit case. It is of interest, not only in its own right, but also as an instrument to approximate the logit case, thanks to the well-known probit-logit ‘conversion.’ Next to the general situation, some important special cases such as exchangeable clustered outcomes receive attention because they produce insightful expressions. The closed-form expressions are contrasted with the generic approximate expressions of Vangeneugden et al. (2010 Vangeneugden, T., Molenberghs, G., Laenen, A., Geys, H., Beunckens, C., Sotto, C. (2010). Marginal correlation in longitudinal binary data based on generalized linear mixed models. Communi. Stati. Theory &; Methods. 39:3542–3557. [Google Scholar]) and with approximations derived for the so-called logistic-beta-normal combined model. A simulation study explores performance of the method proposed. Data from a schizophrenia trial are analyzed and correlation functions derived.     

Nonparametric Testing of an Exclusion Restriction in Quantile Regression

Using the framework proposed by Bickel et al. (2006 Bickel , P. J. , Ritov , Y. , Stoker , T. ( 2006 ). Tailor-made tests for goodness-of-fit to semiparametric hypotheses . Ann. Stat. 34 ( 2 ): 721 – 741 . [Google Scholar]), we provide a score-based testing method to check the exclusion restriction in quantile regression, i.e., H: ν_α(Y|U, V) = ν_α(Y|U) w.p.1, where ν_α denotes the αth (0 < α < 1) quantile. A subsampling method is suggested to acquire the critical values and justified. The tests are all found to be consistent against fixed alternatives and have discriminating power against local alternatives at root-n scale. We address this particular problem as a representative among a wide family of semiparametric model checking problems. The methodology can be carried over to other goodness-of-fit testing of semiparametric models, possibly involve non smooth functions.     

Hartley-Ross Type Estimators for Population Mean Using Known Parameters of Auxiliary Variate

This article proposes Hartley-Ross type unbiased estimators of finite population mean using information on known parameters of auxiliary variate when the study variate and auxiliary variate are positively correlated. The variances of the proposed unbiased estimators are obtained. It has been shown that the proposed estimators are more efficient than the simple mean estimator, usual ratio estimator and estimators proposed by Sisodia and Dwivedi (1981 Sisodia , B. V. S. , Dwivedi , V. K. ( 1981 ). A modified ratio estimator using coefficient of variation of auxiliary variable . J. Indian Soc. Agricultural Statist. 33 ( 1 ): 13 – 18 . [Google Scholar]), Kadilar and Cingi (2006 Kadilar , C. , Cingi , H. ( 2006 ). A new ratio estimator using correlation coefficient . Int. Statist. 1 – 11 . [Google Scholar]), and Kadilar et al. (2007 Kadilar , C. , Candan , M. , Cingi , H. ( 2007 ). Ratio estimators using robust regression . Hacet. J. Math. Statist. 36 ( 2 ): 181 – 188 .[Web of Science ®] , [Google Scholar]) under certain realistic conditions. Empirical studies are also carried out to demonstrate the merits of the proposed unbiased estimators over other estimators considered in this article.     

Random Effects,Fixed Effects and Hausman's Test for the Generalized Mixed Regressive Spatial Autoregressive Panel Data Model

This article suggests random and fixed effects spatial two-stage least squares estimators for the generalized mixed regressive spatial autoregressive panel data model. This extends the generalized spatial panel model of Baltagi et al. (2013 Baltagi, B. H., Egger, P., Pfaffermayr, M. (2013). A generalized spatial panel data model with random effects. Econometric Reviews 32:650–685.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) by the inclusion of a spatial lag term. The estimation method utilizes the Generalized Moments method suggested by Kapoor et al. (2007 Kapoor, M., Kelejian, H. H., Prucha, I. R. (2007). Panel data models with spatially correlated error components. Journal of Econometrics 127(1):97–130.[Crossref], [Web of Science ®] , [Google Scholar]) for a spatial autoregressive panel data model. We derive the asymptotic distributions of these estimators and suggest a Hausman test a la Mutl and Pfaffermayr (2011 Mutl, J., Pfaffermayr, M. (2011). The Hausman test in a Cliff and Ord panel model. Econometrics Journal 14:48–76.[Crossref], [Web of Science ®] , [Google Scholar]) based on the difference between these estimators. Monte Carlo experiments are performed to investigate the performance of these estimators as well as the corresponding Hausman test.     

Methods for Bivariate Survival Data with Mismeasured Covariates Under an Accelerated Failure Time Model

Accelerated failure time models are useful in survival data analysis, but such models have received little attention in the context of measurement error. In this paper we discuss an accelerated failure time model for bivariate survival data with covariates subject to measurement error. In particular, methods based on the marginal and joint models are considered. Consistency and efficiency of the resultant estimators are investigated. Simulation studies are carried out to evaluate the performance of the estimators as well as the impact of ignoring the measurement error of covariates. As an illustration we apply the proposed methods to analyze a data set arising from the Busselton Health Study (Knuiman et al., 1994 Knuiman , M. W. , Cullent , K. J. , Bulsara , M. K. , Welborn , T. A. , Hobbs , M. S. T. ( 1994 ). Mortality trends, 1965 to 1989, in Busselton, the site of repeated health surveys and interventions . Austral. J. Public Health 18 : 129 – 135 . [CSA] [Crossref], [PubMed] , [Google Scholar]).     

More on the Kronecker Structured Covariance Matrix

In this article, the multivariate normal distribution with a Kronecker product structured covariance matrix is studied. Particularly focused is the estimation of a Kronecker structured covariance matrix of order three, the so called double separable covariance matrix. The suggested estimation generalizes the procedure proposed by Srivastava et al. (2008 Srivastava , M. , von Rosen , T. , von Rosen , D. ( 2008 ). Models with a Kronecker product covariance structure: Estimation and testing Mathemat. Meth. Statist. 17 : 357 – 370 .[Crossref] , [Google Scholar]) for a separable covariance matrix. The restrictions imposed by separability and double separability are also discussed.     

Testing for Unit Roots and the Impact of Quadratic Trends,with an Application to Relative Primary Commodity Prices

In practice a degree of uncertainty will always exist concerning what specification to adopt for the deterministic trend function when running unit root tests. While most macroeconomic time series appear to display an underlying trend, it is often far from clear whether this component is best modeled as a simple linear trend (so that long-run growth rates are constant) or by a more complicated nonlinear trend function which may, for instance, allow the deterministic trend component to evolve gradually over time. In this article, we consider the effects on unit root testing of allowing for a local quadratic trend, a simple yet very flexible example of the latter. Where a local quadratic trend is present but not modeled, we show that the quasi-differenced detrended Dickey–Fuller-type test of Elliott et al. (1996 Elliott , G. , Rothenberg , T. J. , Stock , J. H. ( 1996 ). Efficient tests for an autoregressive unit root . Econometrica 64 : 813 – 836 .[Crossref], [Web of Science ®] , [Google Scholar]) has both size and power which tend to zero asymptotically. An extension of the Elliott et al. (1996 Elliott , G. , Rothenberg , T. J. , Stock , J. H. ( 1996 ). Efficient tests for an autoregressive unit root . Econometrica 64 : 813 – 836 .[Crossref], [Web of Science ®] , [Google Scholar]) approach to allow for a quadratic trend resolves this problem but is shown to result in large power losses relative to the standard detrended test when no quadratic trend is present. We consequently propose a simple and practical approach to dealing with this form of uncertainty based on a union of rejections-based decision rule whereby the unit root is rejected whenever either of the detrended or quadratic detrended unit root tests rejects. A modification of this basic strategy is also suggested which further improves on the properties of the procedure. An application to relative primary commodity price data highlights the empirical relevance of the methods outlined in this article. A by-product of our analysis is the development of a test for the presence of a quadratic trend which is robust to whether the data admit a unit root.