Statistical shape analysis via principal factor analysis

Statistical shape analysis techniques commonly employed in the medical imaging community, such as active shape models or active appearance models, rely on principal component analysis (PCA) to decompose shape variability into a reduced set of interpretable components. In this paper we propose principal factor analysis (PFA) as an alternative and complementary tool to PCA providing a decomposition into modes of variation that can be more easily interpretable, while still being a linear efficient technique that performs dimensionality reduction (as opposed to independent component analysis, ICA). The key difference between PFA and PCA is that PFA models covariance between variables, rather than the total variance in the data. The added value of PFA is illustrated on 2D landmark data of corpora callosa outlines. Then, a study of the 3D shape variability of the human left femur is performed. Finally, we report results on vector-valued 3D deformation fields resulting from non-rigid registration of ventricles in MRI of the brain.

Heart rate variability does not discriminate between different levels of haemodynamic responsiveness during surgical anaesthesia

BACKGROUND: Hypnotic depth but not haemodynamic responsiveness is measured with EEG-based monitors. In this study we compared heart rate variability (HRV) in unstimulated patients and stimulation-induced HRV at different levels of anaesthesia. METHODS: A total of 95 ASA I or II patients were randomly assigned to five groups (Group 1: BIS 45(5), remifentanil 1 ng ml(-1); Group 2: BIS 45(5), remifentanil 2 ng ml(-1); Group 3: BIS 45(5), remifentanil 4 ng ml(-1); Group 4: BIS 30(5), remifentanil 2 ng ml(-1); Group 5: BIS 60(5), remifentanil 2 ng ml(-1)). A time- and frequency-domain analysis of the RR interval (RRI) from the electrocardiogram was performed. HRV before induction, before and after a 5 s tetanic stimulus of the ulnar nerve, and before and after tracheal intubation was compared between groups, between stimuli, and between responders to intubation [systolic arterial pressure (SAP) increase >20 mm Hg, a maximal heart rate (HR) after intubation >90 min(-1) or both] and non-responders (anova). RESULTS: Induction of anaesthesia significantly lowered HR and HRV. Mean RRI before stimulation was higher in G3 than in G1, G2, and G4 (P < 0.001), whereas the other HRV parameters were similar. Intubation induced a greater HRV response than tetanic stimulation. The mean RRI after intubation was lower in G3 compared with the other groups and the sd of the RRI after tetanic stimulation was lower in G3 compared with G5. Otherwise, unstimulated HRV and stimulation-induced HRV were similar in responders and non-responders. CONCLUSION: HRV parameters discriminate between awake and general anaesthesia, are different after tracheal intubation and a 5 s ulnar nerve stimulation, but do not discriminate between different levels of haemodynamic responsiveness during surgical anaesthesia.

Fluctuation analysis of lung function as a predictor of long-term response to beta2-agonists

The response to beta(2)-agonists differs between asthmatics and has been linked to subsequent adverse events, even death. Possible determinants include beta(2)-adrenoceptor genotype at position 16, lung function and airway hyperresponsiveness. Fluctuation analysis provides a simple parameter alpha measuring the complex correlation properties of day-to-day peak expiratory flow. The present study investigated whether alpha predicts clinical response to beta(2)-agonist treatment, taking into account other conventional predictors. Analysis was performed on previously published twice-daily peak expiratory flow measurements in 66 asthmatic adults over three 6-month randomised order treatment periods: placebo, salbutamol and salmeterol. Multiple linear regression was used to determine the association between alpha during the placebo period and response to treatment (change in the number of days with symptoms), taking into account other predictors namely beta(2)-adrenoceptor genotype, lung function and its variability, and airway hyperresponsiveness. The current authors found that alpha measured during the placebo period considerably improved the prediction of response to salmeterol treatment, taking into account genotype, lung function or its variability, or airway hyperresponsiveness. The present study provides further evidence that response to beta(2)-agonists is related to the time correlation properties of lung function in asthma. The current authors conclude that fluctuation analysis of lung function offers a novel predictor to identify patients who may respond well or poorly to treatment.

Slope stability analysis of the Pacaya Volcano, Guatemala, using limit equilibrium and finite element method

