CHEMOSTAT, UM SOFTWARE GRATUITO PARA ANÁLISE EXPLORATÓRIA DE DADOS MULTIVARIADOS

The objective of this work was to develop a free access exploratory data analysis software application for academic use that is easy to install and can be handled without user-level programming due to extensive use of chemometrics and its association with applications that require purchased licenses or routines. The developed software, called Chemostat, employs Hierarchical Cluster Analysis (HCA), Principal Component Analysis (PCA), intervals Principal Component Analysis (iPCA), as well as correction methods, data transformation and outlier detection. The data can be imported from the clipboard, text files, ASCII or FT-IR Perkin-Elmer “.sp” files. It generates a variety of charts and tables that allow the analysis of results that can be exported in several formats. The main features of the software were tested using midinfrared and near-infrared spectra in vegetable oils and digital images obtained from different types of commercial diesel. In order to validate the software results, the same sets of data were analyzed using Matlab© and the results in both applications matched in various combinations. In addition to the desktop version, the reuse of algorithms allowed an online version to be provided that offers a unique experience on the web. Both applications are available in English.

Efficiency of distinct data mining algorithms for classifying stress level in piglets from their vocalization

Among the challenges of pig farming in today's competitive market, there is factor of the product traceability that ensures, among many points, animal welfare. Vocalization is a valuable tool to identify situations of stress in pigs, and it can be used in welfare records for traceability. The objective of this work was to identify stress in piglets using vocalization, calling this stress on three levels: no stress, moderate stress, and acute stress. An experiment was conducted on a commercial farm in the municipality of Holambra, São Paulo State , where vocalizations of twenty piglets were recorded during the castration procedure, and separated into two groups: without anesthesia and local anesthesia with lidocaine base. For the recording of acoustic signals, a unidirectional microphone was connected to a digital recorder, in which signals were digitized at a frequency of 44,100 Hz. For evaluation of sound signals, Praat® software was used, and different data mining algorithms were applied using Weka® software. The selection of attributes improved model accuracy, and the best attribute selection was used by applying Wrapper method, while the best classification algorithms were the k-NN and Naive Bayes. According to the results, it was possible to classify the level of stress in pigs through their vocalization.

Spatial variability of leaf nutrient contents in a drip irrigated citrus orchard

This study aimed to evaluate the spatial variability of leaf content of macro and micronutrients. The citrus plants orchard with 5 years of age, planted at regular intervals of 8 x 7 m, was managed under drip irrigation. Leaf samples were collected from each plant to be analyzed in the laboratory. Data were analyzed using the software R, version 2.5.1 Copyright (C) 2007, along with geostatistics package GeoR. All contents of macro and micronutrients studied were adjusted to normal distribution and showed spatial dependence.The best-fit models, based on the likelihood, for the macro and micronutrients were the spherical and matern. It is suggest for the macronutrients nitrogen, phosphorus, potassium, calcium, magnesium and sulfur the minimum distances between samples of 37; 58; 29; 63; 46 and 15 m respectively, while for the micronutrients boron, copper, iron, manganese and zinc, the distances suggests are 29; 9; 113; 35 and 14 m, respectively.

Using data mining to identify factors that influence the degree of leg injuries in broilers

Locomotor problems prevent the bird to move freely, jeopardizing the welfare and productivity, besides generating injuries on the legs of chickens. The objective of this study was to evaluate the influence of age, use of vitamin D, the asymmetry of limbs and gait score, the degree of leg injuries in broilers, using data mining. The analysis was performed on a data set obtained from a field experiment in which it was used two groups of birds with 30 birds each, a control group and one treated with vitamin D. It was evaluated the gait score, the asymmetry between the right and left toes, and the degree of leg injuries. The Weka ® software was used in data mining. In particular, C4.5 algorithm (also known as J48 in Weka environment) was used for the generation of a decision tree. The results showed that age is the factor that most influences the degree of leg injuries and that the data from assessments of gait score were not reliable to estimate leg weakness in broilers.

Software for calculation of reservoir active capacity with the sequent-peak algorithm

It is presented a software developed with Delphi programming language to compute the reservoir's annual regulated active storage, based on the sequent-peak algorithm. Mathematical models used for that purpose generally require extended hydrological series. Usually, the analysis of those series is performed with spreadsheets or graphical representations. Based on that, it was developed a software for calculation of reservoir active capacity. An example calculation is shown by 30-years (from 1977 to 2009) monthly mean flow historical data, from Corrente River, located at São Francisco River Basin, Brazil. As an additional tool, an interface was developed to manage water resources, helping to manipulate data and to point out information that it would be of interest to the user. Moreover, with that interface irrigation districts where water consumption is higher can be analyzed as a function of specific seasonal water demands situations. From a practical application, it is possible to conclude that the program provides the calculation originally proposed. It was designed to keep information organized and retrievable at any time, and to show simulation on seasonal water demands throughout the year, contributing with the elements of study concerning reservoir projects. This program, with its functionality, is an important tool for decision making in the water resources management.

