On Testing Sample Selection Bias Under the Multicollinearity Problem

ABSTRACT

This paper reviews and extends the literature on the finite sample behavior of tests for sample selection bias. Monte Carlo results show that, when the “multicollinearity problem” identified by Nawata (1993 Nawata , K. ( 1993 ). A note on the estimation of models with sample-selection biases . Economics Letters 42 : 15 – 24 . [CSA] [CROSSREF] [Crossref], [Web of Science ®] , [Google Scholar]) is severe, (i) the t-test based on the Heckman–Greene variance estimator can be unreliable, (ii) the Likelihood Ratio test remains powerful, and (iii) nonnormality can be interpreted as severe sample selection bias by Maximum Likelihood methods, leading to negative Wald statistics. We also confirm previous findings (Leung and Yu, 1996 Leung , S. F. , Yu , S. ( 1996 ). On the choice between sample selection and two-part models . Journal of Econometrics 72 : 197 – 229 . [CSA] [CROSSREF] [Crossref], [Web of Science ®] , [Google Scholar]) that the standard regression-based t-test (Heckman, 1979 Heckman , J. J. ( 1979 ). Sample selection bias as a specification error . Econometrica 47 : 153 – 161 . [CSA] [Crossref], [Web of Science ®] , [Google Scholar]) and the asymptotically efficient Lagrange Multiplier test (Melino, 1982 Melino , A. ( 1982 ). Testing for sample selection bias . Review of Economic Studies 49 : 151 – 153 . [CSA] [Crossref], [Web of Science ®] , [Google Scholar]), are robust to nonnormality but have very little power.     

An Extension of the Epsilon-Skew-Normal Distribution

This article is related with the probabilistic and statistical properties of an parametric extension of the so-called epsilon-skew-normal (ESN) distribution introduced by Mudholkar and Hutson (2000 Mudholkar , G. S. , Hutson , A. D. ( 2000 ). The epsilon-skew-normal distribution for analyzing near-normal data . J. Statist. Plann. Infer. 83 : 291 – 309 .[Crossref], [Web of Science ®] , [Google Scholar]), which considers an additional shape parameter in order to increase the flexibility of the ESN distribution. Also, this article concerns likelihood inference about the parameters of the new class. In particular, the information matrix of the maximum likelihood estimators is obtained, showing that it is non singular in the special normal case. Finally, the statistical methods are illustrated with two examples based on real datasets.     

GMM in linear regression for longitudinal data with multiple covariates measured with error

Griliches and Hausman 5 Griliches, Z. and Hausman, J. A. 1986. Errors in variables in panel data. J. Econometrics, 32: 93–118. [Crossref], [Web of Science ®] [Google Scholar] and Wansbeek 11 Wansbeek, T. J. 2001. GMM estimation in panel data models with measurement error. J. Econometrics, 104: 259–268. [Crossref], [Web of Science ®] [Google Scholar] proposed using the generalized method of moments (GMM) to obtain consistent estimators in linear regression models for longitudinal data with measurement error in one covariate, without requiring additional validation or replicate data. For usefulness of this methodology, we must extend it to the more realistic situation where more than one covariate are measured with error. Such an extension is not straightforward, since measurement errors across different covariates may be correlated. By a careful construction of the measurement error correlation structure, we are able to extend Wansbeek's GMM and show that the extended Griliches and Hausman's GMM is equivalent to the extended Wansbeek's GMM. For illustration, we apply the extended GMM to data from two medical studies, and compare it with the naive method and the method assuming only one covariate having measurement error.     

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.     

A Note on the Prior Distributions of Weibull Parameters for the Reliability Function

