Poisson approximation io the distribution of empiy cells in randomized occupancy

Let M_o denote the number of empty cells when n balls are dropped independently and at random in m cells such that each ball stays in its cell with probability p and falls through with probability 1-p. A Poisson limit is known for M_o when (n/m)→∞. We find a corresponding approximation to the distribution of M_o when m is large but finite, Ihe method is elementary, and yields the rate of convergence to the limit law. Ihe results are new for the classical case (p = 1) as well.     

Comparison of approaches for estimating the probability of coastal flooding

Coastal flooding is typically caused by combinations of extreme water-levels and large waves. Two extreme value methods, one univariate and the other multivariate, have been used for estimating the probability of coastal flooding at an existing flood defence structure and for aiding the design of a new structure. The properties of these two methods are compared in terms of extrapolation, sophistication and use of information for a range of extremal dependence structures. We find that, when applied to the assessment of the safety offered by an existing Dutch dike, the multivariate approach provides the more useful and accurate design information and has the substantial benefits of consistency and reduced statistical analysis when applied to several sites along a Dutch coastline.     

On the sensitivity to non-normality of a test for outliers in linear models

In this note we use the class of exponential power distributions to assess the robustness to non-normality of the test for outliers based on the maximum absolute studentized residual. We find that the significance levels can be quite markedly affected by even moderate departures from normality of the error distribution in a regression model when the sample size is moderately large.     

非线性回归模型参数估计方法研究——以C-D生产函数为例

通过理论分析和蒙特卡罗模拟,对C-D生产函数模型参数的估计方法进行比较研究的结果表明：当误差项满足经典假设时,非线性最小二乘估计量具有与线性最小二乘估计类似的、近似BLUE的特性,且当误差项存在异方差时,用加权非线性最小二乘法也能大大改善估计量的性质。     

Estimation of the multi-criteria worths of the alternatives in a hierarchical structure of comparisons

The aggregated worths of the alternatives, when compared with respect to several criteria, are estimated in a hierarchical comparisons model introduced by Saaty (1980). A multiplicative model is used for the paired comparisons data which are collected in a ratio scale in this set-up in any level of this hierarchy. An iterative scheme is found for the maximum likelihood estimation of the worth parameters in this multiplicative model. The iterative values are shown to be convergent monotonically to the estimates. We also obtain the asymptotic dispersion matrix of the maximum likelihood estimates of the relative worths of the alternatives according to a single criterion as well as those according to the over-all suitability when compared under several criteria. A numerical example is presented to illustrate the method developed in this paper. Simulation techniques are employed to find the average number of iterations required for the convergence of the above iterative scheme.     

典型相关分析的延拓研究

在典型相关分析中,求得典型相关变量的表达式并没有全部完成任务,例如需要确定典型相关变量的个数和变量选择。针对典型相关变量的个数问题,发现了常用的卡方检验和冗余分析方法的不足,进而提出了一种新的算法。针对原始变量的选择问题,提出了三种可能的路径。最后利用人体尺寸数据对相关结论进行了验证。     

On the robustness of the extreme deviate test for a single multivariate outlier against heavy-tailed distributions

In this paper we assess the sensitivity of the multivariate extreme deviate test for a single multivariate outlier to non-normality in the form of heavy tails. We find that the empirical significance levels can be markedly affected by even modest departures from multivariate normality. The effects are particularly severe when the sample size is large relative to the dimension. Finally, by way of example we demonstrate that certain graphical techniques may prove useful in identifying the source of rejection for the multivariate extreme deviate test.     

Using a genetic algorithm-based RAROC model for the performance and persistence of the funds

Assisting fund investors in making better investment decisions when faced with market climate change is an important subject. For this purpose, we adopt a genetic algorithm (GA) to search for an optimal decay factor for an exponential weighted moving average model, which is used to calculate the value at risk combined with risk-adjusted return on capital (RAROC). We then propose a GA-based RAROC model. Next, using the model we find the optimal decay factor and investigate the performance and persistence of 31 Taiwanese open-end equity mutual funds over the period from November 2006 to October 2009, divided into three periods: November 2006–October 2007, November 2007–October 2008, and November 2008–October 2009, which includes the global financial crisis. We find that for three periods, the optimal decay factors are 0.999, 0.951, and 0.990, respectively. The rankings of funds between bull and bear markets are quite different. Moreover, the proposed model improves performance persistence. That is, a fund's past performance will continue into the future.     

Naive method to test the convergence of simulation and its applications in the computation of bankruptcy probability

We analyse a naive method using sample mean and sample variance to test the convergence of simulation. We find this method is valid for identically, independently distributed samples, as well as correlated samples with correlation disappearing in long period. Our simulation results on the approximation to bankruptcy probability (BP) show the naive method compares well with the Half-Width, Geweke and CUSUM methods in terms of accuracy and time cost. There are clear evidences of variance reduction from tail-distribution sampling for all convergence test methods when the true BP is very low.     

Bayesian estimation of the Hurst parameter of fractional Brownian motion

The primary purpose of this study was to find Bayesian estimates for the Hurst dimension of a Fractional Brownian motion with a Beta prior when the process is observed at discrete times. Overestimation is observed though the overestimation is less severe as real H goes up. In addition, the estimated H decreases as Beta parameters go up given an Alpha value. In contrast, the estimated H increases as Alpha parameters go up given a Beta value. For the real-world data, the 2011 daily Taiwan stock index was used and the estimated Hurst index was 0.21.     

Comparing multiple treatments to both positive and negative controls

