Support vector machines for detecting age-related changes in running kinematics

Age-related changes in running kinematics have been reported in the literature using classical inferential statistics. However, this approach has been hampered by the increased number of biomechanical gait variables reported and subsequently the lack of differences presented in these studies. Data mining techniques have been applied in recent biomedical studies to solve this problem using a more general approach. In the present work, we re-analyzed lower extremity running kinematic data of 17 young and 17 elderly male runners using the Support Vector Machine (SVM) classification approach. In total, 31 kinematic variables were extracted to train the classification algorithm and test the generalized performance. The results revealed different accuracy rates across three different kernel methods adopted in the classifier, with the linear kernel performing the best. A subsequent forward feature selection algorithm demonstrated that with only six features, the linear kernel SVM achieved 100% classification performance rate, showing that these features provided powerful combined information to distinguish age groups. The results of the present work demonstrate potential in applying this approach to improve knowledge about the age-related differences in running gait biomechanics and encourages the use of the SVM in other clinical contexts. (C) 2010 Elsevier Ltd. All rights reserved.

Kinematic, kinetic and EMG patterns during downward squatting

The aim of this study was to investigate the kinematic, kinetic, and electromyographic pattern before, during and after downward squatting when the trunk movement is restricted in the sagittal plane. Eight healthy subjects performed downward squatting at two different positions, semisquatting (40 degrees knee flexion) and half squatting (70 degrees knee flexion). Electromyographic responses of the vastus medialis oblique, vastus medialis longus, rectus femoris, vastus lateralis, biceps femoris, semitendineous, gastrocnemius lateralis, and tibialis anterior were recorded. The kinematics of the major joints were reconstructed using an optoelectronic system. The center of pressure (COP) was obtained using data collected from one force plate, and the ankle and knee joint torques were calculated using inverse dynamics. In the upright position there were small changes in the COP and in the knee and ankle joint torques. The tibialis anterior provoked the disruption of this upright position initiating the squat. During the acceleration phase of the squat the COP moved posteriorly, the knee joint torque remained in flexion and there was no measurable muscle activation. As the body went into the deceleration phase, the knee joint torque increased towards extension with major muscle activities being observed in the four heads of the quadriceps. Understanding these kinematic, kinetic and EMG strategies before, during and after the squat is expected to be beneficial to practitioners for utilizing squatting as a task for improving motor function. (C) 2006 Elsevier Ltd. All rights reserved.

Human gait recognition using extraction and fusion of global motion features

This paper proposes a novel computer vision approach that processes video sequences of people walking and then recognises those people by their gait. Human motion carries different information that can be analysed in various ways. The skeleton carries motion information about human joints, and the silhouette carries information about boundary motion of the human body. Moreover, binary and gray-level images contain different information about human movements. This work proposes to recover these different kinds of information to interpret the global motion of the human body based on four different segmented image models, using a fusion model to improve classification. Our proposed method considers the set of the segmented frames of each individual as a distinct class and each frame as an object of this class. The methodology applies background extraction using the Gaussian Mixture Model (GMM), a scale reduction based on the Wavelet Transform (WT) and feature extraction by Principal Component Analysis (PCA). We propose four new schemas for motion information capture: the Silhouette-Gray-Wavelet model (SGW) captures motion based on grey level variations; the Silhouette-Binary-Wavelet model (SBW) captures motion based on binary information; the Silhouette-Edge-Binary model (SEW) captures motion based on edge information and the Silhouette Skeleton Wavelet model (SSW) captures motion based on skeleton movement. The classification rates obtained separately from these four different models are then merged using a new proposed fusion technique. The results suggest excellent performance in terms of recognising people by their gait.

Evaluation of flow patterns and elongated bubble characteristics during the flow boiling of halocarbon refrigerants in a micro-scale channel

In the present study, quasi-diabatic two-phase flow pattern visualizations and measurements of elongated bubble velocity, frequency and length were performed. The tests were run for R134a and R245fa evaporating in a stainless steel tube with diameter of 2.32 mm, mass velocities ranging from 50 to 600 kg/m(2) s and saturation temperatures of 22 degrees C, 31 degrees C and 41 degrees C. The tube was heated by applying a direct DC current to its surface. Images from a high-speed video-camera (8000 frames/s) obtained through a transparent tube just downstream the heated sections were used to identify the following flow patterns: bubbly, elongated bubbles, churn and annular flows. The visualized flow patterns were compared against the predictions provided by Barnea et al. (1983) [1], Felcar et al. (2007) [10], Revellin and Thome (2007) [3] and Ong and Thome (2009) [11]. From this comparison, it was found that the methods proposed by Felcar et al. (2007) [10] and Ong and Thome (2009) [1] predicted relatively well the present database. Additionally, elongated bubble velocities, frequencies and lengths were determined based on the analysis of high-speed videos. Results suggested that the elongated bubble velocity depends on mass velocity, vapor quality and saturation temperature. The bubble velocity increases with increasing mass velocity and vapor quality and decreases with increasing saturation temperature. Additionally, bubble velocity was correlated as linear functions of the two-phase superficial velocity. (C) 2010 Elsevier Inc. All rights reserved.