In Bayesian Inference it is often desirable to have a posterior density reflecting mainly the information from sample data. To achieve this purpose it is important to employ prior densities which add little information to the sample. We have in the literature many such prior densities, for example, Jeffreys (1967 Jeffreys , H. ( 1967 ). Theory of Probability , 3rd rev. ed. . London : Oxford University Press . [Google Scholar]), Lindley (1956 Lindley , D. V. ( 1956 ). On a measure of the information provided by an experiment . Ann. Mathemat. Statist. 27 : 986 – 1005 .[Crossref] , [Google Scholar]); (1961 Lindley , D. V. ( 1961 ). The use of prior probability distributions in statistical inference and decisions . In: Neyman , J. , ed. Proceedings of the Fourth Berkeley Symposium on Mathematical Statistics and Probability . Vol. 1. Berkeley : University of California Press , pp. 453 – 468 . [Google Scholar]), Hartigan (1964 Hartigan , J. ( 1964 ). Invariant priors distributions . Ann. Mathemat. Statist. 35 : 836 – 845 .[Crossref] , [Google Scholar]), Bernardo (1979 Bernardo , J. M. ( 1979 ). Reference posterior distributions for Bayesian inference . J. Roy. Statist. Soc. 41 ( 2 ): 113 – 147 . [Google Scholar]), Zellner (1984 Zellner , A. ( 1984 ). Maximal Data Information Prior Distributions, Basic Issues in Econometrics . Chicago : University of Chicago Press . [Google Scholar]), Tibshirani (1989 Tibshirani , R. ( 1989 ). Noninformative priors for one parameters of many . Biometrika 76 : 604 – 608 .[Crossref], [Web of Science ®] , [Google Scholar]), etc. In the present article, we compare the posterior densities of the reliability function by using Jeffreys, the maximal data information (Zellner, 1984 Zellner , A. ( 1984 ). Maximal Data Information Prior Distributions, Basic Issues in Econometrics . Chicago : University of Chicago Press . [Google Scholar]), Tibshirani's, and reference priors for the reliability function R(t) in a Weibull distribution.     

Lower Bounds on the Symmetric L 2-Discrepancy and Their Application

The role of uniformity measured by the symmetric L ₂-discrepancy given in Hickernell (1998 Hickernell , F. J. (1998). A generalized discrepancy and quadrature error bound. Math. Computat. 67:299–322.[Crossref], [Web of Science ®] , [Google Scholar]) has been studied in fractional factorial designs. The issue of lower bounds on the symmetric L ₂-discrepancy is crucial in the construction of uniform designs. This article reports some new lower bounds on the symmetric L ₂-discrepancy for symmetric fractional factorials and for a set of asymmetric fractional factorials. It is valuable to use these lower bounds to measure uniformity of given designs.     

A Variance Estimator for Systematic Sampling from a Deliberately Ordered Population

ABSTRACT

The systematic sampling (SYS) design (Madow and Madow, 1944 Madow , L. H. , Madow , W. G. ( 1944 ). On the theory of systematic sampling . Ann. Math. Statist. 15 : 1 – 24 .[Crossref] , [Google Scholar]) is widely used by statistical offices due to its simplicity and efficiency (e.g., Iachan, 1982 Iachan , R. ( 1982 ). Systematic sampling a critical review . Int. Statist. Rev. 50 : 293 – 303 .[Crossref], [Web of Science ®] , [Google Scholar]). But it suffers from a serious defect, namely, that it is impossible to unbiasedly estimate the sampling variance (Iachan, 1982 Iachan , R. ( 1982 ). Systematic sampling a critical review . Int. Statist. Rev. 50 : 293 – 303 .[Crossref], [Web of Science ®] , [Google Scholar]) and usual variance estimators (Yates and Grundy, 1953 Yates , F. , Grundy , P. M. ( 1953 ). Selection without replacement from within strata with probability proportional to size . J. Roy. Statist. Soc. Series B 1 : 253 – 261 . [Google Scholar]) are inadequate and can overestimate the variance significantly (Särndal et al., 1992 Särndal , C. E. , Swenson , B. , Wretman , J. H. ( 1992 ). Model Assisted Survey Sampling . New York : Springer-Verlag , Ch. 3 .[Crossref] , [Google Scholar]). We propose a novel variance estimator which is less biased and that can be implemented with any given population order. We will justify this estimator theoretically and with a Monte Carlo simulation study.     

Use of Two Decks of Cards in Randomized Response Techniques for Complex Survey Designs

The randomized response (RR) technique with two decks of cards proposed by Odumade and Singh (2009 Odumade , O. , Singh , S. ( 2009 ). Efficient use of two deck of cards in randomized response sampling . Commun. Statist. Theor. Meth. 38 : 439 – 446 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) can always be made more efficient than the RR techniques proposed by Warner (1965 Warner , S. L. ( 1965 ). Randomize response: A survey technique for eliminating evasive answer bias . J. Amer. Statist. Assoc. 60 : 63 – 69 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), Mangat and Singh (1990 Mangat , N. S. , Singh , R. ( 1990 ). An alternative randomized response procedure . Biometrika 77 : 349 – 442 .[Crossref], [Web of Science ®] , [Google Scholar]), and Mangat (1994 Mangat , N. S. ( 1994 ). An improved randomized response strategy . J. Roy. Statist. Soc. B 56 : 93 – 95 . [Google Scholar]) by adjusting the proportion of cards in the decks. The proposed method of Odumade and Singh (2009 Odumade , O. , Singh , S. ( 2009 ). Efficient use of two deck of cards in randomized response sampling . Commun. Statist. Theor. Meth. 38 : 439 – 446 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) is limited to simple random sampling with replacement (SRSWR) sampling only. In this article, generalization of Odumade and Singh strategy is provided for complex survey designs and a wider class of estimators. The results of Odumade and Singh (2009 Odumade , O. , Singh , S. ( 2009 ). Efficient use of two deck of cards in randomized response sampling . Commun. Statist. Theor. Meth. 38 : 439 – 446 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) can be derived from the proposed method as a special case.     

