Real-time GPU color-based segmentation of football players

In this paper, we propose a multi-camera application capable of processing high resolution images and extracting features based on colors patterns over graphic processing units (GPU). The goal is to work in real time under the uncontrolled environment of a sport event like a football match. Since football players are composed for diverse and complex color patterns, a Gaussian Mixture Models (GMM) is applied as segmentation paradigm, in order to analyze sport live images and video. Optimization techniques have also been applied over the C++ implementation using profiling tools focused on high performance. Time consuming tasks were implemented over NVIDIA's CUDA platform, and later restructured and enhanced, speeding up the whole process significantly. Our resulting code is around 4-11 times faster on a low cost GPU than a highly optimized C++ version on a central processing unit (CPU) over the same data. Real time has been obtained processing until 64 frames per second. An important conclusion derived from our study is the scalability of the application to the number of cores on the GPU. © 2011 Springer-Verlag.

Stated choices and benefit estimates in the context of traffic calming schemes:Utility maximization, regret minimization, or both?

This paper proposes a discrete mixture model which assigns individuals, up to a probability, to either a class of random utility (RU) maximizers or a class of random regret (RR) minimizers, on the basis of their sequence of observed choices. Our proposed model advances the state of the art of RU-RR mixture models by (i) adding and simultaneously estimating a membership model which predicts the probability of belonging to a RU or RR class; (ii) adding a layer of random taste heterogeneity within each behavioural class; and (iii) deriving a welfare measure associated with the RU-RR mixture model and consistent with referendum-voting, which is the adequate mechanism of provision for such local public goods. The context of our empirical application is a stated choice experiment concerning traffic calming schemes. We find that the random parameter RU-RR mixture model not only outperforms its fixed coefficient counterpart in terms of fit-as expected-but also in terms of plausibility of membership determinants of behavioural class. In line with psychological theories of regret, we find that, compared to respondents who are familiar with the choice context (i.e. the traffic calming scheme), unfamiliar respondents are more likely to be regret minimizers than utility maximizers. © 2014 Elsevier Ltd.

DNA methylation subgroups and the CpG island methylator phenotype in gastric cancer: a comprehensive profiling approach

BACKGROUND: Methylation-induced silencing of promoter CpG islands in tumor suppressor genes plays an important role in human carcinogenesis. In colorectal cancer, the CpG island methylator phenotype (CIMP) is defined as widespread and elevated levels of DNA methylation and CIMP+ tumors have distinctive clinicopathological and molecular features. In contrast, the existence of a comparable CIMP subtype in gastric cancer (GC) has not been clearly established. To further investigate this issue, in the present study we performed comprehensive DNA methylation profiling of a well-characterised series of primary GC.

METHODS: The methylation status of 1,421 autosomal CpG sites located within 768 cancer-related genes was investigated using the Illumina GoldenGate Methylation Panel I assay on DNA extracted from 60 gastric tumors and matched tumor-adjacent gastric tissue pairs. Methylation data was analysed using a recursively partitioned mixture model and investigated for associations with clinicopathological and molecular features including age, Helicobacter pylori status, tumor site, patient survival, microsatellite instability and BRAF and KRAS mutations.

RESULTS: A total of 147 genes were differentially methylated between tumor and matched tumor-adjacent gastric tissue, with HOXA5 and hedgehog signalling being the top-ranked gene and signalling pathway, respectively. Unsupervised clustering of methylation data revealed the existence of 6 subgroups under two main clusters, referred to as L (low methylation; 28% of cases) and H (high methylation; 72%). Female patients were over-represented in the H tumor group compared to L group (36% vs 6%; P = 0.024), however no other significant differences in clinicopathological or molecular features were apparent. CpG sites that were hypermethylated in group H were more frequently located in CpG islands and marked for polycomb occupancy.

CONCLUSIONS: High-throughput methylation analysis implicates genes involved in embryonic development and hedgehog signaling in gastric tumorigenesis. GC is comprised of two major methylation subtypes, with the highly methylated group showing some features consistent with a CpG island methylator phenotype.