Vocalization data mining for estimating swine stress conditions

This study aimed to identify differences in swine vocalization pattern according to animal gender and different stress conditions. A total of 150 barrow males and 150 females (Dalland® genetic strain), aged 100 days, were used in the experiment. Pigs were exposed to different stressful situations: thirst (no access to water), hunger (no access to food), and thermal stress (THI exceeding 74). For the control treatment, animals were kept under a comfort situation (animals with full access to food and water, with environmental THI lower than 70). Acoustic signals were recorded every 30 minutes, totaling six samples for each stress situation. Afterwards, the audios were analyzed by Praat® 5.1.19 software, generating a sound spectrum. For determination of stress conditions, data were processed by WEKA® 3.5 software, using the decision tree algorithm C4.5, known as J48 in the software environment, considering cross-validation with samples of 10% (10-fold cross-validation). According to the Decision Tree, the acoustic most important attribute for the classification of stress conditions was sound Intensity (root node). It was not possible to identify, using the tested attributes, the animal gender by vocal register. A decision tree was generated for recognition of situations of swine hunger, thirst, and heat stress from records of sound intensity, Pitch frequency, and Formant 1.

Chromosomal abnormalities in couples with recurrent first trimester abortions

PURPOSE: To investigate the prevalence of chromosomal abnormalities in couples with two or more recurrent first trimester miscarriages of unknown cause. METHODS: The study was conducted on 151 women and 94 partners who had an obstetrical history of two or more consecutive first trimester abortions (1-12 weeks of gestation). The controls were 100 healthy women without a history of pregnancy loss. Chromosomal analysis was performed on peripheral blood lymphocytes cultured for 72 hours, using Trypsin-Giemsa (GTG) banding. In all cases, at least 30 metaphases were analyzed and 2 karyotypes were prepared, using light microscopy. The statistical analysis was performed using the Student t-test for normally distributed data and the Mann-Whitney test for non-parametric data. The Kruskal-Wallis test or Analysis of Variance was used to compare the mean values between three or more groups. The software used was Statistical Package for the Social Sciences (SPSS), version 17.0. RESULTS: The frequency of chromosomal abnormalities in women with recurrent miscarriages was 7.3%, including 4.7% with X-chromosome mosaicism, 2% with reciprocal translocations and 0.6% with Robertsonian translocations. A total of 2.1% of the partners of women with recurrent miscarriages had chromosomal abnormalities, including 1% with X-chromosome mosaicism and 1% with inversions. Among the controls, 1% had mosaicism. CONCLUSION: An association between chromosomal abnormalities and recurrent miscarriage in the first trimester of pregnancy (OR=7.7; 95%CI 1.2--170.5) was observed in the present study. Etiologic identification of genetic factors represents important clinical information for genetic counseling and orientation of the couple about the risk for future pregnancies and decreases the number of investigations needed to elucidate the possible causes of miscarriages.

Object-oriented Programming Applied to the Development of Structural Analysis and Optimization Software

This paper presents the development of a two-dimensional interactive software environment for structural analysis and optimization based on object-oriented programming using the C++ language. The main feature of the software is the effective integration of several computational tools into graphical user interfaces implemented in the Windows-98 and Windows-NT operating systems. The interfaces simplify data specification in the simulation and optimization of two-dimensional linear elastic problems. NURBS have been used in the software modules to represent geometric and graphical data. Extensions to the analysis of three-dimensional problems have been implemented and are also discussed in this paper.

Neutron tomography using projection data obtained by Monte Carlo simulation for nondestructive evaluation

This work present the application of a computer package for generating of projection data for neutron computerized tomography, and in second part, discusses an application of neutron tomography, using the projection data obtained by Monte Carlo technique, for the detection and localization of light materials such as those containing hydrogen, concealed by heavy materials such as iron and lead. For tomographic reconstructions of the samples simulated use was made of only six equal projection angles distributed between 0º and 180º, with reconstruction making use of an algorithm (ARIEM), based on the principle of maximum entropy. With the neutron tomography it was possible to detect and locate polyethylene and water hidden by lead and iron (with 1cm-thick). Thus, it is demonstrated that thermal neutrons tomography is a viable test method which can provide important interior information about test components, so, extremely useful in routine industrial applications.

Low-cost automatic activity data recording system

