Calcium carbonate global LGM Burial rate data

Global databases of calcium carbonate concentrations and mass accumulation rates in Holocene and last glacial maximum sediments were used to estimate the deep-sea sedimentary calcium carbonate burial rate during these two time intervals. Sparse calcite mass accumulation rate data were extrapolated across regions of varying calcium carbonate concentration using a gridded map of calcium carbonate concentrations and the assumption that accumulation of noncarbonate material is uncorrelated with calcite concentration within some geographical region. Mean noncarbonate accumulation rates were estimated within each of nine regions, determined by the distribution and nature of the accumulation rate data. For core-top sediments the regions of reasonable data coverage encompass 67% of the high-calcite (>75%) sediments globally, and within these regions we estimate an accumulation rate of 55.9 ± 3.6 x 10**11 mol/yr. The same regions cover 48% of glacial high-CaCO3 sediments (the smaller fraction is due to a shift of calcite deposition to the poorly sampled South Pacific) and total 44.1 ± 6.0 x 10**11 mol/yr. Projecting both estimates to 100 % coverage yields accumulation estimates of 8.3 x 10**12 mol/yr today and 9.2 x 10**12 mol/yr during glacial time. This is little better than a guess given the incomplete data coverage, but it suggests that glacial deep sea calcite burial rate was probably not considerably faster than today in spite of a presumed decrease in shallow water burial during glacial time.

Inferences about Parameters of Trivariate Normal Distribution with Missing Data

Multivariate normal distribution is commonly encountered in any field, a frequent issue is the missing values in practice. The purpose of this research was to estimate the parameters in three-dimensional covariance permutation-symmetric normal distribution with complete data and all possible patterns of incomplete data. In this study, MLE with missing data were derived, and the properties of the MLE as well as the sampling distributions were obtained. A Monte Carlo simulation study was used to evaluate the performance of the considered estimators for both cases when ρ was known and unknown. All results indicated that, compared to estimators in the case of omitting observations with missing data, the estimators derived in this article led to better performance. Furthermore, when ρ was unknown, using the estimate of ρ would lead to the same conclusion.

A hybrid Particle Swarm Optimization - Simplex algorithm (PSOS) for structural damage identification

This study proposes a new PSOS-model based damage identification procedure using frequency domain data. The formulation of the objective function for the minimization problem is based on the Frequency Response Functions (FRFs) of the system. A novel strategy for the control of the Particle Swarm Optimization (PSO) parameters based on the Nelder-Mead algorithm (Simplex method) is presented; consequently, the convergence of the PSOS becomes independent of the heuristic constants and its stability and confidence are enhanced. The formulated hybrid method performs better in different benchmark functions than the Simulated Annealing (SA) and the basic PSO (PSO(b)). Two damage identification problems, taking into consideration the effects of noisy and incomplete data, were studied: first, a 10-bar truss and second, a cracked free-free beam, both modeled with finite elements. In these cases, the damage location and extent were successfully determined. Finally, a non-linear oscillator (Duffing oscillator) was identified by PSOS providing good results. (C) 2009 Elsevier Ltd. All rights reserved

Nonparametric estimation of time-dependent ROC curves conditional on a continuous covariate

The receiver-operating characteristic (ROC) curve is the most widely used measure for evaluating the performance of a diagnostic biomarker when predicting a binary disease outcome. The ROC curve displays the true positive rate (or sensitivity) and the false positive rate (or 1-specificity) for different cut-off values used to classify an individual as healthy or diseased. In time-to-event studies, however, the disease status (e.g. death or alive) of an individual is not a fixed characteristic, and it varies along the study. In such cases, when evaluating the performance of the biomarker, several issues should be taken into account: first, the time-dependent nature of the disease status; and second, the presence of incomplete data (e.g. censored data typically present in survival studies). Accordingly, to assess the discrimination power of continuous biomarkers for time-dependent disease outcomes, time-dependent extensions of true positive rate, false positive rate, and ROC curve have been recently proposed. In this work, we present new nonparametric estimators of the cumulative/dynamic time-dependent ROC curve that allow accounting for the possible modifying effect of current or past covariate measures on the discriminatory power of the biomarker. The proposed estimators can accommodate right-censored data, as well as covariate-dependent censoring. The behavior of the estimators proposed in this study will be explored through simulations and illustrated using data from a cohort of patients who suffered from acute coronary syndrome.