Métodos probabilísticos e estatísticos na gestão de fogos florestais

Tese de doutoramento, Estatística e Investigação Operacional (Probabilidades e Estatística), Universidade de Lisboa, Faculdade de Ciências, 2014

Modeling radionuclides dispersion and deposition downwind of a coal-fired power plant

In this study the inhalation doses and respective risk are calculated for the population living within a 20 km radius of a coal-fired power plant. The dispersion and deposition of natural radionuclides were simulated by a Gaussian dispersion model estimating the ground level activity concentration. The annual effective dose and total risk were 0.03205 mSv/y and 1.25 x 10-8, respectively. The effective dose is lower than the limit established by the ICRP and the risk is lower than the limit proposed by the U.S. EPA, which means that the considered exposure does not pose any risk for the public health.

Modelling HIV/AIDS length of stay in Portugal

Tese apresentada como requisito parcial para obtenção do grau de Doutor em Estatística e Gestão de Informação pelo Instituto Superior de Estatística e Gestão de Informação da Universidade Nova de Lisboa

MODELING HETEROTACHY IN PHYLOGENETICS

Il a été démontré que l’hétérotachie, variation du taux de substitutions au cours du temps et entre les sites, est un phénomène fréquent au sein de données réelles. Échouer à modéliser l’hétérotachie peut potentiellement causer des artéfacts phylogénétiques. Actuellement, plusieurs modèles traitent l’hétérotachie : le modèle à mélange des longueurs de branche (MLB) ainsi que diverses formes du modèle covarion. Dans ce projet, notre but est de trouver un modèle qui prenne efficacement en compte les signaux hétérotaches présents dans les données, et ainsi améliorer l’inférence phylogénétique. Pour parvenir à nos fins, deux études ont été réalisées. Dans la première, nous comparons le modèle MLB avec le modèle covarion et le modèle homogène grâce aux test AIC et BIC, ainsi que par validation croisée. A partir de nos résultats, nous pouvons conclure que le modèle MLB n’est pas nécessaire pour les sites dont les longueurs de branche diffèrent sur l’ensemble de l’arbre, car, dans les données réelles, le signaux hétérotaches qui interfèrent avec l’inférence phylogénétique sont généralement concentrés dans une zone limitée de l’arbre. Dans la seconde étude, nous relaxons l’hypothèse que le modèle covarion est homogène entre les sites, et développons un modèle à mélanges basé sur un processus de Dirichlet. Afin d’évaluer différents modèles hétérogènes, nous définissons plusieurs tests de non-conformité par échantillonnage postérieur prédictif pour étudier divers aspects de l’évolution moléculaire à partir de cartographies stochastiques. Ces tests montrent que le modèle à mélanges covarion utilisé avec une loi gamma est capable de refléter adéquatement les variations de substitutions tant à l’intérieur d’un site qu’entre les sites. Notre recherche permet de décrire de façon détaillée l’hétérotachie dans des données réelles et donne des pistes à suivre pour de futurs modèles hétérotaches. Les tests de non conformité par échantillonnage postérieur prédictif fournissent des outils de diagnostic pour évaluer les modèles en détails. De plus, nos deux études révèlent la non spécificité des modèles hétérogènes et, en conséquence, la présence d’interactions entre différents modèles hétérogènes. Nos études suggèrent fortement que les données contiennent différents caractères hétérogènes qui devraient être pris en compte simultanément dans les analyses phylogénétiques.

Modèle de mélange de lois multinormales appliqué à l'analyse de comportements et d'habiletés cognitives d'enfants.

