Hayter and Tsui's test with double sampling for the vector mean of multivariate normal population

Abstract

In this paper, we introduce a version of Hayter and Tsui's statistical test with double sampling for the vector mean of a population under multivariate normal assumption. A study showed that this new test was more or as efficient than the well-known Hotelling's T² with double sampling. Some nice features of Hayter and Tsui's test are its simplicity of implementation and its capability of identifying the errant variables when the null hypothesis is rejected. Taking that into consideration, a new control chart called HTDS is also introduced as a tool to monitor multivariate process vector mean when using double sampling.     

Inference in log-alpha-power and log-skew-normal multivariate models

ABSTRACT

Random vectors with positive components are common in many applied fields, for example, in meteorology, when daily precipitation is measured through a region Marchenko and Genton (2010 Marchenko, Y., Genton, M. (2010). Multivariate log-skew-elliptical distributions with applications to precipitation data. Environmetrics 21:318–340.[Crossref], [Web of Science ®] , [Google Scholar]). Frequently, the log-normal multivariate distribution is used for modeling this type of data. This modeling approach is not appropriate for data with high asymmetry or kurtosis. Consequently, more flexible multivariate distributions than the log-normal multivariate are required. As an alternative to this distribution, we propose the log-alpha-power multivariate and log-skew-normal multivariate models. The first model is an extension for positive data of the fractional order statistics model Durrans (1992 Durrans, S. (1992). Distributions of fractional order statistics in hydrology. Water Resour. Res. 28:1649–1655.[Crossref], [Web of Science ®] , [Google Scholar]). The second one is an extension of the log-skew-normal model studied by Mateu-Figueras and Pawlowsky-Glahn (2007 Mateu-Figueras, G., Pawlowsky-Glahn, V. (2007). The skew-normal distribution on the simplex. Commun. Stat.-Theory Methods 36:1787–1802.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). We study parameter estimation for these models by means of pseudo-likelihood and maximum likelihood methods. We illustrate the proposal analyzing a real dataset.     

Modeling of covariance structures of random effects and random errors in linear mixed models

Abstract

In this paper, we discuss how to model the mean and covariancestructures in linear mixed models (LMMs) simultaneously. We propose a data-driven method to modelcovariance structures of the random effects and random errors in the LMMs. Parameter estimation in the mean and covariances is considered by using EM algorithm, and standard errors of the parameter estimates are calculated through Louis’ (1982 Louis, T.A. (1982). Finding observed information using the EM algorithm. J. Royal Stat. Soc. B 44:98–130. [Google Scholar]) information principle. Kenward’s (1987 Kenward, M.G. (1987). A method for comparing profiles of repeated measurements. Appl. Stat. 36:296–308.[Crossref], [Web of Science ®] , [Google Scholar]) cattle data sets are analyzed for illustration,and comparison to the literature work is made through simulation studies. Our numerical analysis confirms the superiority of the proposed method to existing approaches in terms of Akaike information criterion.     

PMSE performance of two different types of preliminary test estimators under a multivariate t error term

Abstract

In this paper, assuming that the error terms follow a multivariate t distribution, we derive the exact formula for the predictive mean squared error (PMSE) of two different types of pretest estimators. It is shown analytically that one of the pretest estimator dominates the SR estimator if a critical value of the pretest is chosen appropriately. Also, we compare the PMSE of the pretest estimators with the MMSE, AMMSE, SR and PSR estimators by numerical evaluations. Our results show that the pretest estimators dominate the OLS estimator for all combinations when the degrees of freedom is not more than 5.     

Estimation on semi-functional linear errors-in-variables models

Abstract

Semi-functional linear regression models are important in practice. In this paper, their estimation is discussed when function-valued and real-valued random variables are all measured with additive error. By means of functional principal component analysis and kernel smoothing techniques, the estimators of the slope function and the non parametric component are obtained. To account for errors in variables, deconvolution is involved in the construction of a new class of kernel estimators. The convergence rates of the estimators of the unknown slope function and non parametric component are established under suitable norm and conditions. Simulation studies are conducted to illustrate the finite sample performance of our method.     

An improved comorbidity summary score for measuring disease burden and predicting mortality with applications to two national cohorts

Abstract

Research involving administrative healthcare data to study patient outcomes requires the investigator to account for the patient’s disease burden in order to reduce the potential for biased results. Here we develop a comorbidity summary score based on variable importance measures derived from several statistical and machine learning methods and show it has superior predictive performance to the Elixhauser and Charlson indices when used to predict 1-year, 5-year, and 10-year mortality. We used two large Veterans Administration cohorts to develop and validate the summary score and compared predictive performance using the area under ROC curve (AUC) and the Brier score.     

A detection algorithm for the first jump time in sample trajectories of jump-diffusions driven by α-stable white noise

Abstract

The purpose of this paper is to develop a detection algorithm for the first jump point in sampling trajectories of jump-diffusions which are described as solutions of stochastic differential equations driven by α-stable white noise. This is done by a multivariate Lagrange interpolation approach. To this end, we utilize computer simulation algorithm in MATLAB to visualize the sampling trajectories of the jump-diffusions for various combinations of parameters arising in the modeling structure of stochastic differential equations.     

A two-stage procedure for selecting the best of k exponential distributions

The problem of selecting the best of k exponential distributions with different guarantee times and the same unknown variance is considered. A two-stage procedure, similar to the one considered by Bechhofer, Dunnett and Sobel (1954), is given. Some specific guidelines for selecting the first-stage sample size are also given.     

First and Second Derivatives in Time Series Classification Using DTW

In our previous work, we developed a new distance function based on a derivative and showed that our algorithm is effective. In contrast to well-known measures from the literature, our approach considers the general shape of a time series rather than standard distance of function (value) comparison. The new distance was used in classification with the nearest neighbor rule. Now we improve on our previous technique by adding the second derivative. In order to provide a comprehensive comparison, we conducted a set of experiments, testing effectiveness on 47 time series datasets from a wide variety of application domains. Our experiments show that this new method provides a significantly more accurate classification on the examined datasets.