Access to specialized pediatric cancer care in Switzerland.

BACKGROUND: Specialized pediatric cancer centers (PCCs) are thought to be essential to obtain state-of-the-art care for children and adolescents. We determined the proportion of childhood cancer patients not treated in a PCC, and described their characteristics and place of treatment. PROCEDURE: The Swiss Childhood Cancer Registry (SCCR) registers all children treated in Swiss PCCs. The regional cancer registries (covering 14/26 cantons) register all cancer patients of a region. The children of the SCCR with data from 7 regions (11 cantons) were compared, using specialized software for record linkage. All children <16 years of age at diagnosis with primary malignant tumors, diagnosed between 1990 and 2004, and living in one of these regions were included in the analysis. RESULTS: 22.1% (238/1,077) of patients recorded in regional registries were not registered in the SCCR. Of these, 15.7% (169/1,077) had never been in a PCC while 6.4% (69/1,077) had been in a PCC but were not registered in the SCCR, due to incomplete data flow. In all diagnostic groups and in all age groups, a certain proportion of children was treated outside a PCC, but this proportion was largest in children suffering from malignant bone tumors/soft tissue sarcomas and from malignant epithelial neoplasms, and in older children. The proportion of patients treated in a PCC increased over the study period (P < 0.0001). CONCLUSIONS: One in six childhood cancer patients in Switzerland was not treated in a PCC. Whether these patients have different treatment outcomes remained unclear.

Hindsight judgement on ambiguous episodes of suspected infection in critically ill children: poor consensus amongst experts?

Few episodes of suspected infection observed in paediatric intensive care are classifiable without ambiguity by a priori defined criteria. Most require additional expert judgement. Recently, we observed a high variability in antibiotic prescription rates, not explained by the patients' clinical data or underlying diseases. We hypothesised that the disagreement of experts in adjudication of episodes of suspected infection could be one of the potential causes for this variability. During a 5-month period, we included all patients of a 19-bed multidisciplinary, tertiary, neonatal and paediatric intensive care unit, in whom infection was clinically suspected and antibiotics were prescribed ( n=183). Three experts (two senior ICU physicians and a specialist in infectious diseases) were provided with all patient data, laboratory and microbiological findings. All experts classified episodes according to a priori defined criteria into: proven sepsis, probable sepsis (negative cultures), localised infection and no infection. Episodes of proven viral infection and incomplete data sets were excluded. Of the remaining 167 episodes, 48 were classifiable by a priori criteria ( n=28 proven sepsis, n= 20 no infection). The three experts only achieved limited agreement beyond chance in the remaining 119 episodes (kappa = 0.32, and kappa = 0.19 amongst the ICU physicians). The kappa is a measure of the degree of agreement beyond what would be expected by chance alone, with 0 indicating the chance result and 1 indicating perfect agreement. CONCLUSION: agreement of specialists in hindsight adjudication of episodes of suspected infection is of questionable reliability.

Chylothorax in children: guidelines for diagnosis and management.

STUDY OBJECTIVE: To establish guidelines for the diagnosis and management of chylothorax in children. DESIGN: Retrospective study. PATIENTS: Fifty-one patients with a diagnosis of chylothorax. Twelve patients were excluded because of incomplete data or incorrect diagnosis. The following parameters were analyzed: triglyceride level, total cell number, and lymphocyte percentage; amount of pleural effusion on day of diagnosis, day 5, and day 14; and total time of pleural effusion. Prospectively, the same parameters were analyzed in a control group of 10 patients with pleural drainage. INTERVENTION: Patients with chylothorax were treated primarily with fat-free oral nutrition; if chyle did not stop, total parenteral nutrition with total enteric rest was started. If conservative therapy was not successful, pleurodesis was performed. RESULTS: In children with chylothorax triglyceride, triglyceride content ranged from 0.56 to 26.6 mmol/L; all values except one were > 1.1 mmol/L. In 36 of 39 patients (92%), the cell count was > 1,000 cells/microL. In 33 of 39 patients (85%), lymphocytes were > 90%. In patients without chylothorax triglyceride, triglyceride levels ranged from 0.1 to 0.71 mmol/L (median, 0.38 mmol/L) and cell count was from 20 to 1400 cells/microL (median, 322 cells/microL), with a maximum of 60% lymphocytes. With fat-free nutrition, chyle disappeared in 29 of 39 patients. Five patients died, and five required pleurodesis. CONCLUSIONS: Pleural effusion in children is chyle when it contains > 1.1 mmol/L triglycerides (with oral fat intake) and has a total cell count > or 1,000 cells/microL, with a lymphocyte fraction > 80%. Chylous effusions usually last long; however, after 6 weeks, the majority of the effusions (29 of 39 patients) had ceased. Late surgical interventions reduce the number of thoracotomies substantially, but can lead to very long hospitalization times. Early surgical interventions (after < 3 weeks) lead to a high number of thoracotomies, but certainly reduce hospitalization time.

