Improved P-Values for Testing Marginal Homogeneity in 2 × 2 Contingency Tables

This article considers the problem of testing marginal homogeneity in a 2 × 2 contingency table. We first review some well-known conditional and unconditional p-values appeared in the statistical literature. Then we treat the p-value as the test statistic and use the unconditional approach to obtain the modified p-value, which is shown to be valid. For a given nominal level, the rejection region of the modified p-value test contains that of the original p-value test. Some nice properties of the modified p-value are given. Especially, under mild conditions the rejection region of the modified p-value test is shown to be the Barnard convex set as described by Barnard (1947 Barnard , G. A. ( 1947 ). Significance tests for 2 × 2 tables . Biometrika 34 : 123 – 138 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]). If the one-sided null hypothesis has two nuisance parameters, we show that this result can reduce the dimension of the nuisance parameter space from two to one for computing modified p-values and sizes of tests. Numerical studies including an illustrative example are given. Numerical comparisons show that the sizes of the modified p-value tests are closer to a nominal level than those of the original p-value tests for many cases, especially in the case of small to moderate sample sizes.     

Bayesian P-Values for Testing Independence in 2 × 2 Contingency Tables

This article considers Bayesian p-values for testing independence in 2 × 2 contingency tables with cell counts observed from the two independent binomial sampling scheme and the multinomial sampling scheme. From the frequentist perspective, Fisher's p-value (p _F ) is the most commonly used p-value but it can be conservative for small to moderate sample sizes. On the other hand, from the Bayesian perspective, Bayarri and Berger (2000 Bayarri , M. J. , Berger , J. O. ( 2000 ). P-values for composite null models (with discussion) . J. Amer. Statist. Assoc. 95 : 1127 – 1170 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) first proposed the partial posterior predictive p-value (p _PPOST ), which can avoid the double use of the data that occurs in another Bayesian p-value proposed by Guttman (1967 Guttman , I. ( 1967 ). The use of the concept of a future observation in goodness-of-fit problems . J. Roy. Statist. Soc. Ser. B 29 : 83 – 100 . [Google Scholar]) and Rubin (1984 Rubin , D. B. ( 1984 ). Bayesianly justifiable and relevant frequency calculations for the applied statistician . Ann. Statist. 12 : 1151 – 1172 .[Crossref], [Web of Science ®] , [Google Scholar]), called the posterior predictive p-value (p _POST ). The subjective and objective Bayesian p-values in terms of p _POST and p _PPOST are derived under the beta prior and the (noninformative) Jeffreys prior, respectively. Numerical comparisons among p _F , p _POST , and p _PPOST reveal that p _PPOST performs much better than p _F and p _POST for small to moderate sample sizes from the frequentist perspective.     

Generalized Inferences on the Common Mean in MANOVA Models

This article presents procedures for testing hypothesis and interval estimation of the common mean vector in MANOVA models when the covariance matrices are unknown and unequal. The methods are based on the concepts of generalized p-value and generalized confidence interval. Some important statistical properties of the exact test and confidence region are given. For two multivariate normal populations, a minor modification to the combined tests given by Zhou and Mathew (1994a Zhou , L. P. , Mathew , T. ( 1994a ). Combining independent tests in multivariate linear models . J. Multivariate Anal. 51 : 265 – 276 . [Google Scholar]) is proposed. Some simulation results to compare the performance of the proposed tests with others are reported. The simulation results indicate that new tests have significant gain in the power.     

Improved Likelihood Inference on Testing the Difference of Non Centrality Parameters of Two Independent Non central t Distributions with Identical Degrees of Freedom

Kim and Hayter (2008 Kim , J. , Hayter , A. J. ( 2008 ). Testing the equality of the non centrality parameters of the two non central t-distributions with identical degrees of freedom . Communications in Statistics: Simulation and Computation 37 : 1709 – 1717 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) showed that testing the difference of the non centralities of two non central t distributions, which have identical degrees of freedom, is the equivalent to testing the difference of two signal-to-noise ratios for two independent samples from different normal distributions. Moreover, from a biomedical research point of view, the above test is also equivalent to testing the existence of the effect size when the sample sizes of the control group and treatment group are the same. In this article, two third-order likelihood-based methods are proposed to approximate the p-value for this significance test problem. A real life example is used to illustrate that results from the proposed methods and the existing methods can be quite different. Simulation results show that the proposed methods give remarkable accuracy even when sample size reaches the smallest possible value.     

