LATENT RANK CHANGE DETECTION FOR ANALYSIS OF SPLICE-JUNCTION MICROARRAYS WITH NONLINEAR EFFECTS

Alternative splicing of gene transcripts greatly expands the functional capacity of the genome, and certain splice isoforms may indicate specific disease states such as cancer. Splice junction microarrays interrogate thousands of splice junctions, but data analysis is difficult and error prone because of the increased complexity compared to differential gene expression analysis. We present Rank Change Detection (RCD) as a method to identify differential splicing events based upon a straightforward probabilistic model comparing the over-or underrepresentation of two or more competing isoforms. RCD has advantages over commonly used methods because it is robust to false positive errors due to nonlinear trends in microarray measurements. Further, RCD does not depend on prior knowledge of splice isoforms, yet it takes advantage of the inherent structure of mutually exclusive junctions, and it is conceptually generalizable to other types of splicing arrays or RNA-Seq. RCD specifically identifies the biologically important cases when a splice junction becomes more or less prevalent compared to other mutually exclusive junctions. The example data is from different cell lines of glioblastoma tumors assayed with Agilent microarrays.

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.

The generalized inverse Weibull distribution

The inverse Weibull distribution has the ability to model failure rates which are quite common in reliability and biological studies. A three-parameter generalized inverse Weibull distribution with decreasing and unimodal failure rate is introduced and studied. We provide a comprehensive treatment of the mathematical properties of the new distribution including expressions for the moment generating function and the rth generalized moment. The mixture model of two generalized inverse Weibull distributions is investigated. The identifiability property of the mixture model is demonstrated. For the first time, we propose a location-scale regression model based on the log-generalized inverse Weibull distribution for modeling lifetime data. In addition, we develop some diagnostic tools for sensitivity analysis. Two applications of real data are given to illustrate the potentiality of the proposed regression model.

Residuals for log-Burr XII regression models in survival analysis

In this paper, we compare three residuals to assess departures from the error assumptions as well as to detect outlying observations in log-Burr XII regression models with censored observations. These residuals can also be used for the log-logistic regression model, which is a special case of the log-Burr XII regression model. For different parameter settings, sample sizes and censoring percentages, various simulation studies are performed and the empirical distribution of each residual is displayed and compared with the standard normal distribution. These studies suggest that the residual analysis usually performed in normal linear regression models can be straightforwardly extended to the modified martingale-type residual in log-Burr XII regression models with censored data.

Double generalized linear model for tissue culture proportion data: a Bayesian perspective

Joint generalized linear models and double generalized linear models (DGLMs) were designed to model outcomes for which the variability can be explained using factors and/or covariates. When such factors operate, the usual normal regression models, which inherently exhibit constant variance, will under-represent variation in the data and hence may lead to erroneous inferences. For count and proportion data, such noise factors can generate a so-called overdispersion effect, and the use of binomial and Poisson models underestimates the variability and, consequently, incorrectly indicate significant effects. In this manuscript, we propose a DGLM from a Bayesian perspective, focusing on the case of proportion data, where the overdispersion can be modeled using a random effect that depends on some noise factors. The posterior joint density function was sampled using Monte Carlo Markov Chain algorithms, allowing inferences over the model parameters. An application to a data set on apple tissue culture is presented, for which it is shown that the Bayesian approach is quite feasible, even when limited prior information is available, thereby generating valuable insight for the researcher about its experimental results.

Qualitative assessment of stock prices listed on the Sao Paulo stock exchange: An approach from the perspective of homogeneity analysis

