Modeling in-Hospital Patient Survival During the First 28 Days After Intensive Care Unit Admission: a Prognostic Model for Clinical Trials in General Critically Ill Patients

OBJECTIVE: The objective of the study was to develop a model for estimating patient 28-day in-hospital mortality using 2 different statistical approaches. DESIGN: The study was designed to develop an outcome prediction model for 28-day in-hospital mortality using (a) logistic regression with random effects and (b) a multilevel Cox proportional hazards model. SETTING: The study involved 305 intensive care units (ICUs) from the basic Simplified Acute Physiology Score (SAPS) 3 cohort. PATIENTS AND PARTICIPANTS: Patients (n = 17138) were from the SAPS 3 database with follow-up data pertaining to the first 28 days in hospital after ICU admission. INTERVENTIONS: None. MEASUREMENTS AND RESULTS: The database was divided randomly into 5 roughly equal-sized parts (at the ICU level). It was thus possible to run the model-building procedure 5 times, each time taking four fifths of the sample as a development set and the remaining fifth as the validation set. At 28 days after ICU admission, 19.98% of the patients were still in the hospital. Because of the different sampling space and outcome variables, both models presented a better fit in this sample than did the SAPS 3 admission score calibrated to vital status at hospital discharge, both on the general population and in major subgroups. CONCLUSIONS: Both statistical methods can be used to model the 28-day in-hospital mortality better than the SAPS 3 admission model. However, because the logistic regression approach is specifically designed to forecast 28-day mortality, and given the high uncertainty associated with the assumption of the proportionality of risks in the Cox model, the logistic regression approach proved to be superior.

C-reactive protein in early risk assessment in acute pancreatitis

ABSTRACT:C-reactive protein (CRP) has been widely used in the early risk assessment of patients with acute pancreatitis (AP), but unclear aspects about its prognostic accuracy in this setting persist. This project evaluated first CRP prognostic accuracy for severity, pancreatic necrosis (PNec), and in-hospital mortality (IM) in AP in terms of the best timing for CRP measurement and the optimal CRP cutoff points. Secondly it was evaluated the CRP measured at approximately 24 hours after hospital admission (CRP24) prognostic accuracy for IM in AP individually and in a combined model with a recent developed tool for the early risk assessment of patients with AP, the Bedside Index for Severity in AP (BISAP). Two single-centre retrospective cohort studies were held. The first study included 379 patients and the second study included 134 patients. Statistical methods such as the Hosmer-Lemeshow goodness-of-fit test, the area under the receiver-operating characteristic curve, the net reclassification improvement, and the integrated discrimination improvement were used. It was found that CRP measured at approximately 48 hours after hospital admission (CRP48) had a prognostic accuracy for severity, PNec, and IM in AP better than CRP measured at any other timing. It was observed that the optimal CRP48 cutoff points for severity, PNec, and IM in AP varied from 170mg/l to 190mg/l, values greater than the one most often recommended in the literature – 150mg/l. It was found that CRP24 had a good prognostic accuracy for IM in AP and that the cutoff point of 60mg/l had a negative predictive value of 100%. Finally it was observed that the prognostic accuracy of a combined model including BISAP and CRP24 for IM in AP could perform better than the BISAP alone model. These results might have a direct impact on the early risk assessment of patients with AP in the daily clinical practice.--------- RESUMO: A proteina c-reactiva (CRP) tem sido largamente usada na avaliação precoce do risco em doentes com pancreatite aguda (AP), mas aspectos duvidosos acerca do seu valor prognóstico neste contexto persistem. Este projecto avaliou primeiro o valor prognóstico da CRP para a gravidade, a necrose pancreática (PNec) e a mortalidade intra-hospitalar (IM) na AP em termos do melhor momento para efectuar a sua medição e dos seus pontos-de-corte óptimos. Em segundo lugar foi avaliado o valor prognóstico da proteína c-reactiva medida aproximadamente às 24 horas após a admissão hospitalar (CRP24) para a IM na AP isoladamente e num modelo combinado, que incluiu uma ferramenta de avaliação precoce do risco em doentes com AP recentemente desenvolvida, o Bedside Index for Severity in Acute Pancreatitis (BISAP). Dois estudos unicêntricos de coorte retrospectivo foram realizados. O primeiro estudo incluiu 379 doentes e o segundo estudo incluiu 134 doentes. Metodologias estatísticas como o teste de Hosmer-Lemeshow goodness-of-fit, a area under the receiver-operating characteristic curve, o net reclassification improvement e o integrated discrimination improvement foram usadas. Verificou-se que a CRP medida às 48 horas após a admissão hospitalar (CRP48) teve um valor prognóstico para a gravidade, a PNec e a IM na AP melhor do que a CRP medida em qualquer outro momento. Observou-se que os pontos de corte óptimos da CRP48 para a gravidade, a PNec e a IM na AP variaram entre 170mg/l e 190mg/l, valores acima do valor mais frequentemente recomendado na literatura – 150mg/l. Verificou-se que a CRP medida aproximadamente às 24 horas após a admissão hospitalar (CRP24) teve um bom valor prognóstico para a IM na AP e que o ponto de corte 60mg/l teve um valor preditivo negativo de 100%. Finalmente observou-se que o valor prognóstico de um modelo combinado incluindo o BISAP e a CRP24 para a IM na AP pode ter um desempenho melhor do que o do BISAP isoladamente. Estes resultados podem ter um impacto directo na avaliação precoce do risco em doentes com AP na prática clínica diária.

