Robust gene–environment interaction analysis using penalized trimmed regression

ABSTRACT

In biomedical and epidemiological studies, gene–environment (G–E) interactions have been shown to importantly contribute to the etiology and progression of many complex diseases. Most existing approaches for identifying G–E interactions are limited by the lack of robustness against outliers/contaminations in response and predictor spaces. In this study, we develop a novel robust G–E identification approach using the trimmed regression technique under joint modelling. A robust data-driven criterion and stability selection are adopted to determine the trimmed subset which is free from both vertical outliers and leverage points. An effective penalization approach is developed to identify important G–E interactions, respecting the ‘main effects, interactions’ hierarchical structure. Extensive simulations demonstrate the better performance of the proposed approach compared to multiple alternatives. Interesting findings with superior prediction accuracy and stability are observed in the analysis of The Cancer Genome Atlas data on cutaneous melanoma and breast invasive carcinoma.     

Approximating M/G/1 Waiting Time Tail Probabilities

Abstract

We propose a new approximation formula for the waiting time tail probability of the M/G/1 queue with FIFO discipline and unlimited waiting space. The aim is to address the difficulty of obtaining good estimates when the tail probability has non-exponential asymptotics. We show that the waiting time tail probability can be expressed in terms of the waiting time tail probability of a notional M/G/1 queue with truncated service time distribution plus the tail probability of an extreme order statistic. The Cramér–Lundberg approximation is applied to approximate the tail probability of the notional queue. In essence, our technique extends the applicability of the Cramér–Lundberg approximation to cases where the standard Lundberg condition does not hold. We propose a simple moment-based technique for estimating the parameters of the approximation; numerical results demonstrate that our approximation can yield very good estimates over the whole range of the argument.     

Bivariate exponential-type distributions with linear regression

Let Xbe a bivariate exponential-type random vector (BIDLIKAR, PATIL (1968)), than it is proved:

1. If P(X ≥0) = 1 is valid, then Xhas linear regression to both directions if and only if Xpossesses a symmetric Γ-distribution.

2. Xpossesses linear regression to both directions with constant regression coefficients (independent of the parameter vector ? of the exponential-type distribution (BIDLIKAR, PATIL (1968)) if and only if Xis normal distributed.     

Multilevel regression mixture analysis

Summary.  A two-level regression mixture model is discussed and contrasted with the conventional two-level regression model. Simulated and real data shed light on the modelling alternatives. The real data analyses investigate gender differences in mathematics achievement from the US National Education Longitudinal Survey. The two-level regression mixture analyses show that unobserved heterogeneity should not be presupposed to exist only at level 2 at the expense of level 1. Both the simulated and the real data analyses show that level 1 heterogeneity in the form of latent classes can be mistaken for level 2 heterogeneity in the form of the random effects that are used in conventional two-level regression analysis. Because of this, mixture models have an important role to play in multilevel regression analyses. Mixture models allow heterogeneity to be investigated more fully, more correctly attributing different portions of the heterogeneity to the different levels.     

Estimation of parameters for a Birnbaum–Saunders regression model with censored data

Little work has been published on the analysis of censored data for the Birnbaum–Saunders distribution (BISA). In this article, we implement the EM algorithm to fit a regression model with censored data when the failure times follow the BISA. Three approaches to implement the E-Step of the EM algorithm are considered. In two of these implementations, the M-Step is attained by an iterative least-squares procedure. The algorithm is exemplified with a single explanatory variable in the model.     

Four tests of fit for the beta-binomial distribution

Tests based on the Anderson–Darling statistic, a third moment statistic and the classical Pearson–Fisher X ² statistic, along with its third-order component, are considered. A small critical value and power study are given. Some examples illustrate important applications.     

Bayesian Binomial Regression: Predicting Survival at a Trauma Center

Standard methods for analyzing binomial regression data rely on asymptotic inferences. Bayesian methods can be performed using simple computations, and they apply for any sample size. We provide a relatively complete discussion of Bayesian inferences for binomial regression with emphasis on inferences for the probability of “success.” Furthermore, we illustrate diagnostic tools, perform model selection among nonnested models, and examine the sensitivity of the Bayesian methods.     

Formulating the Matthew Principle: on the role of the middle classes in the welfare state

The Matthew Principle, which formulates the role of the middle classes in the welfare state, is discussed in this article. The middle classes are described as primary beneficiaries of the welfare state. This status is achieved through the ability of this social group to influence the policy-formulating process by way of six different channels of influence. It is also facilitated due to the impact of middle-class clients upon the implementation stage of social policy. A case study of employment policy in Israel illustrates the workings of different aspects of the Matthew Principle.     

Small sample efficiencies of tests based on the method of n rankings

A method is proposed for calculating the small sample powers of rank tests which are based on the method of n rankings. A class of normal shift alternative hypotheses is considered, and Hodges–Lehmann efficiencies are calculated for the Friedman test.     

Tail dependence for skew Laplace distribution and skew Cauchy distribution

ABSTRACT

Coefficient of tail dependence measures the strength of dependence in the tail of a bivariate distribution and it has been found useful in the risk management. In this paper, we derive the upper tail dependence coefficient for a random vector following the skew Laplace distribution and the skew Cauchy distribution, respectively. The result shows that skew Laplace distribution is asymptotically independent in upper tail, however, skew Cauchy distribution has asymptotic upper tail dependence.