Modeling the forensic two-trace problem with Bayesian networks

The forensic two-trace problem is a perplexing inference problem introduced by Evett (J Forensic Sci Soc 27:375-381, 1987). Different possible ways of wording the competing pair of propositions (i.e., one proposition advanced by the prosecution and one proposition advanced by the defence) led to different quantifications of the value of the evidence (Meester and Sjerps in Biometrics 59:727-732, 2003). Here, we re-examine this scenario with the aim of clarifying the interrelationships that exist between the different solutions, and in this way, produce a global vision of the problem. We propose to investigate the different expressions for evaluating the value of the evidence by using a graphical approach, i.e. Bayesian networks, to model the rationale behind each of the proposed solutions and the assumptions made on the unknown parameters in this problem.

Further insights on the French WISC-IV factor structure through Bayesian structural equation modeling

The interpretation of the Wechsler Intelligence Scale for Children-Fourth Edition (WISC-IV) is based on a 4-factor model, which is only partially compatible with the mainstream Cattell-Horn-Carroll (CHC) model of intelligence measurement. The structure of cognitive batteries is frequently analyzed via exploratory factor analysis and/or confirmatory factor analysis. With classical confirmatory factor analysis, almost all crossloadings between latent variables and measures are fixed to zero in order to allow the model to be identified. However, inappropriate zero cross-loadings can contribute to poor model fit, distorted factors, and biased factor correlations; most important, they do not necessarily faithfully reflect theory. To deal with these methodological and theoretical limitations, we used a new statistical approach, Bayesian structural equation modeling (BSEM), among a sample of 249 French-speaking Swiss children (8-12 years). With BSEM, zero-fixed cross-loadings between latent variables and measures are replaced by approximate zeros, based on informative, small-variance priors. Results indicated that a direct hierarchical CHC-based model with 5 factors plus a general intelligence factor better represented the structure of the WISC-IV than did the 4-factor structure and the higher order models. Because a direct hierarchical CHC model was more adequate, it was concluded that the general factor should be considered as a breadth rather than a superordinate factor. Because it was possible for us to estimate the influence of each of the latent variables on the 15 subtest scores, BSEM allowed improvement of the understanding of the structure of intelligence tests and the clinical interpretation of the subtest scores.

VARIABILITY OF PHENOTYPIC, PROTEOMIC AND TRANSCRIPTOMIC EXPRESSION OF STAPHYLOCOCCUS AUREUS SURFACE ADHESINS

