Probabilistic physically-based cloud screening of satellite infra-red imagery for operational sea surface temperature retrieval

We propose and demonstrate a fully probabilistic (Bayesian) approach to the detection of cloudy pixels in thermal infrared (TIR) imagery observed from satellite over oceans. Using this approach, we show how to exploit the prior information and the fast forward modelling capability that are typically available in the operational context to obtain improved cloud detection. The probability of clear sky for each pixel is estimated by applying Bayes' theorem, and we describe how to apply Bayes' theorem to this problem in general terms. Joint probability density functions (PDFs) of the observations in the TIR channels are needed; the PDFs for clear conditions are calculable from forward modelling and those for cloudy conditions have been obtained empirically. Using analysis fields from numerical weather prediction as prior information, we apply the approach to imagery representative of imagers on polar-orbiting platforms. In comparison with the established cloud-screening scheme, the new technique decreases both the rate of failure to detect cloud contamination and the false-alarm rate by one quarter. The rate of occurrence of cloud-screening-related errors of >1 K in area-averaged SSTs is reduced by 83%. Copyright © 2005 Royal Meteorological Society.

Representation errors and retrievals in linear and nonlinear data assimilation

This article shows how one can formulate the representation problem starting from Bayes’ theorem. The purpose of this article is to raise awareness of the formal solutions,so that approximations can be placed in a proper context. The representation errors appear in the likelihood, and the different possibilities for the representation of reality in model and observations are discussed, including nonlinear representation probability density functions. Specifically, the assumptions needed in the usual procedure to add a representation error covariance to the error covariance of the observations are discussed,and it is shown that, when several sub-grid observations are present, their mean still has a representation error ; socalled ‘superobbing’ does not resolve the issue. Connection is made to the off-line or on-line retrieval problem, providing a new simple proof of the equivalence of assimilating linear retrievals and original observations. Furthermore, it is shown how nonlinear retrievals can be assimilated without loss of information. Finally we discuss how errors in the observation operator model can be treated consistently in the Bayesian framework, connecting to previous work in this area.

Methods of Estimation in Multiple Linear Regression: Application to Clinical Data

Nesse artigo, tem-se o interesse em avaliar diferentes estratégias de estimação de parâmetros para um modelo de regressão linear múltipla. Para a estimação dos parâmetros do modelo foram utilizados dados de um ensaio clínico em que o interesse foi verificar se o ensaio mecânico da propriedade de força máxima (EM-FM) está associada com a massa femoral, com o diâmetro femoral e com o grupo experimental de ratas ovariectomizadas da raça Rattus norvegicus albinus, variedade Wistar. Para a estimação dos parâmetros do modelo serão comparadas três metodologias: a metodologia clássica, baseada no método dos mínimos quadrados; a metodologia Bayesiana, baseada no teorema de Bayes; e o método Bootstrap, baseado em processos de reamostragem.

Estudo de caso-controle para desenvolvimento e validação de um questionário de avaliação de abuso sexual em crianças

Objetivos: Desenvolver e validar instrumento que auxilie o pediatra a determinar a probabilidade de ocorrência do abuso sexual em crianças. Métodos: Estudo de caso-controle com 201 crianças que consultaram em ambulatórios de pediatria e locais de referência para vítimas de abuso sexual, entre março e novembro de 2004: grupo caso (com suspeita ou revelação de abuso sexual) e grupo controle (sem suspeita de abuso sexual). Aplicou-se, junto aos responsáveis, um questionário com 18 itens e cinco opções de respostas segundo a escala Likert, abordando comportamento, sintomas físicos e emocionais apresentados pelas crianças. Excluíram-se nove crianças sem controle esfincteriano e um item respondido por poucas pessoas. A validade e consistência interna dos itens foram avaliadas com obtenção de coeficientes de correlação (Pearson, Spearman e Goodman-Kruskal), coeficiente α de Cronbach e cálculo da área da curva ROC. Calculou-se, após, a razão de verossimilhança (RV) e os valores preditivo positivos (VPP) para os cinco itens do questionário que apresentaram os melhores desempenhos. Resultados: Obteve-se um questionário composto pelos cinco itens que melhor discriminaram crianças com e sem abuso sexual em dois contextos. Cada criança recebeu um escore resultante da soma das respostas com pesos de 0 a 4 (amplitude de 0 a 20), o qual, através do teorema de Bayes (RV), indicou sua probabilidade pós-teste (VPP) de abuso sexual. Conclusões: O instrumento proposto é útil por ser de fácil aplicação, auxiliando o pediatra na identificação de crianças vítimas de abuso sexual. Ele fornecerá, conforme o escore obtido, a probabilidade (VPP) de abuso sexual, orientando na conduta de cuidado à criança.

