A note on simplicial depth function

We investigate the behaviour of simplicial depth under the perturbation (1−ε)F+ε δ _z , where F is a p-dimensional probability distribution and δ _z is the point-mass distribution concentrated at the point z. The influence function of simplicial depth at the point x, up to a scalar multiplier, turns out to be the difference between the conditional depth, given that one of the vertices of the random simplex is fixed at the position z, and the unconditional depth. The scalar multiplier is p+1, which suggests that simplicial depth can be more sensitive to perturbations as the dimensionality grows higher. The geometrical properties of the influence function give new insight into the observed behaviour of simplicial depth and its relation with halfspace depth. The behaviour of the perturbed simplicial median is also investigated.     

Two-sample tests for sparse high-dimensional binary data

In this article, we study the methods for two-sample hypothesis testing of high-dimensional data coming from a multivariate binary distribution. We test the random projection method and apply an Edgeworth expansion for improvement. Additionally, we propose new statistics which are especially useful for sparse data. We compare the performance of these tests in various scenarios through simulations run in a parallel computing environment. Additionally, we apply these tests to the 20 Newsgroup data showing that our proposed tests have considerably higher power than the others for differentiating groups of news articles with different topics.     

Multivariate option price models and extremes

We consider the Cox-Ross-Rubinstein model of option prices which is a simple binomial model and deal with its multivariate extensions. The model consists of n independent up or down movements of the (multivariate) price. We discuss the model in the view of the limiting distributions for the price as well for the extreme changes of the prices during a period T which is split up into n small price changes, which depend on n (with nh = T). Interesting is also whether the components of the prices and of the extremes are asymptotically dependent.     

Dominating james-stein positive-part estimator for normal mean with unknown covariance matrix

Assume that we have a random sample of size n from p-variate normal population and we wish to estimate the mean vector under quadratic loss with respect to the inverse of the unknown covariance matrix, A class of superior estimators to James-Stein positive part estimator is given when n>max{9p+10,13p-7}, based on the argument by Shao and Strawderman(1994).     

Optimum allocation in multivariate stratified sampling design in the presence of nonresponse with Gamma cost function

In multivariate stratified sample survey with L strata, let p-characteristics are defined on each unit of the population. To estimate the unknown p-population means of each characteristic, a random sample is taken out from the population. In multivariate stratified sample survey, the optimum allocation of any characteristic may not be optimum for others. Thus the problem arises to find out an allocation which may be optimum for all characteristics in some sense. Therefore a compromise criterion is needed to workout such allocation. In this paper, the procedure of estimation of p-population means is discussed in the presence of nonresponse when the use of linear cost function is not advisable. A solution procedure is suggested by using lexicographic goal programming problem. The numerical illustrations are given for its practical utility.     

Extended lasso-type MARS (LMARS) model in the description of biological network

The multivariate adaptive regression splines (MARS) model is one of the well-known, additive non-parametric models that can deal with highly correlated and nonlinear datasets successfully. From our previous analyses, we have seen that lasso-type MARS (LMARS) can be a strong alternative of the Gaussian graphical model (GGM) which is a well-known probabilistic method to describe the steady-state behaviour of the complex biological systems via the lasso regression. In this study, we extend our original LMARS model by taking into account the second-order interaction effects of genes as the representative of the feed-forward loop in biological networks. By this way, we can describe both linear and nonlinear activations of the genes in the same model. We evaluate the performance of our new model under different dimensional simulated and real systems, and then compare the accuracy of the estimates with GGM and LMARS outputs. The results show the advantage of this new model over its close alternatives.     

Unstructured principal fitted response reduction in multivariate regression

In this paper, an unstructured principal fitted response reduction approach is proposed. The new approach is mainly different from two existing model-based approaches, because a required condition is assumed in a covariance matrix of the responses instead of that of a random error. Also, it is invariant under one of popular ways of standardizing responses with its sample covariance equal to the identity matrix. According to numerical studies, the proposed approach yields more robust estimation than the two existing methods, in the sense that its asymptotic performances are not severely sensitive to various situations. So, it can be recommended that the proposed method should be used as a default model-based method.     

Weighted Estimator for the Linear Transformation Models with Multivariate Failure Time Data

In this article, a simple and efficient weighted method is proposed to improve the estimation efficiency for the linear transformation models with multivariate failure time data. Asymptotic properties of the estimators with a closed-form variance-covariance matrix are established. In addition, a goodness-of-fit test is developed to evaluate the adequacy of the model. The performance of proposed method and the comparison on the efficiency between the proposed method and the working independence method (Lu, 2005) are conducted in finite-sample situation by simulation studies. Finally a real data set from the Busselton Population Health Surveys is illustrated to validate the proposed methodology. The related proofs of the theorems are given in the Appendix.     

Diagnosis Aids in Multivariate Multiple Linear Regression Profiles Monitoring

Diagnosis aids in addition to detecting the out-of-control state is an important issue in multivariate multiple linear regression profiles monitoring; because a large number of parameters and profiles in this structure are involved. In this paper, we specifically concentrate on identification of profile(s) and parameter(s) which have changed during the process in multivariate multiple linear regression profiles structure in Phase II. We demonstrate the effectiveness of our proposed approaches through Monte Carlo simulations and a real case study in terms of accuracy percent.     

An Optimality Probability Index for Correlated Multiple Responses

In the real problems, there are many cases which have correlated quality characteristics so multiple response optimization can be more realistic if we can consider correlation structure of responses. In this study we propose a new method which uses multivariate normal probability to find the optimal treatment in an experimental design. Moreover, a heuristic method is used to find better factors’ level in all possible combinations in the designs with large number of controllable factors and their levels. Some simulated numerical examples and a real case were studied by the proposed approach and the comparison of the results with previous methods show efficiency of the proposed method.