Analysis of nonlinear dynamics of vocal folds using high-speed video observation and biomechanical modeling

Primary voice production occurs in the larynx through vibrational movements carried out by vocal folds. However, many problems can affect this complex system resulting in voice disorders. In this context, time-frequency-shape analysis based on embedding phase space plots and nonlinear dynamics methods have been used to evaluate the vocal fold dynamics during phonation. For this purpose, the present work used high-speed video to record the vocal fold movements of three subjects and extract the glottal area time series using an image segmentation algorithm. This signal is used for an optimization method which combines genetic algorithms and a quasi-Newton method to optimize the parameters of a biomechanical model of vocal folds based on lumped elements (masses, springs and dampers). After optimization, this model is capable of simulating the dynamics of recorded vocal folds and their glottal pulse. Bifurcation diagrams and phase space analysis were used to evaluate the behavior of this deterministic system in different circumstances. The results showed that this methodology can be used to extract some physiological parameters of vocal folds and reproduce some complex behaviors of these structures contributing to the scientific and clinical evaluation of voice production. (C) 2010 Elsevier Inc. All rights reserved.

New Optimization Algorithms for Structural Reliability Analysis

Solution of structural reliability problems by the First Order method require optimization algorithms to find the smallest distance between a limit state function and the origin of standard Gaussian space. The Hassofer-Lind-Rackwitz-Fiessler (HLRF) algorithm, developed specifically for this purpose, has been shown to be efficient but not robust, as it fails to converge for a significant number of problems. On the other hand, recent developments in general (augmented Lagrangian) optimization techniques have not been tested in aplication to structural reliability problems. In the present article, three new optimization algorithms for structural reliability analysis are presented. One algorithm is based on the HLRF, but uses a new differentiable merit function with Wolfe conditions to select step length in linear search. It is shown in the article that, under certain assumptions, the proposed algorithm generates a sequence that converges to the local minimizer of the problem. Two new augmented Lagrangian methods are also presented, which use quadratic penalties to solve nonlinear problems with equality constraints. Performance and robustness of the new algorithms is compared to the classic augmented Lagrangian method, to HLRF and to the improved HLRF (iHLRF) algorithms, in the solution of 25 benchmark problems from the literature. The new proposed HLRF algorithm is shown to be more robust than HLRF or iHLRF, and as efficient as the iHLRF algorithm. The two augmented Lagrangian methods proposed herein are shown to be more robust and more efficient than the classical augmented Lagrangian method.

Space-time cluster analysis of American visceral leishmaniasis in Bauru, Sao Paulo State, Brazil

A space-time analysis of American visceral leishmaniasis (AVL) in humans in the city of Bauru, Sao Paulo State, Brazil was carried out based on 239 cases diagnosed between June 2003 and October 2008. Spatial analysis of the disease showed that cases occurred especially in the city's urban areas. AVL annual incidence rates were calculated, demonstrating that the highest rate occurred in 2006 (19.55/100,000 inhabitants). This finding was confirmed by the time series analysis, which also showed a positive tendency over the period analyzed. The present study allows us to conclude that the disease was clustered in the Southwest side of the city in 2006, suggesting that this area may require special attention with regard to control and prevention measures.

Space-time cluster analysis of American visceral leishmaniasis in Bauru, São Paulo State, Brazil

A space-time analysis of American visceral leishmaniasis (AVL) in humans in the city of Bauru, São Paulo State, Brazil was carried out based on 239 cases diagnosed between June 2003 and October 2008. Spatial analysis of the disease showed that cases occurred especially in the city's urban areas. AVL annual incidence rates were calculated, demonstrating that the highest rate occurred in 2006 (19.55/100,000 inhabitants). This finding was confirmed by the time series analysis, which also showed a positive tendency over the period analyzed. The present study allows us to conclude that the disease was clustered in the Southwest side of the city in 2006, suggesting that this area may require special attention with regard to control and prevention measures.

Crystallization and preliminary X-ray diffraction analysis of endoglucanase III from Trichoderma harzianum