One-Sided Estimating and Testing Problems for Scale Models from Grouped Samples

Abstract

This article mainly analyzes estimating and testing problems for scale models from grouped samples. Suppose the support region of a density function, which does not depend on parameters, is divided into some disjoint intervals, grouped samples are the number of observations falling in each intervals respectively. The studying of grouped samples may be dated back to the beginning of the century, in which only one sample location and/or scale models is considered. (Shi, N.-Z., Gao, W., Zhang, B.-X. (2001 Shi, N.-Z., Gao, W. and Zhang, B.-X. 2001. One-sided estimating and testing problems for location models from grouped samples. Comm. Statist.—Simul. Comput, 30(4): 895–898.  [Google Scholar]). One-sided estimating and testing problems for location models from grouped samples. Comm. Statist.—Simul. Comput. 30(4)) had investigated one-sided problems for location models, this article discusses one-sided estimating and testing problems for scale models. Some algorithms for obtaining the maximum likelihood estimates of the parameters subject to order restrictions are proposed.     

Bootstrap approach to test the homogeneity of order restricted mean vectors when the covariance matrices are unknown

Testing homogeneity of multivariate normal mean vectors under an order restriction when the covariance matrices are unknown, arbitrary positive definite and unequal are considered. This problem of testing has been studied to some extent, for example, by Kulatunga and Sasabuchi (1984 Kulatunga, D. D. S., Sasabuchi, S. (1984). A test of homogeneity of mean vectors against multivariate isotonic alternatives. Mem Fac Sci, Kyushu Univ Ser A Mathemat 38:151–161. [Google Scholar]) when the covariance matrices are known and also Sasabuchi et al. (2003 Sasabuchi, S., Tanaka, K., Tsukamodo, T. (2003). Testing homogeneity of multivariate normal mean vectors under an order restriction when the covariance matrices are common but unknown. Annals of Statistics. 31(5):1517–1536.[Web of Science ®] , [Google Scholar]) and Sasabuchi (2007 Sasabuchi, S. (2007). More powerful tests for homogeneity of multivariate normal mean vectors under an order restriction. Sankhya 69(4):700–716. [Google Scholar]) when the covariance matrices are unknown but common. In this paper, a test statistic is proposed and because of the main advantage of the bootstrap test is that it avoids the derivation of the complex null distribution analytically, a bootstrap test statistic is derived and since the proposed test statistic is location invariance the bootstrap p-value defined logical and some steps are presented to estimate it. Our numerical studies via Monte Carlo simulation show that the proposed bootstrap test can correctly control the type I error rates. The power of the test for some of the p-dimensional normal distributions is computed by Monte Carlo simulation. Also, the null distribution of test statistic is estimated using kernel density. Finally, the bootstrap test is illustrated using a real data.     

Two-sided Confidence Bands with the Inferential Sensitivity of One-Sided Bands

In this article, we consider investigating whether any of k treatments are better than a control under the assumption of each treatment mean being no less than the control mean. A classic problem is to find the simultaneous confidence bounds for the difference between each treatment and the control. Compared with hypothesis testing, confidence bounds have the attractive advantage of telling more information about the effective treatment. Generally, the one-sided lower bounds are provided as it's enough for detecting effective treatment and the one-sided lower bounds has sharper lower bands than two-sided ones. However, a two-sided procedure provides both upper and lower bounds on the differences. In this article, we develop a new procedure which combines the good aspects of both the one-sided and the two-sided procedures. This new procedure has the same inferential sensitivity of the one-sided procedure proposed by Zhao (2007 Zhao , H. B. ( 2007 ). Comparing several treatments with a control . J. Statist. Plann. Infer. 137 : 2996 – 3006 .[Crossref], [Web of Science ®] , [Google Scholar]) while also providing simultaneous two-sided bounds for the differences between treatments and the control. By our computation results, we find the new procedure is better than Hayter, Miwa and Liu's procedure (Hayter et al., 2000 Hayter , A. J. , Miwa , T. , Liu , W. ( 2000 ). Combining the advantages of one-sided and two-sided procedures for comparing several treatments with a control . J. Statist. Plann. Infer. 86 : 81 – 99 .[Crossref], [Web of Science ®] , [Google Scholar]), when the sample size is balanced. We also illustrate the new procedure by an example.     