Staphylococcus aureus is a highly successful pathogen responsible of a wide variety of diseases, from minor skin infection to life-threatening sepsis or infective endocarditis, as well as food poisoning and toxic shock syndrome. This heterogeneity of infections and the ability of S. aureus to develop antibiotic-resistance to virtually any available drugs reflect its extraordinary capacity to adapt and survive in a great variety of environments. The pathogenesis of S. aureus infection involves a wide range of cell wall-associated adhesins and extracellular toxins that promote host colonization and invasion. In addition, S. aureus is extremely well equipped with regulatory systems that sense environmental conditions and respond by fine tuning the expression of metabolic and virulence determinants. Surface adhesins referred to MSCRAMMs - for Microbial Surface Component Recognizing Adherence Matrix Molecules - mediate binding to the host extracellular matrix or serum components, including fibrinogen, fibronectin, collagen and elastin, and promote tissue colonization and invasion. Major MSCRAMMs include a family of surface-attached proteins covalently bound to the cell wall peptidoglycan via a conserved LPXTG motif. Genomic analyses indicate that S. aureus contain up to 22 LPXTG surface proteins, which could potentially act individually or in synergy to promote infection. In the first part of this study we determined the range of adherence phenotypes to fibrinogen and fibronectin among 30 carriage isolates of S. aureus and compared it to the adherence phenotypes of 30 infective endocarditis and 30 blood culture isolates. Overall there were great variations in in vitro adherence, but no differences were observed between carriage and infection strains. We further determined the relation between in vitro adherence and in vivo infectivity in a rat model of experimental endocarditis, using 4 isolates that displayed either extremely low or high adherence phenotypes. Unexpectedly, no differences were observed between the in vivo infectivity of isolates that were poorly and highly adherent in vitro. We concluded that the natural variability of in vitro adherence to fibrinogen and fibronectin did not correlate with in vivo infectivity, and thus that pathogenic differences between various strains might only be expressed in in vivo conditions, but not in vitro. Therefore, considering the importance of adhesins expression for infection, direct measurement of those adhesins present on the bacterial surface were made by proteomic approach. 5 In the second series of experiments we assessed the physical presence of the LPXTG species at the staphylococcal surface, as measured at various time points during growth in different culture media. S. aureus Newman was grown in either tryptic soy broth (TSB) or in Roswell Park Memorial Institute (RPMI) culture medium, and samples were removed from early exponential growth phase to late stationary phase. Experiments were performed with mutants in the global accessory-gene regulator (agr), surface protein A (Spa) and clumping factor A (ClfA). Peptides of surface proteins were recovered by "trypsin-shaving" of live bacteria, and semi-quantitative proteomic analysis was performed by tandem liquid-chromatography and mass-spectrometry (LC-MS). We also determined in parallel the mRNA expression by microarrays analysis, as well as the phenotypic adherence of the bacteria to fibrinogen in vitro. The surface proteome was highly complex and contained numerous proteins theoretically not belonging to the bacterial envelope, including ribosomal proteins and metabolic enzymes. Sixteen of the 21 known LPXTG species were detected, but were differentially expressed. As expected, 9 known agr-regulated proteins (e.g. including Spa, FnBPA, ClfA, IsdA, IsdB, SasH, SasD, SasG and FmtB) increased up to the late exponential growth phase, and were abrogated in agr-negative mutants. However, only Spa and SasH modified their proteomic and mRNA profiles in parallel in the parent and its agr negative mutant, while all other LPXTG proteins modified their proteomic profiles independently of their mRNA. Moreover, ClfA became highly transcribed and active in in vitro fibrinogen adherence tests during late growth (24h), whereas it remained poorly detected by proteomics. Differential expression was also detected in iron-rich TSB versus iron-poor RPMI. Proteins from the iron-regulated surface determinant (isd) system, including IsdA, IsdB and IsdH were barely expressed in iron-rich TSB, whereas they increased their expression by >10 time in iron-poor RPMI. We conclude that semi-quantitative proteomic analysis of specific protein species is feasible in S. aureus and that proteomic, transcriptomic and adherence phenotypes demonstrated differential profiles in S. aureus. Furthermore, peptide signatures released by trypsin shaving suggested differential protein domain exposures in various environments, which might be relevant for antiadhesins vaccines. A comprehensive understanding of the S. aureus physiology should integrate all these approaches.

Groundwater Availability Modeling: Lower Dakota Aquifer, 2008

This monthly report from the Iowa Department of Natural Resources is about the water quality management of Iowa's rivers, streams and lakes.

Molecular modeling study of the binding and signaling properties of the TCRpMHC complex

