Compositional Data Analysis in Forensic Geoscience

Statistics are regularly used to make some form of comparison between trace evidence or deploy the exclusionary principle (Morgan and Bull, 2007) in forensic investigations. Trace evidence are routinely the results of particle size, chemical or modal analyses and as such constitute compositional data. The issue is that compositional data including percentages, parts per million etc. only carry relative information. This may be problematic where a comparison of percentages and other constraint/closed data is deemed a statistically valid and appropriate way to present trace evidence in a court of law. Notwithstanding an awareness of the existence of the constant sum problem since the seminal works of Pearson (1896) and Chayes (1960) and the introduction of the application of log-ratio techniques (Aitchison, 1986; Pawlowsky-Glahn and Egozcue, 2001; Pawlowsky-Glahn and Buccianti, 2011; Tolosana-Delgado and van den Boogaart, 2013) the problem that a constant sum destroys the potential independence of variances and covariances required for correlation regression analysis and empirical multivariate methods (principal component analysis, cluster analysis, discriminant analysis, canonical correlation) is all too often not acknowledged in the statistical treatment of trace evidence. Yet the need for a robust treatment of forensic trace evidence analyses is obvious. This research examines the issues and potential pitfalls for forensic investigators if the constant sum constraint is ignored in the analysis and presentation of forensic trace evidence. Forensic case studies involving particle size and mineral analyses as trace evidence are used to demonstrate the use of a compositional data approach using a centred log-ratio (clr) transformation and multivariate statistical analyses.

Compositional data-driven alternatives to single geochemical component maps

Conventional practice in Regional Geochemistry includes as a final step of any geochemical campaign the generation of a series of maps, to show the spatial distribution of each of the components considered. Such maps, though necessary, do not comply with the compositional, relative nature of the data, which unfortunately make any conclusion based on them sensitive
to spurious correlation problems. This is one of the reasons why these maps are never interpreted isolated. This contribution aims at gathering a series of statistical methods to produce individual maps of multiplicative combinations of components (logcontrasts), much in the flavor of equilibrium constants, which are designed on purpose to capture certain aspects of the data.
We distinguish between supervised and unsupervised methods, where the first require an external, non-compositional variable (besides the compositional geochemical information) available in an analogous training set. This external variable can be a quantity (soil density, collocated magnetics, collocated ratio of Th/U spectral gamma counts, proportion of clay particle fraction, etc) or a category (rock type, land use type, etc). In the supervised methods, a regression-like model between the external variable and the geochemical composition is derived in the training set, and then this model is mapped on the whole region. This case is illustrated with the Tellus dataset, covering Northern Ireland at a density of 1 soil sample per 2 square km, where we map the presence of blanket peat and the underlying geology. The unsupervised methods considered include principal components and principal balances
(Pawlowsky-Glahn et al., CoDaWork2013), i.e. logcontrasts of the data that are devised to capture very large variability or else be quasi-constant. Using the Tellus dataset again, it is found that geological features are highlighted by the quasi-constant ratios Hf/Nb and their ratio against SiO2; Rb/K2O and Zr/Na2O and the balance between these two groups of two variables; the balance of Al2O3 and TiO2 vs. MgO; or the balance of Cr, Ni and Co vs. V and Fe2O3. The largest variability appears to be related to the presence/absence of peat.

Geostatistical Fisher discriminant analysis

A geostatistical version of the classical Fisher rule (linear discriminant analysis) is presented.This method is applicable when a large dataset of multivariate observations is available within a domain split in several known subdomains, and it assumes that the variograms (or covariance functions) are comparable between subdomains, which only differ in the mean values of the available variables. The method consists on finding the eigen-decomposition of the matrix W-1B, where W is the matrix of sills of all direct- and cross-variograms, and B is the covariance matrix of the vectors of weighted means within each subdomain, obtained by generalized least squares. The method is used to map peat blanket occurrence in Northern Ireland, with data from the Tellus
survey, which requires a minimal change to the general recipe: to use compositionally-compliant variogram tools and models, and work with log-ratio transformed data.

What defines a geochemical map? Implications for environmental assessment.