Asymptotics of the Lp-Norms of Density Estimators in the Nonlinear Autoregressive Models

This article investigates the asymptotic behavior of the error density function in nonlinear autoregressive stationary time series regression models. For any 1 ? p < ∞, the kernel density estimator of residuals is shown to be consistent for the error estimator concerning the L_p-distance, which extends the result developed by Cheng and Sun (2008 Cheng, F. X. (2005). Asymptotic distributions of error density estimators in first-order autoregression models. Sankhy Ind. J. Statist. 67:553–567. [Google Scholar]) in L₂-norm. Moreover, the result developed in this article is extended the results of Horváth and Zitikis (2003 Horváth, L., Zitikis, R. (2003). Asymptotics of the L_p-norms of density estimators in the first-order autoregressive models. Statist. Probab. Lett. 65:331342.[Crossref], [Web of Science ®] , [Google Scholar]) to nonlinear autoregressive models.     

Improved Predictive Inference for Responses Based on Location Parameter Median from the Rayleigh Model Given a Doubly Censored Sample

The purpose of this article is to investigate the predictive inference for responses from the location parameter mean as well as from the median given a doubly censored sample from the two-parameter Rayleigh model. The predictive results by Khan et al. (2010 Khan , H. M. R. , Provost , S. B. , Ashima , S. ( 2010 ). Predictive Inference from a Two-Parameter Rayleigh Life Model Given a Doubly Censored Sample . Commun. Statist. Theor. Meth. 39 ( 7 ): 1237 – 1246 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) are used to obtain the predictive inference for responses from the median, where Khan et al. (2010 Khan , H. M. R. , Provost , S. B. , Ashima , S. ( 2010 ). Predictive Inference from a Two-Parameter Rayleigh Life Model Given a Doubly Censored Sample . Commun. Statist. Theor. Meth. 39 ( 7 ): 1237 – 1246 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) obtained the future estimates from the mean. A numerical example representing 66 liver cancer patients is used for predictive analysis. It is concluded that the predictive inference from the median gives precise results as compared with the location parameter mean.     

A Note on Testing Homogeneity of Several Exponential Location Parameters

In this article, the problems of testing homogeneity of several exponential location parameters against simple and tree ordered alternatives are considered separately. Test procedures for both the alternatives are proposed using restricted maximum likelihood estimators (RMLE) of exponential location parameters under the respective orderings. Critical constants for the implementation of the proposed procedures are tabulated. Power comparison of the proposed test procedure under the simple ordered alternative with the procedure of Chen (1982 Chen , H. J. ( 1982 ). A new range statistic for comparisons of several exponential location parameters . Biometrika 69 ( 1 ): 257 – 260 .[Crossref], [Web of Science ®] , [Google Scholar]) and of Dhawan and Gill (1999 Dhawan , A. K. , Gill , A. N. ( 1999 ). A one-sided test for testing homogeneity of scale parameters against ordered alternative . Communication in Statistics—Theory and Methods 28 ( 10 ): 2417 – 2439 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) is carried out using Monte-Carlo simulation.     

Smoothed Mann–Whitney–Wilcoxon Procedure for Two-Sample Location Problem