The activation of the specific immune response against tumor cells is based on the recognition by the CD8+ Cytotoxic Τ Lymphocytes (CTL), of antigenic peptides (p) presented at the surface of the cell by the class I major histocompatibility complex (MHC). The ability of the so-called T-Cell Receptors (TCR) to discriminate between self and non-self peptides constitutes the most important specific control mechanism against infected cells. The TCR/pMHC interaction has been the subject of much attention in cancer therapy since the design of the adoptive transfer approach, in which Τ lymphocytes presenting an interesting response against tumor cells are extracted from the patient, expanded in vitro, and reinfused after immunodepletion, possibly leading to cancer regression. In the last decade, major progress has been achieved by the introduction of engineered lypmhocytes. In the meantime, the understanding of the molecular aspects of the TCRpMHC interaction has become essential to guide in vitro and in vivo studies. In 1996, the determination of the first structure of a TCRpMHC complex by X-ray crystallography revealed the molecular basis of the interaction. Since then, molecular modeling techniques have taken advantage of crystal structures to study the conformational space of the complex, and understand the specificity of the recognition of the pMHC by the TCR. In the meantime, experimental techniques used to determine the sequences of TCR that bind to a pMHC complex have been used intensively, leading to the collection of large repertoires of TCR sequences that are specific for a given pMHC. There is a growing need for computational approaches capable of predicting the molecular interactions that occur upon TCR/pMHC binding without relying on the time consuming resolution of a crystal structure. This work presents new approaches to analyze the molecular principles that govern the recognition of the pMHC by the TCR and the subsequent activation of the T-cell. We first introduce TCRep 3D, a new method to model and study the structural properties of TCR repertoires, based on homology and ab initio modeling. We discuss the methodology in details, and demonstrate that it outperforms state of the art modeling methods in predicting relevant TCR conformations. Two successful applications of TCRep 3D that supported experimental studies on TCR repertoires are presented. Second, we present a rigid body study of TCRpMHC complexes that gives a fair insight on the TCR approach towards pMHC. We show that the binding mode of the TCR is correctly described by long-distance interactions. Finally, the last section is dedicated to a detailed analysis of an experimental hydrogen exchange study, which suggests that some regions of the constant domain of the TCR are subject to conformational changes upon binding to the pMHC. We propose a hypothesis of the structural signaling of TCR molecules leading to the activation of the T-cell. It is based on the analysis of correlated motions in the TCRpMHC structure. - L'activation de la réponse immunitaire spécifique dirigée contre les cellules tumorales est basée sur la reconnaissance par les Lymphocytes Τ Cytotoxiques (CTL), d'un peptide antigénique (p) présenté à la suface de la cellule par le complexe majeur d'histocompatibilité de classe I (MHC). La capacité des récepteurs des lymphocytes (TCR) à distinguer les peptides endogènes des peptides étrangers constitue le mécanisme de contrôle le plus important dirigé contre les cellules infectées. L'interaction entre le TCR et le pMHC est le sujet de beaucoup d'attention dans la thérapie du cancer, depuis la conception de la méthode de transfer adoptif: les lymphocytes capables d'une réponse importante contre les cellules tumorales sont extraits du patient, amplifiés in vitro, et réintroduits après immunosuppression. Il peut en résulter une régression du cancer. Ces dix dernières années, d'importants progrès ont été réalisés grâce à l'introduction de lymphocytes modifiés par génie génétique. En parallèle, la compréhension du TCRpMHC au niveau moléculaire est donc devenue essentielle pour soutenir les études in vitro et in vivo. En 1996, l'obtention de la première structure du complexe TCRpMHC à l'aide de la cristallographie par rayons X a révélé les bases moléculaires de l'interaction. Depuis lors, les techniques de modélisation moléculaire ont exploité les structures expérimentales pour comprendre la spécificité de la reconnaissance du pMHC par le TCR. Dans le même temps, de nouvelles techniques expérimentales permettant de déterminer la séquence de TCR spécifiques envers un pMHC donné, ont été largement exploitées. Ainsi, d'importants répertoires de TCR sont devenus disponibles, et il est plus que jamais nécessaire de développer des approches informatiques capables de prédire les interactions moléculaires qui ont lieu lors de la liaison du TCR au pMHC, et ce sans dépendre systématiquement de la résolution d'une structure cristalline. Ce mémoire présente une nouvelle approche pour analyser les principes moléculaires régissant la reconnaissance du pMHC par le TCR, et l'activation du lymphocyte qui en résulte. Dans un premier temps, nous présentons TCRep 3D, une nouvelle méthode basée sur les modélisations par homologie et ab initio, pour l'étude de propriétés structurales des répertoires de TCR. Le procédé est discuté en détails et comparé à des approches standard. Nous démontrons ainsi que TCRep 3D est le plus performant pour prédire des conformations pertinentes du TCR. Deux applications à des études expérimentales des répertoires TCR sont ensuite présentées. Dans la seconde partie de ce travail nous présentons une étude de complexes TCRpMHC qui donne un aperçu intéressant du mécanisme d'approche du pMHC par le TCR. Finalement, la dernière section se concentre sur l'analyse détaillée d'une étude expérimentale basée sur les échanges deuterium/hydrogène, dont les résultats révèlent que certaines régions clés du domaine constant du TCR sont sujettes à un changement conformationnel lors de la liaison au pMHC. Nous proposons une hypothèse pour la signalisation structurelle des TCR, menant à l'activation du lymphocyte. Celle-ci est basée sur l'analyse des mouvements corrélés observés dans la structure du TCRpMHC.

Intravitreal ranibizumab for age-related macular degeneration : optical coherence tomography guided retreatment intervals and their intra- and inter-individual variability