The stock market suffers uncertain relations throughout the entire negotiation process, with different variables exerting direct and indirect influence on stock prices. This study focuses on the analysis of certain aspects that may influence these values offered by the capital market, based on the Brazil Index of the Sao Paulo Stock Exchange (Bovespa), which selects 100 stocks among the most traded on Bovespa in terms of number of trades and financial volume. The selected variables are characterized by the companies` activity area and the business volume in the month of data collection, i.e. April/2007. This article proposes an analysis that joins the accounting view of the stock price variables that can be influenced with the use of multivariate qualitative data analysis. Data were explored through Correspondence Analysis (Anacor) and Homogeneity Analysis (Homals). According to the research, the selected variables are associated with the values presented by the stocks, which become an internal control instrument and a decision-making tool when it comes to choosing investments.

Neural activity changes to emotional stimuli in healthy individuals under chronic use of clomipramine

Previous functional magnetic resonance imaging (fMRI) studies examined neural activity responses to emotive stimuli in healthy individuals after acute/subacute administration of antidepressants. We now report the effects of repeated use of the antidepressant clomipramine on fMRI data acquired during presentation of emotion-provoking and neutral stimuli on healthy volunteers. A total of 12 volunteers were evaluated with fMRI after receiving low doses of clomipramine for 4 weeks and again after 4 weeks of washout. Fear-, happiness-, anger-provoking and neutral pictures from the International Affective Picture System (IAPS) were used. Data analysis was performed with statistical parametric mapping (P < 0.05). Paired t-test comparisons for each condition between medicated and unmedicated states showed, to negative valence paradigms, decrease in brain activity in the amygdala when participants were medicated. We also demonstrated, across both positive and negative valence paradigms, consistent decreases in brain activity in the medicated state in the anterior cingulate gyrus and insula. This is the first report of modulatory effects of repeated antidepressant use on the central representation of somatic states in response to emotions of both negative and positive valences in healthy individuals. Also, our results corroborate findings of antidepressant-induced temporolimbic activity changes to emotion-provoking stimuli obtained in studies of subjects treated acutely with such agents.

Detection of nosocomial malnutrition is improved in Amazon region by a standard clinical nutrition education program

Background: In Brazil hospital malnutrition is highly prevalent. physician awareness of malnutrition is low, and nutrition therapy is underprescribed. One alternative to approach this problem is to educate health care providers in clinical nutrition. The present study aims to evaluate the effect of an intensive education course given to health care professionals and students on the diagnosis ability concerning to hospital malnutrition. Materials and methods: An intervention study based on a clinical nutrition educational program, offered to medical and nursing students and professionals, was held in a hospital of the Amazon region. Participants were evaluated through improvement of diagnostic ability, according to agreement of malnutrition diagnosis using Subjective Global Assessment before and after the workshop, as compared to independent evaluations (Kappa Index, k). To evaluate the impact of the educational intervention on the hospital malnutrition diagnosis, medical records were reviewed for documentation of parameters associated with nutritional status of in-patients. The SPSS statistical software package was used for data analysis. Results: A total of 165 participants concluded the program. The majority (76.4%) were medical and nursing students. Malnutrition diagnosis improved after the course (before k = 0.5; after k = 0.64; p < 0.05). A reduction of false negatives from 50% to 33.3% was observed. During the course, concern of nutritional diagnosis was increased W = 17.57; p < 0.001) and even after the course, improvement on the height measurement was detected chi(2) 12.87;p < 0.001). Conclusions: Clinical nutrition education improved the ability of diagnosing malnutrition; however the primary impact was on medical and nursing students. To sustain diagnostic capacity a clinical nutrition program should be part of health professional curricula and be coupled with continuing education for health care providers.

Risk of Potential Drug-Drug Interactions among Brazilian Elderly A Population-Based, Cross-Sectional Study

Background: Drug-drug interactions (DDIs) are one of the main causes of adverse reactions related to medications, being responsible for up to 23% of hospital admissions. However, only a few studies have evaluated this problem in elderly Brazilians. Objectives: To determine the prevalence of potential DDIs (PDDIs) in community-dwelling elderly people in Brazil, analyse these interactions with regard to severity and clinical implications, and identify associated factors. Methods: A population-based cross-sectional study was carried out involving 2143 elderly (aged 60 years) residents of the metropolitan area of Sao Paulo, Brazil. Data were obtained from the SABE (Saude, Bem estar e Envelhecimento [Health, Well-Being, and Aging]) survey, which is a multicentre study carried out in seven countries of Latin America and the Caribbean, coordinated by the Pan-American Health Organization. PDDIs were analysed using a computerized program and categorized according to level of severity, onset, mechanism and documentation in the literature. The STATA software statistical package was used for data analysis, and logistic regression was conducted to determine whether variables were associated with PDDIs. Results: Analysis revealed that 568 (26.5%) of the elderly population included in the study were taking medications that could lead to a DDI. Almost two-thirds (64.4%) of the elderly population exposed to PDDIs were women, 50.7% were aged >= 75 years, 71.7% reported having fair or poor health and 65.8% took 2-5 medications. A total of 125 different PDDIs were identified; the treatment combination of an ACE inhibitor with a thiazide or loop diuretic (associated with hypotension) was the most frequent cause of PDDIs (n=322 patients; 56.7% of individuals with PDDIs). Analysis of the PDDIs revealed that 70.4% were of moderate severity, 64.8% were supported by good quality evidence and 56.8% were considered of delayed onset. The multivariate analysis showed that the risk of a PDDI was significantly increased among elderly individuals using six or more medications (odds ratio [OR] 3.37) and in patients with hypertension (OR 2.56), diabetes mellitus (OR 1.73) or heart problems (OR 3.36). Conclusions: Approximately one-quarter of the elderly population living in Sao Paulo could be taking two or more potentially interacting medicines. Polypharmacy predisposes elderly individuals to PDDIs. More than half of these drug combinations (57.6%, n = 72) were part of commonly employed treatment regimens and may be responsible for adverse reactions that compromise the safety of elderly individuals, especially at home. Educational initiatives are needed to avoid unnecessary risks.

Is pain relief equally efficacious and free of side effects with repeated doses of oral sucrose in preterm neonates?

The aim of the present study was to examine the efficacy and potential side effects of repeated doses of oral sucrose for pain relief during procedures in NICU. Thirty-three preterm neonates were randomly allocated in blind fashion into two groups, the sucrose group (SG = 17) and the control group (CG = 16). The responses of neonates to pain and distress were assessed during blood collection on four consecutive assessment (ass.) days. For the first assessment, the neonates did not receive any solution before the blood collection procedure. During the next three days, the SG received oral sucrose (25%; 0.5 ml/kg) and the CG received sterile water, 2 min before each minor acute painful procedure. The neonates were evaluated during blood collection each morning. The assessment was divided into five phases: Baseline (BL), Antisepsis (A), Puncture (P), Dressing (D), and Recovery (R). The neonates` facial activity (NFCS), behavioral state, and heart rate were evaluated. The data analysis used cut-off scores for painful and distressful responses. No side effects of using sucrose were detected. There were significantly fewer SG neonates with facial actions signaling pain than CG neonates in P (ass.2) and in A (ass.3). We found significantly fewer SG neonates in the awake state than CG neonates in P (ass.2 and ass.4). There were significantly fewer SG neonates crying during A (ass.2), P (ass.2 and ass.4), and D (ass.3). There was no statistical difference between-groups for physiological response. The efficacy of sucrose was maintained for pain relief in preterm neonates with no side effects. (C) 2007 International Association for the Study of Pain. Published by Elsevier B.V. All rights reserved.