An exploratory analysis of relationships between socioeconomic, land use, activity participation variables and travel patterns

This paper aims to find relations between the socioeconomic characteristics, activity participation, land use patterns and travel behavior of the residents in the Sao Paulo Metropolitan Area (SPMA) by using Exploratory Multivariate Data Analysis (EMDA) techniques. The variables influencing travel pattern choices are investigated using: (a) Cluster Analysis (CA), grouping and characterizing the Traffic Zones (17), proposing the independent variable called Origin Cluster and, (b) Decision Tree (DT) to find a priori unknown relations among socioeconomic characteristics, land use attributes of the origin TZ and destination choices. The analysis was based on the origin-destination home-interview survey carried out in SPMA in 1997. The DT application revealed the variables of greatest influence on the travel pattern choice. The most important independent variable considered by DT is car ownership, followed by the Use of Transportation ""credits"" for Transit tariff, and, finally, activity participation variables and Origin Cluster. With these results, it was possible to analyze the influence of a family income, car ownership, position of the individual in the family, use of transportation ""credits"" for transit tariff (mainly for travel mode sequence choice), activities participation (activity sequence choice) and Origin Cluster (destination/travel distance choice). (c) 2010 Elsevier Ltd. All rights reserved.

New prediction methods for CO2 evaporation inside tubes: Part II - An updated general flow boiling heat transfer model based on flow patterns

Corresponding to the updated flow pattern map presented in Part I of this study, an updated general flow pattern based flow boiling heat transfer model was developed for CO2 using the Cheng-Ribatski-Wojtan-Thome [L. Cheng, G. Ribatski, L. Wojtan, J.R. Thome, New flow boiling heat transfer model and flow pattern map for carbon dioxide evaporating inside horizontal tubes, Int. J. Heat Mass Transfer 49 (2006) 4082-4094; L. Cheng, G. Ribatski, L. Wojtan, J.R. Thome, Erratum to: ""New flow boiling heat transfer model and flow pattern map for carbon dioxide evaporating inside tubes"" [Heat Mass Transfer 49 (21-22) (2006) 4082-4094], Int. J. Heat Mass Transfer 50 (2007) 391] flow boiling heat transfer model as the starting basis. The flow boiling heat transfer correlation in the dryout region was updated. In addition, a new mist flow heat transfer correlation for CO2 was developed based on the CO2 data and a heat transfer method for bubbly flow was proposed for completeness sake. The updated general flow boiling heat transfer model for CO2 covers all flow regimes and is applicable to a wider range of conditions for horizontal tubes: tube diameters from 0.6 to 10 mm, mass velocities from 50 to 1500 kg/m(2) s, heat fluxes from 1.8 to 46 kW/m(2) and saturation temperatures from -28 to 25 degrees C (reduced pressures from 0.21 to 0.87). The updated general flow boiling heat transfer model was compared to a new experimental database which contains 1124 data points (790 more than that in the previous model [Cheng et al., 2006, 2007]) in this study. Good agreement between the predicted and experimental data was found in general with 71.4% of the entire database and 83.2% of the database without the dryout and mist flow data predicted within +/-30%. However, the predictions for the dryout and mist flow regions were less satisfactory due to the limited number of data points, the higher inaccuracy in such data, scatter in some data sets ranging up to 40%, significant discrepancies from one experimental study to another and the difficulties associated with predicting the inception and completion of dryout around the perimeter of the horizontal tubes. (C) 2007 Elsevier Ltd. All rights reserved.

Registration of temporal sequences of coronal and sagittal MR images through respiratory patterns

