A multivariate solution of the multivariate ranking and selection problem

The problem of selection of the best multivariate population is given a new formulation which does not involve reducing the populations to univariate quantities. This formulation's solution is developed for known, and (using the Heteroscedastic Method) also for unknown, variance-covariance matrices. Preference reversals and arbitrary nonlinear preference functions are explicitly allowed in this new theory     

Nonparametric prediction analysis for binary data

A nonparametric inference algorithm developed by Davis and Geman (1983) is extended problem. The algorithm and applied to a medical prediction employs an estimation procedure for acquiring pairwise statistics among variables of a binary data set, allows for the data-driven creation of interaction terms among the variables, and employs a decision rule which asymptotically gives the minimum expected error. The inference procedure was designed for large data sets but has been extended via the method of cross-validation to encompass smaller data sets.     

Reliability models for categorical data

As assumed hypothetical consensus category corresponding to a case being classified provides a basis for assessment of reliability of judges. Equivalent judges are characterised by the joint probability distribution of the judge assignment and the consensus category. Estimates of the conditional probabilities of judge assignment given consensus category and of consensus category given judge assignments are indices of reliability. All parameters can be estimated if data include classifications of a number of cases by 3 or more judges. Restrictive assumptions are imposed to obtain models for data from classifications by two judges. Maximum likelihood estimation is discussed and illustrated by example for the 3 or more judges case.     

More flexible sequential and non-sequential designs in long-term clinical trial

In this paper, we describe decision making procedures as they exist in most clinical trials,review some recently suggested approaches to monitoring and clarify how these methods allow greater flexibility in monitoring and explicit specification of data monitoring methods in the protocol.     

Clustering Discrete Data Through the Multinomial Mixture Model

In this work, the multinomial mixture model is studied, through a maximum likelihood approach. The convergence of the maximum likelihood estimator to a set with characteristics of interest is shown. A method to select the number of mixture components is developed based on the form of the maximum likelihood estimator. A simulation study is then carried out to verify its behavior. Finally, two applications on real data of multinomial mixtures are presented.     

Bayes Analysis and Comparison of Accelerated Weibull and Accelerated Birnbaum–Saunders Models

Several models are proposed in the literature for modeling fatigue data resulting from materials subject to cyclic stress and strain. Accelerated Weibull and accelerated Birnbaum–Saunders distributions are most commonly used models. Whereas the accelerated Weibull model is easier compared to accelerated Birnbaum–Saunders, it fails to represent the situation equally well. The present article focuses on Bayes analysis of the two models and provides a comparison based on some important Bayesian tools. Model compatibility study using predictive simulation ideas is preceded by the said comparison. Throughout, the posterior simulations are carried out by Markov chain Monte Carlo procedure.     

On the Computation and Finite Sample Behavior of the Constrained M-Estimates for Multivariate Location and Scatter

Abstract

Constrained M (CM) estimates of multivariate location and scatter [Kent, J. T., Tyler, D. E. (1996). Constrained M-estimation for multivariate location and scatter. Ann. Statist. 24:1346–1370] are defined as the global minimum of an objective function subject to a constraint. These estimates combine the good global robustness properties of the S estimates and the good local robustness properties of the redescending M estimates. The CM estimates are not explicitly defined. Numerical methods have to be used to compute the CM estimates. In this paper, we give an algorithm to compute the CM estimates. Using the algorithm, we give a small simulation study to demonstrate the capability of the algorithm finding the CM estimates, and also to explore the finite sample behavior of the CM estimates. We also use the CM estimators to estimate the location and scatter parameters of some multivariate data sets to see the performance of the CM estimates dealing with the real data sets that may contain outliers.     

High-order extensions of the Double Chain Markov Model

The Double Chain Markov Model is a fully Markovian model for the representation of time-series in random environments. In this article, we show that it can handle transitions of high-order between both a set of observations and a set of hidden states. In order to reduce the number of parameters, each transition matrix can be replaced by a Mixture Transition Distribution model. We provide a complete derivation of the algorithms needed to compute the model. Three applications, the analysis of a sequence of DNA, the song of the wood pewee, and the behavior of young monkeys show that this model is of great interest for the representation of data that can be decomposed into a finite set of patterns.     

Correcting for nonlinear measurement errors in the dependent variable in the general linear model

Linear models are considered in which measurement error is present in the dependent variable. Observed values are related to true values via nonlinear regression models with the parameters in the measurement error models being estimated with the use of independent, external data, collected using standards. Pseudo-maximum likelihood estimators and their asymptotic properties are developed under normality assumptions and the common approach of simply analyzing imputed values obtained from the nestimated calibration curves is assessed. A small simulation evaluates the procedures. An example is presented in which urinary neopterin (measured via radioimmunoassay) is nbeing compared between two groups of individuals.     

Kullback-leibler information approach to the optimum measurement point for bayesian estimation

When an appropriate parametric model and a prior distribution of its parameters are given to describe clinical time courses of a dynamic biological process, Bayesian approaches allow us to estimate the entire profiles from a few or even a single observation per subject. The goodness of the estimation depends on the measurement points at which the observations were made. The number of measurement points per subject is generally limited to one or two. The limited measurement points have to be selected carefully. This paper proposes an approach to the selection of the optimum measurement point for Bayesian estimations of clinical time courses. The selection is made among given candidates, based on the goodness of estimation evaluated by the Kullback-Leibler information. This information measures the discrepancy of an estimated time course from the true one specified by a given appropriate model. The proposed approach is applied to a pharmacokinetic analysis, which is a typical clinical example where the selection is required. The results of the present study strongly suggest that the proposed approach is applicable to pharmacokinetic data and has a wide range of clinical applications.