The environmental quality of land is often assessed by the calculation of threshold values which aim to differentiate between concentrations of elements based on whether the soils are in residential or industrial sites. In Europe, for example, soil guideline values exist for agricultural and grazing land. A threshold is often set to differentiate between concentrations of the element that naturally occur in the soil and concentrations that result from diffuse anthropogenic sources. Regional geochemistry and, in particular, single component geochemical maps are increasingly being used to determine these baseline environmental assessments. The key question raised in this paper is whether the geochemical map can provide an accurate interpretation on its own. Implicit is the thought that single component geochemical maps represent absolute abundances. However,because of the compositional (closed) nature of the data univariate geochemical maps cannot be compared directly with one another.. As a result, any interpretation based on them is vulnerable to spurious correlation problems. What does this mean for soil geochemistry mapping, baseline quality documentation, soil resource assessment or risk evaluation? Despite the limitation of relative abundances, individual raw geochemical maps are deemed fundamental to several applications of geochemical maps including environmental assessments. However, element toxicity is related to its bioavailable concentration, which is lowered if its source is mixed with another source. Elements interact, for example under reducing conditions with iron oxides, its solid state is lost and arsenic becomes soluble and mobile. Both of these matters may be more adequately dealt with if a single component map is not interpreted in isolation to determine baseline and threshold assessments. A range of alternative compositionally compliant representations based on log-ratio and log-contrast approaches are explored to supplement the classical single component maps for environmental assessment. Case study examples are shown based on the Tellus soil geochemical dataset, covering Northern Ireland and the results of in vitro oral bioaccessibility testing carried out on a sub-set of archived Tellus Survey shallow soils following the Unified BARGE (Bioaccessibility Research Group of Europe).

Green Space Attachment and Health: A Comparative Study in Two Urban Neighborhoods

The positive relationships between urban green space and health have been well documented. Little is known, however, about the role of residents’ emotional attachment to local green spaces in these relationships, and how attachment to green spaces and health may be promoted by the availability of accessible and usable green spaces. The present research aimed to examine the links between self-reported health, attachment to green space, and the availability of accessible and usable green spaces. Data were collected via paper-mailed surveys in two neighborhoods (n = 223) of a medium-sized Dutch city in the Netherlands. These neighborhoods differ in the perceived and objectively measured accessibility and usability of green spaces, but are matched in the physically available amount of urban green space, as well as in demographic and socio-economic status, and housing conditions. Four dimensions of green space attachment were identified through confirmatory factor analysis: place dependence, affective attachment, place identity and social bonding. The results show greater attachment to local green space and better self-reported mental health in the neighborhood with higher availability of accessible and usable green spaces. The two neighborhoods did not differ, however, in physical and general health. Structural Equation Modelling confirmed the neighborhood differences in green space attachment and mental health, and also revealed a positive path from green space attachment to mental health. These findings convey the message that we should make green places, instead of green spaces.

Visual symptoms and retinal straylight after laser peripheral iridotomy: the Zhongshan Angle-Closure Prevention Trial.

OBJECTIVE: To assess the impact of laser peripheral iridotomy (LPI) on forward-scatter of light and subjective visual symptoms and to identify LPI parameters influencing these phenomena. DESIGN: Cohort study derived from a randomized trial, using an external control group. PARTICIPANTS: Chinese subjects initially aged 50 or older and 70 years or younger with bilateral narrow angles undergoing LPI in 1 eye selected at random, and age- and gender-matched controls. METHODS: Eighteen months after laser, LPI-treated subjects underwent digital iris photography and photogrammetry to characterize the size and location of the LPI, Lens Opacity Classification System III cataract grading, and measurement of retinal straylight (C-Quant; OCULUS, Wetzlar, Germany) in the treated and untreated eyes and completed a visual symptoms questionnaire. Controls answered the questionnaire and underwent straylight measurement and (in a random one-sixth sample) cataract grading. MAIN OUTCOME MEASURES: Retinal straylight levels and subjective visual symptoms. RESULTS: Among 230 LPI-treated subjects (121 [58.8%] with LPI totally covered by the lid, 43 [19.8%] with LPI partly covered by the lid, 53 [24.4%] with LPI uncovered by the lid), 217 (94.3%) completed all testing, as did 250 (93.3%) of 268 controls. Age, gender, and prevalence of visual symptoms did not differ between treated subjects and controls, although nuclear (P<0.01) and cortical (P = 0.03) cataract were less common among controls. Neither presenting visual acuity nor straylight score differed between the treated and untreated eyes among all treated persons, nor among those (n = 96) with LPI partially or totally uncovered. Prevalence of subjective glare did not differ significantly between participants with totally covered LPI (6.61%; 95% confidence interval [CI], 3.39%-12.5%), partially covered LPI (11.6%; 95% CI, 5.07%-24.5%), or totally uncovered LPI (9.43%; 95% CI, 4.10%-10.3%). In regression models, only worse cortical cataract grade (P = 0.01) was associated significantly with straylight score, and no predictors were associated with subjective glare. None of the LPI size or location parameters were associated with straylight or subjective symptoms. CONCLUSIONS: These results suggests that LPI is safe regarding measures of straylight and visual symptoms. This randomized design provides strong evidence that treatment programs for narrow angles would be unlikely to result in important medium-term visual disability.