In the past, most comparison to control problems have dealt with comparing k test treatments to either positive or negative controls. Dasgupta et al. [2006. Using numerical methods to find the least favorable configuration when comparing k test treatments to both positive and negative controls. Journal of Statistical Computation and Simulation 76, 251–265] enumerate situations where it is imperative to compare several test treatments to both a negative as well as a positive control simultaneously. Specifically, the aim is to see if the test treatments are worse than the negative control, or if they are better than the positive control when the two controls are sufficiently apart. To find critical regions for this problem, one needs to find the least favorable configuration (LFC) under the composite null. In their paper, Dasgupta et al. [2006. Using numerical methods to find the least favorable configuration when comparing k test treatments to both positive and negative controls. Journal of Statistical Computation and Simulation 76, 251–265] came up with a numerical technique to find the LFC. In this paper we verify their result analytically. Via Monte Carlo simulation we compare the proposed method to the logical single step alternatives: Dunnett's [1955. A multiple comparison procedure for comparing several treatments with a control. Journal of the American Statistical Association 50, 1096–1121] or the Bonferroni correction. The proposed method is superior in terms of both the Type I error and the marginal power.     

On the behaviour of tests based on sample spacings for moderatesamples

We revisit the question about optimal performance of goodness-of-fit tests based on sample spacings. We reveal the importance of centering of the test-statistic and of the sample size when choosing a suitable test-statistic from a family of statistics based on power transformations of sample spacings. In particular, we find that a test-statistic based on empirical estimation of the Hellinger distance between hypothetical and data-supported distribution does possess some optimality properties for moderate sample sizes. These findings confirm earlier statements about the robust behaviour of the test-statistic based on the Hellinger distance and are in contrast to findings about the asymptotic (when sample size approaches infinity) of statistics such as Moran's and/or Greenwood's statistic. We include simulation results that support our findings.     

Unobserved heterogeneity in Markovian analysis of the size distortion of unit root tests

We confirm that units root tests can exhibit substantial size distortion when breaks in mean are generated by a first-order Markov chain, but unlike previous literature, we find augmentation largely remedies this situation. However, considerable heterogeneity is evident in the size properties of the tests when faced with breaks in mean varying in duration, in number and in position within the sample. This heterogeneity will be hidden when a Markov chain is employed. For instance, when the transition probabilities generate single period outliers, rejection frequencies (RFs) rise substantially with the number of outliers, but augmentation results in approximately nominal RFs. Qualitatively similar results hold when a number of structural breaks are allocated randomly in the central section of the sample. Interestingly, very different behaviour is exposed by a design exploring the impact on the tests of two breaks imposed at a range of fixed intervals, RFs rising when break occur in the extremities of the sample, a situation unaffected by augmentation.     

Block designs robust against the presence of positional effects

The problem considered is to find optimum designs for treatment effects in a block design (BD) setup, when positional effects are also present besides treatment and block effects, but they are ignored while formulating the model. In the class of symmetric balanced incomplete block designs, the Youden square design is shown to be optimal in the sense of minimizing the bias term in the mean squared error (MSE) of the best linear unbiased estimators of the full set of orthonormal treatment contrasts, irrespective of the value of the positional effects.     

Local regression when the responses are Interval-Censored

Conditional expectation imputation and local-likelihood methods are contrasted with a midpoint imputation method for bivariate regression involving interval-censored responses. Although the methods can be extended in principle to higher order polynomials, our focus is on the local constant case. Comparisons are based on simulations of data scattered about three target functions with normally distributed errors. Two censoring mechanisms are considered: the first is analogous to current-status data in which monitoring times occur according to a homogeneous Poisson process; the second is analogous to a coarsening mechanism such as would arise when the response values are binned. We find that, according to a pointwise MSE criterion, no method dominates any other when interval sizes are fixed, but when the intervals have a variable width, the local-likelihood method often performs better than the other methods, and midpoint imputation performs the worst. Several illustrative examples are presented.     

Wald检验不一致性的微分几何解释及改进

Wald检验对于等价的零假设中不同形式的表达式在有限样本的情况下缺乏一致性,而从微分几何的角度来解释这一现象,并发现由于Wald统计量是一个混杂的不恰当的几何量,从而对不同的含参数的等价表达式不具有一致性。同时还展示了芬斯拉（Finsler）测地统计量如何能较为简便的计算出来、它在线性回归模型中的非线性约束条件下如何应用以及两者在什么情况下保持一致,并提出了一种解决Wald检验不一致性的思路。     

Confidence bounds and selection of the t best populations

Confidence statements about location (or scale) parameters associated with K populations, which may be used in making selection decisions about those populations, are investigated. When a subset of fixed size t is selected from the K populations a lower bound is obtained for the minimum selected parameter as a function of the maximum non-selected parameter. Tables are produced for the normal means case when the variance is common but unknown. It is pointed out that these tables may be used to find confidence intervals discussed by Hsu (1984     

Distribution of ranks in a two-sample problem with right censored data

We consider different censoring models for a two-sample and find the joint distribution of the rank vector and number of uncensored observations under each censoring model when the distributions of life times and/or distributions of censoring of censoring variables satisfy the condition for the Lehmann type of alternatives.     

Optimal two-point designs for the michaelis-menten model with heteroscedastic errors

We construct D-optimal designs for the Michaelis-Menten model when the variance of the response depends on the independent variable. However, this dependence is only partially known. A Bayesian approacn is used to find an optimal design by incorporating the prior lnformation about the variance structure. We demonstrate the method for a class of error variance structures and present efficiencies of these optimal designs under prior mis-specifications. In particular, we show that an erroneous assumption on the variance structure for the Michaelis-Menten model can have serious consequences.     

Statistical planning when the parameters in the model are of unequal importance

In this paper, we consider statistical planning of experiments when the parameters in the linear model assumed are divided into disjoint sets; the parameters in one set are more influential than the parameters in the other set and require more precise estimation. Me characterize the plans (or, designs) and discuss some symmetric designs which are easier to find.