An R Package for Value at Risk and Expected Shortfall

Value at risk and expected shortfall are the two most popular measures of financial risk. But the available R packages for their computation are limited. Here, we introduce an R contributed package written by the authors. It computes the two measures for over 100 parametric distributions, including all commonly known distributions. We expect that the R package could be useful to researchers and to the financial community.     

Nonparametric Estimation of the Probability of Detection of Flaws in an Industrial Component,from Destructive and Nondestructive Testing Data,Using Approximate Bayesian Computation

We consider the problem of estimating the probability of detection (POD) of flaws in an industrial steel component. Modeled as an increasing function of the flaw height, the POD characterizes the detection process; it is also involved in the estimation of the flaw size distribution, a key input parameter of physical models describing the behavior of the steel component when submitted to extreme thermodynamic loads. Such models are used to assess the resistance of highly reliable systems whose failures are seldom observed in practice. We develop a Bayesian method to estimate the flaw size distribution and the POD function, using flaw height measures from periodic in‐service inspections conducted with an ultrasonic detection device, together with measures from destructive lab experiments. Our approach, based on approximate Bayesian computation (ABC) techniques, is applied to a real data set and compared to maximum likelihood estimation (MLE) and a more classical approach based on Markov Chain Monte Carlo (MCMC) techniques. In particular, we show that the parametric model describing the POD as the cumulative distribution function (cdf) of a log‐normal distribution, though often used in this context, can be invalidated by the data at hand. We propose an alternative nonparametric model, which assumes no predefined shape, and extend the ABC framework to this setting. Experimental results demonstrate the ability of this method to provide a flexible estimation of the POD function and describe its uncertainty accurately.     

Écrire la loi. Un travail de bureau pour hauts fonctionnaires du ministère de l’Écologie

This article analyses the process of the drafting of the Bachelot Bill on technological risks, passed on 30 July 2003 but initiated by Lionel Jospin's government immediately following the AZF factory disaster in September 2001. It focuses on the practical work of the civil servants at the Ministry of Ecology responsible for setting the orientations of the reform, transcribing them into legal provisions and ensuring that they would be supported and passed in Parliament. The analysis approaches the task by tracking both the successive versions of the bill and the multiple documents used in the preparation of these versions. This immersion in the administrative production of law provides a better understanding of the relations between the civil service and the political authorities, and shows how the boundary between these two worlds of practice is difficult to draw in the legal drafting process.     

SELECTION OF WEIGHTS FOR WEIGHTED MODEL AVERAGING

We address the task of choosing prior weights for models that are to be used for weighted model averaging. Models that are very similar should usually be given smaller weights than models that are quite distinct. Otherwise, the importance of a model in the weighted average could be increased by augmenting the set of models with duplicates of the model or virtual duplicates of it. Similarly, the importance of a particular model feature (a certain covariate, say) could be exaggerated by including many models with that feature. Ways of forming a correlation matrix that reflects the similarity between models are suggested. Then, weighting schemes are proposed that assign prior weights to models on the basis of this matrix. The weighting schemes give smaller weights to models that are more highly correlated. Other desirable properties of a weighting scheme are identified, and we examine the extent to which these properties are held by the proposed methods. The weighting schemes are applied to real data, and prior weights, posterior weights and Bayesian model averages are determined. For these data, empirical Bayes methods were used to form the correlation matrices that yield the prior weights. Predictive variances are examined, as empirical Bayes methods can result in unrealistically small variances.     

Deux définitions de la collégialité et leur articulation: le cas d’un diocèse catholique

Against the background of the conflict observed between managers and professionals, two definitions of collegiality are identified: on the one hand, as a specific organizational form (bottom up) and, on the other, as a procedure of bureaucratic management (top down). A study of networks of priests in a catholic diocese in France is used to explore how these two definitions are related. Questions are raised about the effects of a too narrow organizational rationalization that uses collegiality only as a top down, bureaucratic managerial procedure. This always entails the risk of making the work done by experts sterile because it overlooks the first type of collegiality, which is based on the nature of the nonroutine tasks that members perform jointly thanks to an endogenous organizational structure of a bottom up type.     

Some Epidemiological Characteristics of Perpetrators and Victims of Incest in Contemporary Ghana: Analysis of Media Reports

ABSTRACT