This study is mainly concerned with estimating a shift parameter in the two-sample location problem. The proposed Smoothed Mann–Whitney–Wilcoxon method smooths the empirical distribution functions of each sample by using convolution technique, and it replaces unknown distribution functions F(x) and G(x ? Δ₀) with the new smoothed distribution functions F _s (x) and G _s (x ? Δ₀), respectively. The unknown shift parameter Δ₀ is estimated by solving the gradient function S _n (Δ) with respect to an arbitrary variable Δ. The asymptotic properties of the new estimator are established under some conditions that are similar to the Generalized Wilcoxon procedure proposed by Anderson and Hettmansperger (1996 Anderson , G. F. , Hettmansperger , T. P. ( 1996 ). Generalized Wilcoxon methods for the one and two-sample location models . In: Brunner , E. , Denker , M. , eds. Research Developments in Probability and Statistics: Festschrift in Honor of Madan L. Puri on the Occasion of his 65th Birthday . Zeist, The Netherlands : VSP BV , pp. 303 – 317 . [Google Scholar]). Some of these properties are asymptotic normality, asymptotic level confidence interval, and hypothesis testing for Δ₀. Asymptotic relative efficiency of the proposed method with respect to the least squares, Generalized Wilcoxon and Hodges and Lehmann (1963 Hodges , J. L. , Lehmann , E. L. ( 1963 ). Estimates of location based on rank tests . Ann. Mathemat. Statist. 34 : 598 – 611 .[Crossref] , [Google Scholar]) procedures are also calculated under the contaminated normal model.     

Generalizing the Marshall's Inequality

In the article, we generalizes the Marshall's inequality in Marshall (1960 Marshall , A. W. ( 1960 ). A one-sided analog of Kolmogorov's inequality . Ann. Math. Statist. 31 : 483 – 487 .[Crossref] , [Google Scholar]) under the condition that the p-order moment exists by a new approach. Some properties of the inequality and further conclusions are given.     

Is the Selected Population the Best?—Location and Scale Parameter Cases

We observe X ₁,…,X _k , where X _i has density f(x,θ _i ) possessing monotone likelihood ratio. The best population corresponds to the largest θ _i . We select the population corresponding to the largest X _i . The goal is to attach the best possible p-value to the inference: the selected population has the uniquely largest θ _i . Gutmann and Maymin (1987 Gutmann , S. , Maymin , Z. ( 1987 ). Is the selected population the best? Ann. Statist . 15 : 456 – 461 .[Crossref], [Web of Science ®] , [Google Scholar]) considered the location parameter case and derived the supremum of the error probability by conditioning on S, the index of the largest X _i . Using this conditioning approach, Kannan and Panchapakesan (2009 Kannan , N. , Panchapakesan , S. ( 2009 ). Does the selected normal population have the smallest variance? Amer. J. Math. Management Sci . 29 : To appear . [Google Scholar]) considered the problem for the gamma family. We consider here a unified approach to both the location and scale parameter cases, and obtain the supremum of the error probability without using conditioning.     

A Simple Method for Testing Variance Components in Unbalanced Nested Model

In this article, we consider the unbalanced case of the three fold nested random effects model under partial balance. The distributions of unweighted sums of squares are obtained first. Using the method of generalized p value introduced in Tsui and Weerahandi (1989 Tsui , K. , Weerahandi , S. ( 1989 ). Generalized p-values in significance testing of hypotheses in the presence of nuisance parameters . Journal of the American Statistical Association 84 : 602 – 607 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), a new method is proposed for hypothesis tests involving functions of variance components. To evaluate the sizes of the generalized p value, a simulation study is conducted. The results indicate that the proposed method performs well under all examined conditions.     

D-Optimal Designs for Covariate Models

