On the Bias of the Maximum Likelihood Estimator for the Two-Parameter Lomax Distribution

The Lomax (Pareto II) distribution has found wide application in a variety of fields. We analyze the second-order bias of the maximum likelihood estimators of its parameters for finite sample sizes, and show that this bias is positive. We derive an analytic bias correction which reduces the percentage bias of these estimators by one or two orders of magnitude, while simultaneously reducing relative mean squared error. Our simulations show that this performance is very similar to that of a parametric bootstrap correction based on a linear bias function. Three examples with actual data illustrate the application of our bias correction.     

Generalized Jackknife-Based Estimators for Univariate Extreme-Value Modeling

In this article, we revisit the importance of the generalized jackknife in the construction of reliable semi-parametric estimates of some parameters of extreme or even rare events. The generalized jackknife statistic is applied to a minimum-variance reduced-bias estimator of a positive extreme value index—a primary parameter in statistics of extremes. A couple of refinements are proposed and a simulation study shows that these are able to achieve a lower mean square error. A real data illustration is also provided.     

A Note on Wald's Type Chi-squared Tests of Fit

The Wald's method for constructing chi-squared tests of fit has been formulated more accurately. It is shown that Wald's type statistics will follow the central chi-squared distribution if and only if the limit covariance matrix of standardized frequencies will not depend on unknown parameters. Several examples that illustrate this important fact are presented. In particular, it is shown that the goodness-of-fit statistic developed by Moore and Stubblebine does not follow the chi-squared limit distribution, and, hence, cannot be used for testing multivariate normality.     

A Compressive Sensing Based Analysis of Anomalies in Generalized Linear Models

In this article, we present a compressive sensing based framework for generalized linear model regression that employs a two-component noise model and convex optimization techniques to simultaneously detect outliers and determine optimally sparse representations of noisy data from arbitrary sets of basis functions. We then extend our model to include model order reduction capabilities that can uncover inherent sparsity in regression coefficients and achieve simple, superior fits. Second, we use the mixed ?₂/?₁ norm to develop another model that can efficiently uncover block-sparsity in regression coefficients. By performing model order reduction over all independent variables and basis functions, our algorithms successfully deemphasize the effect of independent variables that become uncorrelated with dependent variables. This desirable property has various applications in real-time anomaly detection, such as faulty sensor detection and sensor jamming in wireless sensor networks. After developing our framework and inheriting a stable recovery theorem from compressive sensing theory, we present two simulation studies on sparse or block-sparse problems that demonstrate the superior performance of our algorithms with respect to (1) classic outlier-invariant regression techniques like least absolute value and iteratively reweighted least-squares and (2) classic sparse-regularized regression techniques like LASSO.     

The Generalized Pascal Triangle and the Matrix Variate Jensen-Logistic Distribution

This article defines the so called Generalized Matrix Variate Jensen-Logistic distribution. The relevant applications of this class of distributions in Configuration Shape Theory consist of a more efficient computation, supported by the corresponding inference. This demands the solution of two important problems: (1) the development of analytical and efficient formulae for their k-th derivatives and (2) the use of the derivatives to transform the configuration density into a polynomial density under some special matrix Kummer relation, indexed in this case by the Jensen-Logistic kernel. In this article, we solve these problems by deriving a simple formula for the k-th derivative of the density function, avoiding the usual partition theory framework and using a generalization of Pascal triangles. Then we apply the results by obtaining the associated Jensen-Logistic Kummer relations and the configuration polynomial density in the setting of Statistical Shape Theory.     

Bonferroni-type Plug-in Procedure Controlling Generalized Familywise Error Rate

Consider the multiple hypotheses testing problem controlling the generalized familywise error rate k-FWER, the probability of at least k false rejections. We propose a plug-in procedure based on the estimation of the number of true null hypotheses. Under the independence assumption of the p-values corresponding to the true null hypotheses, we first introduce the least favorable configuration (LFC) of k-FWER for Bonferroni-type plug-in procedure, then we construct a plug-in k-FWER-controlled procedure based on LFC. For dependent p-values, we establish the asymptotic k-FWER control under some mild conditions. Simulation studies suggest great improvement over generalized Bonferroni test and generalized Holm test.     

Generalized Analysis-of-variance-type Test for the Single-index Quantile Model

In this paper, we are concerned with a test for the index parameter and index function in the single-index model. Based on the estimates obtained by the quantile regression, we extend the generalized analysis-of-variance-type test to the single-index model. We investigate the asymptotic behavior of the proposed test and demonstrate that its limiting null distribution follows an asymptotically χ²-distribution. The simulation studies and real data applications are conducted to illustrate the finite sample performance of the proposed methods.     

Bayesian,MLE, and GMM Estimation of a Spot Rate Model

ABSTRACT

We develop Markov chain Monte Carlo algorithms for estimating the parameters of the short-term interest rate model. Using Monte Carlo experiments we compare the Bayes estimators with the maximum likelihood and generalized method of moments estimators. We estimate the model using the Japanese overnight call rate data.     

From Ordinary to Shrinkage Square-Root Estimators

Consider a skewed population. Suppose an intelligent guess could be made about an interval that contains the population mean. There may exist biased estimators with smaller mean squared error than the arithmetic mean within such an interval. This article indicates when it is advisable to shrink the arithmetic mean towards a guessed interval using root estimators. The goal is to obtain an estimator that is better near the average of natural origins. An estimator proposed. This estimator contains the Thompson (1968 Thompson , J. R. ( 1968 ). Accuracy borrowing in the estimation of the mean by shrinkage towards an interval . J. Amer. Statist. Assoc. 63 : 953 – 963 . [CSA] [CROSSREF] [Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) ordinary shrinkage estimator, the Jenkins et al. (1973 Jenkins , O. C. , Ringer , L. J. , Hartley , H. O. ( 1973 ). Root estimators . J Amer. Statist. Assoc. 68 : 414 – 419 . [CSA] [CROSSREF] [Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) square-root estimator, and the arithmetic sample mean as special cases. The bias and the mean squared error of the proposed more general estimator is compared with the three special cases. Shrinkage coefficients that yield minimum mean squared error estimators are obtained. The proposed estimator is considerably more efficient than the three special cases. This remains true for highly skewed populations. The merits of the proposed shrinkage square-root estimator are supported by the results of numerical and simulation studies.     

Inference for a Linear Functional of Cumulative Hazard Function Via Empirical Likelihood

This article discusses sampling plans, that is, the allocation of sampling units, for computing tolerance limits in a balanced one--way random-effects model. The expected width of the tolerance interval is derived and used as the basis for comparing different sampling plans. A well-known cost function and examples are used to facilitate the discussion.