Estimation and testing in generalized partial linear models—A comparative study

A particular semiparametric model of interest is the generalized partial linear model (GPLM) which extends the generalized linear model (GLM) by a nonparametric component.The paper reviews different estimation procedures based on kernel methods as well as test procedures on the correct specification of this model (vs. a parametric generalized linear model). Simulations and an application to a data set on East–West German migration illustrate similarities and dissimilarities of the estimators and test statistics.     

A Composite Likelihood Approach to Multivariate Survival Data

This paper is about the statistical analysis of multivariate survival data. We discuss the additive and multiplicative frailty models which have been the most popular models for multivariate survival data. As an alternative to the additive and multiplicative frailty models, we propose basing inference on a composite likelihood function that only requires modelling of the marginal distribution of pairs of failure times. Each marginal distribution of a pair of failure times is here assumed to follow a shared frailty model. The method is illustrated with a real-life example.     

A multivariate von mises distribution with applications to bioinformatics

Motivated by problems of modelling torsional angles in molecules, Singh, Hnizdo & Demchuk (2002) proposed a bivariate circular model which is a natural torus analogue of the bivariate normal distribution and a natural extension of the univariate von Mises distribution to the bivariate case. The authors present here a multivariate extension of the bivariate model of Singh, Hnizdo & Demchuk (2002). They study the conditional distributions and investigate the shapes of marginal distributions for a special case. The methods of moments and pseudo‐likelihood are considered for the estimation of parameters of the new distribution. The authors investigate the efficiency of the pseudo‐likelihood approach in three dimensions. They illustrate their methods with protein data of conformational angles     

Selection of Prior for the Variance Component and Approximations for Posterior Moments in the Fay-Herriot Model

In the Fay-Herriot model, a prior distribution for the variance component allows posterior moments to be approximated with the Laplace method, avoiding computer intensive Monte Carlo Markov chains. The extremely skewed posterior distribution of the variance component results from the asymmetry of the parameter space with variance parameters constrained to be positive. The prior avoids the extreme skewness of the posterior in contrast to the commonly used uniform prior. With this prior, the mean squared error and coverage in the approximate hierarchical Bayes method are satisfactory when used to estimate small area means. Computation time is shorter than with Monte Carlo Markov chains. The approximations give easy interpretations of Bayesian methods and highlight frequentist properties of the parameters.     

Testing treatment‐by‐period interaction in four‐period crossover trials

Statistical analyses of crossover clinical trials have mainly focused on assessing the treatment effect, carryover effect, and period effect. When a treatment‐by‐period interaction is plausible, it is important to test such interaction first before making inferences on differences among individual treatments. Considerably less attention has been paid to the treatment‐by‐period interaction, which has historically been aliased with the carryover effect in two‐period or three‐period designs. In this article, from the data of a newly developed four‐period crossover design, we propose a statistical method to compare the effects of two active drugs with respect to two response variables. We study estimation and hypothesis testing considering the treatment‐by‐period interaction. Constrained least squares is used to estimate the treatment effect, period effect, and treatment‐by‐period interaction. For hypothesis testing, we extend a general multivariate method for analyzing the crossover design with multiple responses. Results from simulation studies have shown that this method performs very well. We also illustrate how to apply our method to the real data problem.     

A new useful three-parameter extension of the exponential distribution

A three-parameter extension of the exponential distribution is introduced and studied in this paper. The new distribution is quite flexible and can be used effectively in modelling survival data, reliability problems, fatigue life studies and hydrological data. It can have constant, decreasing, increasing, upside-down bathtub (unimodal), bathtub-shaped and decreasing–increasing–decreasing hazard rate functions. We provide a comprehensive account of the mathematical properties of the new distribution and various structural quantities are derived. We discuss maximum likelihood estimation of the model parameters for complete sample and for censored sample. An empirical application of the new model to real data is presented for illustrative purposes. We hope that the new distribution will serve as an alternative model to other models available in the literature for modelling real data in many areas.     

On Nonsmooth Estimating Functions via Jackknife Empirical Likelihood

In many applications, the parameters of interest are estimated by solving non‐smooth estimating functions with U‐statistic structure. Because the asymptotic covariances matrix of the estimator generally involves the underlying density function, resampling methods are often used to bypass the difficulty of non‐parametric density estimation. Despite its simplicity, the resultant‐covariance matrix estimator depends on the nature of resampling, and the method can be time‐consuming when the number of replications is large. Furthermore, the inferences are based on the normal approximation that may not be accurate for practical sample sizes. In this paper, we propose a jackknife empirical likelihood‐based inferential procedure for non‐smooth estimating functions. Standard chi‐square distributions are used to calculate the p‐value and to construct confidence intervals. Extensive simulation studies and two real examples are provided to illustrate its practical utilities.     

Bayes estimation for exponential distributions with common location parameter and applications to multi-state reliability models

This paper considers the estimation of the stress–strength reliability of a multi-state component or of a multi-state system where its states depend on the ratio of the strength and stress variables through a kernel function. The article presents a Bayesian approach assuming the stress and strength as exponentially distributed with a common location parameter but different scale parameters. We show that the limits of the Bayes estimators of both location and scale parameters under suitable priors are the maximum likelihood estimators as given by Ghosh and Razmpour [15 M. Ghosh and A. Razmpour, Estimation of the common location parameter of several exponentials, Sankhyā, Ser. A 46 (1984), pp. 383–394. [Google Scholar]]. We use the Bayes estimators to determine the multi-state stress–strength reliability of a system having states between 0 and 1. We derive the uniformly minimum variance unbiased estimators of the reliability function. Interval estimation using the bootstrap method is also considered. Under the squared error loss function and linex loss function, risk comparison of the reliability estimators is carried out using extensive simulations.     

Second‐order Accurate Confidence Regions Based on Members of the Generalized Power Divergence Family

Recently, a technique based on pseudo‐observations has been proposed to tackle the so‐called convex hull problem for the empirical likelihood statistic. The resulting adjusted empirical likelihood also achieves the high‐order precision of the Bartlett correction. Nevertheless, the technique induces an upper bound on the resulting statistic that may lead, in certain circumstances, to worthless confidence regions equal to the whole parameter space. In this paper, we show that suitable pseudo‐observations can be deployed to make each element of the generalized power divergence family Bartlett‐correctable and released from the convex hull problem. Our approach is conceived to achieve this goal by means of two distinct sets of pseudo‐observations with different tasks. An important effect of our formulation is to provide a solution that permits to overcome the problem of the upper bound. The proposal, which effectiveness is confirmed by simulation results, gives back attractiveness to a broad class of statistics that potentially contains good alternatives to the empirical likelihood.