New perspectives in the use of ink evidence in forensic science

Resume : L'utilisation de l'encre comme indice en sciences forensiques est décrite et encadrée par une littérature abondante, comprenant entre autres deux standards de l'American Society for Testing and Materials (ASTM). La grande majorité de cette littérature se préoccupe de l'analyse des caractéristiques physiques ou chimiques des encres. Les standards ASTM proposent quelques principes de base qui concernent la comparaison et l'interprétation de la valeur d'indice des encres en sciences forensiques. L'étude de cette littérature et plus particulièrement des standards ASTM, en ayant a l'esprit les développements intervenus dans le domaine de l'interprétation de l'indice forensique, montre qu'il existe un potentiel certain pour l'amélioration de l'utilisation de l'indice encre et de son impact dans l'enquête criminelle. Cette thèse propose d'interpréter l'indice encre en se basant sur le cadre défini par le théorème de Bayes. Cette proposition a nécessité le développement d'un système d'assurance qualité pour l'analyse et la comparaison d'échantillons d'encre. Ce système d'assurance qualité tire parti d'un cadre théorique nouvellement défini. La méthodologie qui est proposée dans ce travail a été testée de manière compréhensive, en tirant parti d'un set de données spécialement créer pour l'occasion et d'outils importés de la biométrie. Cette recherche répond de manière convaincante à un problème concret généralement rencontré en sciences forensiques. L'information fournie par le criminaliste, lors de l'examen de traces, est souvent bridée, car celui-ci essaie de répondre à la mauvaise question. L'utilisation d'un cadre théorique explicite qui définit et formalise le goal de l'examen criminaliste, permet de déterminer les besoins technologiques et en matière de données. Le développement de cette technologie et la collection des données pertinentes peut être justifiées économiquement et achevée de manière scientifique. Abstract : The contribution of ink evidence to forensic science is described and supported by an abundant literature and by two standards from the American Society for Testing and Materials (ASTM). The vast majority of the available literature is concerned with the physical and chemical analysis of ink evidence. The relevant ASTM standards mention some principles regarding the comparison of pairs of ink samples and the evaluation of their evidential value. The review of this literature and, more specifically, of the ASTM standards in the light of recent developments in the interpretation of forensic evidence has shown some potential improvements, which would maximise the benefits of the use of ink evidence in forensic science. This thesis proposes to interpret ink evidence using the widely accepted and recommended Bayesian theorem. This proposition has required the development of a new quality assurance process for the analysis and comparison of ink samples, as well as of the definition of a theoretical framework for ink evidence. The proposed technology has been extensively tested using a large dataset of ink samples and state of the art tools, commonly used in biometry. Overall, this research successfully answers to a concrete problem generally encountered in forensic science, where scientists tend to self-limit the usefulness of the information that is present in various types of evidence, by trying to answer to the wrong questions. The declaration of an explicit framework, which defines and formalises their goals and expected contributions to the criminal and civil justice system, enables the determination of their needs in terms of technology and data. The development of this technology and the collection of the data is then justified economically, structured scientifically and can be proceeded efficiently.

Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study.

The identification of genetically homogeneous groups of individuals is a long standing issue in population genetics. A recent Bayesian algorithm implemented in the software STRUCTURE allows the identification of such groups. However, the ability of this algorithm to detect the true number of clusters (K) in a sample of individuals when patterns of dispersal among populations are not homogeneous has not been tested. The goal of this study is to carry out such tests, using various dispersal scenarios from data generated with an individual-based model. We found that in most cases the estimated 'log probability of data' does not provide a correct estimation of the number of clusters, K. However, using an ad hoc statistic DeltaK based on the rate of change in the log probability of data between successive K values, we found that STRUCTURE accurately detects the uppermost hierarchical level of structure for the scenarios we tested. As might be expected, the results are sensitive to the type of genetic marker used (AFLP vs. microsatellite), the number of loci scored, the number of populations sampled, and the number of individuals typed in each sample.

An empirical hierarchical Bayesian unification of occupational exposure assessment methods.