Purpose: To determine whether the need for retreatment after an initial phase of 3 monthly intravitreal injections of ranibizumab shows an intra-individual regular rhythm and to what degree it varies between different patients. Methods: Prospective study with 42 patients with exudative AMD, treatment naïve. Loading dose of 3 monthly doses of ranibizumab (0,5 mg), followed by a 12 months pro re nata (PRN) regimen according to early exudative signs on HD-OCT Cirrus, Zeiss. The follow-up visits were intensified (week 4, 5, 6, 7, 8, 10, 12, 14, 16, 20, etc after each injection) in order to detect recurrences early, and injection followed within 3 days in cases of subretinal fluid, cysts, or central thickness increase of>50microns. Intervals were calculated between injections for the 12 month follow-up with PRN treatment. Variability was expressed as standard deviation (SD). Results: Visual acuity (VA) improved from a mean ETDRS score of 61.6 (SD 10.8) at baseline to 68.0 (SD 10.2) at month 3 and to 74.7(SD 9.0) at month 12. The 15 patients who have already completed the study showed maintenance of the VA improvement. Central foveal thickness improved from a mean value of 366 microns (baseline) to 253 microns (month 3), well maintained thereafter. Mean number of injections was 8.8 (SD 3.5,range 0-12) per 12 months of follow-up (after 3 doses), with mean individual treatment-recurrence (TR) intervals ranging from 28->365 days (mean 58). Intraindividual variability of TR intervals (SD) was 7.1 days as a mean value (range 1.7¡V22.6). It ranged within 20% of the mean intra-individual interval for 30 (91%) and within 15% for 21 patients (64%). The first interval was within 1 week of the mean intra-individual interval in 64% and within 2 weeks in 89% of patients. Conclusions: The majority of AMD patients showed a relatively stable rhythm for PRN injections of ranibizumab after initial loading phase, associated with excellent functional/anatomical results. The initial interval last loading dose-first recurrence may have a predictive value for further need of treatment, potentially facilitating follow-up and patient care.

On the relevance of modeling volatility for pricing purposes

This paper presents a two-factor (Vasicek-CIR) model of the term structure of interest rates and develops its pricing and empirical properties. We assume that default free discount bond prices are determined by the time to maturity and two factors, the long-term interest rate and the spread. Assuming a certain process for both factors, a general bond pricing equation is derived and a closed-form expression for bond prices is obtained. Empirical evidence of the model's performance in comparisson with a double Vasicek model is presented. The main conclusion is that the modeling of the volatility in the long-term rate process can help (in a large amount) to fit the observed data can improve - in a reasonable quantity - the prediction of the future movements in the medium- and long-term interest rates. However, for shorter maturities, it is shown that the pricing errors are, basically, negligible and it is not so clear which is the best model to be used.

Aplicação dos Sistemas de Informação Geográfica em Estudos de Caracterização Biofísica da Ilha de Santo Antão em Cabo Verde