Analgesic exposure in pregnant rats affects fetal germ cell development with inter-generational reproductive consequences

Analgesics which affect prostaglandin (PG) pathways are used by most pregnant women. As germ cells (GC) undergo developmental and epigenetic changes in fetal life and are PG targets, we investigated if exposure of pregnant rats to analgesics (indomethacin or acetaminophen) affected GC development and reproductive function in resulting offspring (F1) or in the F2 generation. Exposure to either analgesic reduced F1 fetal GC number in both sexes and altered the tempo of fetal GC development sex-dependently, with delayed meiotic entry in oogonia but accelerated GC differentiation in males. These effects persisted in adult F1 females as reduced ovarian and litter size, whereas F1 males recovered normal GC numbers and fertility by adulthood. F2 offspring deriving from an analgesic-exposed F1 parent also exhibited sex-specific changes. F2 males exhibited normal reproductive development whereas F2 females had smaller ovaries and reduced follicle numbers during puberty/adulthood; as similar changes were found for F2 offspring of analgesic-exposed F1 fathers or mothers, we interpret this as potentially indicating an analgesic-induced change to GC in F1. Assuming our results are translatable to humans, they raise concerns that analgesic use in pregnancy could potentially affect fertility of resulting daughters and grand-daughters.

The Top 10 Oomycete Pathogens in Molecular Plant Pathology

Oomycetes form a deep lineage of eukaryotic organisms that includes a large number of plant pathogens that threaten natural and managed ecosystems. We undertook a survey to query the community for their ranking of plant pathogenic oomycete species based on scientific and economic importance. In total, we received 263 votes from 62 scientists in 15 countries for a total of 33 species. The Top 10 species and their ranking are: (1) Phytophthora infestans; (2, tied) Hyaloperonospora arabidopsidis; (2, tied) Phytophthora ramorum; (4) Phytophthora sojae; (5) Phytophthora capsici; (6) Plasmopara viticola; (7) Phytophthora cinnamomi; (8, tied) Phytophthora parasitica; (8, tied) Pythium ultimum; and (10) Albugo candida. The article provides an introduction to these 10 taxa and a snapshot of current research. We hope that the list will serve as a benchmark for future trends in oomycete research.

Prediction of breast cancer risk based on profiling with common genetic variants

