A ratio-of-uniforms method for generating exponential power variates

A ratio-of-uniforms method of generating exponential power variates is presented. It is compared to an established generalized rejection method developed by Tadikamalla (1980) and shown to be faster and more easily implemented.     

An error-bounded polynomial approximation for bivariate normal probabilities

Although the bivariate normal distribution is frequently employed in the development of screening models, the formulae for computing bivariate normal probabilities are quite complicated. A simple and accurate error-bounded, noniterative approximation for bivariate normal probabilities based on a simple univariate normal quadratic or cubic approximation is developed for use in screening applications. The approximation, which is most accurate for large absolute correlation coefficients, is especially suitable for screening applications (e.g., in quality control), where large absolute correlations between performance and screening variables are desired. A special approximation for conditional bivariate normal probabilities is also provided which in quality control screening applications improves the accuracy of estimating the average outgoing product quality. Some anomalies in computing conditional bivariate normal probabilities using BNRDF and NORDF in IMSL are also discussed.     

Inequalities for tail probabilities for the multivariate normal distribution

Inequalities for tail probabilities of the multivariate normal distribution are obtained, as a generalization of those given by Feller (1966). Upper and lower bounds are given in the equi-correlated case. For an arbitrary correlation matrix R, an upper bound is obtained, using a result of Slepian (1962) which asserts that certain multivariate normal probabilities are a non-decreasing function of correlations.     

On the estimation of normal copula discrete regression models using the continuous extension and simulated likelihood

The continuous extension of a discrete random variable is amongst the computational methods used for estimation of multivariate normal copula-based models with discrete margins. Its advantage is that the likelihood can be derived conveniently under the theory for copula models with continuous margins, but there has not been a clear analysis of the adequacy of this method. We investigate the asymptotic and small-sample efficiency of two variants of the method for estimating the multivariate normal copula with univariate binary, Poisson, and negative binomial regressions, and show that they lead to biased estimates for the latent correlations, and the univariate marginal parameters that are not regression coefficients. We implement a maximum simulated likelihood method, which is based on evaluating the multidimensional integrals of the likelihood with randomized quasi-Monte Carlo methods. Asymptotic and small-sample efficiency calculations show that our method is nearly as efficient as maximum likelihood for fully specified multivariate normal copula-based models. An illustrative example is given to show the use of our simulated likelihood method.     

Dropping out of university: A statistical analysis of the probability of withdrawal for UK university students

From individual level data for an entire cohort of undergraduate students in the 'old' universities in the UK, we use a binomial probit model to estimate the probability that an individual will 'drop out' of university before the completion of their degree course. We examine the cohort of students enrolling full time for a 3- or 4-year degree in the academic year 1989–1990. We find evidence to support both the hypothesis that the completion of courses by students is influenced by the extent of prior academic preparedness and the hypothesis that social integration at university is important. We also find an influence of unemployment in the county of prior residence, especially for poorer male students. Finally, we draw conclusions regarding the public policy of constructing university performance indicators in this area.     

Simple detection of outlying short time series

Sets of relatively short time series arise in many situations. One aspect of their analysis may be the detection of outlying series. We examine the performance of standard normal outlier tests applied to the means, or to simple functions of the means, of AR(1) series, not necessarily of equal lengths. Although unequal lengths of series implies that the means have unequal variances, that are only known approximately, it is shown that nominal significance levels hold good under most circumstances. Thus a standard outlier test can usefully be applied, avoiding the complication of estimating the time series' parameters. The test's power is affected by unequal lengths, being higher when the slippage occurs in one of the longer series     

THE EFFECTS OF CORRELATION AMONG OBSERVATIONS ON THE ACCURACY OF APPROXIMATING THE DISTRIBUTION OF SAMPLE MEAN BY ITS ASYMPTOTIC DISTRIBUTION1

In this paper, we investigate the effects of correlation among observations on the accuracy of approximating the distribution of sample mean by its asymptotic distribution. The accuracy is investigated by the Berry-Esseen bound (BEB), which gives an upper bound on the error of approximation of the distribution function of the sample mean from its asymptotic distribution for independent observations. For a given sample size (n₀) the BEB is obtained when the observations are independent. Let this be BEB. We then find the sample size (n_*) required to have BEB below BEB₀, when the observations are dependent. Comparison of n_* with n₀ reveals the effects of correlation among observations on the accuracy of the asymptotic distribution as an approximation. It is shown that the effects of correlation among observations are not appreciable if the correlation is moderate to small but it can be severe for extreme correlations.     

A hazard model of the probability of medical school drop-out in the UK

Summary.  From individual level longitudinal data for two entire cohorts of medical students in UK universities, we use multilevel models to analyse the probability that an individual student will drop out of medical school. We find that academic preparedness—both in terms of previous subjects studied and levels of attainment therein—is the major influence on withdrawal by medical students. Additionally, males and more mature students are more likely to withdraw than females or younger students respectively. We find evidence that the factors influencing the decision to transfer course differ from those affecting the decision to drop out for other reasons.     

Testing linear inequality constraints in the standard linear model

In this paper we are concerned with the problem of testing whether the â-parameters of the standard linear model satisfy the linear equality constraints R = r when they are known to satisfy the corresponding linear inequality constraints Râ ? r. In particular we will show that the exact finite sample null distributions of the Likelihood Ratio, Wald and Kuhn-Tucker