Development of an open case-based decision-support system for diagnosis in oral pathology

Making diagnoses in oral pathology are often difficult and confusing in dental practice, especially for the lessexperienced dental student. One of the most promising areas in bioinformatics is computer-aided diagnosis, where a computer system is capable of imitating human reasoning ability and provides diagnoses with an accuracy approaching that of expert professionals. This type of system could be an alternative tool for assisting dental students to overcome the difficulties of the oral pathology learning process. This could allow students to define variables and information, important to improving the decision-making performance. However, no current open data management system has been integrated with an artificial intelligence system in a user-friendly environment. Such a system could also be used as an education tool to help students perform diagnoses. The aim of the present study was to develop and test an open case-based decisionsupport system.Methods: An open decision-support system based on Bayes' theorem connected to a relational database was developed using the C++ programming language. The software was tested in the computerisation of a surgical pathology service and in simulating the diagnosis of 43 known cases of oral bone disease. The simulation was performed after the system was initially filled with data from 401 cases of oral bone disease.Results: the system allowed the authors to construct and to manage a pathology database, and to simulate diagnoses using the variables from the database.Conclusion: Combining a relational database and an open decision-support system in the same user-friendly environment proved effective in simulating diagnoses based on information from an updated database.

A Bayesian approach for constructing genetic maps when markers are miscoded

The advent of molecular markers has created opportunities for a better understanding of quantitative inheritance and for developing novel strategies for genetic improvement of agricultural species, using information on quantitative trait loci (QTL). A QTL analysis relies on accurate genetic marker maps. At present, most statistical methods used for map construction ignore the fact that molecular data may be read with error. Often, however, there is ambiguity about some marker genotypes. A Bayesian MCMC approach for inferences about a genetic marker map when random miscoding of genotypes occurs is presented, and simulated and real data sets are analyzed. The results suggest that unless there is strong reason to believe that genotypes are ascertained without error, the proposed approach provides more reliable inference on the genetic map.

A Bayesian approach to estimate the accuracy of in-house ELISA assay to measure rabies antibodies from compulsory vaccinated dogs and cattle

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Molecular phylogenetic relationships and phenotypic diversity in miniaturized toadlets, genus Brachycephalus (Amphibia: Anura: Brachycephalidae)

Toadlets of the genus Brachycephalus are endemic to the Atlantic rainforests of southeastern and southern Brazil. The 14 species currently described have snout-vent lengths less than 18. mm and are thought to have evolved through miniaturization: an evolutionary process leading to an extremely small adult body size. Here, we present the first comprehensive phylogenetic analysis for Brachycephalus, using a multilocus approach based on two nuclear (Rag-1 and Tyr) and three mitochondrial (Cyt b, 12S, and 16S rRNA) gene regions. Phylogenetic relationships were inferred using a partitioned Bayesian analysis of concatenated sequences and the hierarchical Bayesian method (BEST) that estimates species trees based on the multispecies coalescent model. Individual gene trees showed conflict and also varied in resolution. With the exception of the mitochondrial gene tree, no gene tree was completely resolved. The concatenated gene tree was completely resolved and is identical in topology and degree of statistical support to the individual mtDNA gene tree. On the other hand, the BEST species tree showed reduced significant node support relative to the concatenate tree and recovered a basal trichotomy, although some bipartitions were significantly supported at the tips of the species tree. Comparison of the log likelihoods for the concatenated and BEST trees suggests that the method implemented in BEST explains the multilocus data for Brachycephalus better than the Bayesian analysis of concatenated data. Landmark-based geometric morphometrics revealed marked variation in cranial shape between the species of Brachycephalus. In addition, a statistically significant association was demonstrated between variation in cranial shape and genetic distances estimated from the mtDNA and nuclear loci. Notably, B. ephippium and B. garbeana that are predicted to be sister-species in the individual and concatenated gene trees and the BEST species tree share an evolutionary novelty, the hyperossified dorsal plate. © 2011 Elsevier Inc.

Random regression models using Legendre polynomials or linear splines for test-day milk yield of dairy Gyr (Bos indicus) cattle