Endoglucanases are enzymes that hydrolyze cellulose and are important components of the cellulolytic complex. In contrast to other members of the complex, they cleave internal beta-1,4-glycosidic bonds in the cellulose polymer, allowing cellulose to be used as an energy source. Since biomass is an important renewable source of energy, the structural and functional characterization of these enzymes is of interest. In this study, endoglucanase III from Trichoderma harzianum was produced in Pichia pastoris and purified. Crystals belonging to the orthorhombic space group P212121, with unit-cell parameters a = 47.54, b = 55.57, c = 157.3 angstrom, were obtained by the sitting-drop vapour-diffusion method and an X-ray diffraction data set was collected to 2.07 angstrom resolution.

Bayesian analysis of weak gravitational lensing and Sunyaev-Zel'dovich data for six galaxy clusters

We present an analysis of observations made with the Arcminute Microkelvin Imager (AMI) and the CanadaFranceHawaii Telescope (CFHT) of six galaxy clusters in a redshift range of 0.160.41. The cluster gas is modelled using the SunyaevZeldovich (SZ) data provided by AMI, while the total mass is modelled using the lensing data from the CFHT. In this paper, we (i) find very good agreement between SZ measurements (assuming large-scale virialization and a gas-fraction prior) and lensing measurements of the total cluster masses out to r200; (ii) perform the first multiple-component weak-lensing analysis of A115; (iii) confirm the unusual separation between the gas and mass components in A1914 and (iv) jointly analyse the SZ and lensing data for the relaxed cluster A611, confirming our use of a simulation-derived masstemperature relation for parametrizing measurements of the SZ effect.

Identifying regulational alterations in gene regulatory networks by state space representation of vector autoregressive models and variational annealing

Background: In the analysis of effects by cell treatment such as drug dosing, identifying changes on gene network structures between normal and treated cells is a key task. A possible way for identifying the changes is to compare structures of networks estimated from data on normal and treated cells separately. However, this approach usually fails to estimate accurate gene networks due to the limited length of time series data and measurement noise. Thus, approaches that identify changes on regulations by using time series data on both conditions in an efficient manner are demanded. Methods: We propose a new statistical approach that is based on the state space representation of the vector autoregressive model and estimates gene networks on two different conditions in order to identify changes on regulations between the conditions. In the mathematical model of our approach, hidden binary variables are newly introduced to indicate the presence of regulations on each condition. The use of the hidden binary variables enables an efficient data usage; data on both conditions are used for commonly existing regulations, while for condition specific regulations corresponding data are only applied. Also, the similarity of networks on two conditions is automatically considered from the design of the potential function for the hidden binary variables. For the estimation of the hidden binary variables, we derive a new variational annealing method that searches the configuration of the binary variables maximizing the marginal likelihood. Results: For the performance evaluation, we use time series data from two topologically similar synthetic networks, and confirm that our proposed approach estimates commonly existing regulations as well as changes on regulations with higher coverage and precision than other existing approaches in almost all the experimental settings. For a real data application, our proposed approach is applied to time series data from normal Human lung cells and Human lung cells treated by stimulating EGF-receptors and dosing an anticancer drug termed Gefitinib. In the treated lung cells, a cancer cell condition is simulated by the stimulation of EGF-receptors, but the effect would be counteracted due to the selective inhibition of EGF-receptors by Gefitinib. However, gene expression profiles are actually different between the conditions, and the genes related to the identified changes are considered as possible off-targets of Gefitinib. Conclusions: From the synthetically generated time series data, our proposed approach can identify changes on regulations more accurately than existing methods. By applying the proposed approach to the time series data on normal and treated Human lung cells, candidates of off-target genes of Gefitinib are found. According to the published clinical information, one of the genes can be related to a factor of interstitial pneumonia, which is known as a side effect of Gefitinib.

Shape analysis using fractal dimension: A curvature based approach

The present work shows a novel fractal dimension method for shape analysis. The proposed technique extracts descriptors from a shape by applying a multi-scale approach to the calculus of the fractal dimension. The fractal dimension is estimated by applying the curvature scale-space technique to the original shape. By applying a multi-scale transform to the calculus, we obtain a set of descriptors which is capable of describing the shape under investigation with high precision. We validate the computed descriptors in a classification process. The results demonstrate that the novel technique provides highly reliable descriptors, confirming the efficiency of the proposed method. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4757226]

Probing the anomalous extinction of four young star clusters: the use of colour-excess, main-sequence fitting and fractal analysis