In occupational exposure assessment of airborne contaminants, exposure levels can either be estimated through repeated measurements of the pollutant concentration in air, expert judgment or through exposure models that use information on the conditions of exposure as input. In this report, we propose an empirical hierarchical Bayesian model to unify these approaches. Prior to any measurement, the hygienist conducts an assessment to generate prior distributions of exposure determinants. Monte-Carlo samples from these distributions feed two level-2 models: a physical, two-compartment model, and a non-parametric, neural network model trained with existing exposure data. The outputs of these two models are weighted according to the expert's assessment of their relevance to yield predictive distributions of the long-term geometric mean and geometric standard deviation of the worker's exposure profile (level-1 model). Bayesian inferences are then drawn iteratively from subsequent measurements of worker exposure. Any traditional decision strategy based on a comparison with occupational exposure limits (e.g. mean exposure, exceedance strategies) can then be applied. Data on 82 workers exposed to 18 contaminants in 14 companies were used to validate the model with cross-validation techniques. A user-friendly program running the model is available upon request.

Multi-ethnic fine-mapping of 14 central adiposity loci.

The Genetic Investigation of Anthropometric Traits (GIANT) consortium identified 14 loci in European Ancestry (EA) individuals associated with waist-to-hip ratio (WHR) adjusted for body mass index. These loci are wide and narrowing the signals remains necessary. Twelve of 14 loci identified in GIANT EA samples retained strong associations with WHR in our joint EA/individuals of African Ancestry (AA) analysis (log-Bayes factor >6.1). Trans-ethnic analyses at five loci (TBX15-WARS2, LYPLAL1, ADAMTS9, LY86 and ITPR2-SSPN) substantially narrowed the signals to smaller sets of variants, some of which are in regions that have evidence of regulatory activity. By leveraging varying linkage disequilibrium structures across different populations, single-nucleotide polymorphisms (SNPs) with strong signals and narrower credible sets from trans-ethnic meta-analysis of central obesity provide more precise localizations of potential functional variants and suggest a possible regulatory role. Meta-analysis results for WHR were obtained from 77 167 EA participants from GIANT and 23 564 AA participants from the African Ancestry Anthropometry Genetics Consortium. For fine mapping we interrogated SNPs within ± 250 kb flanking regions of 14 previously reported index SNPs from loci discovered in EA populations by performing trans-ethnic meta-analysis of results from the EA and AA meta-analyses. We applied a Bayesian approach that leverages allelic heterogeneity across populations to combine meta-analysis results and aids in fine-mapping shared variants at these locations. We annotated variants using information from the ENCODE Consortium and Roadmap Epigenomics Project to prioritize variants for possible functionality.

Associer et évaluer

Le forensicien joue fréquemment le rôle d'expert judiciaire. Dans cette responsabilité, il doit par exemple évaluer le lien qu'est susceptible d'entretenir une trace avec une source putative (une personne, un objet) donnée. L'expert examine aussi parfois la relation entre deux traces (proviennent-elles de la même source ?), combine les indices, ou interprète les activités possibles en regard des traces recueillies. Un mouvement dominant et intégrateur préconise d'effectuer cette évaluation au moyen de probabilités qui expriment dans un cadre formel les incertitudes inhérentes à ce genre de questions. Basé sur le théorème découvert par le révérend Bayes sur la probabilité des causes, le rapport de vraisemblance équilibre le point de vue de l'accusation (par exemple, probabilité de la trace, sachant que x en est à l'origine) et de la défense (la probabilité de la trace si la source n'est pas x). L'exploitation de ce modèle fait appel à l'expérience de l'expert, mais aussi à de grandes quantités de données représentatives qui doivent encore être rassemblées. Une tendance consiste aussi à exploiter pour cela les données stockées par l'utilisation de systèmes automatisés (empreintes digitales, projectiles, etc.). De nouvelles possibilités d'extraire des caractéristiques encore inexplorées à partir des traces sont susceptibles d'augmenter considérablement les moyens pour les experts judiciaires d'apporter des réponses probabilistes dans des circonstances toujours plus variées aux autorités judiciaires qui les mandatent.

Synchronization of EEG: bivariate and multivariate measures.