The Stability of the Blood Oxygenation Level-Dependent Functional MRI Response to Motor Tasks Is Altered in Patients With Chronic Ischemic Stroke

Background and Purpose-Functional MRI is a powerful tool to investigate recovery of brain function in patients with stroke. An inherent assumption in functional MRI data analysis is that the blood oxygenation level-dependent (BOLD) signal is stable over the course of the examination. In this study, we evaluated the validity of such assumption in patients with chronic stroke. Methods-Fifteen patients performed a simple motor task with repeated epochs using the paretic and the unaffected hand in separate runs. The corresponding BOLD signal time courses were extracted from the primary and supplementary motor areas of both hemispheres. Statistical maps were obtained by the conventional General Linear Model and by a parametric General Linear Model. Results-Stable BOLD amplitude was observed when the task was executed with the unaffected hand. Conversely, the BOLD signal amplitude in both primary and supplementary motor areas was progressively attenuated in every patient when the task was executed with the paretic hand. The conventional General Linear Model analysis failed to detect brain activation during movement of the paretic hand. However, the proposed parametric General Linear Model corrected the misdetection problem and showed robust activation in both primary and supplementary motor areas. Conclusions-The use of data analysis tools that are built on the premise of a stable BOLD signal may lead to misdetection of functional regions and underestimation of brain activity in patients with stroke. The present data urge the use of caution when relying on the BOLD response as a marker of brain reorganization in patients with stroke. (Stroke. 2010; 41:1921-1926.)

Diagnostic tool to analyse colour-magnitude diagrams of poorly populated stellar concentrations

The dynamical processes that lead to open cluster disruption cause its mass to decrease. To investigate such processes from the observational point of view, it is important to identify open cluster remnants (OCRs), which are intrinsically poorly populated. Due to their nature, distinguishing them from field star fluctuations is still an unresolved issue. In this work, we developed a statistical diagnostic tool to distinguish poorly populated star concentrations from background field fluctuations. We use 2MASS photometry to explore one of the conditions required for a stellar group to be a physical group: to produce distinct sequences in a colour-magnitude diagram (CMD). We use automated tools to (i) derive the limiting radius; (ii) decontaminate the field and assign membership probabilities; (iii) fit isochrones; and (iv) compare object and field CMDs, considering the isochrone solution, in order to verify the similarity. If the object cannot be statistically considered as a field fluctuation, we derive its probable age, distance modulus, reddening and uncertainties in a self-consistent way. As a test, we apply the tool to open clusters and comparison fields. Finally, we study the OCR candidates DoDz 6, NGC 272, ESO 435 SC48 and ESO 325 SC15. The tool is optimized to treat these low-statistic objects and to separate the best OCR candidates for studies on kinematics and chemical composition. The study of the possible OCRs will certainly provide a deep understanding of OCR properties and constraints for theoretical models, including insights into the evolution of open clusters and dissolution rates.