Studies investigating the use of random regression models for genetic evaluation of milk production in Zebu cattle are scarce. In this study, 59,744 test-day milk yield records from 7,810 first lactations of purebred dairy Gyr (Bos indicus) and crossbred (dairy Gyr × Holstein) cows were used to compare random regression models in which additive genetic and permanent environmental effects were modeled using orthogonal Legendre polynomials or linear spline functions. Residual variances were modeled considering 1, 5, or 10 classes of days in milk. Five classes fitted the changes in residual variances over the lactation adequately and were used for model comparison. The model that fitted linear spline functions with 6 knots provided the lowest sum of residual variances across lactation. On the other hand, according to the deviance information criterion (DIC) and Bayesian information criterion (BIC), a model using third-order and fourth-order Legendre polynomials for additive genetic and permanent environmental effects, respectively, provided the best fit. However, the high rank correlation (0.998) between this model and that applying third-order Legendre polynomials for additive genetic and permanent environmental effects, indicates that, in practice, the same bulls would be selected by both models. The last model, which is less parameterized, is a parsimonious option for fitting dairy Gyr breed test-day milk yield records. © 2013 American Dairy Science Association.

Analysis of genetic correlations of hip height with selection indices and mature weight in Nelore cattle

Body size is directly related to the productive and reproductive performance of beef cattle raised under free-range conditions. In an attempt to better plan selection criteria, avoiding extremes in body size, this study estimated the heritabilities and genetic correlations of yearling hip height (YH) and mature hip height (MH) with selection indices obtained at weaning (WI) and yearling (YI) and mature weight (MW). Data from 102,373 Nelore animals born between 1984 and 2010, which belong to 263 farms that participate in genetic evaluation programmes of beef cattle conducted in Brazil and Paraguay, were used. The (co)variance components and genetic parameters were estimated by Bayesian inference in multi-trait analysis using an animal model. The mean heritabilities for YH, MH and MW were 0. 56 ± 0. 06, 0. 47 ± 0. 02 and 0. 42 ± 0. 02, respectively. The genetic correlation of YH with WI (0. 13 ± 0. 01) and YI (0. 11 ± 0. 01) was practically zero, whereas a higher correlation was observed with MW (0. 22 ± 0. 03). Positive genetic correlations of medium magnitude were estimated between MH and WI and YI (0. 23 ± 0. 01 and 0. 43 ± 0. 02, respectively). On the other hand, a high genetic correlation (0. 68 ± 0. 03) was observed between the indicator traits of mature body size (MH and MW). Considering the top 20 % of sire (896 sires) in terms of breeding values for the yearling index, the rank sire correlations between breeding values for MH and MW was 0. 62. In general, the results indicate that selection based on WI and YI should not lead to important changes in YH. However, an undesired correlated response in mature cow height is expected, particularly when selection is performed using YI. Therefore, changes in the body structure of Nelore females can be obtained when MH and MW is used as a selection criterion for cows. © 2012 Institute of Plant Genetics, Polish Academy of Sciences, Poznan.

Comparative cytogenetic analysis of some species of the Dendropsophus microcephalus group (Anura, Hylidae) in the light of phylogenetic inferences

Background: Dendropsophus is a monophyletic anuran genus with a diploid number of 30 chromosomes as an important synapomorphy. However, the internal phylogenetic relationships of this genus are poorly understood. Interestingly, an intriguing interspecific variation in the telocentric chromosome number has been useful in species identification. To address certain uncertainties related to one of the species groups of Dendropsophus, the D. microcephalus group, we carried out a cytogenetic analysis combined with phylogenetic inferences based on mitochondrial sequences, which aimed to aid in the analysis of chromosomal characters. Populations of Dendropsophus nanus, Dendropsophus walfordi, Dendropsophus sanborni, Dendropsophus jimi and Dendropsophus elianeae, ranging from the extreme south to the north of Brazil, were cytogenetically compared. A mitochondrial region of the ribosomal 12S gene from these populations, as well as from 30 other species of Dendropsophus, was used for the phylogenetic inferences. Phylogenetic relationships were inferred using maximum parsimony and Bayesian analyses.Results: The species D. nanus and D. walfordi exhibited identical karyotypes (2n = 30; FN = 52), with four pairs of telocentric chromosomes and a NOR located on metacentric chromosome pair 13. In all of the phylogenetic hypotheses, the paraphyly of D. nanus and D. walfordi was inferred. D. sanborni from Botucatu-SP and Torres-RS showed the same karyotype as D. jimi, with 5 pairs of telocentric chromosomes (2n = 30; FN = 50) and a terminal NOR in the long arm of the telocentric chromosome pair 12. Despite their karyotypic similarity, these species were not found to compose a monophyletic group. Finally, the phylogenetic and cytogenetic analyses did not cluster the specimens of D. elianeae according to their geographical occurrence or recognized morphotypes.Conclusions: We suggest that a taxonomic revision of the taxa D. nanus and D. walfordi is quite necessary. We also observe that the number of telocentric chromosomes is useful to distinguish among valid species in some cases, although it is unchanged in species that are not necessarily closely related phylogenetically. Therefore, inferences based on this chromosomal character must be made with caution; a proper evolutionary analysis of the karyotypic variation in Dendropsophus depends on further characterization of the telocentric chromosomes found in this group. © 2013 Medeiros et al.; licensee BioMed Central Ltd.