In Ghana, incest is considered sinful, taboo, and illegal. However, recent media reports show that incest has become a daily reality in Ghana. This study is a situational analysis of the pattern of incest in Ghana as reported in the media from January 2008 through July 2015. Qualitative content analysis was conducted on 48 incest news reports in Ghana. The findings showed that father–daughter incest was most frequent across the study period. Forty-seven females aged 3 to 25 years and a male aged 3 years were identified as victims. Generally, the incest lasted between 1 day and 13 years before disclosure. Perpetrators employed psychological and/or physical methods to coerce their victims. Marital difficulties, diabolical control, and seduction by victim featured prominently as alleged motives behind the abuse. The study observes that the recent increase in father–daughter incest warrants an immediate shift of research attention onto men’s mental health in Ghana.     

On the effect of the number of quadrature points in a logistic random effects model: an example

Although generalized linear mixed models are recognized to be of major practical importance, it is also known that they can be computationally demanding. The problem is the evaluation of the integral in calculating the marginalized likelihood. The straightforward method is based on the Gauss–Hermite technique, based on Gaussian quadrature points. Another approach is provided by the class of penalized quasi-likelihood methods. It is commonly believed that the Gauss–Hermite method works relatively well in simple situations but fails in more complicated structures. However, we present here a strikingly simple example of a logistic random-intercepts model in the context of a longitudinal clinical trial where the method gives valid results only for a high number of quadrature points ( Q ). As a consequence, this result warns the practitioner to examine routinely the dependence of the results on Q . The adaptive Gaussian quadrature, as implemented in the new SAS procedure NLMIXED, offered the solution to our problem. However, even the adaptive version of Gaussian quadrature needs careful handling to ensure convergence.     

Screening for prostate cancer by using random-effects models

Summary. Random-effects models are used to screen male participants in a long-term longitudinal study for prostate cancer. By using posterior probabilities, each male can be classified into one of four diagnostic states for prostate disease: normal, benign prostatic hyperplasia, local cancer and metastatic cancer. Repeated measurements of prostate-specific antigen, collected when there was no clinical evidence of prostate disease, are used in the classification process. An individual's screening data are considered one repeated measurement at a time as his data are collected longitudinally over time. Posterior probabilities are calculated on the basis of data from other individuals with confirmed diagnoses of each of the four diagnostic states.     

A Theory of Optimal Inheritance Taxation

This paper derives optimal inheritance tax formulas that capture the key equity‐efficiency trade‐off, are expressed in terms of estimable sufficient statistics, and are robust to the underlying structure of preferences. We consider dynamic stochastic models with general and heterogeneous bequest tastes and labor productivities. We limit ourselves to simple but realistic linear or two‐bracket tax structures to obtain tractable formulas. We show that long‐run optimal inheritance tax rates can always be expressed in terms of aggregate earnings and bequest elasticities with respect to tax rates, distributional parameters, and social preferences for redistribution. Those results carry over with tractable modifications to (a) the case with social discounting (instead of steady‐state welfare maximization), (b) the case with partly accidental bequests, (c) the standard Barro–Becker dynastic model. The optimal tax rate is positive and quantitatively large if the elasticity of bequests to the tax rate is low, bequest concentration is high, and society cares mostly about those receiving little inheritance. We propose a calibration using micro‐data for France and the United States. We find that, for realistic parameters, the optimal inheritance tax rate might be as large as 50%–60%—or even higher for top bequests, in line with historical experience.     

Investigation of mixed model repeated measures analyses and non‐linear random coefficient models in the context of long‐term efficacy data

The longitudinal data from 2 published clinical trials in adult subjects with upper limb spasticity (a randomized placebo‐controlled study [NCT01313299] and its long‐term open‐label extension [NCT01313312]) were combined. Their study designs involved repeat intramuscular injections of abobotulinumtoxinA (Dysport®), and efficacy endpoints were collected accordingly. With the objective of characterizing the pattern of response across cycles, Mixed Model Repeated Measures analyses and Non‐Linear Random Coefficient (NLRC) analyses were performed and their results compared. The Mixed Model Repeated Measures analyses, commonly used in the context of repeated measures with missing dependent data, did not involve any parametric shape for the curve of changes over time. Based on clinical expectations, the NLRC included a negative exponential function of the number of treatment cycles, with its asymptote and rate included as random coefficients in the model. Our analysis focused on 2 specific efficacy parameters reflecting complementary aspects of efficacy in the study population. A simulation study based on a similar study design was also performed to further assess the performance of each method under different patterns of response over time. This highlighted a gain of precision with the NLRC model, and most importantly the need for its assumptions to be verified to avoid potentially biased estimates. These analyses describe a typical situation and the conditions under which non‐linear mixed modeling can provide additional insights on the behavior of efficacy parameters over time. Indeed, the resulting estimates from the negative exponential NLRC can help determine the expected maximal effect and the treatment duration required to reach it.