Aims. We studied four young star clusters to characterise their anomalous extinction or variable reddening and asses whether they could be due to contamination by either dense clouds or circumstellar effects. Methods. We evaluated the extinction law (R-V) by adopting two methods: (i) the use of theoretical expressions based on the colour-excess of stars with known spectral type; and (ii) the analysis of two-colour diagrams, where the slope of the observed colour distribution was compared to the normal distribution. An algorithm to reproduce the zero-age main-sequence (ZAMS) reddened colours was developed to derive the average visual extinction (A(V)) that provides the closest fit to the observational data. The structure of the clouds was evaluated by means of a statistical fractal analysis, designed to compare their geometric structure with the spatial distribution of the cluster members. Results. The cluster NGC 6530 is the only object of our sample affected by anomalous extinction. On average, the other clusters suffer normal extinction, but several of their members, mainly in NGC 2264, seem to have high R-V, probably because of circumstellar effects. The ZAMS fitting provides A(V) values that are in good agreement with those found in the literature. The fractal analysis shows that NGC 6530 has a centrally concentrated distribution of stars that differs from the substructures found in the density distribution of the cloud projected in the A(V) map, suggesting that the original cloud was changed by the cluster formation. However, the fractal dimension and statistical parameters of Berkeley 86, NGC 2244, and NGC 2264 indicate that there is a good cloud-cluster correlation, when compared to other works based on an artificial distribution of points.

Combinatorial-topological framework for the analysis of global dynamics

We discuss an algorithmic framework based on efficient graph algorithms and algebraic-topological computational tools. The framework is aimed at automatic computation of a database of global dynamics of a given m-parameter semidynamical system with discrete time on a bounded subset of the n-dimensional phase space. We introduce the mathematical background, which is based upon Conley's topological approach to dynamics, describe the algorithms for the analysis of the dynamics using rectangular grids both in phase space and parameter space, and show two sample applications. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4767672]

Morphological Homogeneity of Neurons: Searching for Outlier Neuronal Cells

We report a morphology-based approach for the automatic identification of outlier neurons, as well as its application to the NeuroMorpho.org database, with more than 5,000 neurons. Each neuron in a given analysis is represented by a feature vector composed of 20 measurements, which are then projected into a two-dimensional space by applying principal component analysis. Bivariate kernel density estimation is then used to obtain the probability distribution for the group of cells, so that the cells with highest probabilities are understood as archetypes while those with the smallest probabilities are classified as outliers. The potential of the methodology is illustrated in several cases involving uniform cell types as well as cell types for specific animal species. The results provide insights regarding the distribution of cells, yielding single and multi-variate clusters, and they suggest that outlier cells tend to be more planar and tortuous. The proposed methodology can be used in several situations involving one or more categories of cells, as well as for detection of new categories and possible artifacts.

Transiting exoplanets from the CoRoT space mission XXI. CoRoT-19b: a low density planet orbiting an old inactive F9V-star

Context. Observations of transiting extrasolar planets are of key importance to our understanding of planets because their mass, radius, and mass density can be determined. These measurements indicate that planets of similar mass can have very different radii. For low-density planets, it is generally assumed that they are inflated owing to their proximity to the host-star. To determine the causes of this inflation, it is necessary to obtain a statistically significant sample of planets with precisely measured masses and radii. Aims. The CoRoT space mission allows us to achieve a very high photometric accuracy. By combining CoRoT data with high-precision radial velocity measurements, we derive precise planetary radii and masses. We report the discovery of CoRoT-19b, a gas-giant planet transiting an old, inactive F9V-type star with a period of four days. Methods. After excluding alternative physical configurations mimicking a planetary transit signal, we determine the radius and mass of the planet by combining CoRoT photometry with high-resolution spectroscopy obtained with the echelle spectrographs SOPHIE, HARPS, FIES, and SANDIFORD. To improve the precision of its ephemeris and the epoch, we observed additional transits with the TRAPPIST and Euler telescopes. Using HARPS spectra obtained during the transit, we then determine the projected angle between the spin of the star and the orbit of the planet. Results. We find that the host star of CoRoT-19b is an inactive F9V-type star close to the end of its main-sequence life. The host star has a mass M-* = 1.21 +/- 0.05 M-circle dot and radius R-* = 1.65 +/- 0.04 R-circle dot. The planet has a mass of M-P = 1.11 +/- 0.06 M-Jup and radius of R-P = 1.29 +/- 0.03 R-Jup. The resulting bulk density is only rho = 0.71 +/- 0.06 g cm (3), which is much lower than that for Jupiter. Conclusions. The exoplanet CoRoT-19b is an example of a giant planet of almost the same mass as Jupiter but a approximate to 30% larger radius.

