Source inference of exogenous gamma-hydroxybutyric acid (GHB) administered to humans by means of carbon isotopic ratio analysis: novel perspectives regarding forensic investigation and intelligence issues.

γ-Hydroxybutyric acid (GHB) is an endogenous short-chain fatty acid popular as a recreational drug due to sedative and euphoric effects, but also often implicated in drug-facilitated sexual assaults owing to disinhibition and amnesic properties. Whilst discrimination between endogenous and exogenous GHB as required in intoxication cases may be achieved by the determination of the carbon isotope content, such information has not yet been exploited to answer source inference questions of forensic investigation and intelligence interests. However, potential isotopic fractionation effects occurring through the whole metabolism of GHB may be a major concern in this regard. Thus, urine specimens from six healthy male volunteers who ingested prescription GHB sodium salt, marketed as Xyrem(®), were analysed by means of gas chromatography/combustion/isotope ratio mass spectrometry to assess this particular topic. A very narrow range of δ(13)C values, spreading from -24.810/00 to -25.060/00, was observed, whilst mean δ(13)C value of Xyrem(®) corresponded to -24.990/00. Since urine samples and prescription drug could not be distinguished by means of statistical analysis, carbon isotopic effects and subsequent influence on δ(13)C values through GHB metabolism as a whole could be ruled out. Thus, a link between GHB as a raw matrix and found in a biological fluid may be established, bringing relevant information regarding source inference evaluation. Therefore, this study supports a diversified scope of exploitation for stable isotopes characterized in biological matrices from investigations on intoxication cases to drug intelligence programmes.

The wage effect of an indexation clause: Evidence from Spanish manufacturing firms

An endogenous switching model of ex-ante wage changes under indexed and non-indexed settlements is estimated for the Spanish manufacturing sector using collective bargaining firm data for the 1984-1991 period. The likelihood of indexing the settlement is higher for nationwide unions than for other union groups within the works council and increases with the expected level of inflation. For wage change equations, a common structure for indexed and non-indexed settlements is strongly rejected, showing a source of nominal rigidity. For indexed contracts, the expected ex-ante total inflation coverage is nearly complete. It is also shown that workers pay a significant ex-ante wage change premium (differential) to obtain a cost of living allowance clause. However, the realised contingent compensation exceeds such a premium for all industries. Finally, important spillover effects in wage setting and the decision to index the settlement have been detected.

A unifying approach to the empirical evaluation of asset pricing models

Two main approaches are commonly used to empirically evaluate linear factor pricingmodels: regression and SDF methods, with centred and uncentred versions of the latter.We show that unlike standard two-step or iterated GMM procedures, single-step estimatorssuch as continuously updated GMM yield numerically identical values for prices of risk,pricing errors, Jensen s alphas and overidentifying restrictions tests irrespective of the modelvalidity. Therefore, there is arguably a single approach regardless of the factors being tradedor not, or the use of excess or gross returns. We illustrate our results by revisiting Lustigand Verdelhan s (2007) empirical analysis of currency returns.

Variable Kernel estimates: On the impossibility of tuning the parameters

For the standard kernel density estimate, it is known that one can tune the bandwidth such that the expected L1 error is within a constant factor of the optimal L1 error (obtained when one is allowed to choose the bandwidth with knowledge of the density). In this paper, we pose the same problem for variable bandwidth kernel estimates where the bandwidths are allowed to depend upon the location. We show in particular that for positive kernels on the real line, for any data-based bandwidth, there exists a densityfor which the ratio of expected L1 error over optimal L1 error tends to infinity. Thus, the problem of tuning the variable bandwidth in an optimal manner is ``too hard''. Moreover, from the class of counterexamples exhibited in the paper, it appears thatplacing conditions on the densities (monotonicity, convexity, smoothness) does not help.

Power transformations in correspondence analysis

Power transformations of positive data tables, prior to applying the correspondence analysis algorithm, are shown to open up a family of methods with direct connections to the analysis of log-ratios. Two variations of this idea are illustrated. The first approach is simply to power the original data and perform a correspondence analysis this method is shown to converge to unweighted log-ratio analysis as the power parameter tends to zero. The second approach is to apply the power transformation to thecontingency ratios, that is the values in the table relative to expected values based on the marginals this method converges to weighted log-ratio analysis, or the spectral map. Two applications are described: first, a matrix of population genetic data which is inherently two-dimensional, and second, a larger cross-tabulation with higher dimensionality, from a linguistic analysis of several books.

Assessment of a new MRI classification of spinal stenosis based on cross-sectional morphology of the dural sac : FM2