Cette étude aborde le thème de l’utilisation des modèles de mélange de lois pour analyser des données de comportements et d’habiletés cognitives mesurées à plusieurs moments au cours du développement des enfants. L’estimation des mélanges de lois multinormales en utilisant l’algorithme EM est expliquée en détail. Cet algorithme simplifie beaucoup les calculs, car il permet d’estimer les paramètres de chaque groupe séparément, permettant ainsi de modéliser plus facilement la covariance des observations à travers le temps. Ce dernier point est souvent mis de côté dans les analyses de mélanges. Cette étude porte sur les conséquences d’une mauvaise spécification de la covariance sur l’estimation du nombre de groupes formant un mélange. La conséquence principale est la surestimation du nombre de groupes, c’est-à-dire qu’on estime des groupes qui n’existent pas. En particulier, l’hypothèse d’indépendance des observations à travers le temps lorsque ces dernières étaient corrélées résultait en l’estimation de plusieurs groupes qui n’existaient pas. Cette surestimation du nombre de groupes entraîne aussi une surparamétrisation, c’est-à-dire qu’on utilise plus de paramètres qu’il n’est nécessaire pour modéliser les données. Finalement, des modèles de mélanges ont été estimés sur des données de comportements et d’habiletés cognitives. Nous avons estimé les mélanges en supposant d’abord une structure de covariance puis l’indépendance. On se rend compte que dans la plupart des cas l’ajout d’une structure de covariance a pour conséquence d’estimer moins de groupes et les résultats sont plus simples et plus clairs à interpréter.

Modélisation statistique de l’érosion de cavitation d’une turbine hydraulique selon les paramètres d’opération

Dans une turbine hydraulique, la rotation des aubes dans l’eau crée une zone de basse pression, amenant l’eau à passer de l’état liquide à l’état gazeux. Ce phénomène de changement de phase est appelé cavitation et est similaire à l’ébullition. Lorsque les cavités de vapeur formées implosent près des parois, il en résulte une érosion sévère des matériaux, accélérant de façon importante la dégradation de la turbine. Un système de détection de l’érosion de cavitation à l’aide de mesures vibratoires, employable sur les turbines en opération, a donc été installé sur quatre groupes turbine-alternateur d’une centrale et permet d’estimer précisément le taux d’érosion en kg/ 10 000 h. Le présent projet vise à répondre à deux objectifs principaux. Premièrement, étudier le comportement de la cavitation sur un groupe turbine-alternateur cible et construire un modèle statistique, dans le but de prédire la variable cavitation en fonction des variables opératoires (tels l’ouverture de vannage, le débit, les niveaux amont et aval, etc.). Deuxièmement, élaborer une méthodologie permettant la reproductibilité de l’étude à d’autres sites. Une étude rétrospective sera effectuée et on se concentrera sur les données disponibles depuis la mise à jour du système en 2010. Des résultats préliminaires ont mis en évidence l’hétérogénéité du comportement de cavitation ainsi que des changements entre la relation entre la cavitation et diverses variables opératoires. Nous nous proposons de développer un modèle probabiliste adapté, en utilisant notamment le regroupement hiérarchique et des modèles de régression linéaire multiple.

Studies on air pollution meteorology with special reference to Madras city

The main aim of the study is to give a clear picture of various meteorological factors affecting the dispersal of pollutants. One such important developing metropolis, namely Madras, is chosen for the present study. The study throws light into the occurrence of inversions, isothermals and lapse conditions and the vertical and horizontal extent of mixing of pollutants. The thesis also aims to study the wind climatology and atmospheric stability .The study gives a insight to the spatial distribution of sulphudioxide concentration using the Gaussian plume model, which accounts for various industrial sources. The researcher suggests optimum locations for industries and various steps to reduce air pollution.

Studies on the meteriological aspects of air pollution over Cochin