Behavioural assessment impacts on the performance results

A Work Project, presented as part of the requirements for the Award of a Masters Degree in Management from the NOVA – School of Business and Economics

Spatial Patterns and Irregularities of the electoral data: general elections in Canada

Dissertation submitted in partial fulfillment of the requirements for the Degree of Master of Science in Geospatial Technologies.

Evolution of methodology for the study of adaptability and stability in cultivated species.

The GxE interaction only became widely discussed from evolutionary studies and evaluations of the causes of behavioral changes of species cultivated in environments. In the last 60 years, several methodologies for the study of adaptability and stability of genotypes in multiple environments trials were developed in order to assist the breeder's choice regarding which genotypes are more stable and which are the most suitable for the crops in the most diverse environments. The methods that use linear regression analysis were the first to be used in a general way by breeders, followed by multivariate analysis methods and mixed models. The need to identify the genetic and environmental causes that are behind the GxE interaction led to the development of new models that include the use of covariates and which can also include both multivariate methods and mixed modeling. However, further studies are needed to identify the causes of GxE interaction as well as for the more accurate measurement of its effects on phenotypic expression of varieties in competition trials carried out in genetic breeding programs.

Politicization of non-political events

Spatial analysis and social network analysis typically take into consideration social processes in specific contexts of geographical or network space. The research in political science increasingly strives to model heterogeneity and spatial dependence. To better understand and geographically model the relationship between “non-political” events, streaming data from social networks, and political climate was the primary objective of the current study. Geographic information systems (GIS) are useful tools in the organization and analysis of streaming data from social networks. In this study, geographical and statistical analysis were combined in order to define the temporal and spatial nature of the data eminating from the popular social network Twitter during the 2014 FIFA World Cup. The study spans the entire globe because Twitter’s geotagging function, the fundamental data that makes this study possible, is not limited to a geographic area. By examining the public reactions to an inherenlty non-political event, this study serves to illuminate broader questions about social behavior and spatial dependence. From a practical perspective, the analyses demonstrate how the discussion of political topics fluсtuate according to football matches. Tableau and Rapidminer, in addition to a set basic statistical methods, were applied to find patterns in the social behavior in space and time in different geographic regions. It was found some insight into the relationship between an ostensibly non-political event – the World Cup - and public opinion transmitted by social media. The methodology could serve as a prototype for future studies and guide policy makers in governmental and non-governmental organizations in gauging the public opinion in certain geographic locations.

Experimental and modeling studies on reverse electrodialysis for sustainable power generation