Introduction: Quantitative measures of degree of lumbar spinal stenosis (LSS) such as antero-posterior diameter of the canal or dural sac cross sectional area vary widely and do not correlate with clinical symptoms or results of surgical decompression. In an effort to improve quantification of stenosis we have developed a grading system based on the morphology of the dural sac and its contents as seen on T2 axial images. The grading comprises seven categories ranging form normal to the most severe stenosis and takes into account the ratio of rootlet/CSF content. Material and methods: Fifty T2 axial MRI images taken at disc level from twenty seven symptomatic lumbar spinal stenosis patients who underwent decompressive surgery were classified into seven categories by five observers and reclassified 2 weeks later by the same investigators. Intra- and inter-observer reliability of the classification were assessed using Cohen's and Fleiss' kappa statistics, respectively. Results: Generally, the morphology grading system itself was well adopted by the observers. Its success in application is strongly influenced by the identification of the dural sac. The average intraobserver Cohen's kappa was 0.53 ± 0.2. The inter-observer Fleiss' kappa was 0.38 ± 0.02 in the first rating and 0.3 ± 0.03 in the second rating repeated after two weeks. Discussion: In this attempt, the teaching of the observers was limited to an introduction to the general idea of the morphology grading system and one example MRI image per category. The identification of the dimension of the dural sac may be a difficult issue in absence of complete T1 T2 MRI image series as it was the case here. The similarity of the CSF to possibly present fat on T2 images was the main reason of mismatch in the assignment of the cases to a category. The Fleiss correlation factors of the five observers are fair and the proposed morphology grading system is promising.

Combined analysis of grey matter voxel-based morphometry and white matter tract-based spatial statistics in late-life bipolar disorder.

Background: Previous magnetic resonance imaging (MRI) studies in young patients with bipolar disorder indicated the presence of grey matter concentration changes as well as microstructural alterations in white matter in various neocortical areas and the corpus callosum. Whether these structural changes are also present in elderly patients with bipolar disorder with long-lasting clinical evolution remains unclear. Methods: We performed a prospective MRI study of consecutive elderly, euthymic patients with bipolar disorder and healthy, elderly controls. We conducted a voxel-based morphometry (VBM) analysis and a tract-based spatial statistics (TBSS) analysis to assess fractional anisotropy and longitudinal, radial and mean diffusivity derived by diffusion tensor imaging (DTI). Results: We included 19 patients with bipolar disorder and 47 controls in our study. Fractional anisotropy was the most sensitive DTI marker and decreased significantly in the ventral part of the corpus callosum in patients with bipolar disorder. Longitudinal, radial and mean diffusivity showed no significant between-group differences. Grey matter concentration was reduced in patients with bipolar disorder in the right anterior insula, head of the caudate nucleus, nucleus accumbens, ventral putamen and frontal orbital cortex. Conversely, there was no grey matter concentration or fractional anisotropy increase in any brain region in patients with bipolar disorder compared with controls. Limitations: The major limitation of our study is the small number of patients with bipolar disorder. Conclusion: Our data document the concomitant presence of grey matter concentration decreases in the anterior limbic areas and the reduced fibre tract coherence in the corpus callosum of elderly patients with long-lasting bipolar disorder.

HACEK infective endocarditis: characteristics and outcomes from a large, multi-national cohort.

The HACEK organisms (Haemophilus species, Aggregatibacter species, Cardiobacterium hominis, Eikenella corrodens, and Kingella species) are rare causes of infective endocarditis (IE). The objective of this study is to describe the clinical characteristics and outcomes of patients with HACEK endocarditis (HE) in a large multi-national cohort. Patients hospitalized with definite or possible infective endocarditis by the International Collaboration on Endocarditis Prospective Cohort Study in 64 hospitals from 28 countries were included and characteristics of HE patients compared with IE due to other pathogens. Of 5591 patients enrolled, 77 (1.4%) had HE. HE was associated with a younger age (47 vs. 61 years; p<0.001), a higher prevalence of immunologic/vascular manifestations (32% vs. 20%; p<0.008) and stroke (25% vs. 17% p = 0.05) but a lower prevalence of congestive heart failure (15% vs. 30%; p = 0.004), death in-hospital (4% vs. 18%; p = 0.001) or after 1 year follow-up (6% vs. 20%; p = 0.01) than IE due to other pathogens (n = 5514). On multivariable analysis, stroke was associated with mitral valve vegetations (OR 3.60; CI 1.34-9.65; p<0.01) and younger age (OR 0.62; CI 0.49-0.90; p<0.01). The overall outcome of HE was excellent with the in-hospital mortality (4%) significantly better than for non-HE (18%; p<0.001). Prosthetic valve endocarditis was more common in HE (35%) than non-HE (24%). The outcome of prosthetic valve and native valve HE was excellent whether treated medically or with surgery. Current treatment is very successful for the management of both native valve prosthetic valve HE but further studies are needed to determine why HE has a predilection for younger people and to cause stroke. The small number of patients and observational design limit inferences on treatment strategies. Self selection of study sites limits epidemiological inferences.

Contribution biplots

