Beyond tail median and conditional tail expectation: Extreme risk estimation using tail Lp-optimization

The conditional tail expectation (CTE) is an indicator of tail behavior that takes into account both the frequency and magnitude of a tail event. However, the asymptotic normality of its empirical estimator requires that the underlying distribution possess a finite variance; this can be a strong restriction in actuarial and financial applications. A valuable alternative is the median shortfall (MS), although it only gives information about the frequency of a tail event. We construct a class of tail L^p-medians encompassing the MS and CTE. For p in (1,2), a tail L^p-median depends on both the frequency and magnitude of tail events, and its empirical estimator is, within the range of the data, asymptotically normal under a condition weaker than a finite variance. We extrapolate this estimator and another technique to extreme levels using the heavy-tailed framework. The estimators are showcased on a simulation study and on real fire insurance data.     

A class of symmetric bivariate uniform distributions

A class of symmetric bivariate uniform distributions is proposed for use in statistical modeling. The distributions may be constructed to be absolutely continuous with correlations as close to±1 as desired. Expressions for the correlations, regressions and copulas are found. An extension to three dimensions is proposed.     

Measuring agreement between two statistics with applications to age standardization

Summary.  We detail a general method for measuring agreement between two statistics. An application is two ratios of directly standardized rates which differ only by the choice of the standard. If the statistics have a high value for the coefficient of agreement then the expected squared difference between the statistics is small relative to the variance of the average of the two statistics, and inferences vary little by changing statistics. The estimation of a coefficient of agreement between two statistics is not straightforward because there is only one pair of observed values, each statistic calculated from the data. We introduce estimators of the coefficient of agreement for two statistics and discuss their use, especially as applied to functions of standardized rates.     

On Optimality of Bayesian Wavelet Estimators

Abstract.  We investigate the asymptotic optimality of several Bayesian wavelet estimators, namely, posterior mean, posterior median and Bayes Factor, where the prior imposed on wavelet coefficients is a mixture of a mass function at zero and a Gaussian density. We show that in terms of the mean squared error, for the properly chosen hyperparameters of the prior, all the three resulting Bayesian wavelet estimators achieve optimal minimax rates within any prescribed Besov space     for p  ≥ 2. For 1 ≤  p  < 2, the Bayes Factor is still optimal for (2 s +2)/(2 s +1) ≤  p  < 2 and always outperforms the posterior mean and the posterior median that can achieve only the best possible rates for linear estimators in this case.     

A Comparison of Partially Adaptive and Reweighted Least Squares Estimation

The small sample performance of least median of squares, reweighted least squares, least squares, least absolute deviations, and three partially adaptive estimators are compared using Monte Carlo simulations. Two data problems are addressed in the paper: (1) data generated from non-normal error distributions and (2) contaminated data. Breakdown plots are used to investigate the sensitivity of partially adaptive estimators to data contamination relative to RLS. One partially adaptive estimator performs especially well when the errors are skewed, while another partially adaptive estimator and RLS perform particularly well when the errors are extremely leptokur-totic. In comparison with RLS, partially adaptive estimators are only moderately effective in resisting data contamination; however, they outperform least squares and least absolute deviation estimators.     

La médiane simpliciale d'Oja: Existence,unicité et stabilité

Oja (1983) examined various ways of measuring location, scatter, skewness, and kurtosis for multivariate distributions. Among other measures of location, he introduced a generalised median known in this paper under the name of the Oja median. In our study of the existence of that median, we show that Oja's definition can only be applied to distributions having a mean. In dimension d θ 2, we establish that the usual method of extension breaks down, which raises the question of the validity of the concept as a notion of median. Two fundamental theoretical properties of that median are also considered: uniqueness and consistency.     

Estimating the variance of the sample median

The small-sample bias and root mean squared error of several distribution-free estimators of the variance of the sample median are examined. A new estimator is proposed that is easy to compute and tends to have the smallest bias and root mean squared error.     

On Pitman's measure of closeness of k-records

Pitman closeness of both the upper and lower k-record statistics to the population quantiles of a location–scale family of distributions is studied. For the population median, the Pitman-closest k-record is also determined. In the case of symmetric distributions, the Pitman closeness probabilities of k-record statistics are shown to be distribution-free, and explicit expressions are also derived for these probabilities. Exact expressions are derived for the required probabilities for uniform and exponential distributions. Numerical results are given for these families and also the Pitman-closest k-record is determined.     

Bootstrap methods for bias correction and confidence interval estimation for nonlinear quantile regression of longitudinal data

This paper examines the use of bootstrapping for bias correction and calculation of confidence intervals (CIs) for a weighted nonlinear quantile regression estimator adjusted to the case of longitudinal data. Different weights and types of CIs are used and compared by computer simulation using a logistic growth function and error terms following an AR(1) model. The results indicate that bias correction reduces the bias of a point estimator but fails for CI calculations. A bootstrap percentile method and a normal approximation method perform well for two weights when used without bias correction. Taking both coverage and lengths of CIs into consideration, a non-bias-corrected percentile method with an unweighted estimator performs best.     

Effect of tuberculosis on immune restoration among HIV-infected patients receiving antiretroviral therapy

In this article, time to immune recovery during antiretroviral therapy was estimated and compared between HIV-infected children with and without tuberculosis (TB). CD4?T-cell restoration was used as a criterion for determining immune recovery. The median residual lifetime function, which is more intuitive and robust compared to the frequently used measures of lifetime data, was used to estimate time to CD4?T-cell restoration. The median residual lifetime is not influenced by extreme observations and heavy-tailed distributions which are commonly encountered in clinical studies. Permutation-based methods were used to compare the CD4?T-cell restoration times between the two groups of patients. Our results indicate that children with TB had uniformly higher median residual lifetimes to immune recovery compared to those without TB. Although TB was associated with slower CD4?T-cell restoration, the differences between the restoration times of the two groups were not statistically significant.