Breast-conserving surgery for non-palpable breast cancer: relationship of lumpectomy resection margins measurements between remote perioperative ultrasound and postoperative histopathology

Purpose: Tumour-free resection margins (RMs) are mandatory in breast-conserving surgery. On-site intraoperative ultrasound (US)-guided tumour resection with extemporaneous histopathological assessment of RMs has been described. Remote intraoperative US assessment of RMs is an alternative. The purpose of this study was to evaluate the relationship of lumpectomy RMs measurements between remote intraoperative US and postoperative histopathology.Methods and Materials: In a retrospective IRB-approved review of 100 consecutive lumpectomies performed between October 2009 and April 2011 for presumed non-palpable breast cancer, 71 women (mean age 63.8years) were included. Twenty-nine patients were excluded because of absence of cancer at histopathology and/or incomplete data. Measurements of lumpectomy minimal RMs and tumour maximal diameter obtained on remote intraoperative US and postoperative histopathology were compared.Results: Minimal RMs were 0.35±0.32 (mean±SD) and 0.35±0.32cm on remote intraoperative US and postoperative histopathology, respectively. No significant difference was found between these measurements (p=0.37). Tumour maximal diameter was 1.02±0.51 (mean±SD) and 1.33±0.74cm on remote intraoperative US and postoperative histopathology, respectively. US measurements were significantly smaller (p<0.001). The 71 breast carcinoma (CA) consisted of: invasive canalar (n=49), invasive lobular (n=11), in situ (n=3) and other types of CA (n=8). Twenty-nine patients had intraoperative re-excision (24 without residual CA), while 16 patients were re-operated due to insufficient histopathological RMs (12 without residual CA).Conclusion: Good correlation of minimal RMs between remote intraoperative US and postoperative histopathology warrants use of both techniques in a complementary manner. Remote intraoperative US is helpful in taking rapid decision of re-excision and maintaining low re-operation rate after breast-conserving surgery for non-palpable cancer.

Does this child have a foreign body aspiration?

Context: Foreign body aspiration (FbA) is a serious problem in children. Accurate clinical and radiographic diagnosis is important because missed or delayed diagnosis can result in respiratory difficulties ranging from life-treatening airway obstruction to chronic wheezing or recurrent pneumonia. Bronchoscopy also has risks and accurate clinical and radiographc diagnosis can support the decision of bronchoscopy. Objective: To rewiev the diagnostic accuracy of clinical presentation (CP) and pulmonary radiograph (PR) for the diagnosis of FbA. There is no previous rewievMethods: A search of Medline is conducted for articles containing data regarding CP and PR signes of FbA. Calculation of likelihood ratios (LR) and pre and post test probability using Bayes theorem were performed for all signs of CP and PR. Inclusion criteria: Articles containing prospective data regarding CP and PR of FbA. Exclusion criteria: Retrospectives studies. Articles containing incomplete data for calculation of LR. Results: Five prospectives studies are included with a total of 585 patients. Prevalence of FbA is 63% in children suspected of FbA. If CP is normal, probability of FbA is 25% and if PR is normal, probability is 14%. If CP is pathologic, probability of FbA is 69-76% with presence of cough (LR = 1.32) or dyspnea (LR = 1.84) or localized crackles (LR = 1.5). Probability is 81-88% if cyanosis (LR = 4.8) or decreased breaths sounds (LR = 4.3) or asymetric auscultation (LR = 2.9) or localized wheezing (LR = 2.5) are present. When CP is anormal and PR show mediatinal shift (LR = 100), pneumomediatin (LR = 100), radio opaque foreign body (LR = 100), lobar distention (LR = 4), atelectasis (LR = 2.5), inspiratory/expiratory abnormal (LR = 7), the probability of FbA is 96-100%. If CP is normal and PR is abnormal the probability is 40-100%. If CP is abnormal and PR is normal the probability is 55-75%. Conclusions: This rewiev of prospective studies demonstrates the importance of CP and PR and an algorithm can be proposed. When CP is abnormal with or without PR pathologic, the probability of FbA is high and bronchoscopy is indicated. When CP and PR are normal the probability of FbA is low and bronchoscopy is not necessary immediatly, observation should be proposed. This approach should be validated with prospective study.