A potentially renewable and sustainable source of energy is the chemical energy associated with solvation of salts. Mixing of two aqueous streams with different saline concentrations is spontaneous and releases energy. The global theoretically obtainable power from salinity gradient energy due to World’s rivers discharge into the oceans has been estimated to be within the range of 1.4-2.6 TW. Reverse electrodialysis (RED) is one of the emerging, membrane-based, technologies for harvesting the salinity gradient energy. A common RED stack is composed by alternately-arranged cation- and anion-exchange membranes, stacked between two electrodes. The compartments between the membranes are alternately fed with concentrated (e.g., sea water) and dilute (e.g., river water) saline solutions. Migration of the respective counter-ions through the membranes leads to ionic current between the electrodes, where an appropriate redox pair converts the chemical salinity gradient energy into electrical energy. Given the importance of the need for new sources of energy for power generation, the present study aims at better understanding and solving current challenges, associated with the RED stack design, fluid dynamics, ionic mass transfer and long-term RED stack performance with natural saline solutions as feedwaters. Chronopotentiometry was used to determinate diffusion boundary layer (DBL) thickness from diffusion relaxation data and the flow entrance effects on mass transfer were found to avail a power generation increase in RED stacks. Increasing the linear flow velocity also leads to a decrease of DBL thickness but on the cost of a higher pressure drop. Pressure drop inside RED stacks was successfully simulated by the developed mathematical model, in which contribution of several pressure drops, that until now have not been considered, was included. The effect of each pressure drop on the RED stack performance was identified and rationalized and guidelines for planning and/or optimization of RED stacks were derived. The design of new profiled membranes, with a chevron corrugation structure, was proposed using computational fluid dynamics (CFD) modeling. The performance of the suggested corrugation geometry was compared with the already existing ones, as well as with the use of conductive and non-conductive spacers. According to the estimations, use of chevron structures grants the highest net power density values, at the best compromise between the mass transfer coefficient and the pressure drop values. Finally, long-term experiments with natural waters were performed, during which fouling was experienced. For the first time, 2D fluorescence spectroscopy was used to monitor RED stack performance, with a dedicated focus on following fouling on ion-exchange membrane surfaces. To extract relevant information from fluorescence spectra, parallel factor analysis (PARAFAC) was performed. Moreover, the information obtained was then used to predict net power density, stack electric resistance and pressure drop by multivariate statistical models based on projection to latent structures (PLS) modeling. The use in such models of 2D fluorescence data, containing hidden, but extractable by PARAFAC, information about fouling on membrane surfaces, considerably improved the models fitting to the experimental data.

Changes in prescribing patterns of benzodiazepines after training of general practitioners