Cross-entropy optimizer: a new tool to study precession in astrophysical jets

Evidence of jet precession in many galactic and extragalactic sources has been reported in the literature. Much of this evidence is based on studies of the kinematics of the jet knots, which depends on the correct identification of the components to determine their respective proper motions and position angles on the plane of the sky. Identification problems related to fitting procedures, as well as observations poorly sampled in time, may influence the follow-up of the components in time, which consequently might contribute to a misinterpretation of the data. In order to deal with these limitations, we introduce a very powerful statistical tool to analyse jet precession: the cross-entropy method for continuous multi-extremal optimization. Only based on the raw data of the jet components (right ascension and declination offsets from the core), the cross-entropy method searches for the precession model parameters that better represent the data. In this work we present a large number of tests to validate this technique, using synthetic precessing jets built from a given set of precession parameters. With the aim of recovering these parameters, we applied the cross-entropy method to our precession model, varying exhaustively the quantities associated with the method. Our results have shown that even in the most challenging tests, the cross-entropy method was able to find the correct parameters within a 1 per cent level. Even for a non-precessing jet, our optimization method could point out successfully the lack of precession.

MODELING VERY LONG BASELINE INTERFEROMETRIC IMAGES WITH THE CROSS-ENTROPY GLOBAL OPTIMIZATION TECHNIQUE

We present a new technique for obtaining model fittings to very long baseline interferometric images of astrophysical jets. The method minimizes a performance function proportional to the sum of the squared difference between the model and observed images. The model image is constructed by summing N(s) elliptical Gaussian sources characterized by six parameters: two-dimensional peak position, peak intensity, eccentricity, amplitude, and orientation angle of the major axis. We present results for the fitting of two main benchmark jets: the first constructed from three individual Gaussian sources, the second formed by five Gaussian sources. Both jets were analyzed by our cross-entropy technique in finite and infinite signal-to-noise regimes, the background noise chosen to mimic that found in interferometric radio maps. Those images were constructed to simulate most of the conditions encountered in interferometric images of active galactic nuclei. We show that the cross-entropy technique is capable of recovering the parameters of the sources with a similar accuracy to that obtained from the very traditional Astronomical Image Processing System Package task IMFIT when the image is relatively simple (e. g., few components). For more complex interferometric maps, our method displays superior performance in recovering the parameters of the jet components. Our methodology is also able to show quantitatively the number of individual components present in an image. An additional application of the cross-entropy technique to a real image of a BL Lac object is shown and discussed. Our results indicate that our cross-entropy model-fitting technique must be used in situations involving the analysis of complex emission regions having more than three sources, even though it is substantially slower than current model-fitting tasks (at least 10,000 times slower for a single processor, depending on the number of sources to be optimized). As in the case of any model fitting performed in the image plane, caution is required in analyzing images constructed from a poorly sampled (u, v) plane.

Consistency in approximate entropy given by a volumetric estimate

Non-linear methods for estimating variability in time-series are currently of widespread use. Among such methods are approximate entropy (ApEn) and sample approximate entropy (SampEn). The applicability of ApEn and SampEn in analyzing data is evident and their use is increasing. However, consistency is a point of concern in these tools, i.e., the classification of the temporal organization of a data set might indicate a relative less ordered series in relation to another when the opposite is true. As highlighted by their proponents themselves, ApEn and SampEn might present incorrect results due to this lack of consistency. In this study, we present a method which gains consistency by using ApEn repeatedly in a wide range of combinations of window lengths and matching error tolerance. The tool is called volumetric approximate entropy, vApEn. We analyze nine artificially generated prototypical time-series with different degrees of temporal order (combinations of sine waves, logistic maps with different control parameter values, random noises). While ApEn/SampEn clearly fail to consistently identify the temporal order of the sequences, vApEn correctly do. In order to validate the tool we performed shuffled and surrogate data analysis. Statistical analysis confirmed the consistency of the method. (C) 2008 Elsevier Ltd. All rights reserved.