Objective Bayesian Analysis for Linear Degradation Models

We propose an objective Bayesian approach to analyze degradation models. For the linear degradation models, two reference priors are derived, and based on this we show the posterior distributions are proper. Since the lifetime of the product is of interest in practice, a transformation is introduced to obtain the reference priors of the medium lifetime. In the posterior analysis, we explore two sampling procedures: Monte Carlo (MC) procedure and Monte Carlo Markov Chain (MCMC) procedure. A real data from Takeda and Suzuki (1983 Takeda , E. , Suzuki , N. ( 1983 ). An empirical model for device degradation due to hot-carrier injection . IEEE Electron Dev. Lett. 4 : 111 – 113 .[Crossref], [Web of Science ®] , [Google Scholar]) is analyzed, and we find the results obtained by both procedures are close to the given literature.     

A New Variant of the Bivariate Setting the Clock Back to Zero Property

Raja Rao et al. (1993 Raja Rao , B. , Damaraju , C. V. , Alhumoud , J. M. ( 1993 ). Setting the clock back to zero property of a class of bivariate life distributions . Commun. Statist. Theor. Meth. 22 ( 7 ): 2067 – 2080 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) introduced the bivariate setting the clock back to zero property. A new variant of this property is introduced that is appropriate for analysing a broader area of practical situations. Some distributions possessing the proposed property are presented. Applications of this property for simplifying the computation of the bivariate mean residual life function and the bivariate percentile residual life function are studied. The relation between the proposed property with the one studied by Raja Rao and Talwalker (1990 Raja Rao , B. , Talwalker , S. (1990). Setting the clock back to zero property of a family of life distributions. J. Statist. Plann. Infer. 24:347–352. [Google Scholar]) and the bivariate lack of memory property is studied.     

A Measure of Multivariate Kurtosis with Principal Components

Srivastava (1984 Srivastava , M. S. ( 1984 ). A measure of skewness and kurtosis and a graphical method for assessing multivariate normality . Statist. Probab. Lett. 2 ( 5 ): 263 – 267 .[Crossref], [Web of Science ®] , [Google Scholar]) defined a measure of multivariate kurtosis and derived its asymptotic distribution for samples from a multivariate normal population. Some new results are obtained by generalizing Srivastava's theorem to an asymptotic expansion up to higher order. Finally, two numerical examples are presented.     

Bayesian Robustness of the Positive False Discovery Rate

The positive false discovery rate was introduced by Storey (2003 Storey , J. D. (2003). The positive false discovery rate: a Bayesian interpretation and the q-value. Ann. Statist. 31:2013–2035.[Crossref], [Web of Science ®] , [Google Scholar]) as an alternative to the family wise error rate for the case in which we are simultaneously testing a large amount of hypotheses. The positive false discovery rate has a very nice Bayesian interpretation (as it was shown by Storey, 2003 Storey , J. D. (2003). The positive false discovery rate: a Bayesian interpretation and the q-value. Ann. Statist. 31:2013–2035.[Crossref], [Web of Science ®] , [Google Scholar]) and its robustness is analyzed. The emphasis is on the ε-contamination class (one of the most used classes of priors for Bayesian robustness) and it is shown that robustness is not obtained when the basic prior concentrates the probability on the null hypothesis.     

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.     

A Class of Improved Estimators for the Scale Parameter of a Mixture Model of Exponential Distribution with Unknown Location

Under Stein's loss, a class of improved estimators for the scale parameter of a mixture of exponential distribution with unknown location is constructed. The method is analogous to Maruyama's (1998 Maruyama , Y. ( 1998 ). Minimax estimators of a normal variance . Metrika 48 : 209 – 214 .[Crossref], [Web of Science ®] , [Google Scholar]) construction for the variance of a normal distribution and also an extension of the result produced in Petropoulos and Kourouklis (2002 Petropoulos , C. , Kourouklis , S. ( 2002 ). A class of improved estimators for the scale parameter of an exponential distribution with unknown location . Commun. Statist. Theor. Meth. 31 : 325 – 335 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). Also, robustness properties are considered.     

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 Study of Expansions of Posterior Distributions

