Variable Selection in Joint Location and Scale Models of the Skew-t-Normal Distribution

Variable selection is an important issue in all regression analysis, and in this article, we investigate the simultaneous variable selection in joint location and scale models of the skew-t-normal distribution when the dataset under consideration involves heavy tail and asymmetric outcomes. We propose a unified penalized likelihood method which can simultaneously select significant variables in the location and scale models. Furthermore, the proposed variable selection method can simultaneously perform parameter estimation and variable selection in the location and scale models. With appropriate selection of the tuning parameters, we establish the consistency and the oracle property of the regularized estimators. These estimators are compared by simulation studies.     

Computation of the central and noncentral f distributions

Recursion relations suitable for rapid computation are derived for the cumulative distribution of F′ = (X/m)/(Y/n) where X is χ²(λ, m) and Y is independently χ²(n). When n is even no complicated function evaluations are needed. For n odd, a special doubly noncentral t distribution is needed to start the computation. Series representations for this t distribution are given with rigorous bounds on truncation errors. Proper recursion techniques for numerical evaluation of the special functions are given.     

Tests of hypotheses based on ranks in the general linear model

A unified approach is developed for testing hypotheses in the general linear model based on the ranks of the residuals. It complements the nonparametric estimation procedures recently reported in the literature. The testing and estimation procedures together provide a robust alternative to least squares. The methods are similar in spirit to least squares so that results are simple to interpret. Hypotheses concerning a subset of specified parameters can be tested, while the remaining parameters are treated as nuisance parameters. Asymptotically, the test statistic is shown to have a chi-square distribution under the null hypothesis. This result is then extended to cover a sequence of contiguous alternatives from which the Pitman efficacy is derived. The general application of the test requires the consistent estimation of a functional of the underlying distribution and one such estimate is furnished.     

On a characterization of the dispersion matrix based on the properties of regression

This paper provides a partial solution to a problem posed by J. Neyman (1965) regarding the characterization of multivariate negative binomial distribution based on the properties of regression. It is shown that some of the properties of regression characterize the form of the nonsingular dispersion matrix of the parent distribution, which, interestingly enough, corresponds to only two types viz. those of positive and negative multivariate binomial distributions.     

The rank distribution of the maximin in a random matrix

Let S be a set of tm distinct real numbers and R a random t × m matrix of these tm numbers with rows {r_i} and columns (c_i}. Define b = Max Min x. l≤i≤t x?r_i. Let c be the event Max Min x = Min Max x. l≤i≤t x?r_i l≤i≤m x?c_i. This paper derives the probability distribution of the rank of b in S, as well as the same distribution conditional on c.     

A Case Study in Nonresponse: Plaintiff vs. California State Board of Equalization

In any sample survey, nonresponse bias is a potential issue. Even with a moderately high nonresponse rate, however, covariates can sometimes be used to show that the nonresponse bias is likely to be small. This note presents such an argument, which was used by the winning side in a tax case.     

An EM algorithm for the estimation of parameters of bivariate generalized exponential distribution under random left censoring

In this paper, we consider the four-parameter bivariate generalized exponential distribution proposed by Kundu and Gupta [Bivariate generalized exponential distribution, J. Multivariate Anal. 100 (2009), pp. 581–593] and propose an expectation–maximization algorithm to find the maximum-likelihood estimators of the four parameters under random left censoring. A numerical experiment is carried out to discuss the properties of the estimators obtained iteratively.     

Construction of a Markov Chain with Given Stationary Distribution

Given a rational, finite probability vector, a Markov chain is constructed having the given vector as its stationary distribution.     

An Improved Boxplot for Univariate Data

The boxplot is an effective data-visualization tool useful in diverse applications and disciplines. Although more sophisticated graphical methods exist, the boxplot remains relevant due to its simplicity, interpretability, and usefulness, even in the age of big data. This article highlights the origins and developments of the boxplot that is now widely viewed as an industry standard as well as its inherent limitations when dealing with data from skewed distributions, particularly when detecting outliers. The proposed Ratio-Skewed boxplot is shown to be practical and suitable for outlier labeling across several parametric distributions.