辅以排序集样本的分层抽样方法研究

分层抽样中,样本在各层中的不同获取方式会对估计量的精度和试验费用产生一定的影响,而已有的理论方法大多不能在提高精度的同时降低调查费用。为此,将排序抽样与分层抽样方法相结合,提出了辅以排序集样本的分层抽样方案,并得到了总体均值的估计量以及这一估计量的良好性质。这些结果表明,与单一的分层随机抽样相比,这种抽样设计的估计量具有更高的精度,同时也节约了各层抽样调查的费用。     

Nearly exact tests in factorial experiments using the aligned rank transform

A procedure is studied that uses rank-transformed data to perform exact and estimated exact tests, which is an alternative to the commonly used F-ratio test procedure. First, a common parametric test statistic is computed using rank-transformed data, where two methods of ranking-ranks taken for the original observations and ranks taken after aligning the observations-are studied. Significance is then determined using either the exact permutation distribution of the statistic or an estimate of this distribution based on a random sample of all possible permutations. Simulation studies compare the performance of this method with the normal theory parametric F-test and the traditional rank transform procedure. Power and nominal type I error rates are compared under conditions when normal theory assumptions are satisfied, as well as when these assumptions are violated. The method is studied for a two-factor factorial arrangement of treatments in a completely randomized design and for a split-unit experiment. The power of the tests rivals the parametric F-test when normal theory assumptions are satisfied, and is usually superior when normal theory assumptions are not satisfied. Based on the evidence of this study, the exact aligned rank procedure appears to be the overall best choice for performing tests in a general factorial experiment.     

Simultaneous rank tests for detecting differentially expressed genes

Rank tests are known to be robust to outliers and violation of distributional assumptions. Two major issues besetting microarray data are violation of the normality assumption and contamination by outliers. In this article, we formulate the normal theory simultaneous tests and their aligned rank transformation (ART) analog for detecting differentially expressed genes. These tests are based on the least-squares estimates of the effects when data follow a linear model. Application of the two methods are then demonstrated on a real data set. To evaluate the performance of the aligned rank transform method with the corresponding normal theory method, data were simulated according to the characteristics of a real gene expression data. These simulated data are then used to compare the two methods with respect to their sensitivity to the distributional assumption and to outliers for controlling the family-wise Type I error rate, power, and false discovery rate. It is demonstrated that the ART generally possesses the robustness of validity property even for microarray data with small number of replications. Although these methods can be applied to more general designs, in this article the simulation study is carried out for a dye-swap design since this design is broadly used in cDNA microarray experiments.     

Estimation of a population size through capture-mark-recapture method: a comparison of various point and interval estimators

This article deals with the estimation of a fixed population size through capture-mark-recapture method that gives rise to hypergeometric distribution. There are a few well-known and popular point estimators available in the literature, but no good comprehensive comparison is available about their merits. Apart from the available estimators, an empirical Bayes (EB) estimator of the population size is proposed. We compare all the point estimators in terms of relative bias and relative mean squared error. Next, two new interval estimators – (a) an EB highest posterior distribution interval and (b) a frequentist interval estimator based on a parametric bootstrap method, are proposed. The comparison is then carried among the two proposed interval estimators and interval estimators derived from the currently available estimators in terms of coverage probability and average length (AL). Based on comprehensive numerical results, we rank and recommend the point estimators as well as interval estimators for practical use. Finally, a real-life data set for a green treefrog population is used as a demonstration for all the methods discussed.     

A test for bivariate two sample location problem

A consistent test for difference in locations between two bivariate populations is proposed, The test is similar as the Mann-Whitney test and depends on the exceedances of slopes of the two samples where slope for each sample observation is computed by taking the ratios of the observed values. In terms of the slopes, it reduces to a univariate problem, The power of the test has been compared with those of various existing tests by simulation. The proposed test statistic is compared with Mardia's(1967) test statistics, Peters-Randies(1991) test statistic, Wilcoxon's rank sum test. statistic and Hotelling' T² test statistic using Monte Carlo technique. It performs better than other statistics compared for small differences in locations between two populations when underlying population is population 7(light tailed population) and sample size 15 and 18 respectively. When underlying population is population 6(heavy tailed population) and sample sizes are 15 and 18 it performas better than other statistic compared except Wilcoxon's rank sum test statistics for small differences in location between two populations. It performs better than Mardia's(1967) test statistic for large differences in location between two population when underlying population is bivariate normal mixture with probability p=0.5, population 6, Pearson type II population and Pearson type VII population for sample size 15 and 18 .Under bivariate normal population it performs as good as Mardia' (1967) test statistic for small differences in locations between two populations and sample sizes 15 and 18. For sample sizes 25 and 28 respectively it performs better than Mardia's (1967) test statistic when underlying population is population 6, Pearson type II population and Pearson type VII population     