statistics are known when R is of full row rank but not known when R has less than full row rank. The less than full row rank problem has not been discussed previously but it is of considerable potential importance.

This paper contains several simple numerical examples which illustrate the computational details of the tests     

SOME SIMPLE APPROXIMATIONS FOR THE DOUBLY NONCENTRAL z DISTRIBUTION

We give two simple approximations for evaluating the cumulative probabilities of the doubly noncentral z distribution. These can easily be used for evaluating the cumulative probabilities of the doubly noncentral F distribution as well. We compare our results with those obtained by Tiku (1965) using series expansion. An industrial situation where a quality characteristic of interest follows the doubly noncentral z distribution is also cited. However, in this case the exact probabilities could be calculated using results on the ratio of two normal variables.     

Application of box-cox transformations to discrimination for the two-class problem

The performance of Box-Cox power transformations in classification using Hinkley's (1975) method is studied. Misclassification probabilities before and after transformation are compared. It is found that the use of Box-Cox transformations can sometimes substantially reduce the error probabilities. Estimates of error probabilities are obtained and certain properties are derived. Examples for a number of distributions are given.     

Measure of Asymmetry for Square Contingency Tables Having Ordered Categories

For the analysis of square contingency tables with nominal categories, Tomizawa and coworkers have considered measures that represent the degree of departure from symmetry. This paper proposes a measure that represents the degree of asymmetry for square contingency tables with ordered categories (instead of those with nominal categories). The measure proposed is expressed using the Cressie–Read power-divergence or Patil–Taillie diversity index, defined for the cumulative probabilities that an observation falls in row (column) category i or below and column (row) category j (> i ) or above. The measure depends on the order of listing the categories. It should be useful for comparing the degree of asymmetry in several tables with ordered categories. The relationship between the measure and the normal distribution is shown.     

Testing for or against a trend in odds ratios

Consider data arranged into k × 2 × 2 contingency tables. The principal result of this paper is the derivation of the likelihood ratio test and its asymptotic distribution for testing for or against an order restriction placed upon the odds ratios. We will show that the limiting distributions are of chi-bar square type and give the expression of the weighting values.     

非农业户籍会使人更幸福吗

中国的户籍制度相对比较特殊,它曾衍生了大量“农转非”等户口买卖现象,这使人不禁产生疑问：非农业户籍会使人更幸福吗？本文基于北京大学CFPS数据,构建了一个影响幸福程度的包括户籍因素在内的有序因变量模型。我们不但详述了户籍类型与居民幸福程度的内在联系机制,同时还分析众多主客观因素对幸福的影响。研究发现,户籍类型转变的实质是社会行动者基于户籍制度功能而做出的理性选择行为,不同时间和空间的户籍制度功能会引起幸福程度的差异;生活态度、人际和教育等因素也会显著影响幸福程度。因此,户籍制度改革过程中,重视户籍制度的功能而非形式是明智之举;同时,完善农村就业机制、增强个人自信心是提升社会行动者幸福程度的有效手段。     

Exact tests for equality of two proportions: fisher v. bayes

Yates (1984) using theoretical and philosophical arguments claims to have proved that the Fisher exact test for comparing the proportions of two binomial experiments is the best exact test. The present article uses objective and practical arguments to confront the Fisher exact test with a Bayes exact test. Using simulated samples we claim to have proved here the inferiority of the Fisher exact test in relation to a Bayes exact test. The comparison is based on the quality concept of Dawid (1982).     

Unified Conditional Frequentist and Bayesian Testing of Composite Hypotheses

Testing of a composite null hypothesis versus a composite alternative is considered when both have a related invariance structure. The goal is to develop conditional frequentist tests that allow the reporting of data-dependent error probabilities, error probabilities that have a strict frequentist interpretation and that reflect the actual amount of evidence in the data. The resulting tests are also seen to be Bayesian tests, in the strong sense that the reported frequentist error probabilities are also the posterior probabilities of the hypotheses under default choices of the prior distribution. The new procedures are illustrated in a variety of applications to model selection and multivariate hypothesis testing.     

中国股票市场的长期记忆性与趋势预测研究

本文首次在国内利用较新的精准局部似然函数法（Exact Local Whittle）,以上证指数为对象,估计了ARFIMA(p,d,q)模型的长期记忆参数d,并分析了上证指数的趋势性变化。估计结果和稳健性检验均表明,上证指数具有长期记忆性,以上证指数为代表的股票市场并非有效;模拟结果显示,当滚动窗口n=260,带宽m=[n0.65]时,长期记忆参数即估计量d既具备一致性,又具有渐进正态性。在2004年10月8日至2015年11月13日期间,模型给出了8次上涨或下跌的趋势转换信号,其中7次信号是正确的,仅有1次给出了错误信号;股票价格由下跌趋势转为上涨趋势、由上涨趋势转为下跌趋势两种情况相比,记忆参数d在前一种情况时下跌幅度更大,给出的信号更明显。     

The weighted compounding of probabilities from three independent hypothesis tests

Tail probabilities from three independent hypothesis tests can be combined to form a test statistic of the form P₁,P₂ ^θ2,P₃ ^θ3.The null distribution of the combined test statistic is presented and critical values for α=0.01 and 0.05 are provided.The power of this test is discussed for the special case ofthree independent F-tests.