Are patients at nutritional risk more prone to complications after major urological surgery?

PURPOSE: The nutritional risk score is a recommended screening tool for malnutrition. While a nutritional risk score of 3 or greater predicts adverse outcomes after digestive surgery, to our knowledge its predictive value for morbidity after urological interventions is unknown. We determined whether urological patients at nutritional risk are at higher risk for complications after major surgery than patients not at nutritional risk. MATERIALS AND METHODS: We performed a prospective observational study in consecutive patients undergoing major surgery. A priori sample calculation resulted in a study cohort of 220 patients. Interim analysis was planned after 110 patients. The nutritional risk score was assessed preoperatively by a specialized study nurse. Nutritional care was standardized in all patients. Postoperative complications were defined previously using the standardized Dindo-Clavien classification. The primary end point was 30-day morbidity. Univariate and multivariate analysis was performed to identify predictors of complications. RESULTS: The study was discontinued due to significant results after interim analysis. A total of 125 patients were included in analysis from June 2011 to June 2012 and 15 were excluded because of incomplete data. Of 51 patients at nutritional risk 38 (74%) presented with at least 1 complication compared to 28 of 59 controls (47%). Patients at nutritional risk were at threefold risk for complications on univariate and multivariate analysis (OR 3.3, 95% CI 1.3-8.0). Cystectomy was the only other predictor of morbidity (OR 10, 95% CI 2-48). CONCLUSIONS: Patients at nutritional risk are more prone to complications after major urological procedures. Whether this increased morbidity can be reversed by perioperative nutritional support should be studied.

Algorithmic Techniques in Gene Expression Processing. From Imputation to Visualization

The amount of biological data has grown exponentially in recent decades. Modern biotechnologies, such as microarrays and next-generation sequencing, are capable to produce massive amounts of biomedical data in a single experiment. As the amount of the data is rapidly growing there is an urgent need for reliable computational methods for analyzing and visualizing it. This thesis addresses this need by studying how to efficiently and reliably analyze and visualize high-dimensional data, especially that obtained from gene expression microarray experiments. First, we will study the ways to improve the quality of microarray data by replacing (imputing) the missing data entries with the estimated values for these entries. Missing value imputation is a method which is commonly used to make the original incomplete data complete, thus making it easier to be analyzed with statistical and computational methods. Our novel approach was to use curated external biological information as a guide for the missing value imputation. Secondly, we studied the effect of missing value imputation on the downstream data analysis methods like clustering. We compared multiple recent imputation algorithms against 8 publicly available microarray data sets. It was observed that the missing value imputation indeed is a rational way to improve the quality of biological data. The research revealed differences between the clustering results obtained with different imputation methods. On most data sets, the simple and fast k-NN imputation was good enough, but there were also needs for more advanced imputation methods, such as Bayesian Principal Component Algorithm (BPCA). Finally, we studied the visualization of biological network data. Biological interaction networks are examples of the outcome of multiple biological experiments such as using the gene microarray techniques. Such networks are typically very large and highly connected, thus there is a need for fast algorithms for producing visually pleasant layouts. A computationally efficient way to produce layouts of large biological interaction networks was developed. The algorithm uses multilevel optimization within the regular force directed graph layout algorithm.

Génération efficace de graphes d’appels dynamiques complets

