Dirichlet Processes in Nonlinear Mixed Effects Models

In this article, we use two efficient approaches to deal with the difficulty in computing the intractable integrals when implementing Gibbs sampling in the nonlinear mixed effects model (NLMM) based on Dirichlet processes (DP). In the first approach, we compute the Laplace's approximation to the integral for its high accuracy, low cost, and ease of implementation. The second approach uses the no-gaps algorithm of MacEachern and Müller (1998 MacEachern , S. , Müller , P. ( 1998 ). Estimating mixtures of Dirichlet process models . Journal of Computational and Graphical Statistics 7 : 223 – 238 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) to perform Gibbs sampling without evaluating the difficult integral. We apply both approaches to real problems and simulations. Results show that both approaches perform well in density estimation and prediction and are superior to the parametric analysis in that they can detect important model features, such as skewness, long tails, and multimodality, whereas the parametric analysis cannot.     

New Prediction Interval and Band in the Nonlinear Regression Model: Application to Predictive Modeling in Foods

This article is concerned with the proposal of a new prediction interval and band for the nonlinear regression model. The construction principle of this interval and band is based on an exact (the meaning of the term “exact” will be given later) confidence region for parameters of the nonlinear regression model. This region, fully described in Vila and Gauchi (2007 Vila , J.-P. , Gauchi , J.-P. ( 2007 ). Optimal designs based on exact confidence regions for parameter estimation of a nonlinear regression model . Journal of Statistical Planning and Inference 137 ( 9 ): 2935 – 2953 .[Crossref], [Web of Science ®] , [Google Scholar]), provides a rigorous justification for the new prediction interval and band that we propose. This new band is then compared to the classical bands (which are asymptotic and thus approximate for small n), and also to the band based on the bootstrap resampling method. The comparison of these bands is undertaken with simulated and real data from predictive modeling in food science.     

Discriminant Analysis for the von Mises-Fisher Distribution

The von Mises-Fisher distribution is widely used for modeling directional data. In this article, we derive the discriminant rules based on this distribution to assign objects into pre-existing classes. We determine a distance between two von Mises-Fisher populations and we calculate estimates of the misclassification probabilities. We also analyze the behavior of the distance between two von Mises-Fisher populations and of the estimates of the misclassification probabilities when we modify the parameters of the populations or the samples size or the dimension of the sphere. Finally, we present an example with real spherical data available in the literature.     

Variable Selection for Partially Linear Models with Randomly Censored Data

This article proposes a variable selection procedure for partially linear models with right-censored data via penalized least squares. We apply the SCAD penalty to select significant variables and estimate unknown parameters simultaneously. The sampling properties for the proposed procedure are investigated. The rate of convergence and the asymptotic normality of the proposed estimators are established. Furthermore, the SCAD-penalized estimators of the nonzero coefficients are shown to have the asymptotic oracle property. In addition, an iterative algorithm is proposed to find the solution of the penalized least squares. Simulation studies are conducted to examine the finite sample performance of the proposed method.     

The computational aspects of multiple attribute sampling procedures

One method of controlling the quality of incoming lots is through attribute sampling. To simultaneously control several (possibly dependent) attributes, properly chosen single attribute sampling plans can be merged into a multiple attribute sampling plan. The general form of such a plan is given and various alternatives are discussed. The multinomial distribution is used to develop formulae necessary for an analysis of a multiple attribute plan. Due to the lengthy nature of the calculations involved, a computer algorithm is outlined.     

Comparison of conditional and unconditional confidence intervals for the double exponential distribution

Conditional confidence intervals for the location parameter of the double exponential distribution based on maximum likelihood estimators conditioned on a set of ancillary statistics and the corresponding unconditional confidence intervals based on the maximum likelihood estimators alone are compared in two ways. Monte Carlo techniques are used and the conditional approach appears to give slightly better results although agreement as n becomes larger is noted     

Goodness-of-Fit Test and Parameter Estimation for a Proportional Odds Model of Random Censorship

Censored data arise naturally in a number of fields, particularly in problems of reliability and survival analysis. There are several types of censoring, in this article, we will confine ourselves to the right randomly censoring type. Recently, Ahmadi et al. (2010 Ahmadi , J. , Doostparast , M. , Parsian , A. ( 2010 ). Bayes estimation based on random censored data for some life time models under symmetric and asymmetric loss functions . Communcations in Statistics-Theory and Methods , 39 : 3058 – 3071 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) considered the problem of estimating unknown parameters in a general framework based on the right randomly censored data. They assumed that the survival function of the censoring time is free of the unknown parameter. This assumption is sometimes inappropriate. In such cases, a proportional odds (PO) model may be more appropriate (Lam and Leung, 2001 Lam , K. F. , Leung , T. L. ( 2001 ). Marginal likelihood estimation for proportional odds models with right censored data . Lifetime Data Analysis 7 : 39 – 54 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]). Under this model, in this article, point and interval estimations for the unknown parameters are obtained. Since it is important to check the adequacy of models upon which inferences are based (Lawless, 2003 Lawless , J. F. (2003). Statistical Models and Methods for Lifetime Data. , 2nd ed. New York : John Wiley & Sons. [Google Scholar], p. 465), two new goodness-of-fit tests for PO model based on right randomly censored data are proposed. The proposed procedures are applied to two real data sets due to Smith (2002 Smith , P. J. ( 2002 ). Analysis of Failure and Survival Data . London : Chapman & Hall, CRC . [Google Scholar]). A Monte Carlo simulation study is conducted to carry out the behavior of the estimators obtained.     

A Comparison of Approaches to Interval Estimation of Mediated Effects in Linear Regression Models

A novel approach based on the concepts of a generalized pivotal quantity (GPQ) is developed to construct confidence intervals for the mediated effect. Thereafter, its performance is compared with six interval estimation approaches in terms of empirical coverage probability and expected length via simulation and two real examples. The results show that the GPQ-based and bootstrap percentile methods outperform other methods when mediated effects exist in small and medium samples. Moreover, the GPQ-based method exhibits a more stable performance in small and non-normal samples. A discussion on how to choose the best interval estimation method for mediated effects is presented.     

Infinite Variation Tempered Stable Ornstein–Uhlenbeck Processes with Discrete Observations

We investigate transition law between consecutive observations of Ornstein–Uhlenbeck processes of infinite variation with tempered stable stationary distribution. Thanks to the Markov autoregressive structure, the transition law can be written in the exact sense as a convolution of three random components; a compound Poisson distribution and two independent tempered stable distributions, one with stability index in (0, 1) and the other with index in (1, 2). We discuss simulation techniques for those three random elements. With the exact transition law and proposed simulation techniques, sample paths simulation proves significantly more efficient, relative to the known approximative technique based on infinite shot noise series representation of tempered stable Lévy processes.