We describe a low-cost, high quality device capable of monitoring indirect activity by detecting touch-release events on a conducting surface, i.e., the animal's cage cover. In addition to the detecting sensor itself, the system includes an IBM PC interface for prompt data storage. The hardware/software design, while serving for other purposes, is used to record the circadian activity rhythm pattern of rats with time in an automated computerized fashion using minimal cost computer equipment (IBM PC XT). Once the sensor detects a touch-release action of the rat in the upper portion of the cage, the interface sends a command to the PC which records the time (hours-minutes-seconds) when the activity occurred. As a result, the computer builds up several files (one per detector/sensor) containing a time list of all recorded events. Data can be visualized in terms of actograms, indicating the number of detections per hour, and analyzed by mathematical tools such as Fast Fourier Transform (FFT) or cosinor. In order to demonstrate method validation, an experiment was conducted on 8 Wistar rats under 12/12-h light/dark cycle conditions (lights on at 7:00 a.m.). Results show a biological validation of the method since it detected the presence of circadian activity rhythm patterns in the behavior of the rats

Thermodynamic evaluation and modeling of proton and water exchange associated with benzamidine and berenil binding to ß-trypsin

Serine-proteases are involved in vital processes in virtually all species. They are important targets for researchers studying the relationships between protein structure and activity, for the rational design of new pharmaceuticals. Trypsin was used as a model to assess a possible differential contribution of hydration water to the binding of two synthetic inhibitors. Thermodynamic parameters for the association of bovine ß-trypsin (homogeneous material, observed 23,294.4 ± 0.2 Da, theoretical 23,292.5 Da) with the inhibitors benzamidine and berenil at pH 8.0, 25ºC and with 25 mM CaCl2, were determined using isothermal titration calorimetry and the osmotic stress method. The association constant for berenil was about 12 times higher compared to the one for benzamidine (binding constants are K = 596,599 ± 25,057 and 49,513 ± 2,732 M-1, respectively; the number of binding sites is the same for both ligands, N = 0.99 ± 0.05). Apparently the driving force responsible for this large difference of affinity is not due to hydrophobic interactions because the variation in heat capacity (DCp), a characteristic signature of these interactions, was similar in both systems tested (-464.7 ± 23.9 and -477.1 ± 86.8 J K-1 mol-1 for berenil and benzamidine, respectively). The results also indicated that the enzyme has a net gain of about 21 water molecules regardless of the inhibitor tested. It was shown that the difference in affinity could be due to a larger number of interactions between berenil and the enzyme based on computational modeling. The data support the view that pharmaceuticals derived from benzamidine that enable hydrogen bond formation outside the catalytic binding pocket of ß-trypsin may result in more effective inhibitors.

Mechanisms underlying gas exchange alterations in an experimental model of pulmonary embolism

The aim of the present study was to determine the ventilation/perfusion ratio that contributes to hypoxemia in pulmonary embolism by analyzing blood gases and volumetric capnography in a model of experimental acute pulmonary embolism. Pulmonary embolization with autologous blood clots was induced in seven pigs weighing 24.00 ± 0.6 kg, anesthetized and mechanically ventilated. Significant changes occurred from baseline to 20 min after embolization, such as reduction in oxygen partial pressures in arterial blood (from 87.71 ± 8.64 to 39.14 ± 6.77 mmHg) and alveolar air (from 92.97 ± 2.14 to 63.91 ± 8.27 mmHg). The effective alveolar ventilation exhibited a significant reduction (from 199.62 ± 42.01 to 84.34 ± 44.13) consistent with the fall in alveolar gas volume that effectively participated in gas exchange. The relation between the alveolar ventilation that effectively participated in gas exchange and cardiac output (V Aeff/Q ratio) also presented a significant reduction after embolization (from 0.96 ± 0.34 to 0.33 ± 0.17 fraction). The carbon dioxide partial pressure increased significantly in arterial blood (from 37.51 ± 1.71 to 60.76 ± 6.62 mmHg), but decreased significantly in exhaled air at the end of the respiratory cycle (from 35.57 ± 1.22 to 23.15 ± 8.24 mmHg). Exhaled air at the end of the respiratory cycle returned to baseline values 40 min after embolism. The arterial to alveolar carbon dioxide gradient increased significantly (from 1.94 ± 1.36 to 37.61 ± 12.79 mmHg), as also did the calculated alveolar (from 56.38 ± 22.47 to 178.09 ± 37.46 mL) and physiological (from 0.37 ± 0.05 to 0.75 ± 0.10 fraction) dead spaces. Based on our data, we conclude that the severe arterial hypoxemia observed in this experimental model may be attributed to the reduction of the V Aeff/Q ratio. We were also able to demonstrate that V Aeff/Q progressively improves after embolization, a fact attributed to the alveolar ventilation redistribution induced by hypocapnic bronchoconstriction.