RESUMO: Introdução: As benzodiazepinas são os fármacos ansiolíticos e hipnóticos mais utilizados. O elevado consumo destes fármacos tem representado uma preocupação devido aos efeitos secundários do seu uso prolongado e dependência. Portugal tem a maior utilização de benzodiazepinas na Europa. Este estudo pretende analisar a alteração do padrão de prescrição de benzodiazepinas após uma intervenção com clínicos gerais. Métodos: A intervenção consistiu numa sessão educacional a um grupo de clínicos gerais. Foi comparado o padrão de prescrição de benzodiazepinas dos médicos intervencionados com o de um grupo de médicos não intervencionado da mesma região e com o de um grupo de médicos não intervencionados de outra região. Analisaram-­‐se as prescrições de 12 meses antes e depois da intervenção. A análise do padrão de prescrição utilizou como metodologia a Dose Diária Definida (DDD) e a Dose Diária Definida por 1000 pacientes por dia (DHD). A análise estatística recorreu a métodos de regressão segmentada. Resultados: Houve uma diminuição no padrão de prescrição de benzodiazepinas no grupo intervencionado após a intervenção (p=0.005). Houve também uma redução no padrão de prescrição no grupo não intervencionada da mesma região (p=0.037) e no grupo não-intervencionado da região diferente (p=0.010). Analisando por género, prescritores do género feminino prescrevem uma quantidade maior de benzodiazepinas. Os clínicos gerais do género feminino intervencionados tiveram a maior redução na prescrição após a intervenção (p=0.008). Discussão: Os dados demonstraram que a intervenção reduziu a prescrição de benzodiazepinas após a intervenção. A diminuição geral do padrão de prescrição poderá ser explicada pelo efeito de Hawthorne ou pela contaminação entre os três grupos de clínicos gerais. Os dados disponíveis não explicam as diferenças nos padrões de prescrição por género. Conclusão: Este estudo demonstra como uma única intervenção tem um impacto positivo na melhoria dos padrões de prescrição. A replicação desta intervenção poderá representar uma oportunidade para alterar a prescrição de benzodiazepinas em Portugal. -----------------------------ABSTRACT: Introduction: Benzodiazepines are the most utilized anxiolytic and hypnotic drugs. The high consumption of benzodiazepines has been a concern due to the reported side effects of long-­‐term use and dependence. Portugal has the highest benzodiazepine utilisation in Europe. This study aims to analyse the change in General Practitioners’ (GPs) benzodiazepine prescription pattern after na intervention period. Methods: An educational session was delivered to a group of intervened GPs. The benzodiazepine prescription pattern of the intervened group was compared to the pattern of a non-­‐intervened matched group from the same region, and to the pattern of another non-­‐intervened matched group from a diferente region. The research time frame was 12 month before and after intervention. The analysis of the prescription trends used the Defined Daily Dose (DDD) and Defined Daily Dose per 1000 patients per day (DHD) methodology. The statistical methods consisted of segmented regression analysis. Results: There was a decrease in benzodiazepine prescription pattern of intervened GPs after intervention (p=0.005). There was also a decrease in benzodiazepine prescription pattern for the non-­‐intervened group from the same region (p=0.037) and for the non-­‐ intervened group from a diferente region (p=0.010). Concerningthe analysis by gender, female gender prescribed a higher amount of benzodiazepines. The intervened female gender prescribers presented the highest decrease in prescription trend after intervention (p=0.008). Discussion: The data demonstrated that the intervention was effective in reducing benzodiazepine prescription after intervention. The general decrease in prescription trend might be explained by a Hawthorne effect or a contamination effect between the three groups of GPs. The available data couldn´t explain the diferences in prescription patterns by gender. Conclusion: This study demonstrates how a single intervention has a positive impact on improving prescription trends. The replication of this intervention might be an opportunity to changing the worrying benzodiazepine utilisation in Portugal.

Análise de dados e aplicação de modelos estatísticos para apoio à decisão de contratação na Petróleo Brasileiro

Dissertação de mestrado em Engenharia Industrial

Assistive locomotion strategies for active lower limb devices

Dissertação de mestrado integrado em Engenharia Biomédica (área de especialização em Eletrónica Médica)

A gestão de recursos humanos nas pequenas e médias empresas em Portugal: Definição de um modelo de atuação

Tese de Doutoramento em Ciências Empresariais

Skull morphometrics of adult male Antartic fur seal, Arctocephalus gazella, and South American fur seal A. australis

The skull morphometrics of adult male Antarctic fur seal, Arctocephalus gazella (Peters, 1875) and South American fur seal, A. australis (Zimmermann, 1783) were investigated using a collection of 45 and 38 skulls, respectively. Eighteen measurements were taken for each specimen. Comparative univariate and multivariate statistical analyses included standard statistics, one-way analysis of variance, principal component analysis and discriminant analysis. Individual variation was relatively high for some variables, as expressed by the coefficient of variation. Skulls of A. gazella were larger than those of A. australis for all but two variables: squamosal jugal suture and rostral length. Both species differed significantly as shown by both univariate and multivariate analyses. The discriminant function correctly classified all specimens. The standardized canonical coefficients showed that the variables which most contribute to the differentiation between species were, in decreasing order, the rostral length, palatal length, palatal width at postcanine 5 and braincase width. The present study corroborates that A. gazella and A. australis are phenotipically distinct species.

Cranial sexual discrimination in hatchling broad-snouted caiman (Caiman latirostris)