This work discusses the determination of the breathing patterns in time sequence of images obtained from magnetic resonance (MR) and their use in the temporal registration of coronal and sagittal images. The registration is made without the use of any triggering information and any special gas to enhance the contrast. The temporal sequences of images are acquired in free breathing. The real movement of the lung has never been seen directly, as it is totally dependent on its surrounding muscles and collapses without them. The visualization of the lung in motion is an actual topic of research in medicine. The lung movement is not periodic and it is susceptible to variations in the degree of respiration. Compared to computerized tomography (CT), MR imaging involves longer acquisition times and it is preferable because it does not involve radiation. As coronal and sagittal sequences of images are orthogonal to each other, their intersection corresponds to a segment in the three-dimensional space. The registration is based on the analysis of this intersection segment. A time sequence of this intersection segment can be stacked, defining a two-dimension spatio-temporal (2DST) image. The algorithm proposed in this work can detect asynchronous movements of the internal lung structures and lung surrounding organs. It is assumed that the diaphragmatic movement is the principal movement and all the lung structures move almost synchronously. The synchronization is performed through a pattern named respiratory function. This pattern is obtained by processing a 2DST image. An interval Hough transform algorithm searches for synchronized movements with the respiratory function. A greedy active contour algorithm adjusts small discrepancies originated by asynchronous movements in the respiratory patterns. The output is a set of respiratory patterns. Finally, the composition of coronal and sagittal image pairs that are in the same breathing phase is realized by comparing of respiratory patterns originated from diaphragmatic and upper boundary surfaces. When available, the respiratory patterns associated to lung internal structures are also used. The results of the proposed method are compared with the pixel-by-pixel comparison method. The proposed method increases the number of registered pairs representing composed images and allows an easy check of the breathing phase. (C) 2010 Elsevier Ltd. All rights reserved.

Ciratefi: An RST-invariant template matching with extension to color images

Template matching is a technique widely used for finding patterns in digital images. A good template matching should be able to detect template instances that have undergone geometric transformations. In this paper, we proposed a grayscale template matching algorithm named Ciratefi, invariant to rotation, scale, translation, brightness and contrast and its extension to color images. We introduce CSSIM (color structural similarity) for comparing the similarity of two color image patches and use it in our algorithm. We also describe a scheme to determine automatically the appropriate parameters of our algorithm and use pyramidal structure to improve the scale invariance. We conducted several experiments to compare grayscale and color Ciratefis with SIFT, C-color-SIFT and EasyMatch algorithms in many different situations. The results attest that grayscale and color Ciratefis are more accurate than the compared algorithms and that color-Ciratefi outperforms grayscale Ciratefi most of the time. However, Ciratefi is slower than the other algorithms.

Recognition of noisy images by PLL networks

We use networks composed of three phase-locked loops (PLLs), where one of them is the master, for recognizing noisy images. The values of the coupling weights among the PLLs control the noise level which does not affect the successful identification of the input image. Analytical results and numerical tests are presented concerning the scheme performance. (c) 2008 Elsevier B.V. All rights reserved.

Synaptic compensation on Hopfield network: implications for memory rehabilitation

The discrete-time neural network proposed by Hopfield can be used for storing and recognizing binary patterns. Here, we investigate how the performance of this network on pattern recognition task is altered when neurons are removed and the weights of the synapses corresponding to these deleted neurons are divided among the remaining synapses. Five distinct ways of distributing such weights are evaluated. We speculate how this numerical work about synaptic compensation may help to guide experimental studies on memory rehabilitation interventions.

Identification of protein coding regions using the modified Gabor-wavelet transform

An important topic in genomic sequence analysis is the identification of protein coding regions. In this context, several coding DNA model-independent methods based on the occurrence of specific patterns of nucleotides at coding regions have been proposed. Nonetheless, these methods have not been completely suitable due to their dependence on an empirically predefined window length required for a local analysis of a DNA region. We introduce a method based on a modified Gabor-wavelet transform (MGWT) for the identification of protein coding regions. This novel transform is tuned to analyze periodic signal components and presents the advantage of being independent of the window length. We compared the performance of the MGWT with other methods by using eukaryote data sets. The results show that MGWT outperforms all assessed model-independent methods with respect to identification accuracy. These results indicate that the source of at least part of the identification errors produced by the previous methods is the fixed working scale. The new method not only avoids this source of errors but also makes a tool available for detailed exploration of the nucleotide occurrence.

Watershed scale temporal stability of soil water content

The recognition of temporally stable locations with respect to soil water content is of importance for soil water management decisions, especially in sloping land of watersheds. Neutron probe soil water content (0 to 0.8 m), evaluated at 20 dates during a year in the Loess Plateau of China, in a 20 ha watershed dominated by Ust-Sandiic Entisols and Aeolian sandy soils, were used to define their temporal stability through two indices: the standard deviation of relative difference (SDRD) and the mean absolute bias error (MABE). Specific concerns were (a) the relationship of temporal stability with soil depth, (b) the effects of soil texture and land use on temporal stability, and (c) the spatial pattern of the temporal stability. Results showed that temporal stability of soil water content at 0.2 m was significantly weaker than those at the soil depths of 0.6 and 0.8 m. Soil texture can significantly (P<0.05) affect the stability of soil water content except for the existence of an insignificant difference between sandy loam and silt loam textures, while temporal stability of areas covered by bunge needlegrass land was not significantly different from those covered by korshinsk peashrub. Geostatistical analysis showed that the temporal stability was spatially variable in an organized way as inferred by the degree of spatial dependence index. With increasing soil depth, the range of both temporal stability indices showed an increasing trend, being 65.8-120.5 m for SDRD and 148.8-214.1 m for MABE, respectively. This study provides a valuable support for soil water content measurements for soil water management and hydrological applications on sloping land areas. (C) 2010 Elsevier B.V. All rights reserved.