Effect of pentoxifylline on lung inflammation and gas exchange in a sepsis-induced acute lung injury model

Experimental models of sepsis-induced pulmonary alterations are important for the study of pathogenesis and for potential intervention therapies. The objective of the present study was to characterize lung dysfunction (low PaO2 and high PaCO2, and increased cellular infiltration, protein extravasation, and malondialdehyde (MDA) production assessed in bronchoalveolar lavage) in a sepsis model consisting of intraperitoneal (ip) injection of Escherichia coli and the protective effects of pentoxifylline (PTX). Male Wistar rats (weighing between 270 and 350 g) were injected ip with 10(7) or 10(9) CFU/100 g body weight or saline and samples were collected 2, 6, 12, and 24 h later (N = 5 each group). PaO2, PaCO2 and pH were measured in blood, and cellular influx, protein extravasation and MDA concentration were measured in bronchoalveolar lavage. In a second set of experiments either PTX or saline was administered 1 h prior to E. coli ip injection (N = 5 each group) and the animals were observed for 6 h. Injection of 10(7) or 10(9) CFU/100 g body weight of E. coli induced acidosis, hypoxemia, and hypercapnia. An increased (P < 0.05) cell influx was observed in bronchoalveolar lavage, with a predominance of neutrophils. Total protein and MDA concentrations were also higher (P < 0.05) in the septic groups compared to control. A higher tumor necrosis factor-alpha (P < 0.05) concentration was also found in these animals. Changes in all parameters were more pronounced with the higher bacterial inoculum. PTX administered prior to sepsis reduced (P < 0.05) most functional alterations. These data show that an E. coli ip inoculum is a good model for the induction of lung dysfunction in sepsis, and suitable for studies of therapeutic interventions.

Pipeline for macro- and microarray analyses

The pipeline for macro- and microarray analyses (PMmA) is a set of scripts with a web interface developed to analyze DNA array data generated by array image quantification software. PMmA is designed for use with single- or double-color array data and to work as a pipeline in five classes (data format, normalization, data analysis, clustering, and array maps). It can also be used as a plugin in the BioArray Software Environment, an open-source database for array analysis, or used in a local version of the web service. All scripts in PMmA were developed in the PERL programming language and statistical analysis functions were implemented in the R statistical language. Consequently, our package is a platform-independent software. Our algorithms can correctly select almost 90% of the differentially expressed genes, showing a superior performance compared to other methods of analysis. The pipeline software has been applied to 1536 expressed sequence tags macroarray public data of sugarcane exposed to cold for 3 to 48 h. PMmA identified thirty cold-responsive genes previously unidentified in this public dataset. Fourteen genes were up-regulated, two had a variable expression and the other fourteen were down-regulated in the treatments. These new findings certainly were a consequence of using a superior statistical analysis approach, since the original study did not take into account the dependence of data variability on the average signal intensity of each gene. The web interface, supplementary information, and the package source code are available, free, to non-commercial users at http://ipe.cbmeg.unicamp.br/pub/PMmA.

Manual and semi-automatic quantification of in vivo ¹H-MRS data for the classification of human primary brain tumors

In vivo proton magnetic resonance spectroscopy (¹H-MRS) is a technique capable of assessing biochemical content and pathways in normal and pathological tissue. In the brain, ¹H-MRS complements the information given by magnetic resonance images. The main goal of the present study was to assess the accuracy of ¹H-MRS for the classification of brain tumors in a pilot study comparing results obtained by manual and semi-automatic quantification of metabolites. In vivo single-voxel ¹H-MRS was performed in 24 control subjects and 26 patients with brain neoplasms that included meningiomas, high-grade neuroglial tumors and pilocytic astrocytomas. Seven metabolite groups (lactate, lipids, N-acetyl-aspartate, glutamate and glutamine group, total creatine, total choline, myo-inositol) were evaluated in all spectra by two methods: a manual one consisting of integration of manually defined peak areas, and the advanced method for accurate, robust and efficient spectral fitting (AMARES), a semi-automatic quantification method implemented in the jMRUI software. Statistical methods included discriminant analysis and the leave-one-out cross-validation method. Both manual and semi-automatic analyses detected differences in metabolite content between tumor groups and controls (P < 0.005). The classification accuracy obtained with the manual method was 75% for high-grade neuroglial tumors, 55% for meningiomas and 56% for pilocytic astrocytomas, while for the semi-automatic method it was 78, 70, and 98%, respectively. Both methods classified all control subjects correctly. The study demonstrated that ¹H-MRS accurately differentiated normal from tumoral brain tissue and confirmed the superiority of the semi-automatic quantification method.