BACKGROUND: Data for multiple common susceptibility alleles for breast cancer may be combined to identify women at different levels of breast cancer risk. Such stratification could guide preventive and screening strategies. However, empirical evidence for genetic risk stratification is lacking. METHODS: We investigated the value of using 77 breast cancer-associated single nucleotide polymorphisms (SNPs) for risk stratification, in a study of 33 673 breast cancer cases and 33 381 control women of European origin. We tested all possible pair-wise multiplicative interactions and constructed a 77-SNP polygenic risk score (PRS) for breast cancer overall and by estrogen receptor (ER) status. Absolute risks of breast cancer by PRS were derived from relative risk estimates and UK incidence and mortality rates. RESULTS: There was no strong evidence for departure from a multiplicative model for any SNP pair. Women in the highest 1% of the PRS had a three-fold increased risk of developing breast cancer compared with women in the middle quintile (odds ratio [OR] = 3.36, 95% confidence interval [CI] = 2.95 to 3.83). The ORs for ER-positive and ER-negative disease were 3.73 (95% CI = 3.24 to 4.30) and 2.80 (95% CI = 2.26 to 3.46), respectively. Lifetime risk of breast cancer for women in the lowest and highest quintiles of the PRS were 5.2% and 16.6% for a woman without family history, and 8.6% and 24.4% for a woman with a first-degree family history of breast cancer. CONCLUSIONS: The PRS stratifies breast cancer risk in women both with and without a family history of breast cancer. The observed level of risk discrimination could inform targeted screening and prevention strategies. Further discrimination may be achievable through combining the PRS with lifestyle/environmental factors, although these were not considered in this report.

Quality assurance in the 22991 EORTC ROG trial in localized prostate cancer: dummy run and individual case review.

INTRODUCTION: EORTC trial 22991 was designed to evaluate the addition of concomitant and adjuvant short-term hormonal treatments to curative radiotherapy in terms of disease-free survival for patients with intermediate risk localized prostate cancer. In order to assess the compliance to the 3D conformal radiotherapy protocol guidelines, all participating centres were requested to participate in a dummy run procedure. An individual case review was performed for the largest recruiting centres as well. MATERIALS AND METHODS: CT-data of an eligible prostate cancer patient were sent to 30 centres including a description of the clinical case. The investigator was requested to delineate the volumes of interest and to perform treatment planning according to the protocol. Thereafter, the investigators of the 12 most actively recruiting centres were requested to provide data on five randomly selected patients for an individual case review. RESULTS: Volume delineation varied significantly between investigators. Dose constraints for organs at risk (rectum, bladder, hips) were difficult to meet. In the individual case review, no major protocol deviations were observed, but a number of dose reporting problems were documented for centres using IMRT. CONCLUSIONS: Overall, results of this quality assurance program were satisfactory. The efficacy of the combination of a dummy run procedure with an individual case review is confirmed in this study, as none of the evaluated patient files harboured a major protocol deviation. Quality assurance remains a very important tool in radiotherapy to increase the reliability of the trial results. Special attention should be given when designing quality assurance programs for more complex irradiation techniques.

Low-Dose Vascular Photodynamic Therapy Decreases Tumor Interstitial Fluid Pressure, which Promotes Liposomal Doxorubicin Distribution in a Murine Sarcoma Metastasis Model.

INTRODUCTION: Solid tumors are known to have an abnormal vasculature that limits the distribution of chemotherapy. We have recently shown that tumor vessel modulation by low-dose photodynamic therapy (L-PDT) could improve the uptake of macromolecular chemotherapeutic agents such as liposomal doxorubicin (Liporubicin) administered subsequently. However, how this occurs is unknown. Convection, the main mechanism for drug transport between the intravascular and extravascular spaces, is mostly related to interstitial fluid pressure (IFP) and tumor blood flow (TBF). Here, we determined the changes of tumor and surrounding lung IFP and TBF before, during, and after vascular L-PDT. We also evaluated the effect of these changes on the distribution of Liporubicin administered intravenously (IV) in a lung sarcoma metastasis model. MATERIALS AND METHODS: A syngeneic methylcholanthrene-induced sarcoma cell line was implanted subpleurally in the lung of Fischer rats. Tumor/surrounding lung IFP and TBF changes induced by L-PDT were determined using the wick-in-needle technique and laser Doppler flowmetry, respectively. The spatial distribution of Liporubicin in tumor and lung tissues following IV drug administration was then assessed in L-PDT-pretreated animals and controls (no L-PDT) by epifluorescence microscopy. RESULTS: L-PDT significantly decreased tumor but not lung IFP compared to controls (no L-PDT) without affecting TBF. These conditions were associated with a significant improvement in Liporubicin distribution in tumor tissues compared to controls (P < .05). DISCUSSION: L-PDT specifically enhanced convection in blood vessels of tumor but not of normal lung tissue, which was associated with a significant improvement of Liporubicin distribution in tumors compared to controls.