Methylation patterns revealed by MSAP profiling in genetically stable somatic embryogenic cultures of Ocotea catharinensis (Lauraceae)

Ocotea catharinensis is a rare tree species indigenous to the Atlantic rainforest of South America. In spite of its value as a hardwood species, it is in danger of extinction. The species erratically produces seeds showing irregular flowering and slow growth. Therefore, plants are not easily replaced. Tissue culture-based techniques are commonly used for obtaining living material for tree propagation and in vitro preservation. Therefore, a high-frequency somatic embryogenic system was developed for the species. In the present work, the genetic fidelity of cell aggregates and somatic embryos at various stages of in vitro development of O. catharinensis was investigated using RAPD and AFLP markers. Both analyses confirmed the absence of genetic variation in all developmental stages of O. catharinensis embryogenic cultures, verifying that the in vitro system is genetically stable. The cultures were also analyzed for their methylation profiles at 5`-CCGG-3` sites by identifying methylation-sensitive amplification polymorphisms. Some of these markers differentiated cell aggregates from embryo bodies. The sequencing of ten MSAP markers revealed that four sequences showed significant similarity to genes encoding plant proteins. Particularly, the predicted amino acid sequence of the fragment designated as OcEaggHMttc155 was similar to the enzyme 1-aminocyclopropane-1-carboxylate oxidase (ACO), which is involved in the biosynthesis of ethylene, and its expression was reported to occur from the beginning to the intermediate stages of plant embryo development. Here, we suggest that this enzyme is possibly involved in the control of the earliest stages of somatic embryogenesis of O. catharinensis, and an approach to study ACO expression during somatic embryogenesis is proposed.

Spatial patterns of soil attributes and components in a mangrove system in Southeast Brazil (Sao Paulo)

Purpose Among environmental factors governing innumerous processes that are active in estuarine environments, those of edaphic character have received special attention in recent studies. With the objectives of determining the spatial patterns of soil attributes and components across different mangrove forest landscapes and obtaining additional information on the cause-effect relationships between these variables and position within the estuary, we analyzed several soil attributes in 31 mangrove soil profiles from the state of So Paulo (Guaruja, Brazil). Materials and methods Soil samples were collected at low tide along two transects within the CrumahA(0) mangrove forest. Samples were analyzed to determine pH, Eh, salinity, and the percentages of sand, silt, clay, total organic carbon (TOC), and total S. Mineralogy of the clay fraction (< 2 mm) was also studied by X-ray diffraction analysis, and partitioning of solid-phase Fe was performed by sequential extraction. Results and discussion The results obtained indicate important differences in soil composition at different depths and landscape positions, causing variations in physicochemical parameters, clay mineralogy, TOC contents, and iron geochemistry. The results also indicate that physicochemical conditions may vary in terms of different local microtopographies. Soil salinity was determined by relative position in relation to flood tide and transition areas with highlands. The proportions of TOC and total S are conditioned by the sedimentation of organic matter derived from vegetation and by the prevailing redox conditions, which clearly favored intense sulfate reduction in the soils (similar to 80% of the total Fe is Fe-pyrite). Particle-size distribution is conditioned by erosive/deposition processes (present and past) and probably by the positioning of ancient and reworked sandy ridges. The existing physicochemical conditions appear to contribute to the synthesis (smectite) and transformation (kaolinite) of clay minerals. Conclusions The results demonstrate that the position of soils in the estuary greatly affects soil attributes. Differences occur even at small scales (meters), indicating that both edaphic (soil classification, soil mineralogy, and soil genesis) and environmental (contamination and carbon stock) studies should take such variability into account.

Out of the Andes: patterns of diversification in clearwing butterflies

Global biodiversity peaks in the tropical forests of the Andes, a striking geological feature that has likely been instrumental in generating biodiversity by providing opportunities for both vicariant and ecological speciation. However, the role of these mountains in the diversification of insects, which dominate biodiversity, has been poorly explored using phylogenetic methods. Here we study the role of the Andes in the evolution of a diverse Neotropical insect group, the clearwing butterflies. We used dated species-level phylogenies to investigate the time course of speciation and to infer ancestral elevation ranges for two diverse genera. We show that both genera likely originated at middle elevations in the Andes in the Middle Miocene, contrasting with most published results in vertebrates that point to a lowland origin. Although we detected a signature of vicariance caused by the uplift of the Andes at the Miocene-Pliocene boundary, most sister species were parapatric without any obvious vicariant barrier. Combined with an overall decelerating speciation rate, these results suggest an important role for ecological speciation and adaptive radiation, rather than simple vicariance.