Maximum likelihood estimation of bivariate circular hidden Markov models from incomplete data

In this paper, we propose a hidden Markov model for the analysis of the time series of bivariate circular observations, by assuming that the data are sampled from bivariate circular densities, whose parameters are driven by the evolution of a latent Markov chain. The model segments the data by accounting for redundancies due to correlations along time and across variables. A computationally feasible expectation maximization (EM) algorithm is provided for the maximum likelihood estimation of the model from incomplete data, by treating the missing values and the states of the latent chain as two different sources of incomplete information. Importance-sampling methods facilitate the computation of bootstrap standard errors of the estimates. The methodology is illustrated on a bivariate time series of wind and wave directions and compared with popular segmentation models for bivariate circular data, which ignore correlations across variables and/or along time.     

More on progressively Type‐II censored order statistics for multivariate observations

In this paper, we consider some results on distribution theory of multivariate progressively Type‐II censored order statistics. We also establish some characterizations of Freund's bivariate exponential distribution based on the lack of memory property.     

A test for bivariate exponentiality against bifr alternative

A test is proposed for testing bivariate exponentiality against the Bivariate Increasing Failure Rate (BIFR) class of alternatives. The test statistic is a function of U-statistics and hence asymptotically normally distributed and consistent     

Probabilistic tolerance design for a subsystem under burr distribution using taguchi's loss functions

Taguchi (1986) has derived tolerances for subcomponents, subsystems, parts and materials in which the relationship between a higher-level (Y) and a lower-level (X) quality characteristic is assumed to be deterministic and linear, namely, Y=_α+βX, without an error term. Tsai (1990) developed a probabilistic tolerance design for a subsystem in which a bivariate normal distribution between the above two quality characteristics as well as Taguchi's quadratic loss function were considered together to develop a closed form solution of the tolerance design for a subsystem. The Burr family is very rich for fitting sample data, and has positive domain. A bivariate Burr distribution can describe a nonlinear relationship between two quality characteristics, hence, it is adopted instead of a bivariate normal distribution and the simple solutions of three probabilistic tolerance desings for a subsystem are obtained for three cases of “nominal-is-best”, “smaller-is-berrer”, and “larger-is-beter” quality characteristics, by using Taguchi’ los functions, respectively.     

A test for the composite hypothesis that a population has a gamma distribution

A test of the composite hypothesis that a population has a gamma distribution is presented. The test is conducted by using a rank test of bivariate independence, such as the one .based on Kendallfs sample tau coefficient. The performance of the test is examined by means of a Monte Carlo study.     

Estimating the precision of a bivariate population

Quick efficient estimates are proposed for estimating the standard deviation of a circular bivariate population. Two procedures based on extreme observations are considered. The first of these employs the 100 p percent largest observations, while the second utilizes the extreme observations in k radial sectors.     

On the behaviour of some transforms of the sample correlation coefficient in samples from the bivariate t and the bivariate X2distribution

The present paper studies the normality of five transformations suggested in the literature to normalize the sample correlation coefficient. The parent populations are the bivariate t and the bivariate X ²The results in the previous work of Subrahmaniam and Gajjar are exploited to assess their performance. The density estimation procedure of Tarter and Kronmal is used to provide empiric support to the asymptotic results     

Evaluation and prediction of polygon approximations of planar contours for shape analysis

Contours may be viewed as the 2D outline of the image of an object. This type of data arises in medical imaging as well as in computer vision and can be modeled as data on a manifold and can be studied using statistical shape analysis. Practically speaking, each observed contour, while theoretically infinite dimensional, must be discretized for computations. As such, the coordinates for each contour as obtained at k sampling times, resulting in the contour being represented as a k-dimensional complex vector. While choosing large values of k will result in closer approximations to the original contour, this will also result in higher computational costs in the subsequent analysis. The goal of this study is to determine reasonable values for k so as to keep the computational cost low while maintaining accuracy. To do this, we consider two methods for selecting sample points and determine lower bounds for k for obtaining a desired level of approximation error using two different criteria. Because this process is computationally inefficient to perform on a large scale, we then develop models for predicting the lower bounds for k based on simple characteristics of the contours.     

Adaptive robust estimation in joint mean–covariance regression model for bivariate longitudinal data

The estimation of the covariance matrix is important in the analysis of bivariate longitudinal data. A good estimator for the covariance matrix can improve the efficiency of the estimators of the mean regression coefficients. Furthermore, the covariance estimation itself is also of interest, but it is a challenging job to model the covariance matrix of bivariate longitudinal data due to the complex structure and positive definite constraint. In addition, most of existing approaches are based on the maximum likelihood, which is very sensitive to outliers or heavy-tail error distributions. In this article, an adaptive robust estimation method is proposed for bivariate longitudinal data. Unlike the existing likelihood-based methods, the proposed method can adapt to different error distributions. Specifically, at first, we utilize the modified Cholesky block decomposition to parameterize the covariance matrices. Secondly, we apply the bounded Huber's score function to develop a set of robust generalized estimating equations to estimate the parameters both in the mean and the covariance models simultaneously. A data-driven approach is presented to select the parameter c in the Huber's score function, which can ensure that the proposed method is robust and efficient. A simulation study and a real data analysis are conducted to illustrate the robustness and efficiency of the proposed approach.     

Estimation of monotone bivariate quantile residual life

The bivariate quantile residual life function can play an important role in statistical reliability and survival analysis. In many situations assuming a decreasing form for it is recommended. Here, we propose a new non-parametric estimator of this measure under such restriction. It has been shown that the new estimator is consistent and, with proper normalization, weakly converges to a bivariate Gaussian process. A simulation study shows that the proposed estimator is an alternative to the unrestricted estimator when the bivariate quantile residual life is decreasing. Finally, the new estimators are applied to two real data sets.