Synchronization behavior of electroencephalographic (EEG) signals is important for decoding information processing in the human brain. Modern multichannel EEG allows a transition from traditional measurements of synchronization in pairs of EEG signals to whole-brain synchronization maps. The latter can be based on bivariate measures (BM) via averaging over pair-wise values or, alternatively, on multivariate measures (MM), which directly ascribe a single value to the synchronization in a group. In order to compare BM versus MM, we applied nine different estimators to simulated multivariate time series with known parameters and to real EEGs.We found widespread correlations between BM and MM, which were almost frequency-independent for all the measures except coherence. The analysis of the behavior of synchronization measures in simulated settings with variable coupling strength, connection probability, and parameter mismatch showed that some of them, including S-estimator, S-Renyi, omega, and coherence, aremore sensitive to linear interdependences,while others, like mutual information and phase locking value, are more responsive to nonlinear effects. Onemust consider these properties together with the fact thatMM are computationally less expensive and, therefore, more efficient for the large-scale data sets than BM while choosing a synchronization measure for EEG analysis.

Parametric estimation of pulse arrival time: a robust approach to pulse wave velocity.

Pulse wave velocity (PWV) is a surrogate of arterial stiffness and represents a non-invasive marker of cardiovascular risk. The non-invasive measurement of PWV requires tracking the arrival time of pressure pulses recorded in vivo, commonly referred to as pulse arrival time (PAT). In the state of the art, PAT is estimated by identifying a characteristic point of the pressure pulse waveform. This paper demonstrates that for ambulatory scenarios, where signal-to-noise ratios are below 10 dB, the performance in terms of repeatability of PAT measurements through characteristic points identification degrades drastically. Hence, we introduce a novel family of PAT estimators based on the parametric modeling of the anacrotic phase of a pressure pulse. In particular, we propose a parametric PAT estimator (TANH) that depicts high correlation with the Complior(R) characteristic point D1 (CC = 0.99), increases noise robustness and reduces by a five-fold factor the number of heartbeats required to obtain reliable PAT measurements.

Abnormal EEG Synchronization Correlates with Demyelination in Alzheimer's Disease

Introduction : The pathological processes caused by Alzheimer's disease (AD) supposedly disrupt communication between and within the distributed cortical networks due to the dysfunction/loss of synapses and myelination breakdown. Indeed, recently (Knyazeva et al. 2008), we have revealed the whole-head topography of EEG synchronization specific to AD. Here we analyze whether and how these abnormalities of synchronization are related to the demyelination of cortico-cortical fibers. Methods : Fifteen newly diagnosed AD patients (CDR 0.5-1) and 15 controls matched for age, participated in the study. Their multichannel (128) EEGs were recorded during 3-5 min at rest. They were submitted to the multivariate phase synchronization (MPS) analysis for mapping regional synchronization. To obtain individual whole-head maps, the MPS was computed for each sensor considering its 2nd nearest topographical neighbors. Separate calculations were performed for the delta, theta, alpha-1/−2, and beta-1/−2 EEG bands. The same subjects were scanned on a 3 Tesla Philips scanner. The protocol included a high-resolution T1-weighted sequence and a Magnetization Transfer Imaging (MTI) acquisition. For each subject, we defined a 3mm thick layer of white matter exactly below the cortical gray matter. The magnetization transfer ratio (MTR) - an estimator of myelination - was calculated for this layer in 39 Brodmann-defined ROIs per hemisphere. To assess the between-group differences, we used a permutation version of Hotelling's T2 test or two-sample T-test (Pcorrected <0.05). For correlation analysis, Spearman Rank Correlation was calculated. Results : In AD patients, we have found an abnormal landscape of synchronization characterized by a decrease in MPS over the fronto-temporal region of the left hemisphere and an increase over the temporo-parieto-occipital regions bilaterally. Also, we have shown a widespread decrease in regional MTR in the AD patients for all the areas excluding motor, premotor, and primary sensory ones. Assuming that AD-related changes in synchronization are associated with demyelination, we hypothesized a correlation between the regional MTR values and MPS values in the hypo- and hyper-synchronized clusters. We found that MPS in the left fronto-temporal hypo-synchronized cluster directly correlates with myelination in BA42-46 of the left hemisphere: the lower the myelination in individual patients, the lower the EEG synchronization. By contrast, in the posterior hyper-synchronized cluster, MPS inversely correlated with myelination, i.e., the lower the myelination, the higher the synchronization. This posterior hyper-synchronization, more characteristic for early-onset AD, probably, results from the initial effect of the disease on cortical inhibition, reducing cortical capacity for decoupling irrelevant connections. Remarkably, it showed different topography of correlations in early- vs. late-onset patients. In the early-onset patients, hyper-synchronization was mainly related to demyelination in posterior BAs, the effect being significant in all the EEG frequency bands. In the late-onset patients, widely distributed correlations were significant for the EEG delta band, suggesting an interaction between the cerebral manifestations of AD and the age of its onset, i.e., topographically selective impairment of cortical inhibition in early-onset AD vs. its wide-spread weakening in old age. Conclusions : Overall, our results document that the degradation of white matter is a significant factor of AD pathogenesis leading to functional dysconnection, the latter being reflected in EEG synchronization abnormalities.