Pollutants that once enter into the earth’s atmosphere become part of the atmosphere and hence their dispersion, dilution, direction of transportation etc. are governed by the meteorological conditions. The thesis deals with the study of the atmospheric dispersion capacity, wind climatology, atmospheric stability, pollutant distribution by means of a model and the suggestions for a comprehensive planning for the industrially developing city, Cochin. The definition, sources, types and effects of air pollution have been dealt with briefly. The influence of various meteorological parameters such as vector wind, temperature and its vertical structure and atmospheric stability in relation to pollutant dispersal have been studied. The importance of inversions, mixing heights, ventilation coefficients were brought out. The spatial variation of mixing heights studies for the first time on a microscale region, serves to delineate the regions of good and poor dispersal capacity. A study of wind direction fluctuation, σθ and its relation to stability and mixing heights were shown to be much useful. It was shown that there is a necessity to look into the method of σθ computation. The development of Gausssian Plume Model along with the application for multiple sources was presented. The pollutant chosen was sulphur dioxide and industrial sources alone were considered. The percentage frequency of occurrence of inversions and isothermals are found to be low in all months during the year. The spatial variation of mixing heights revealed that a single mixing height cannot be taken as a representative for the whole city have low mixing heights and monsoonal months showed lowest mixing heights. The study of ventilation co-efficients showed values less than the required optimum value 6000m2/5. However, the low values may be due to the consideration of surface wind alone instead of the vertically averaged wind. Relatively more calm conditions and light winds during night and strong winds during day time were observed. During the most of the year westerlies during day time and northeasterlies during night time are the dominant winds. Unstable conditions with high values of σθ during day time and stable conditions with lower values of σθ during night time are the prominent features. Monsoonal months showed neutral stability for most of the time. A study σθ of and Pasquill Stability category has revealed the difficulty in giving a unique value of for each stability category. For the first time regression equations have been developed relating mixing heights and σθ. A closer examination of σθ revealed that half of the range of wind direction fluctuations is to be taken, instead of one by sixth, to compute σθ. The spatial distribution of SO2 showed a more or less uniform distribution with a slight intrusion towards south. Winter months showed low concentrations contrary to the expectations. The variations of the concentration is found to be influenced more by the mixing height and the stack height rather than wind speed. In the densely populated areas the concentration is more than the threshold limit value. However, the values reported appear to be high, because no depletion of the material is assumed through dry or wet depositions and also because of the inclusion of calm conditions with a very light wind speed. A reduction of emission during night time with a consequent rise during day time would bring down the levels of pollution. The probable locations for the new industries could be the extreme southeast parts because the concentration towards the north falls off very quickly resulting low concentrations. In such a case pollutant spread would be towards south and west, thus keeping the city interior relatively free from pollution. A more detailed examination of the pollutant spread by means of models that would take the dry and wet depositions may be necessary. Nevertheless, the present model serves to give the trend of the distribution of pollutant concentration with which one can suggest the optimum locations for the new industries

Studis on the climatological aspects of air pollution potential over trivandrum

The deteriorating air quality especially in urban environments is a cause of serious concern. In spite of being an effective sink, the atmosphere also has its own limitations in effectively dispersing the pollutants being dumped into it continuously by various sources, mainly industries. Many a time, it is not the higher emissions that cause alarming level of pollutants but the unfavourable atmospheric conditions under which the atmosphere is not able to disperse them effectively, leading to accumulation of pollutants near the ground. Hence, it is imperative to have an estimate of the atmospheric potential for dispersal of the substances emitted into it. This requires a knowledge of mixing height, ventilation coefficient, wind and stability of the region under study. Mere estimation of such pollution potential is not adequate, unless the probable distribution of concentration of pollutants is known. This can be obtained by means of mathematical models. The pollution potential coupled with the distribution of concentration provides a good basis for initiating steps to mitigate air pollution in any developing urban area. In this thesis, a fast developing industrial city, namely, Trivandrum is chosen for estimating the pollution potential and determining the spatial distribution of sulphur dioxide concentration. Each of the parameters required for pollution potential is discussed in detail separately. The thesis is divided into nine chapters.

Organic field-effect transistors and phototransistors based on amorphous materials