Com características morfológicas e edafo-climáticas extremamente diversificadas, a ilha de Santo Antão em Cabo Verde apresenta uma reconhecida vulnerabilidade ambiental a par de uma elevada carência de estudos científicos que incidam sobre essa realidade e sirvam de base à uma compreensão integrada dos fenómenos. A cartografia digital e as tecnologias de informação geográfica vêm proporcionando um avanço tecnológico na colecção, armazenamento e processamento de dados espaciais. Várias ferramentas actualmente disponíveis permitem modelar uma multiplicidade de factores, localizar e quantificar os fenómenos bem como e definir os níveis de contribuição de diferentes factores no resultado final. No presente estudo, desenvolvido no âmbito do curso de pós-graduação e mestrado em sistemas de Informação geográfica realizado pela Universidade de Trás-os-Montes e Alto Douro, pretende-se contribuir para a minimização do deficit de informação relativa às características biofísicas da citada ilha, recorrendo-se à aplicação de tecnologias de informação geográfica e detecção remota, associadas à análise estatística multivariada. Nesse âmbito, foram produzidas e analisadas cartas temáticas e desenvolvido um modelo de análise integrada de dados. Com efeito, a multiplicidade de variáveis espaciais produzidas, de entre elas 29 variáveis com variação contínua passíveis de influenciar as características biofísicas da região e, possíveis ocorrências de efeitos mútuos antagónicos ou sinergéticos, condicionam uma relativa complexidade à interpretação a partir dos dados originais. Visando contornar este problema, recorre-se a uma rede de amostragem sistemática, totalizando 921 pontos ou repetições, para extrair os dados correspondentes às 29 variáveis nos pontos de amostragem e, subsequente desenvolvimento de técnicas de análise estatística multivariada, nomeadamente a análise em componentes principais. A aplicação destas técnicas permitiu simplificar e interpretar as variáreis originais, normalizando-as e resumindo a informação contida na diversidade de variáveis originais, correlacionadas entre si, num conjunto de variáveis ortogonais (não correlacionadas), e com níveis de importância decrescente, as componentes principais. Fixou-se como meta a concentração de 75% da variância dos dados originais explicadas pelas primeiras 3 componentes principais e, desenvolveu-se um processo interactivo em diferentes etapas, eliminando sucessivamente as variáveis menos representativas. Na última etapa do processo as 3 primeiras CP resultaram em 74,54% da variância dos dados originais explicadas mas, que vieram a demonstrar na fase posterior, serem insuficientes para retratar a realidade. Optou-se pela inclusão da 4ª CP (CP4), com a qual 84% da referida variância era explicada e, representando oito variáveis biofísicas: a altitude, a densidade hidrográfica, a densidade de fracturação geológica, a precipitação, o índice de vegetação, a temperatura, os recursos hídricos e a distância à rede hidrográfica. A subsequente interpolação da 1ª componente principal (CP1) e, das principais variáveis associadas as componentes CP2, CP3 e CP4 como variáveis auxiliares, recorrendo a técnicas geoestatística em ambiente ArcGIS permitiu a obtenção de uma carta representando 84% da variação das características biofísicas no território. A análise em clusters validada pelo teste “t de Student” permitiu reclassificar o território em 6 unidades biofísicas homogéneas. Conclui-se que, as tecnologias de informação geográfica actualmente disponíveis a par de facilitar análises interactivas e flexíveis, possibilitando que se faça variar temas e critérios, integrar novas informações e introduzir melhorias em modelos construídos com bases em informações disponíveis num determinado contexto, associadas a técnicas de análise estatística multivariada, possibilitam, com base em critérios científicos, desenvolver a análise integrada de múltiplas variáveis biofísicas cuja correlação entre si, torna complexa a compreensão integrada dos fenómenos.

Pitfalls in the modeling of forward-looking price setting and investment decisions

The present paper makes progress in explaining the role of capital for inflation and output dynamics. We followWoodford (2003, Ch. 5) in assuming Calvo pricing combined with a convex capital adjustment cost at the firm level. Our main result is that capital accumulation affects inflation dynamics primarily through its impact on the marginal cost. This mechanism is much simpler than the one implied by the analysis in Woodford's text. The reason is that his analysis suffers from a conceptual mistake, as we show. The latter obscures the economic mechanism through which capital affects inflation and output dynamics in the Calvo model, as discussed in Woodford (2004).

Flexible multivariate GARCH modeling with an application to international stock markets

The goal of this paper is to estimate time-varying covariance matrices.Since the covariance matrix of financial returns is known to changethrough time and is an essential ingredient in risk measurement, portfolioselection, and tests of asset pricing models, this is a very importantproblem in practice. Our model of choice is the Diagonal-Vech version ofthe Multivariate GARCH(1,1) model. The problem is that the estimation ofthe general Diagonal-Vech model model is numerically infeasible indimensions higher than 5. The common approach is to estimate more restrictive models which are tractable but may not conform to the data. Our contributionis to propose an alternative estimation method that is numerically feasible,produces positive semi-definite conditional covariance matrices, and doesnot impose unrealistic a priori restrictions. We provide an empiricalapplication in the context of international stock markets, comparing thenew estimator to a number of existing ones.

Macroeconomic modeling for monetary policy evaluation

We describe some of the main features of the recent vintage macroeconomic models used for monetary policy evaluation. We point to some of the key differences with respect to the earlier generation ofmacro models, and highlight the insights for policy that these new frameworks have to offer. Our discussion emphasizes two key aspects of the new models: the significant role of expectations of future policy actions in the monetary transmission mechanism, and the importance for the central bank of tracking of the flexible price equilibrium values of the natural levels of output and the real interest rate. We argue that both features have important implications for the conduct of monetary policy.