Spatio-temporal trends of mortality in small areas of Southern Spain

Background: Most mortality atlases show static maps from count data aggregated over time. This procedure has several methodological problems and serious limitations for decision making in Public Health. The evaluation of health outcomes, including mortality, should be approached from a dynamic time perspective that is specific for each gender and age group. At the moment, researches in Spain do not provide a dynamic image of the population’s mortality status from a spatio-temporal point of view. The aim of this paper is to describe the spatial distribution of mortality from all causes in small areas of Andalusia (Southern Spain) and evolution over time from 1981 to 2006. Methods: A small-area ecological study was devised using the municipality as the unit for analysis. Two spatiotemporal hierarchical Bayesian models were estimated for each age group and gender. One of these was used to estimate the specific mortality rate, together with its time trends, and the other to estimate the specific rate ratio for each municipality compared with Spain as a whole. Results: More than 97% of the municipalities showed a diminishing or flat mortality trend in all gender and age groups. In 2006, over 95% of municipalities showed male and female mortality specific rates similar or significantly lower than Spanish rates for all age groups below 65. Systematically, municipalities in Western Andalusia showed significant male and female mortality excess from 1981 to 2006 only in age groups over 65. Conclusions: The study shows a dynamic geographical distribution of mortality, with a different pattern for each year, gender and age group. This information will contribute towards a reflection on the past, present and future of mortality in Andalusia.

Area deprivation and mortality in the provincial capital cities of Andalusia and Catalonia (Spain)

Objective: To study the linkage between material deprivation and mortality from all causes, for men and women separately, in the capital cities of the provinces in Andalusia and Catalonia (Spain). Methods: A small-area ecological study was devised using the census section as the unit for analysis. 188 983 Deaths occurring in the capital cities of the Andalusian provinces and 109 478 deaths recorded in the Catalan capital cities were examined. Principal components factorial analysis was used to devise a material deprivation index comprising the percentage of manual labourers, unemployment and illiteracy. A hierarchical Bayesian model was used to study the relationship between mortality and area deprivation. Main results: In most cities, results show an increased male mortality risk in the most deprived areas in relation to the least depressed. In Andalusia, the relative risks between the highest and lowest deprivation decile ranged from 1.24 (Malaga) to 1.40 (Granada), with 95% credibility intervals showing a significant excess risk. In Catalonia, relative risks ranged between 1.08 (Girona) and 1.50 (Tarragona). No evidence was found for an excess of female mortality in most deprived areas in either of the autonomous communities. Conclusions: Within cities, gender-related differences were revealed when deprivation was correlated geographically with mortality rates. These differences were found from an ecological perspective. Further research is needed in order to validate these results from an individual approach. The idea to be analysed is to identify those factors that explain these differences at an individual level.

Cancer mortality inequalities in urban areas: a Bayesian small area analysis in Spanish cities