A comparison between the charge transport properties in low molecular amorphous thin films of spiro-linked compound and their corresponding parent compound has been demonstrated. The field-effect transistor method is used for extracting physical parameters such as field-effect mobility of charge carriers, ON/OFF ratios, and stability. In addition, phototransistors have been fabricated and demonstrated for the first time by using organic materials. In this case, asymmetrically spiro-linked compounds are used as active materials. The active materials used in this study can be divided into three classes, namely Spiro-linked compounds (symmetrically spiro-linked compounds), the corresponding parent-compounds, and photosensitive spiro-linked compounds (asymmetrically spiro-linked com-pounds). Some of symmetrically spiro-linked compounds used in this study were 2,2',7,7'-Tetrakis-(di-phenylamino)-9,9'-spirobifluorene (Spiro-TAD),2,2',7,7'-Tetrakis-(N,N'-di-p-methylphenylamino)-9,9'-spirobifluorene (Spiro-TTB), 2,2',7,7'-Tetra-(m-tolyl-phenylamino)-9,9'-spirobifluorene (Spiro-TPD), and 2,2Ž,7,7Ž-Tetra-(N-phenyl-1-naphtylamine)-9,9Ž-spirobifluorene (Spiro alpha-NPB). Related parent compounds of the symmetrically spiro-linked compound used in this study were N,N,N',N'-Tetraphenylbenzidine (TAD), N,N,N',N'-Tetrakis(4-methylphenyl)benzidine (TTB), N,N'-Bis(3-methylphenyl)-(1,1'-biphenyl)-4,4'-diamine (TPD), and N,N'-Diphenyl-N,N'-bis(1-naphthyl)-1,1'-biphenyl-4,4'-diamine (alpha-NPB). The photosensitive asymmetrically spiro-linked compounds used in this study were 2,7-bis-(N,N'-diphenylamino)-2',7'-bis(biphenyl-4-yl)-9,9'-spirobifluorene (Spiro-DPSP), and 2,7-bis-(N,N'-diphenylamino)-2',7'-bis(spirobifluorene-2-yl)-9,9'-spirobifluorene (Spiro-DPSP^2). It was found that the field-effect mobilities of charge carriers in thin films of symmetrically spiro-linked compounds and their corresponding parent compounds are in the same order of magnitude (~10^-5 cm^2/Vs). However, the thin films of the parent compounds were easily crystallized after the samples have been exposed in ambient atmosphere and at room temperature for three days. In contrast, the thin films and the transistor characteristics of symmetrically spiro-linked compound did not change significantly after the samples have been stored in ambient atmosphere and at room temperature for several months. Furthermore, temperature dependence of the mobility was analyzed in two models, namely the Arrhenius model and the Gaussian Disorder model. The Arrhenius model tends to give a high value of the prefactor mobility. However, it is difficult to distinguish whether the temperature behaviors of the material under consideration follows the Arrhenius model or the Gaussian Disorder model due to the narrow accessible range of the temperatures. For the first time, phototransistors have been fabricated and demonstrated by using organic materials. In this case, asymmetrically spiro-linked compounds are used as active materials. Intramolecular charge transfer between a bis(diphenylamino)biphenyl unit and a sexiphenyl unit leads to an increase in charge carrier density, providing the amplification effect. The operational responsivity of better than 1 A/W can be obtained for ultraviolet light at 370 nm, making the device interesting for sensor applications. This result offers a new potential application of organic thin film phototransistors as low-light level and low-cost visible blind ultraviolet photodetectors.

Hierarchical Mixtures of Experts and the EM Algorithm

We present a tree-structured architecture for supervised learning. The statistical model underlying the architecture is a hierarchical mixture model in which both the mixture coefficients and the mixture components are generalized linear models (GLIM's). Learning is treated as a maximum likelihood problem; in particular, we present an Expectation-Maximization (EM) algorithm for adjusting the parameters of the architecture. We also develop an on-line learning algorithm in which the parameters are updated incrementally. Comparative simulation results are presented in the robot dynamics domain.

Multivariate Density Estimation: An SVM Approach

We formulate density estimation as an inverse operator problem. We then use convergence results of empirical distribution functions to true distribution functions to develop an algorithm for multivariate density estimation. The algorithm is based upon a Support Vector Machine (SVM) approach to solving inverse operator problems. The algorithm is implemented and tested on simulated data from different distributions and different dimensionalities, gaussians and laplacians in $R^2$ and $R^{12}$. A comparison in performance is made with Gaussian Mixture Models (GMMs). Our algorithm does as well or better than the GMMs for the simulations tested and has the added advantage of being automated with respect to parameters.