Broad-snouted caiman (Caiman latirostris) hatchlings present a consistent sexual dimorphism in their cranium shape and size. Male hatchlings have smaller crania than females. Using multivariate statistical analyses it is possible to discriminate sex in broad-snouted caiman hatchlings by their cranial shape with a reasonable efficiency. The understanding of sexual dimorphism of crocodilian hatchlings might be possibly improved by experimental approach considering, genetic and phenotypic variables such as incubation temperature and clutch of origin.

Les traces de pneumatiques en sciences forensiques : de l'analyse chimique à l'évaluation de la force probante

This work is focused on the development of a methodology for the use of chemical characteristic of tire traces to help answer the following question: "Is the offending tire at the origin of the trace found on the crime scene?". This methodology goes from the trace sampling on the road to statistical analysis of its chemical characteristics. Knowledge about the composition and manufacture of tread tires as well as a review of instrumental techniques used for the analysis of polymeric materials were studied to select, as an ansi vi cal technique for this research, pyrolysis coupled to a gas Chromatograph with a mass spectrometry detector (Py-GC/MS). An analytical method was developed and optimized to obtain the lowest variability between replicates of the same sample. Within-variability of the tread was evaluated regarding width and circumference with several samples taken from twelve tires of different brands and/or models. The variability within each of the treads (within-variability) and between the treads (between-variability) could be quantified. Different statistical methods have shown that within-variability is lower than between-variability, which helped differentiate these tires. Ten tire traces were produced with tires of different brands and/or models by braking tests. These traces have been adequately sampled using sheets of gelatine. Particles of each trace were analysed using the same methodology as for the tires at their origin. The general chemical profile of a trace or of a tire has been characterized by eighty-six compounds. Based on a statistical comparison of the chemical profiles obtained, it has been shown that a tire trace is not differentiable from the tire at its origin but is generally differentiable from tires that are not at its origin. Thereafter, a sample containing sixty tires was analysed to assess the discrimination potential of the developed methodology. The statistical results showed that most of the tires of different brands and models are differentiable. However, tires of the same brand and model with identical characteristics, such as country of manufacture, size and DOT number, are not differentiable. A model, based on a likelihood ratio approach, was chosen to evaluate the results of the comparisons between the chemical profiles of the traces and tires. The methodology developed was finally blindly tested using three simulated scenarios. Each scenario involved a trace of an unknown tire as well as two tires possibly at its origin. The correct results for the three scenarios were used to validate the developed methodology. The different steps of this work were useful to collect the required information to test and validate the underlying assumption that it is possible to help determine if an offending tire » or is not at the origin of a trace, by means of a statistical comparison of their chemical profile. This aid was formalized by a measure of the probative value of the evidence, which is represented by the chemical profile of the trace of the tire. - Ce travail s'est proposé de développer une méthodologie pour l'exploitation des caractéristiques chimiques des traces de pneumatiques dans le but d'aider à répondre à la question suivante : «Est-ce que le pneumatique incriminé est ou n'est pas à l'origine de la trace relevée sur les lieux ? ». Cette méthodologie s'est intéressée du prélèvement de la trace de pneumatique sur la chaussée à l'exploitation statistique de ses caractéristiques chimiques. L'acquisition de connaissances sur la composition et la fabrication de la bande de roulement des pneumatiques ainsi que la revue de techniques instrumentales utilisées pour l'analyse de matériaux polymériques ont permis de choisir, comme technique analytique pour la présente recherche, la pyrolyse couplée à un chromatographe en phase gazeuse avec un détecteur de spectrométrie de masse (Py-GC/MS). Une méthode analytique a été développée et optimisée afin d'obtenir la plus faible variabilité entre les réplicas d'un même échantillon. L'évaluation de l'intravariabilité de la bande de roulement a été entreprise dans sa largeur et sa circonférence à l'aide de plusieurs prélèvements effectués sur douze pneumatiques de marques et/ou modèles différents. La variabilité au sein de chacune des bandes de roulement (intravariabilité) ainsi qu'entre les bandes de roulement considérées (intervariabilité) a pu être quantifiée. Les différentes méthodes statistiques appliquées ont montré que l'intravariabilité est plus faible que l'intervariabilité, ce qui a permis de différencier ces pneumatiques. Dix traces de pneumatiques ont été produites à l'aide de pneumatiques de marques et/ou modèles différents en effectuant des tests de freinage. Ces traces ont pu être adéquatement prélevées à l'aide de feuilles de gélatine. Des particules de chaque trace ont été analysées selon la même méthodologie que pour les pneumatiques à leur origine. Le profil chimique général d'une trace de pneumatique ou d'un pneumatique a été caractérisé à l'aide de huitante-six composés. Sur la base de la comparaison statistique des profils chimiques obtenus, il a pu être montré qu'une trace de pneumatique n'est pas différenciable du pneumatique à son origine mais est, généralement, différenciable des pneumatiques qui ne sont pas à son origine. Par la suite, un échantillonnage comprenant soixante pneumatiques a été analysé afin d'évaluer le potentiel de discrimination de la méthodologie développée. Les méthodes statistiques appliquées ont mis en évidence que des pneumatiques de marques et modèles différents sont, majoritairement, différenciables entre eux. La méthodologie développée présente ainsi un bon potentiel de discrimination. Toutefois, des pneumatiques de la même marque et du même modèle qui présentent des caractéristiques PTD (i.e. pays de fabrication, taille et numéro DOT) identiques ne sont pas différenciables. Un modèle d'évaluation, basé sur une approche dite du likelihood ratio, a été adopté pour apporter une signification au résultat des comparaisons entre les profils chimiques des traces et des pneumatiques. La méthodologie mise en place a finalement été testée à l'aveugle à l'aide de la simulation de trois scénarios. Chaque scénario impliquait une trace de pneumatique inconnue et deux pneumatiques suspectés d'être à l'origine de cette trace. Les résultats corrects obtenus pour les trois scénarios ont permis de valider la méthodologie développée. Les différentes étapes de ce travail ont permis d'acquérir les informations nécessaires au test et à la validation de l'hypothèse fondamentale selon laquelle il est possible d'aider à déterminer si un pneumatique incriminé est ou n'est pas à l'origine d'une trace, par le biais d'une comparaison statistique de leur profil chimique. Cette aide a été formalisée par une mesure de la force probante de l'indice, qui est représenté par le profil chimique de la trace de pneumatique.