Microsatellite marker-based assessment of the biodiversity of native bioethanol yeast strains

Although many Brazilian sugar mills initiate the fermentation process by inoculating selected commercial Saccharomyces cerevisiae strains, the unsterile conditions of the industrial sugar cane ethanol fermentation process permit the constant entry of native yeast strains. Certain of those native strains are better adapted and tend to predominate over the initial strain, which may cause problems during fermentation. In the industrial fermentation process, yeast cells are often exposed to stressful environmental conditions, including prolonged cell recycling, ethanol toxicity and osmotic, oxidative or temperature stress. Little is known about these S. cerevisiae strains, although recent studies have demonstrated that heterogeneous genome architecture is exhibited by some selected well-adapted Brazilian indigenous yeast strains that display high performance in bioethanol fermentation. In this study, 11 microsatellite markers were used to assess the genetic diversity and population structure of the native autochthonous S. cerevisiae strains in various Brazilian sugar mills. The resulting multilocus data were used to build a similarity-based phenetic tree and to perform a Bayesian population structure analysis. The tree revealed the presence of great genetic diversity among the strains, which were arranged according to the place of origin and the collection year. The population structure analysis revealed genotypic differences among populations; in certain populations, these genotypic differences are combined to yield notably genotypically diverse individuals. The high yeast diversity observed among native S. cerevisiae strains provides new insights on the use of autochthonous high-fitness strains with industrial characteristics as starter cultures at bioethanol plants. © 2013 John Wiley & Sons, Ltd.

Clinical judgment to estimate pretest probability in the diagnosis of Cushing's syndrome under a Bayesian perspective

OBJECTIVE: To estimate the pretest probability of Cushing's syndrome (CS) diagnosis by a Bayesian approach using intuitive clinical judgment. MATERIALS AND METHODS: Physicians were requested, in seven endocrinology meetings, to answer three questions: "Based on your personal expertise, after obtaining clinical history and physical examination, without using laboratorial tests, what is your probability of diagnosing Cushing's Syndrome?"; "For how long have you been practicing Endocrinology?"; and "Where do you work?". A Bayesian beta regression, using the WinBugs software was employed. RESULTS: We obtained 294 questionnaires. The mean pretest probability of CS diagnosis was 51.6% (95%CI: 48.7-54.3). The probability was directly related to experience in endocrinology, but not with the place of work. CONCLUSION: Pretest probability of CS diagnosis was estimated using a Bayesian methodology. Although pretest likelihood can be context-dependent, experience based on years of practice may help the practitioner to diagnosis CS. Arq Bras Endocrinol Metab. 2012;56(9):633-7

Der g-Faktor des Protons

In dieser Arbeit wird die bisher präziseste und erste direkte Hochpräzisionsmessung des g-Faktors eines einzelnen Protons präsentiert. Die Messung beruht auf der nicht-destruktiven Bestimmung der Zyklotronfrequenz und der Larmorfrequenz eines in einer Penning-Falle gespeicherten Protons. Zur Bestimmung der Larmorfrequenz wird die Spin-Flip-Wahrscheinlichkeit als Funktion einer externen Spin-Flip-Anregung aufgenommen. Zu diesem Zweck wird der kontinuierliche Stern-Gerlach Effekt verwendet, welcher zu einer Kopplung des Spin-Moments an die axiale Bewegung des Protons führt. Ein Spin-Flip zeigt sich dabei in einem Sprung der axialen Bewegungsfrequenz. Die Schwierigkeit besteht darin, diesen Frequenzsprung auf einem Hintergrund axialer Frequenzfluktuationen zu detektieren. Um diese Herausforderung zu bewältigen, wurden neuartige Methoden und Techniken angewandt. Zum einen wurden supraleitende Nachweise mit höchster Empfindlichkeit entwickelt, welche schnelle und damit präzise Frequenzmessungen erlauben. Zum anderen wurde eine auf dem statistischen Bayes Theorem basierende Spin-Flip-Analyse-Methode angewandt. Mit diesen Verbesserungen war es möglich, einzelne Spin-Flips eines einzelnen Protons zu beobachten. Dies wiederum ermöglichte die Anwendung der sogenannten Doppelfallen-Methode, und damit die eingangs erwähnte Messung des g-Faktors mit einer Präzision von 4.3 10^-9.