Analyser le code permet de vérifier ses fonctionnalités, détecter des bogues ou améliorer sa performance. L’analyse du code peut être statique ou dynamique. Des approches combinants les deux analyses sont plus appropriées pour les applications de taille industrielle où l’utilisation individuelle de chaque approche ne peut fournir les résultats souhaités. Les approches combinées appliquent l’analyse dynamique pour déterminer les portions à problèmes dans le code et effectuent par la suite une analyse statique concentrée sur les parties identifiées. Toutefois les outils d’analyse dynamique existants génèrent des données imprécises ou incomplètes, ou aboutissent en un ralentissement inacceptable du temps d’exécution. Lors de ce travail, nous nous intéressons à la génération de graphes d’appels dynamiques complets ainsi que d’autres informations nécessaires à la détection des portions à problèmes dans le code. Pour ceci, nous faisons usage de la technique d’instrumentation dynamique du bytecode Java pour extraire l’information sur les sites d’appels, les sites de création d’objets et construire le graphe d’appel dynamique du programme. Nous démontrons qu’il est possible de profiler dynamiquement une exécution complète d’une application à temps d’exécution non triviale, et d’extraire la totalité de l’information à un coup raisonnable. Des mesures de performance de notre profileur sur trois séries de benchmarks à charges de travail diverses nous ont permis de constater que la moyenne du coût de profilage se situe entre 2.01 et 6.42. Notre outil de génération de graphes dynamiques complets, nommé dyko, constitue également une plateforme extensible pour l’ajout de nouvelles approches d’instrumentation. Nous avons testé une nouvelle technique d’instrumentation des sites de création d’objets qui consiste à adapter les modifications apportées par l’instrumentation au bytecode de chaque méthode. Nous avons aussi testé l’impact de la résolution des sites d’appels sur la performance générale du profileur.

The significance of volcanic eruption strength and frequency for climate

A simple physical model of the atmospheric effects of large explosive volcanic eruptions is developed. Using only one input parameter - the initial amount of sulphur dioxide injected into the stratosphere - the global-average stratospheric optical-depth perturbation and surface temperature response are modelled. The simplicity of this model avoids issues of incomplete data (applicable to more comprehensive models), making it a powerful and useful tool for atmospheric diagnostics of this climate forcing mechanism. It may also provide a computationally inexpensive and accurate way of introducing volcanic activity into larger climate models. The modelled surface temperature response for an initial sulphur-dioxide injection, coupled with emission-history statistics, is used to demonstrate that the most climatically significant volcanic eruptions are those of sufficient explosivity to just reach into the stratosphere (and achieve longevity). This study also highlights the fact that this measure of significance is highly sensitive to the representation of the climatic response and the frequency data used, and that we are far from producing a definitive history of explosive volcanism for at least the past 1000 years. Given this high degree of uncertainty, these results suggest that eruptions that release around and above 0.1 Mt SO2 into the stratosphere have the maximum climatic impact.

Local influence for Student-t partially linear models

In this paper we extend partial linear models with normal errors to Student-t errors Penalized likelihood equations are applied to derive the maximum likelihood estimates which appear to be robust against outlying observations in the sense of the Mahalanobis distance In order to study the sensitivity of the penalized estimates under some usual perturbation schemes in the model or data the local influence curvatures are derived and some diagnostic graphics are proposed A motivating example preliminary analyzed under normal errors is reanalyzed under Student-t errors The local influence approach is used to compare the sensitivity of the model estimates (C) 2010 Elsevier B V All rights reserved

On estimation and influence diagnostics for the Grubbs` model under heavy-tailed distributions

The Grubbs` measurement model is frequently used to compare several measuring devices. It is common to assume that the random terms have a normal distribution. However, such assumption makes the inference vulnerable to outlying observations, whereas scale mixtures of normal distributions have been an interesting alternative to produce robust estimates, keeping the elegancy and simplicity of the maximum likelihood theory. The aim of this paper is to develop an EM-type algorithm for the parameter estimation, and to use the local influence method to assess the robustness aspects of these parameter estimates under some usual perturbation schemes, In order to identify outliers and to criticize the model building we use the local influence procedure in a Study to compare the precision of several thermocouples. (C) 2008 Elsevier B.V. All rights reserved.