On the asymptotic normality of linear rank statistics

A method for proving the asymptotic normality of linear rank statistics is presented. Use is made of the fact that the variances of many rank statistics met in practice admit a factorization expressed in terms of the sample size. The method of proof consists of two steps corresponding to the limiting behavior of the modulus and argument of the characteristic function of the rank statistic.     

Adjusted Wilcoxon signed rank test tables for ratio of percentiles

The ordinary Wilcoxon signed rank test table provides confidence intervals for the median of one population. Adjusted Wilcoxon signed rank test tables which can provide confidence intervals for the median and the 10th percentile of one population are created in this paper. Base-(n + 1) number system and theorems about property of symmetry of the adjusted Wilcoxon signed rank test statistic are derived for programming. Theorem 1 states that the adjusted Wilcoxon signed rank test statistic are symmetric around n(n + 1)/4. Theorem 2 states that the adjusted Wilcoxon signed rank test statistic with the same number of negative ranks m are symmetric around m(n+1)/2. 87.5% and 85% confidence intervals of the median are given in the table for n = 12, 13,…, 29 to create approximated 95% confidence intervals of the ratio of medians for two independent populations. 95% and 92.5% confidence intervals of the 10th percentile are given in the table for n = 26, 27, 28, 29 to create approximated 95% confidence regions of the ratio of the 10th percentiles for two independent populations. Finally two large datasets from wood industry will be partitioned to verify the correctness of adjusted Wilcoxon signed rank test tables for small samples.     

Linear rank tests for homogeneity against ordered alternatives in anova and anocova

A unified method of constructing rank tests for homogeneity against ordered alternatives in unbalanced analysis of variance and analysis of covariance is considered. The relationship between these tests with some of the existing methods are studied. The normal theory likelihood ratio tests are also derived and the asymptotic relative efficiency comparisons, in Pitman sense, of the rank tests with respect to the likelihood ratio tests are carried out.     

Affine Invariant Multivariate Sign and Rank Tests and Corresponding Estimates: a Review

The paper reviews recent contributions to the statistical inference methods, tests and estimates, based on the generalized median of Oja. Multivariate analogues of sign and rank concepts, affine invariant one-sample and two-sample sign tests and rank tests, affine equivariant median and Hodges–Lehmann-type estimates are reviewed and discussed. Some comparisons are made to other generalizations. The theory is illustrated by two examples.     

A robust test for non-nested hypotheses

This paper uses a modified rank score test for non-nested linear regression models. The modified rank score test is robust with respect to models with non-normal distributions and can be viewed as a robust version of the J test of Davidson and MacKinnon (Econometrica 49:781–793, 1981). Therefore, this test does not require a specification of error density function and is easy to implement. Also, a modified rank score test for multiple non-nested models is provided. Monte Carlo simulation results show that the test has good finite sample performances. Financial applications for two competing theories, the capital asset pricing model and the arbitrage pricing theory, are considered herein. Empirical evidence from the modified rank score test shows that the former is a better model for asset pricing.     

Nonparametric repeated significance tests for one-way anova with adoptation to multiple comparisons

Based on two-sample rank order statistics, a repeated significance testing procedure for a multi-sample location problem is considered. The asymptotic distribution theory of the proposed tests is given under the null hypothesis as well as under local alternatives. A Bahadur efficiency result of the repeated significance test relative to the terminal test based solely on the target sample size is presented. In the adaptation of the proposed tests to multiple comparisons, an asymptotically equivalent test statistic in terms of the rank estimators of the location parameters is derived from which the Scheffé method of multiple comparisons can be obtained in a convinient way.     

Nonparametric estimation in one-way random effects models

We present a rank based method for obtaining point and interval estimates of a scale version of the intraclass correlation coefficient in a one-way random effects model. When compared to the method of Arvesen and Schmitz (1970), the new method is not only applicable to a broader class of situations, but also much easier to implement. Results of a simulation study indicate that the new procedure compares favorably with the Arvesen-Schmitz procedure and the classical normal theory procedure especially If the random components have heavy tailed distributions.     

