Stochastic compartmental mooeld with branching particles and disasters: sojourn time and related characteristics

A semi-Markov cornpartmental model with branching particies is considered. The notion of disaster is incorporated into the structure. The means of (i) the total sojourn time, (ii) the number of deaths, (iii)the number of births and (iv)the number of emigrant particles in the system are analysed. Some interesting relations connecting these means are established. A few special cases are discussed in detail.     

Multivariate analysis of variance with additional data from an unknown subset of the populations

In this paper we consider the problem of testing the means of k multivariate normal populations with additional data from an unknown subset of the k populations. The purpose of this research is to offer test procedures utilizing all the available data for the multivariate analysis of variance problem because the additional data may contain valuable information about the parameters of the k populations. The standard procedure uses only the data from identified populations. We provide a test using all available data based upon Hotelling' s generalized T²statistic. The power of this test is computed using Betz's approximation of Hotelling' s generalized T²statistic by an F-distribution. A comparison of the power of the test and the standard test procedure is also given.     

Reference regions for beat‐to‐beat ECG data

The QT interval is regarded as an important biomarker for the assessment of arrhythmia liability, and evidence of QT prolongation has led to the withdrawal and relabeling of numerous compounds. Traditional methods of assessing QT prolongation correct the QT interval for the length of the RR interval (which varies inversely with heart‐rate) in a variety of ways. These methods often disagree with each other and do not take into account changes in autonomic state. Correcting the QT interval for RR reduces a bivariate observation (RR, QT) to a univariate observation (QTc). The development of automatic electrocardiogram (ECG) signal acquisition systems has made it possible to collect continuous (so called ‘beat‐to‐beat’) ECG data. ECG data collected prior to administration of a compound allow us to define a region for (RR, QT) values that encompasses typical activity. Such reference regions are used in clinical applications to define the ‘normal’ region of clinical or laboratory measurements. This paper motivates the need for reference regions of (RR, QT) values from beat‐to‐beat ECG data, and describes a way of constructing these. We introduce a measure of agreement between two reference regions that points to the reliability of 12‐lead digital Holter data. We discuss the use of reference regions in establishing baselines for ECG parameters to assist in the evaluation of cardiac risk and illustrate using data from two methodological studies. Copyright © 2010 John Wiley & Sons, Ltd.     

Multivariate discrete exponential family of distributions and their properties

The paper generalizes the univariate discrete exponential family of distributions to the multivariate situation, and this generalization includes the multivariate power series distributions, the multivariate Lagrangian distributions, and the modified multivariate power-series distributions. This provides a unified approach for the study of these three classes of distributions. We obtain recurrence relations for moments and cumulants, and the maximum likelihood estimation for the discrete exponential family. These results are applied to some multivariate discrete distributions like the Lagrangian Poisson, Lagrangian (negative) multinomial, logarithmic series distributions and multivariate Lagrangian negative binomial distribution.     

Multivariate two‐sample permutation tests for trials with multiple time‐to‐event outcomes

Clinical trials involving multiple time‐to‐event outcomes are increasingly common. In this paper, permutation tests for testing for group differences in multivariate time‐to‐event data are proposed. Unlike other two‐sample tests for multivariate survival data, the proposed tests attain the nominal type I error rate. A simulation study shows that the proposed tests outperform their competitors when the degree of censored observations is sufficiently high. When the degree of censoring is low, it is seen that naive tests such as Hotelling's T² outperform tests tailored to survival data. Computational and practical aspects of the proposed tests are discussed, and their use is illustrated by analyses of three publicly available datasets. Implementations of the proposed tests are available in an accompanying R package.     

New control charts for simultaneous monitoring of the mean vector and covariance matrix of multivariate multiple linear profiles

Simultaneous monitoring of the mean vector and covariance matrix in multivariate processes allows practitioners to avoid the inflated false alarm rate that results from using two independent control charts. In this paper, we extend exponentially weighted moving average semicircle and generally weighted moving average semicircle control charts to monitor the mean vector and covariance matrix of multivariate multiple linear regression profiles in Phase II simultaneously. These new control charts are compared with the existing control charts in the literature in terms of the average run length criterion. Finally, a case is considered to show the application of the proposed charts.     

On the property of multivariate generalized hyperbolic distribution and the Stein-type inequality

This note consists of two parts . In the first part, we provide a pedagogic review on the multivariate generalized hyperbolic (MGH) distribution. We show that this probability family is close under margining, conditioning, and linear transforms; however, such property does not hold for its subclasses. In the second part, we obtain the Stein-type inequality in the context of MGH distribution. Moreover, we apply the Stein-type inequality to prove a lower bound for Var[h(X)]. Particularly, we present examples when X belongs to some well-known subclasses in MGH family.     

Bootstrap maximum likelihood for quasi-stationary distributions

Quasi-stationary distributions have many applications in diverse research fields. We develop a bootstrap-based maximum likelihood (BML) method to deal with quasi-stationary distributions in statistical inference. To efficiently implement a bootstrap procedure that can handle the dependence among observations and speed up the computation, a novel block bootstrap algorithm is proposed to accommodate parallel bootstrap. In particular, we select a suitable block length for use with the parallel bootstrap. The estimation error is investigated to show its convergence. The proposed BML is shown to be asymptotically unbiased. Some numerical studies are given to examine the performance of the new algorithm. The advantages are evidenced through a comparison with some competitors and some examples are analysed for illustration.     

Response dimension reduction: model-based approach

In this paper, a model-based approach to reduce the dimension of response variables in multivariate regression is newly proposed, following the existing context of the response dimension reduction developed by Yoo and Cook [Response dimension reduction for the conditional mean in multivariate regression. Comput Statist Data Anal. 2008;53:334–343]. The related dimension reduction subspace is estimated by maximum likelihood, assuming an additive error. In the new approach, the linearity condition, which is assumed for the methodological development in Yoo and Cook (2008), is understood through the covariance matrix of the random error. Numerical studies show potential advantages of the proposed approach over Yoo and Cook (2008). A real data example is presented for illustration.     

A study of geometric and statistical curvatures for the skew-normal family of distributions

With special reference to the family of skew-normal distributions, we consider geometric curvature of a probability density function as a means to define and identify rare or catastrophic events—a phenomenon common in studying the financial instruments. Further, we study the statistical curvature properties of this family of distributions and discuss the sample size issue, to assess, to what extent the linear and likelihood-based inference of exponential family of distribution can be applicable for the skew-normal family.