Urinary oxalate and urate to creatinine ratios in a healthy pediatric population.

The purpose of the study was to determine reference percentiles for the urinary (U) oxalate (Ox) and urate (Ura) to creatinine (Cr) concentration ratios in the second morning urine of healthy infants, children, and adolescents. The urinary oxalate and urate to creatinine ratios were determined in the spontaneously voided second morning urine sample. To test reproducibility, two urine samples were analyzed on 2 consecutive weeks in 63% of the subjects. Three hundred eighty-four healthy children (181 girls, 203 boys), aged 1 month to 17 years, from nurseries, kindergartens, and schools of Lausanne, Switzerland, were studied. The 5th and 95th percentiles were determined from the total number of urine samples (627) after confirmation that there was no order effect between repeated measurements and there were no significant sex differences. A nonlinear regression analysis in terms of age was used to smooth the calculated percentiles. In this manner, curves were obtained from which the reference values can be read at any given age. The 95th percentiles decreased with age: for UOx/Cr from 0.175 mg/mg (0.22 mol/mol) at 1 to 6 months to 0.048 mg/mg (0.06 mol/mol) from 7 years and beyond; and UUra/Cr from 2.378 mg/mg (1.6 mol/mol) at 1 to 6 months to 0.594 mg/mg (0.4 mol/mol) in adolescence. We provide 5th and 95th percentile curves for the UOx/Cr and UUra/Cr ratios determined from the second morning urine samples in a large cohort of healthy infants, children, and adolescents. Values were determined by standard analytical chemical techniques and were analyzed by powerful statistical methods. The calculated 95th percentile for the UOx/Cr values fell rather rapidly and reached normal adult values by the age of 7 years, whereas for UUra/Cr, the 95th percentile decreased slowly and stabilized in adolescence.