Nonparametric Two-Sample Tests of Longitudinal Data with Termination Events

In longitudinal data, observations of response variables occur at fixed or random time points, and can be stopped by a termination event. When comparing longitudinal data for two groups, such irregular observation behavior must be considered to yield suitable results. In this article, we propose the use of nonparametric tests based on the difference between weighted cumulative mean functions for comparing two mean functions with an adjustment for difference in the timing of termination events. We also derive the asymptotic null distributions of the test statistics and examine their small sample properties through simulations. We apply our method to data from a study of liver cirrhosis.     

On conditional correlations defined by a class of nonlinearly truncated distributions and their applications

This article derives and studies several types of conditional correlations. The correlations are obtained by a class of two-piece scale mixture skew-normal distributions. The class is obtained by applying a set of nonlinear constraints to the bivariate scale mixture of normal distributions. The correlations of the class are invariant with respect to the choice of the scale mixing function, however, they are dependent upon the type of the nonlinear truncation. Moreover, their respective upper and lower limits are no longer 1.00 and?1.00. They are useful for the truncated data analysis, the multivariate interdependence methods (such as the principal component analysis and the factor analysis), and the random truncation modelling. Some distributional properties and the Bayesian computation of the correlations are considered when developing necessary theories and providing illustrative examples, respectively. Two applications are also given to demonstrate the usefulness of the conditional correlations in a multivariate analysis.     

Logrank test for bivariate survival data

The logrank test procedure for testing bivariate symmetry against asymmetry in matched-pair data is proposed. The presented test statistic is based on Mantel-Haenszel type statistics evaluated at diagonal grid points on the plane obtained from distinct uncensored failure times. The asymptotic results of the proposed test are derived and an example is shown to illustrate the methodology.     

Simulation-Extrapolation via the Bezier Curve in Measurement Error Models

Simulation-extrapolation (SIMEX) is a method for correcting for bias in measurement error models, and parametric SIMEX estimates are often used. In this paper, we propose a nonparametric method for computing the SIMEX estimate via the Bezier curve, which is a popular smoothing technique in the computer graphics area. Comparisons are done for the bias of the limit values of parametric SIMEX estimates and the Bezier estimate in the various nonlinear measurement error models.     

Analysis of Panel Count Data with Dependent Observation Times

In this article, a semiparametric approach is proposed for the regression analysis of panel count data. Panel count data commonly arise in clinical trials and demographical studies where the response variable is the number of multiple recurrences of the event of interest and observation times are not fixed, varying from subject to subject. It is assumed that two processes exist in this data: the first is for a recurrent event and the second is for observation time. Many studies have been done to estimate mean function and regression parameters under the independency between recurrent event process and observation time process. In this article, the same statistical inference is studied, but the situation where these two processes may be related is also considered. The mixed Poisson process is applied for the recurrent event processes, and a frailty intensity function for the observation time is also used, respectively. Simulation studies are conducted to study the performance of the suggested methods. The bladder tumor data are applied to compare previous studie' results.     

Hierarchical Bayesian Approach to a Multi-Site Hidden Markov Model

Multivariate data with a sequential or temporal structure occur in various fields of study. The hidden Markov model (HMM) provides an attractive framework for modeling long-term persistence in areas of pattern recognition through the extension of independent and identically distributed mixture models. Unlike in typical mixture models, the heterogeneity of data is represented by hidden Markov states. This article extends the HMM to a multi-site or multivariate case by taking a hierarchical Bayesian approach. This extension has many advantages over a single-site HMM. For example, it can provide more information for identifying the structure of the HMM than a single-site analysis. We evaluate the proposed approach by exploiting a spatial correlation that depends on the distance between sites.     

Simulation-Based Inspection Policies for a One-Shot System in Storage Over a Finite Time Span

In this paper, we consider an inspection policy problem for a one-shot system with two types of units over a finite time span and want to determine inspection intervals optimally with given replacement points of Type 2 units. The interval availability and life cycle cost are used as optimization criteria and estimated by simulation. Two optimization models are proposed to find the optimal inspection intervals for the exponential and general distributions. A heuristic method and a genetic algorithm are proposed to find the near-optimal inspection intervals, to satisfy the target interval availability and minimize the life-cycle cost. We study numerical examples to compare the heuristic method with the genetic algorithm and investigate the effect of model parameters to the optimal solutions.     

Influential subsets on the variable selection

When one or few observations are deleted in the multiple linear regression model, they can affect the variable selection. In this paper we derived the formula for the Mallows C_p criterion when k observations are deleted and express it as a functionn of basic building blocks such as residuals and leverages. Also, two real date sets are used to see how the selected model changes as few observations re deleted.     

Asymptotic test of mixture model and its applications to QTL interval mapping

Quantitative trait loci (QTL) mapping has been a standard means in identifying genetic regions harboring potential genes underlying complex traits. Likelihood ratio test (LRT) has been commonly applied to assess the significance of a genetic locus in a mixture model content. Given the time constraint in commonly used permutation tests to assess the significance of LRT in QTL mapping, we study the behavior of the LRT statistic in mixture model when the proportions of the distributions are unknown. We found that the asymptotic null distribution is stationary Gaussian process after suitable transformation. The result can be applied to one-parameter exponential family mixture model. Under certain condition, such as in a backcross mapping model, the tail probability of the supremum of the process is calculated and the threshold values can be determined by solving the distribution function. Simulation studies were performed to evaluate the asymptotic results.     

On the problem of selecting good populations

The problem of selecting good populations out of k normal populations is considered in a Bayesian framework under exchangeable normal priors and additive loss functions. Some basic approximations to the Bayes rules are discussed. These approximations suggest that some well-known classical rules are "approximate" Bayes rules. Especially, it is shown that Gupta-type rules are extended Bayes with respect to a family of the exchangeable normal priors for any bounded and additive loss function. Furthermore, for a simple loss function, the results of a Monte Carlo comparison of Gupta-type rules and Seal-type rules are presented. They indicate that, in general, Gupta-type rules perform better than Seal-type rules