Transiting exoplanets from the CoRoT space mission XXII. CoRoT-16b: a hot Jupiter with a hint of eccentricity around a faint solar-like star

Aims. We report the discovery of CoRoT-16b, a low density hot jupiter that orbits a faint G5V star (mV = 15.63) in 5.3523 +/- 0.0002 days with slight eccentricity. A fit of the data with no a priori assumptions on the orbit leads to an eccentricity of 0.33 +/- 0.1. We discuss this value and also derive the mass and radius of the planet. Methods. We analyse the photometric transit curve of CoRoT-16 given by the CoRoT satellite, and radial velocity data from the HARPS and HIRES spectrometers. A combined analysis using a Markov chain Monte Carlo algorithm is used to get the system parameters. Results. CoRoT-16b is a 0.535 -0.083/+0.085 M-J, 1.17 -0.14/+0.16 R-J hot Jupiter with a density of 0.44 -0.14/+0.21 g cm(-3). Despite its short orbital distance (0.0618 +/- 0.0015 AU) and the age of the parent star (6.73 +/- 2.8 Gyr), the planet orbit exhibits significantly non-zero eccentricity. This is very uncommon for this type of objects as tidal effects tend to circularise the orbit. This value is discussed taking into account the characteristics of the star and the observation accuracy.

An anatomical analysis of the mini-modified orbitozygomatic and supra-orbital approaches

Seven sides of cadaver heads were used to compare the surgical exposures provided by the mini-modified orbitozygomatic (MOz) and supra-orbital (SO) approaches. The Optotrak 3020 computerized tracking system (Northern Digital, Waterloo, ON, Canada) was utilized to evaluate the area of anatomical exposure defined by six points: (1) ipsilateral sphenoid ridge; (2) most distal point of the ipsilateral middle cerebral artery (MCA); (3) most distal point of the ipsilateral posterior cerebral artery (PCA); (4) most distal point of the contralateral PCA; (5) most distal point of the contralateral MCA; and (6) contralateral sphenoid ridge. Additionally, angles of approach for the ipsilateral MCA bifurcation, ipsilateral ICA bifurcation, basilar artery tip, contralateral MCA and ICA bifurcation and anterior communicating artery (AcomA) were evaluated, first for SO and then for MOz. An image guidance system was used to evaluate the limits of surgical exposure. No differences in the area of surgical exposure were noted (p > 0.05). Vertical angles were significantly wider for the ipsilateral and contralateral ICA bifurcation, AcomA, contralateral MCA and basilar tip (p < 0.05) for MOz. No differences in horizontal angles were observed between the approaches for the six targets (p > 0.05). There were no differences in the limits of exposure. MOz affords no additional surgical working space. However, our results demonstrate systematically that vertical exposure is improved. The MOz should be performed while planning an approach to these regions and a wider exposure in the vertical axis is needed. (C) 2012 Elsevier Ltd. All rights reserved.

Interpolation estimate for a finite-element space with embedded discontinuities

We consider a recently proposed finite-element space that consists of piecewise affine functions with discontinuities across a smooth given interface Γ (a curve in two dimensions, a surface in three dimensions). Contrary to existing extended finite element methodologies, the space is a variant of the standard conforming Formula space that can be implemented element by element. Further, it neither introduces new unknowns nor deteriorates the sparsity structure. It is proved that, for u arbitrary in Formula, the interpolant Formula defined by this new space satisfies Graphic where h is the mesh size, Formula is the domain, Formula, Formula, Formula and standard notation has been adopted for the function spaces. This result proves the good approximation properties of the finite-element space as compared to any space consisting of functions that are continuous across Γ, which would yield an error in the Formula-norm of order Graphic. These properties make this space especially attractive for approximating the pressure in problems with surface tension or other immersed interfaces that lead to discontinuities in the pressure field. Furthermore, the result still holds for interfaces that end within the domain, as happens for example in cracked domains.