The Pacaya volcanic complex is part of the Central American volcanic arc, which is associated with the subduction of the Cocos tectonic plate under the Caribbean plate. Located 30 km south of Guatemala City, Pacaya is situated on the southern rim of the Amatitlan Caldera. It is the largest post-caldera volcano, and has been one of Central America’s most active volcanoes over the last 500 years. Between 400 and 2000 years B.P, the Pacaya volcano had experienced a huge collapse, which resulted in the formation of horseshoe-shaped scarp that is still visible. In the recent years, several smaller collapses have been associated with the activity of the volcano (in 1961 and 2010) affecting its northwestern flanks, which are likely to be induced by the local and regional stress changes. The similar orientation of dry and volcanic fissures and the distribution of new vents would likely explain the reactivation of the pre-existing stress configuration responsible for the old-collapse. This paper presents the first stability analysis of the Pacaya volcanic flank. The inputs for the geological and geotechnical models were defined based on the stratigraphical, lithological, structural data, and material properties obtained from field survey and lab tests. According to the mechanical characteristics, three lithotechnical units were defined: Lava, Lava-Breccia and Breccia-Lava. The Hoek and Brown’s failure criterion was applied for each lithotechnical unit and the rock mass friction angle, apparent cohesion, and strength and deformation characteristics were computed in a specified stress range. Further, the stability of the volcano was evaluated by two-dimensional analysis performed by Limit Equilibrium (LEM, ROCSCIENCE) and Finite Element Method (FEM, PHASE 2 7.0). The stability analysis mainly focused on the modern Pacaya volcano built inside the collapse amphitheatre of “Old Pacaya”. The volcanic instability was assessed based on the variability of safety factor using deterministic, sensitivity, and probabilistic analysis considering the gravitational instability and the effects of external forces such as magma pressure and seismicity as potential triggering mechanisms of lateral collapse. The preliminary results from the analysis provide two insights: first, the least stable sector is on the south-western flank of the volcano; second, the lowest safety factor value suggests that the edifice is stable under gravity alone, and the external triggering mechanism can represent a likely destabilizing factor.

Development of record-breaking statistics for climatological time-series analysis

The number of record-breaking events expected to occur in a strictly stationary time-series depends only on the number of values in the time-series, regardless of distribution. This holds whether the events are record-breaking highs or lows and whether we count from past to present or present to past. However, these symmetries are broken in distinct ways by trends in the mean and variance. We define indices that capture this information and use them to detect weak trends from multiple time-series. Here, we use these methods to answer the following questions: (1) Is there a variability trend among globally distributed surface temperature time-series? We find a significant decreasing variability over the past century for the Global Historical Climatology Network (GHCN). This corresponds to about a 10% change in the standard deviation of inter-annual monthly mean temperature distributions. (2) How are record-breaking high and low surface temperatures in the United States affected by time period? We investigate the United States Historical Climatology Network (USHCN) and find that the ratio of record-breaking highs to lows in 2006 increases as the time-series extend further into the past. When we consider the ratio as it evolves with respect to a fixed start year, we find it is strongly correlated with the ensemble mean. We also compare the ratios for USHCN and GHCN (minus USHCN stations). We find the ratios grow monotonically in the GHCN data set, but not in the USHCN data set. (3) Do we detect either mean or variance trends in annual precipitation within the United States? We find that the total annual and monthly precipitation in the United States (USHCN) has increased over the past century. Evidence for a trend in variance is inconclusive.

Analysis of in-cylinder pressure transducer data quality utilizing a SIDI turbocharged engine

In-cylinder pressure transducers have been used for decades to record combustion pressure inside a running engine. However, due to the extreme operating environment, transducer design and installation must be considered in order to minimize measurement error. One such error is caused by thermal shock, where the pressure transducer experiences a high heat flux that can distort the pressure transducer diaphragm and also change the crystal sensitivity. This research focused on investigating the effects of thermal shock on in-cylinder pressure transducer data quality using a 2.0L, four-cylinder, spark-ignited, direct-injected, turbo-charged GM engine. Cylinder four was modified with five ports to accommodate pressure transducers of different manufacturers. They included an AVL GH14D, an AVL GH15D, a Kistler 6125C, and a Kistler 6054AR. The GH14D, GH15D, and 6054AR were M5 size transducers. The 6125C was a larger, 6.2mm transducer. Note that both of the AVL pressure transducers utilized a PH03 flame arrestor. Sweeps of ignition timing (spark sweep), engine speed, and engine load were performed to study the effects of thermal shock on each pressure transducer. The project consisted of two distinct phases which included experimental engine testing as well as simulation using a commercially available software package. A comparison was performed to characterize the quality of the data between the actual cylinder pressure and the simulated results. This comparison was valuable because the simulation results did not include thermal shock effects. All three sets of tests showed the peak cylinder pressure was basically unaffected by thermal shock. Comparison of the experimental data with the simulated results showed very good correlation. The spark sweep was performed at 1300 RPM and 3.3 bar NMEP and showed that the differences between the simulated results (no thermal shock) and the experimental data for the indicated mean effective pressure (IMEP) and the pumping mean effective pressure (PMEP) were significantly less than the published accuracies. All transducers had an IMEP percent difference less than 0.038% and less than 0.32% for PMEP. Kistler and AVL publish that the accuracy of their pressure transducers are within plus or minus 1% for the IMEP (AVL 2011; Kistler 2011). In addition, the difference in average exhaust absolute pressure between the simulated results and experimental data was the greatest for the two Kistler pressure transducers. The location and lack of flame arrestor are believed to be the cause of the increased error. For the engine speed sweep, the torque output was held constant at 203 Nm (150 ft-lbf) from 1500 to 4000 RPM. The difference in IMEP was less than 0.01% and the PMEP was less than 1%, except for the AVL GH14D which was 5% and the AVL GH15DK which was 2.25%. A noticeable error in PMEP appeared as the load increased during the engine speed sweeps, as expected. The load sweep was conducted at 2000 RPM over a range of NMEP from 1.1 to 14 bar. The difference in IMEP values were less 0.08% while the PMEP values were below 1% except for the AVL GH14D which was 1.8% and the AVL GH15DK which was at 1.25%. In-cylinder pressure transducer data quality was effectively analyzed using a combination of experimental data and simulation results. Several criteria can be used to investigate the impact of thermal shock on data quality as well as determine the best location and thermal protection for various transducers.