复杂抽样推断方法体系的比较研究

对复杂样本进行推断通常有两种体系,一种是传统的基于随机化理论的统计推断,另一种是基于模型的统计推断。传统的抽样理论以随机化理论为基础,将总体取值视为固定,随机性仅体现在样本的选取上,对总体的推断依赖于抽样设计。该方法在大样本情况下具有稳健估计量,但在小样本、数据缺失等情况下失效。基于模型的抽样推断认为总体是超总体模型中抽取的一个随机样本,对总体的推断取决于模型的建立,但在不可忽略抽样设计下估计量是有偏估计。在对这两类推断方法分析的基础上,提出抽样设计辅助的模型推断,并指出该方法在复杂抽样中具有重要的应用价值。     

Quad folding: a simple idea for the subjective property characterization of large sample sets

Quad analysis has proven useful for characterizing subjective properties, primarily properties such as carpet hand and body in the textile industry. In essence, it provides an efficient method for conducting paired comparisons, the preferred method for quantifying such properties. An extension to quad analysis, quad folding of one quad design into another, is likewise an efficient method to rank order the subjective properties of larger data sets. A rank ordering of 62 carpets by their body is provided as an example of folding six groups of eight carpets (two replicated) into one another.     

Sample Size Calculation for Rank Tests Comparing K Survival Distributions

Rank tests, such as logrank or Wilcoxon rank sum tests, have been popularly used to compare survival distributions of two or more groups in the presence of right censoring. However, there has been little research on sample size calculation methods for rank tests to compare more than two groups. An existing method is based on a crude approximation, which tends to underestimate sample size, i.e., the calculated sample size has lower power than projected. In this paper we propose an asymptotically correct method and an approximate method for sample size calculation. The proposed methods are compared to other methods through simulation studies.     

Multivariate outlier detection in incomplete survey data: the epidemic algorithm and transformed rank correlations

Summary.  As a part of the EUREDIT project new methods to detect multivariate outliers in incomplete survey data have been developed. These methods are the first to work with sampling weights and to be able to cope with missing values. Two of these methods are presented here. The epidemic algorithm simulates the propagation of a disease through a population and uses extreme infection times to find outlying observations. Transformed rank correlations are robust estimates of the centre and the scatter of the data. They use a geometric transformation that is based on the rank correlation matrix. The estimates are used to define a Mahalanobis distance that reveals outliers. The two methods are applied to a small data set and to one of the evaluation data sets of the EUREDIT project.     

模糊Borda法的缺陷分析及其改进思路

 通过对模糊Borda法的分析，笔者认为该方法在技术处理上偏离了评价的基本思想，根本无法实现将评价值和评价序相结合的初衷。为了实现将两者相结合的目的，本文提出了两个解决思路：利用被评价单位在评价值的差距计算相对位次；或计算名次转换值与评价值的交叉乘积。     

Union-intersection rank tests for ordered alternatives in some simple linear models

The theory of locally most powerful rank tests and the union intersection principle are incorporated in the formulation of some distribution-free rank tests for ordered alternatives in some simple linear models (including the classical one-way layout as a special case). Some approximations for the null-hypothesis distributions of the test statistics ( for finite sample sizes) are considered, and the theory is supplemented by some numerical and simulation results.     

Rank estimation of regression coefficients using iterated reweighted least squares

This paper is concerned with the rank estimator for the parameter vector β in a linear model which is obtained by the minimization of a rank dispersion function. The rank estimator has many advantages over the regular least squares estimator, but the inaccessibility of software to carry out its computation has limited its use. An iterated reweighted least squares algorithm is presented for the computation of the rank estimator. The method is simple in concept and can be carried out readily with a wide variety of statistical software. Details of the method are discussed along with some results on its asymptotic distribution and numerical stability. Some examples are presented to show advantages of the rank method.     

Understanding Canonical Correlation through the General Linear Model and Principal Components

Canonical correlation has been little used and little understood, even by otherwise sophisticated analysts. An alternative approach to canonical correlation, based on a general linear multivariate model, is presented. Properties of principal component analysis are used to help explain the method. Standard computational methods for full rank canonical correlation, techniques for canonical correlation on component scores, and canonical correlation with less than full rank are discussed. They are seen to be essentially equivalent when the model equation for canonical correlation on component scores is presented. The two approaches to less than full rank situations are equivalent in some senses, but quite different in usefulness, depending on the application. An example dataset is analyzed in detail to help demonstrate the conclusions.