Background Intra-urban inequalities in mortality have been infrequently analysed in European contexts. The aim of the present study was to analyse patterns of cancer mortality and their relationship with socioeconomic deprivation in small areas in 11 Spanish cities. Methods It is a cross-sectional ecological design using mortality data (years 1996-2003). Units of analysis were the census tracts. A deprivation index was calculated for each census tract. In order to control the variability in estimating the risk of dying we used Bayesian models. We present the RR of the census tract with the highest deprivation vs. the census tract with the lowest deprivation. Results In the case of men, socioeconomic inequalities are observed in total cancer mortality in all cities, except in Castellon, Cordoba and Vigo, while Barcelona (RR = 1.53 95%CI 1.42-1.67), Madrid (RR = 1.57 95%CI 1.49-1.65) and Seville (RR = 1.53 95%CI 1.36-1.74) present the greatest inequalities. In general Barcelona and Madrid, present inequalities for most types of cancer. Among women for total cancer mortality, inequalities have only been found in Barcelona and Zaragoza. The excess number of cancer deaths due to socioeconomic deprivation was 16,413 for men and 1,142 for women. Conclusion This study has analysed inequalities in cancer mortality in small areas of cities in Spain, not only relating this mortality with socioeconomic deprivation, but also calculating the excess mortality which may be attributed to such deprivation. This knowledge is particularly useful to determine which geographical areas in each city need intersectorial policies in order to promote a healthy environment.

Atlas interactivo de mortalidad de Andalucía (AIMA)

Until now, mortality atlases have been static. Most of them describe the geographical distribution of mortality using count data aggregated over time and standardized mortality rates. However, this methodology has several limitations. Count data aggregated over time produce a bias in the estimation of death rates. Moreover, this practice difficult the study of temporal changes in geographical distribution of mortality. On the other hand, using standardized mortality hamper to check differences in mortality among groups. The Interactive Mortality Atlas in Andalusia (AIMA) is an alternative to conventional static atlases. It is a dynamic Geographical Information System that allows visualizing in web-site more than 12.000 maps and 338.00 graphics related to the spatio-temporal distribution of the main death causes in Andalusia by age and sex groups from 1981. The objective of this paper is to describe the methods used for AIMA development, to show technical specifications and to present their interactivity. The system is available from the link products in www.demap.es. AIMA is the first interactive GIS that have been developed in Spain with these characteristics. Spatio-temporal Hierarchical Bayesian Models were used for statistical data analysis. The results were integrated into web-site using a PHP environment and a dynamic cartography in Flash. Thematic maps in AIMA demonstrate that the geographical distribution of mortality is dynamic, with differences among year, age and sex groups. The information nowadays provided by AIMA and the future updating will contribute to reflect on the past, the present and the future of population health in Andalusia.

Population pharmacokinetics and effects of efavirenz in patients with human immunodeficiency virus infection.

OBJECTIVE: The reverse transcriptase inhibitor efavirenz is currently used at a fixed dose of 600 mg/d. However, dosage individualization based on plasma concentration monitoring might be indicated. This study aimed to assess the efavirenz pharmacokinetic profile and interpatient versus intrapatient variability in patients who are positive for human immunodeficiency virus, to explore the relationship between drug exposure, efficacy, and central nervous system toxicity and to build up a Bayesian approach for dosage adaptation. METHODS: The population pharmacokinetic analysis was performed by use of NONMEM based on plasma samples from a cohort of unselected patients receiving efavirenz. With the use of a 1-compartment model with first-order absorption, the influence of demographic and clinical characteristics on oral clearance and oral volume of distribution was examined. The average drug exposure during 1 dosing interval was estimated for each patient and correlated with markers of efficacy and toxicity. The population kinetic parameters and the variabilities were integrated into a Bayesian equation for dosage adaptation based on a single plasma sample. RESULTS: Data from 235 patients with a total of 719 efavirenz concentrations were collected. Oral clearance was 9.4 L/h, oral volume of distribution was 252 L, and the absorption rate constant was 0.3 h(-1). Neither the demographic covariates evaluated nor the comedications showed a clinically significant influence on efavirenz pharmacokinetics. A large interpatient variability was found to affect efavirenz relative bioavailability (coefficient of variation, 54.6%), whereas the intrapatient variability was small (coefficient of variation, 26%). An inverse correlation between average drug exposure and viral load and a trend with central nervous system toxicity were detected. This enabled the derivation of a dosing adaptation strategy suitable to bring the average concentration into a therapeutic target from 1000 to 4000 microg/L to optimize viral load suppression and to minimize central nervous system toxicity. CONCLUSIONS: The high interpatient and low intrapatient variability values, as well as the potential relationship with markers of efficacy and toxicity, support the therapeutic drug monitoring of efavirenz. However, further evaluation is needed before individualization of an efavirenz dosage regimen based on routine drug level monitoring should be recommended for optimal patient management.