Relationship between heart rate variability, interleukin-6, and soluble tissue factor in healthy subjects

Decreased heart rate variability (HRV) has been associated with an increased risk of atherosclerosis. We hypothesized that a decrease in frequency domains of resting HRV would be associated with elevated plasma levels of interleukin (IL)-6 and soluble tissue factor (sTF) both previously shown to prospectively predict atherothrombotic events in healthy subjects. Subjects were 102 healthy and unmedicated black and white middle-aged men and women. We determined IL-6 and sTF antigen in plasma and HRV measures from surface electrocardiogram data using spectral analysis. All statistical analyses controlled for age, gender, ethnicity, smoking status, blood pressure, and body mass index. Low amounts of low frequency (LF) power (beta=-0.31, p=0.007) and high frequency (HF) power (beta=-0.36, p=0.002) were associated with increased amounts of IL-6, explaining 7% and 9% of the variance, respectively. Interactions between LF power and IL-6 (p=0.002) and between HF power and IL-6 (p=0.012) explained 8% and 5%, respectively, of the variance in sTF. Post hoc analyses showed associations between IL-6 and sTF when LF power (beta=0.51, p<0.001) and HF power (beta=0.48, p<0.001) were low but not when LF power and high HF power were high. The findings suggest that systemic low-grade inflammatory activity is associated with a decrease in HRV. Furthermore, there was a positive relationship between plasma levels of IL-6 and sTF antigen when HRV was low. Inflammation and related hypercoagulability might particularly contribute to atherothrombotic events in a setting of decreased HRV.

Risk of severe asthma episodes predicted from fluctuation analysis of airway function

Asthma is an increasing health problem worldwide, but the long-term temporal pattern of clinical symptoms is not understood and predicting asthma episodes is not generally possible. We analyse the time series of peak expiratory flows, a standard measurement of airway function that has been assessed twice daily in a large asthmatic population during a long-term crossover clinical trial. Here we introduce an approach to predict the risk of worsening airflow obstruction by calculating the conditional probability that, given the current airway condition, a severe obstruction will occur within 30 days. We find that, compared with a placebo, a regular long-acting bronchodilator (salmeterol) that is widely used to improve asthma control decreases the risk of airway obstruction. Unexpectedly, however, a regular short-acting beta2-agonist bronchodilator (albuterol) increases this risk. Furthermore, we find that the time series of peak expiratory flows show long-range correlations that change significantly with disease severity, approaching a random process with increased variability in the most severe cases. Using a nonlinear stochastic model, we show that both the increased variability and the loss of correlations augment the risk of unstable airway function. The characterization of fluctuations in airway function provides a quantitative basis for objective risk prediction of asthma episodes and for evaluating the effectiveness of therapy.

In vivo analysis of efavirenz metabolism in individuals with impaired CYP2A6 function