In order to interpret the biplot it is necessary to know which points usually variables are the ones that are important contributors to the solution, and this information is available separately as part of the biplot s numerical results. We propose a new scaling of the display, called the contribution biplot, which incorporates this diagnostic directly into the graphical display, showing visually the important contributors and thus facilitating the biplot interpretation and often simplifying the graphical representation considerably. The contribution biplot can be applied to a wide variety of analyses such as correspondence analysis, principal component analysis, log-ratio analysis and the graphical results of a discriminant analysis/MANOVA, in fact to any method based on the singular-value decomposition. In the contribution biplot one set of points, usually the rows of the data matrix, optimally represent the spatial positions of the cases or sample units, according to some distance measure that usually incorporates some form of standardization unless all data are comparable in scale. The other set of points, usually the columns, is represented by vectors that are related to their contributions to the low-dimensional solution. A fringe benefit is that usually only one common scale for row and column points is needed on the principal axes, thus avoiding the problem of enlarging or contracting the scale of one set of points to make the biplot legible. Furthermore, this version of the biplot also solves the problem in correspondence analysis of low-frequency categories that are located on the periphery of the map, giving the false impression that they are important, when they are in fact contributing minimally to the solution.

Some hypothesis tests for the covariance matrix when the dimension is large compared to the sample size

This paper analyzes whether standard covariance matrix tests work whendimensionality is large, and in particular larger than sample size. Inthe latter case, the singularity of the sample covariance matrix makeslikelihood ratio tests degenerate, but other tests based on quadraticforms of sample covariance matrix eigenvalues remain well-defined. Westudy the consistency property and limiting distribution of these testsas dimensionality and sample size go to infinity together, with theirratio converging to a finite non-zero limit. We find that the existingtest for sphericity is robust against high dimensionality, but not thetest for equality of the covariance matrix to a given matrix. For thelatter test, we develop a new correction to the existing test statisticthat makes it robust against high dimensionality.

Pregnancy outcome following maternal exposure to statins: a multicentre prospective study.

OBJECTIVE: This contribution addresses the risk associated with exposure to statins during pregnancy. DESIGN: Multicentre observational prospective controlled study. SETTING: European Network of Teratology Information Services. POPULATION: Pregnant women who contacted one of 11 participating centres, seeking advice about exposure to statins during pregnancy, or to agents known to be nonteratogenic. METHODS: Pregnancies exposed during first trimester to statins were followed up prospectively, and their outcomes were compared with a matched control group. MAIN OUTCOME MEASURES: Rates of major birth defects, live births, miscarriages, elective terminations, preterm deliveries and gestational age and birthweight at delivery. RESULTS: We collected observations from 249 exposed pregnancies and 249 controls. The difference in the rate of major birth defects between the statin-exposed and the control groups was small and statistically nonsignificant (4.1% versus 2.7% odds ratio [OR] 1.5; 95% confidence interval [95% CI] 0.5-4.5, P = 0.43). In an adjusted Cox model, the difference between miscarriage rates was also small and not significant (hazard ratio 1.36, 95% CI 0.63-2.93, P = 0.43). Premature birth was more frequent in exposed pregnancies (16.1% versus 8.5%; OR 2.1, 95% CI 1.1-3.8, P = 0.019). Nonetheless, median gestational age at birth (39 weeks, interquartile range [IQR] 37-40 versus 39 weeks, IQR 38-40, P = 0.27) and birth weight (3280 g, IQR 2835-3590 versus 3250 g, IQR 2880-3630, P = 0.95) did not differ between exposed and non-exposed pregnancies. CONCLUSIONS: This study did not detect a teratogenic effect of statins. Its statistical power remains insufficient to challenge current recommendations of treatment discontinuation during pregnancy.

Vital Statistics of Iowa in Brief, 2006

This publication is an historical recording of the most requested statistics on vital events and is a source of information that can be used in further analysis.

Statistics in archaeology: New directions

Connections between Statistics and Archaeology have always appeared veryfruitful. The objective of this paper is to offer an outlook of somestatistical techniques that are being developed in the most recentyears and that can be of interest for archaeologists in the short run.

A singular function and its relation with the number systems involved in its definition

Minkowski's ?(x) function can be seen as the confrontation of two number systems: regular continued fractions and the alternated dyadic system. This way of looking at it permits us to prove that its derivative, as it also happens for many other non-decreasing singular functions from [0,1] to [0,1], when it exists can only attain two values: zero and infinity. It is also proved that if the average of the partial quotients in the continued fraction expansion of x is greater than k* =5.31972, and ?'(x) exists then ?'(x)=0. In the same way, if the same average is less than k**=2 log2(F), where F is the golden ratio, then ?'(x)=infinity. Finally some results are presented concerning metric properties of continued fraction and alternated dyadic expansions.

Correspondence analysis of raw data

Correspondence analysis has found extensive use in ecology, archeology, linguisticsand the social sciences as a method for visualizing the patterns of association in a table offrequencies or nonnegative ratio-scale data. Inherent to the method is the expression of the datain each row or each column relative to their respective totals, and it is these sets of relativevalues (called profiles) that are visualized. This relativization of the data makes perfect sensewhen the margins of the table represent samples from sub-populations of inherently differentsizes. But in some ecological applications sampling is performed on equal areas or equalvolumes so that the absolute levels of the observed occurrences may be of relevance, in whichcase relativization may not be required. In this paper we define the correspondence analysis ofthe raw unrelativized data and discuss its properties, comparing this new method to regularcorrespondence analysis and to a related variant of non-symmetric correspondence analysis.