Municipal distribution of bladder cancer mortality in Spain: Possible role of mining and industry

BACKGROUND Spain shows the highest bladder cancer incidence rates in men among European countries. The most important risk factors are tobacco smoking and occupational exposure to a range of different chemical substances, such as aromatic amines. METHODS This paper describes the municipal distribution of bladder cancer mortality and attempts to "adjust" this spatial pattern for the prevalence of smokers, using the autoregressive spatial model proposed by Besag, York and Molliè, with relative risk of lung cancer mortality as a surrogate. RESULTS It has been possible to compile and ascertain the posterior distribution of relative risk for bladder cancer adjusted for lung cancer mortality, on the basis of a single Bayesian spatial model covering all of Spain's 8077 towns. Maps were plotted depicting smoothed relative risk (RR) estimates, and the distribution of the posterior probability of RR>1 by sex. Towns that registered the highest relative risks for both sexes were mostly located in the Provinces of Cadiz, Seville, Huelva, Barcelona and Almería. The highest-risk area in Barcelona Province corresponded to very specific municipal areas in the Bages district, e.g., Suría, Sallent, Balsareny, Manresa and Cardona. CONCLUSION Mining/industrial pollution and the risk entailed in certain occupational exposures could in part be dictating the pattern of municipal bladder cancer mortality in Spain. Population exposure to arsenic is a matter that calls for attention. It would be of great interest if the relationship between the chemical quality of drinking water and the frequency of bladder cancer could be studied.

Inflammatory Markers and Incident Heart Failure Risk in Older Adults: The Health, Aging, and Body Composition Study

Background: Inflammation is associated with heart failure (HF) risk factors and also directly affects myocardial function. However, the association between inflammation and HF risk in older adults has not been adequately evaluated. Methods: The association of baseline serum concentrations of interleukin-6 (IL-6), tumor necrosis factor alpha (TNF- ), and C-reactive protein (CRP) with incident HF was assessed with Cox proportional hazards models among 2610 older persons without prevalent HF enrolled in the Health, Aging, and Body Composition (Health ABC) Study (age, 73.6±2.9 years; 48.3% men; 59.6% white). Results: Median (interquartile range) baseline concentrations of IL-6, TNF- , and CRP were 1.80 (1.23, 2.76) pg/mL, 3.14 (2.41, 4.06) pg/mL, and 1.64 (0.99, 3.04) µg/mL, respectively. On follow-up (median, 9.4 years), 311 participants (11.9%) developed HF. In models controlling for clinical predictors of HF and incident coronary heart disease, doubling of IL-6, TNF- , and CRP concentrations was associated with 34% (95% CI, 18 -52%; P<.001), 33% (95% CI, 9 - 63%; P=.006), and 13% (95% CI, 3-24%; P=.01) increase in HF risk, respectively. In models including all 3 markers, IL-6 and TNF- , but not CRP, remained significant. Findings were similar across sex and race. Post-HF ejection fraction (EF) was available in 239 (76.8%) cases. When only cases with preserved EF were considered (n=105), IL-6 (HR per doubling, 1.57; 95% CI, 1.28 -1.94; P<.001), TNF- (HR per doubling, 1.59; 95% CI, 1.12-2.26; P=.01), and CRP (HR per doubling, 1.23; 95% CI, 1.05-1.44; P=.01) were all associated with HF risk in adjusted models. In contrast, when only cases with reduced EF (n=134) were considered, only IL-6 attained marginal significance in adjusted models (HR per doubling, 1.20; 95% CI, 0.99 -1.46; P=.06). Participants with 2 or 3 markers above median had pronounced HF risk in adjusted models (HR, 1.66; 95% CI, 1.12-2.46; P=.01; and HR, 1.76; 95% CI, 1.16 -2.65; P=.007, respectively). Addition of IL-6 to the clinical Health ABC HF model improved discrimination (C index from 0.717 to 0.734; P=.001) and fit (decreased Bayes information criterion by 17.8; P<.001). Conclusions: Inflammatory markers are associated with HF risk among older adults and may improve HF risk stratification.