INTRODUCTION: The antiretroviral drug efavirenz (EFV) is extensively metabolized into three primary metabolites: 8-hydroxy-EFV, 7-hydroxy-EFV and N-glucuronide-EFV. There is a wide interindividual variability in EFV plasma exposure, explained to a great extent by cytochrome P450 2B6 (CYP2B6), the main isoenzyme responsible for EFV metabolism and involved in the major metabolic pathway (8-hydroxylation) and to a lesser extent in 7-hydroxylation. When CYP2B6 function is impaired, the relevance of CYP2A6, the main isoenzyme responsible for 7-hydroxylation may increase. We hypothesize that genetic variability in this gene may contribute to the particularly high, unexplained variability in EFV exposure in individuals with limited CYP2B6 function. METHODS: This study characterized CYP2A6 variation (14 alleles) in individuals (N=169) previously characterized for functional variants in CYP2B6 (18 alleles). Plasma concentrations of EFV and its primary metabolites (8-hydroxy-EFV, 7-hydroxy-EFV and N-glucuronide-EFV) were measured in different genetic backgrounds in vivo. RESULTS: The accessory metabolic pathway CYP2A6 has a critical role in limiting drug accumulation in individuals characterized as CYP2B6 slow metabolizers. CONCLUSION: Dual CYP2B6 and CYP2A6 slow metabolism occurs at significant frequency in various human populations, leading to extremely high EFV exposure.

Climate Variability and Extremes during the Past 100 Years

Large progress has been made in the past few years towards quantifying and understanding climate variability during past centuries. At the same time, present-day climate has been studied using state-of-the-art data sets and tools with respect to the physical and chemical mechanisms governing climate variability. Both the understanding of the past and the knowledge of the processes are important for assessing and attributing the anthropogenic effect on present and future climate. The most important time period in this context is the past approximately 100 years, which comprises large natural variations and extremes (such as long droughts) as well as anthropogenic influences, most pronounced in the past few decades. Recent and ongoing research efforts steadily improve the observational record of the 20th century, while atmospheric circulation models are used to underpin the mechanisms behind large climatic variations. Atmospheric chemistry and composition are important for understanding climate variability and change, and considerable progress has been made in the past few years in this field. The evolving integration of these research areas in a more comprehensive analysis of recent climate variability was reflected in the organisation of a workshop “Climate variability and extremes in the past 100 years” in Gwatt near Thun (Switzerland), 24–26 July 2006. The aim of this workshop was to bring together scientists working on data issues together with statistical climatologists, modellers, and atmospheric chemists to discuss gaps in our understanding of climate variability during the past approximately 100 years.

Tree phenology and carbon dioxide fluxes - use of digital photography for process-based interpretation at the ecosystem scale

Vegetation phenology is an important indicator of climate change and climate variability and it is strongly connected to biospheric–atmospheric gas exchange. We aimed to evaluate the applicability of phenological information derived from digital imagery for the interpretation of CO2 exchange measurements. For the years 2005–2007 we analyzed seasonal phenological development of 2 temperate mixed forests using tower-based imagery from standard RGB cameras. Phenological information was jointly analyzed with gross primary productivity (GPP) derived from net ecosystem exchange data. Automated image analysis provided reliable information on vegetation developmental stages of beech and ash trees covering all seasons. A phenological index derived from image color values was strongly correlated with GPP, with a significant mean time lag of several days for ash trees and several weeks for beech trees in early summer (May to mid-July). Leaf emergence dates for the dominant tree species partly explained temporal behaviour of spring GPP but were also masked by local meteorological conditions. We conclude that digital cameras at flux measurement sites not only provide an objective measure of the physiological state of a forest canopy at high temporal and spatial resolutions, but also complement CO2 and water exchange measurements, improving our knowledge of ecosystem processes.

Multi-locus microsatellite analysis supports the hypothesis of an autochthonous focus of Echinococcus multilocularis in northern Italy

Echinococcus multilocularis is characterised by a wide geographical distribution, encompassing three continents (North America, Asia and Europe) yet very low genetic variability is documented. Recently, this parasite has been detected in red foxes (Vulpes vulpes) circulating in an Alpine region of Italy, close to Austria. This finding raised the question as to whether an autochthonous cycle exists in Italy or whether the infected foxes originated from the neighbouring regions of Austria. Studies have shown that multi-locus microsatellite analysis can identify genomic regions carrying mutations that result in a local adaptation. We used a tandem repeated multi-locus microsatellite (EmsB) to evaluate the genetic differences amongst adult worms of E. multilocularis collected in Italy, worms from neighbouring Austria and from other European and extra-European countries. Fluorescent PCR was performed on a panel of E. multilocularis samples to assess intra-specific polymorphism. The analysis revealed four closed genotypes for Italian samples of E. multilocularis which were unique compared with the other 25 genotypes from Europe and the five genotypes from Alaska. An analysis in the Alpine watershed, comparing Italian adult worms with those from neighbouring areas in Austria, showed a unique cluster for Italian samples. This result supports the hypothesis of the presence of an autochthonous cycle of E. multilocularis in Italy. EmsB can be useful for 'tracking' the source of infection of this zoonotic parasite and developing appropriate measures for preventing or reducing the risk of human alveolar echinococcosis.