Johnson (1970 Johnson , R. ( 1970 ). Asymptotic expansions associated with posterior distributions . Ann. Math. Statist. 41 : 851 – 864 .[Crossref] , [Google Scholar]) obtained expansions for marginal posterior distributions through Taylor expansions. Here, the posterior expansion is expressed in terms of the likelihood and the prior together with their derivatives. Recently, Weng (2010 Weng , R. C. ( 2010 ). A Bayesian Edgeworth expansion by Stein's Identity . Bayesian Anal. 5 ( 4 ): 741 – 764 .[Crossref], [Web of Science ®] , [Google Scholar]) used a version of Stein's identity to derive a Bayesian Edgeworth expansion, expressed by posterior moments. Since the pivots used in these two articles are the same, it is of interest to compare these two expansions.

We found that our O(t ^?1/2) term agrees with Johnson's arithmetically, but the O(t ^?1) term does not. The simulations confirmed this finding and revealed that our O(t ^?1) term gives better performance than Johnson's.     

Improved strategy to collect sensitive data by using geometric distribution as a randomization device

Abstract

In this article, we improvise Singh and Grewal (2013 Singh, S., and I. S. Grewal. 2013. Geometric distribution as a randomization device implemented in the Kuk’s model. International Journal of Contemporary Mathematical Sciences 8:243–8.[Crossref] , [Google Scholar]) and Hussain et al. (2016 Hussain, Z., J. Shabbir, Z. Pervez, S. F. Shah, and M. Khan. 2016. Generalized geometric distribution of order k: A flexible choice to randomize the response. Communications in Statistics: Simulation and Computation 46:4708–21.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) techniques by introducing a new two-stage randomization response process. Using the proposed new technique, we achieve better efficiency and increasing protection of privacy of respondents than the Kuk (1990 Kuk, A. Y. C. 1990. Asking sensitive questions indirectly. Biometrika 77 (2):436–8.[Crossref], [Web of Science ®] , [Google Scholar]), Singh and Grewal (2013 Singh, S., and I. S. Grewal. 2013. Geometric distribution as a randomization device implemented in the Kuk’s model. International Journal of Contemporary Mathematical Sciences 8:243–8.[Crossref] , [Google Scholar]) and Hussain et al. (2016 Hussain, Z., J. Shabbir, Z. Pervez, S. F. Shah, and M. Khan. 2016. Generalized geometric distribution of order k: A flexible choice to randomize the response. Communications in Statistics: Simulation and Computation 46:4708–21.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) models. The relative efficiency and protection of the respondents of the proposed two-stage randomization device have been investigated through simulation study, and the situations are reported where the proposed estimator performs better than its competitors. The SAS code used to investigate the performance of the proposed strategy are also provided.     

IPPS Sampling Plans Excluding Adjacent Units

The concept of inclusion probability proportional to size sampling plans excluding adjacent units separated by at most a distance of m (≥ 1) units {IPPSEA plans} is introduced. IPPSEA plans ensure that the first-order inclusion probabilities of units are proportional to size measures of the units, while the second-order inclusion probabilities are zero for pairs of adjacent units separated by a distance of m units or less. IPPSEA plans have been obtained by making use of binary, proper, and unequireplicated block designs and linear programing approach. The performance of IPPSEA plans using Horvitz–Thompson estimator of population total has been compared with existing sampling plans such as simple random sampling without replacement (SRSWOR), balanced sampling plans excluding adjacent units {BSA (m) plans}, probability proportional to size with replacement, Hartley and Rao's plan (1962 Hartley , H. O. , Rao , J. N. K. ( 1962 ). Sampling with unequal probabilities and without replacement . Ann. Math. Statist. 33 : 350 – 374 .[Crossref] , [Google Scholar]), Rao et al.'s strategy (1962 Rao , J. N. K. , Hartley , H. O. , Cochran , W. G. ( 1962 ). On a simple procedure of unequal probability sampling without replacement . J. Roy. Statist. Soc. B 24 : 482 – 491 . [Google Scholar]), and Sampford's IPPS plan (1967 Sampford , M. R. ( 1967 ). On sampling without replacement with unequal probabilities of selection . Biometrika 54 ( 3 ): 499 – 513 .[Crossref], [PubMed] , [Google Scholar]) using a real life population. Unbiased estimation of Horvitz–Thompson estimator of population total is not possible in these types of plans because some of the second-order inclusion probabilities are zero. To resolve this problem, one approximate variance estimation technique has been suggested.     

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.     

Note on posterior inference for the Bingham distribution

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