The problem of finding D-optimal designs in the presence of a number of covariates has been considered in the one-way set-up. This is an extension of Dey and Mukerjee (2006 Dey , A. , Mukerjee , R. ( 2006 ). D-optimal designs for covariate models . Statistics 40 : 297 – 305 .[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) in the sense that for fixed replication numbers of each treatment, an alternative upper bound to the determinant of the information matrix has been found through completely symmetric C-matrices for the regression coefficients; this upper bound includes the upper bound given in Dey and Mukerjee (2006 Dey , A. , Mukerjee , R. ( 2006 ). D-optimal designs for covariate models . Statistics 40 : 297 – 305 .[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) obtained through diagonal C-matrices. Because of the fact that a smaller class of C-matrices was used at the intermediate stage where the replication numbers were fixed, ultimately some optimal designs remained unidentified there. These designs have been identified here and thereby the conjecture made in Dey and Mukerjee (2006 Dey , A. , Mukerjee , R. ( 2006 ). D-optimal designs for covariate models . Statistics 40 : 297 – 305 .[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) has been settled.     

Paired-Samples Tests of Equivalence

Equivalence tests are used when the objective is to find that two or more groups are nearly equivalent on some outcome, such that any difference is inconsequential. Equivalence tests are available for several research designs, however, paired-samples equivalence tests that are accessible and relevant to the research performed by psychologists have been understudied. This study evaluated parametric and nonparametric two one-sided paired-samples equivalence tests and a standardized paired-samples equivalence test developed by Wellek (2003 Wellek , S. (2003). Testing Statistical Hypotheses of Equivalence . New York : Chapman & Hall/CRC. [Google Scholar]). The two one-sided procedures had better Type I error control and greater power than Wellek's test, with the nonparametric procedure having increased power with non normal distributions.     

The p-values for one-sided hypothesis testing in univariate linear calibration

In this article, we focus on the one-sided hypothesis testing for the univariate linear calibration, where a normally distributed response variable and an explanatory variable are involved. The observations of the response variable corresponding to known values of the explanatory variable are used to make inferences on a single unknown value of the explanatory variable. We apply the generalized inference to the calibration problem, and take the generalized p-value as the test statistic to develop a new p-value for one-sided hypothesis testing, which we refer to as the one-sided posterior predictive p-value. The behavior of the one-sided posterior predictive p-value is numerically compared with that of the generalized p-value, and simulations show that the proposed p-value is quite satisfactory in the frequentist performance.     

Two Statistics for Unconditional Exact Tests in Matched-Pairs Design

Abstract

In a recent article Hsueh et al. (Hsueh, H.-M., Liu, J.-P., Chen, J. J. (2001 Hsueh, H.-M., Liu, J.-P. and Chen, J. J. 2001. Unconditional exact tests for equivalence or noninferiority for paired binary endpoints. Biometrics, 57: 478–483. [Crossref], [PubMed], [Web of Science ®] , [Google Scholar]). Unconditional exact tests for equivalence or noninferiority for paired binary endpoints. Biometrics 57:478–483.) considered unconditional exact tests for paired binary endpoints. They suggested two statistics one of which is based on the restricted maximum-likelihood estimator. Properties of these statistics and the related tests are treated in this article.     

A Note on Criterion-Robust Optimal Designs for Model Discrimination and Parameter Estimation in Polynomial Regression Models

Consider the problem of discriminating between the polynomial regression models on [?1, 1] and estimating parameters in the models. Zen and Tsai (2002 Zen , M. M. , Tsai , M. H. ( 2002 ). Some criterion-robust optimal designs for the dual problem of model discrimination and parameter estimation . Sankhya Ind. J. Statist. 64 : (Series B, Pt. 3) : 322 – 338 . [Google Scholar]) proposed a multiple-objective optimality criterion, M _γ-criterion, which uses weight γ (0 ≤ γ ≤ 1) for model discrimination and α = β = (1 ? γ)/2 for parameter estimation in each model. In this article, we generalize it to a wider setup with different values of α and β. For instance, α = 2 β suggests that the “smaller” model is more likely to be the true model. Using similar techniques, the corresponding criterion-robust optimal design is investigated. A study for the original criterion-robust optimal design with α = β, through M-efficiency, shows that it is good enough for any wider setup.     

Normality of Posterior Distribution Under Misspecification and Nonsmoothness,and Bayes Factor for Davies' Problem

We examine the large sample properties of Bayes procedures in a general framework, where data may be dependent and models may be misspecified and nonsmooth. The posterior distribution of parameters is shown to be asymptotically normal, centered at the quasi maximum likelihood estimator, under mild conditions. In this framework, the Bayes factor for the test problem of Davies (1997, 1987 Davies , R. B. ( 1987 ). Hypothesis testing when a nuisance parameter is present only under the alternative . Biometrika 74 : 33 – 43 .[Web of Science ®] , [Google Scholar]), where a parameter is unidentified under the null hypothesis, is analyzed. The probability that the Bayes factor leads to a correct conclusion about the hypotheses in Davies’ problem is shown to approach to one.