Variability of in vivo fish acute toxicity data

The variability of toxicity data contained within databases was investigated using the widely used US EPA ECOTOX database as an example. Fish acute lethality (LC50) values for 44 compounds (for which at least 10 data entries existed) were extracted from the ECOTOX database yielding a total of 4654 test records. Significant variability of LC50 test results was observed, exceeding several orders of magnitude. In an attempt to systematically explore potential causes of the data variability, the influence of biological factors (such as test species or life stages) and physical factors (such as water temperature, pH or water hardness) were examined. Even after eliminating the influence of these inherent factors, considerable data variability remained, suggesting an important role of factors relating to technical and measurement procedures. The analysis, however, was limited by pronounced gaps in the test documentation. Of the 4654 extracted test reports, 66.5% provided no information on the fish life stage used for testing. Likewise, water temperature, hardness or pH were not recorded in 19.6%, 48.2% and 41.2% of the data entries, respectively. From these findings, we recommend the rigorous control of data entries ensuring complete recording of testing conditions. A more consistent database will help to better discriminate between technical and natural variability of the test data, which is of importance in ecological risk assessment for extrapolation from laboratory tests to the field, and also might help to develop correction factors that account for systematic differences in test results caused by species, life stage or test conditions.

1200 years of decadal-scale variability of Mediterranean vegetation and climate at Pantelleria Island, Italy

A new sedimentary sequence from Lago di Venere on Pantelleria Island, located in the Strait of Sicily between Tunisia and Sicily was recovered. The lake is located in the coastal infra-Mediterranean vegetation belt at 2 m a.s.l. Pollen, charcoal and sedimentological analyses are used to explore linkages among vegetation, fire and climate at a decadal scale over the past 1200 years. A dry period from ad 800 to 1000 that corresponds to the ‘Medieval Warm Period’ (WMP) is inferred from sedimentological analysis. The high content of carbonate recorded in this period suggests a dry phase, when the ratio of evaporation/precipitation was high. During this period the island was dominated by thermophilous and drought-tolerant taxa, such as Quercus ilex, Olea, Pistacia and Juniperus. A marked shift in the sediment properties is recorded at ad 1000, when carbonate content became very low suggesting wetter conditions until ad 1850–1900. Broadly, this period coincides with the ‘Little Ice Age’ (LIA), which was characterized by wetter and colder conditions in Europe. During this time rather mesic conifers (i.e. Pinus pinaster), shrubs and herbs (e.g. Erica arborea and Selaginella denticulata) expanded, whereas more drought-adapted species (e.g. Q. ilex) declined. Charcoal data suggest enhanced fire activity during the LIA probably as a consequence of anthropogenic burning and/or more flammable fuel (e.g. resinous Pinus biomass). The last century was characterized by a shift to high carbonate content, indicating a change towards drier conditions, and re-expansion of Q. ilex and Olea. The post-LIA warming is in agreement with historical documents and meteorological time series. Vegetation dynamics were co-determined by agricultural activities on the island. Anthropogenic indicators (e.g. Cerealia-type, Sporormiella) reveal the importance of crops and grazing on the island. Our pollen data suggest that extensive logging caused the local extinction of deciduous Quercus pubescens around ad1750.

Phylogenetic analysis of Prevotella nigrescens, Prevotella intermedia and Porphyromonas gingivalis clinical strains reveals a clear species clustering

Prevotella nigrescens, Prevotella intermedia and Porphyromonas gingivalis are oral pathogens from the family Bacteroidaceae, regularly isolated from cases of gingivitis and periodontitis. In this study, the phylogenetic variability of these three bacterial species was investigated by means of 16S rRNA (rrs) gene sequence comparisons of a set of epidemiologically and geographically diverse isolates. For each of the three species, the rrs gene sequences of 11 clinical isolates as well as the corresponding type strains was determined. Comparison of all rrs sequences obtained with those of closely related species revealed a clear clustering of species, with only a little intraspecies variability but a clear difference in the rrs gene with respect to the next related taxon. The results indicate that the three species form stable, homogeneous genetic groups, which favours an rrs-based species identification of these oral pathogens. This is especially useful given the 7% sequence divergence between Prevotella intermedia and Prevotella nigrescens, since phenotypic distinction between the two Prevotella species is inconsistent or involves techniques not applicable in routine identification.