Robust integration of real time optimization with linear model predictive control

Here, we study the stable integration of real time optimization (RTO) with model predictive control (MPC) in a three layer structure. The intermediate layer is a quadratic programming whose objective is to compute reachable targets to the MPC layer that lie at the minimum distance to the optimum set points that are produced by the RTO layer. The lower layer is an infinite horizon MPC with guaranteed stability with additional constraints that force the feasibility and convergence of the target calculation layer. It is also considered the case in which there is polytopic uncertainty in the steady state model considered in the target calculation. The dynamic part of the MPC model is also considered unknown but it is assumed to be represented by one of the models of a discrete set of models. The efficiency of the methods presented here is illustrated with the simulation of a low order system. (C) 2010 Elsevier Ltd. All rights reserved.

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.

A general hazard model for lifetime data in the presence of cure rate

Historically, the cure rate model has been used for modeling time-to-event data within which a significant proportion of patients are assumed to be cured of illnesses, including breast cancer, non-Hodgkin lymphoma, leukemia, prostate cancer, melanoma, and head and neck cancer. Perhaps the most popular type of cure rate model is the mixture model introduced by Berkson and Gage [1]. In this model, it is assumed that a certain proportion of the patients are cured, in the sense that they do not present the event of interest during a long period of time and can found to be immune to the cause of failure under study. In this paper, we propose a general hazard model which accommodates comprehensive families of cure rate models as particular cases, including the model proposed by Berkson and Gage. The maximum-likelihood-estimation procedure is discussed. A simulation study analyzes the coverage probabilities of the asymptotic confidence intervals for the parameters. A real data set on children exposed to HIV by vertical transmission illustrates the methodology.

Bovine gene polymorphisms related to fat deposition and meat tenderness

Leptin, thyroglobulin and diacylglycerol O-acyltransferase play important roles in fat metabolism. Fat deposition has an influence on meat quality and consumers' choice. The aim of this study was to determine allele and genotype frequencies of polymorphisms of the bovine genes, which encode leptin (LEP), thyroglobulin (TG) and diacylglycerol O-acyltransferase (DGAT1). A further objective was to establish the effects of these polymorphisms on meat characteristics. We genotyped 147 animals belonging to the Nelore (Bos indicus), Canchim (5/8 Bos taurus + 3/8 Bos indicus), Rubia Gallega X Nelore (1/2 Bos taurus + 1/2 Bos indicus), Brangus Three-way cross (9/16 Bos taurus + 7/16 Bos indicus) and Braunvieh Three-way cross (3/4 Bos taurus + 1/4 Bos indicus) breeds. Backfat thickness, total lipids, marbling score, ribeye area and shear force were fitted, using the General Linear Model (GLM) procedure of the SAS software. The least square means of genotypes and genetic groups were compared using Tukey's test. Allele frequencies vary among the genetic groups, depending on Bos indicus versus Bos taurus influence. The LEP polymorphism segregates in pure Bos indicus Nelore animals, which is a new finding. The T allele of TG is fixed in Nelore, and DGAT1 segregates in all groups, but the frequency of allele A is lower in Nelore animals. The results showed no association between the genotypes and traits studied, but a genetic group effect on these traits was found. So, the genetic background remains relevant for fat deposition and meat tenderness, but the gene markers developed for Bos taurus may be insufficient for Bos indicus.

Viability of fibroblasts cultured under nutritional stress irradiated with red laser, infrared laser, and red light-emitting diode

Phototherapy is noninvasive, painless and has no known side effect. However, for its incorporation into clinical practice, more well-designed studies are necessary to define optimal parameters for its application. The viability of fibroblasts cultured under nutritional stress irradiated with either a red laser, an infrared laser, or a red light-emitting diode (LED) was analyzed. Irradiation parameters were: red laser (660 nm, 40 mW, 1 W/cm(2)), infrared laser (780 nm, 40 mW, 1 W/cm(2)), and red LED (637 +/- 15 nm, 40 mW, 1 W/cm(2)). All applications were punctual and performed with a spot with 0.4 mm(2) of diameter for 4 or 8 s. The Kruskal-Wallis test and analysis of variance of the general linear model (p <= 0.05) were used for statistical analysis. After 72 h, phototherapy with low-intensity laser and LED showed no toxicity at the cellular level. It even stimulated methylthiazol tetrazolium assay (MTT) conversion and neutral red uptake of fibroblasts cultured under nutritional stress, especially in the group irradiated with infrared laser (p = 0.004 for MTT conversion and p < 0.001 for neutral red uptake). Considering the parameters and protocol of phototherapy used, it can be concluded that phototherapy stimulated the viability of fibroblasts cultured under nutritional deficit resembling those found in traumatized tissue in which cell viability is reduced. (C) 2011 Society of Photo-Optical Instrumentation Engineers (SPIE). [DOI: 10.1117/1.3602850]

Prevalence, magnitude, and methods of rapid weight loss among judo competitors

ARTIOLI, G. G., B. GUALANO, E. FRANCHINI, F. B. SCAGLIUSI, M. TAKESIAN, M. FUCHS, and A. H. LANCHA. Prevalence, Magnitude, and Methods of Rapid Weight Loss among Judo Competitors. Med. Sci. Sports Exerc., Vol. 42, No. 3, pp. 436-442, 2010. Purpose: To identify the prevalence, magnitude, and methods of rapid weight loss among judo competitors. Methods: Athletes (607 males and 215 females; age = 19.3 +/- 5.3 yr, weight = 70 +/- 7.5 kg, height = 170.6 +/- 9.8 cm) completed a previously validated questionnaire developed to evaluate rapid weight loss in judo athletes, which provides a score. The higher the score obtained, the more aggressive the weight loss behaviors. Data were analyzed using descriptive statistics and frequency analyses. Mean scores obtained in the questionnaire were used to compare specific groups of athletes using, when appropriate, Mann-Whitney U-test or general linear model one-way ANOVA followed by Tamhane post hoc test. Results: Eighty-six percent of athletes reported that have already lost weight to compete. When heavyweights are excluded, this percentage rises to 89%. Most athletes reported reductions of up to 5% of body weight (mean +/- SD: 2.5 +/- 2.3%). The most weight ever lost was 2%-5%, whereas a great part of athletes reported reductions of 5%-10% (mean +/- SD: 6 +/- 4%). The number of reductions underwent in a season was 3 +/- 5. The reductions usually occurred within 7 +/- 7 d. Athletes began cutting weight at 12.6 +/- 6.1 yr. No significant differences were found in the score obtained by male versus female athletes as well as by athletes from different weight classes. Elite athletes scored significantly higher in the questionnaire than nonelite. Athletes who began cutting weight earlier also scored higher than those who began later. Conclusions: Rapid weight loss is highly prevalent in judo competitors. The level of aggressiveness in weight management behaviors seems to not be influenced by the gender or by the weight class, but it seems to be influenced by competitive level and by the age at which athletes began cutting weight.

Wavelet correlation between subjects: A time-scale data driven analysis for brain mapping using fMRI

Functional magnetic resonance imaging (fMRI) based on BOLD signal has been used to indirectly measure the local neural activity induced by cognitive tasks or stimulation. Most fMRI data analysis is carried out using the general linear model (GLM), a statistical approach which predicts the changes in the observed BOLD response based on an expected hemodynamic response function (HRF). In cases when the task is cognitively complex or in cases of diseases, variations in shape and/or delay may reduce the reliability of results. A novel exploratory method using fMRI data, which attempts to discriminate between neurophysiological signals induced by the stimulation protocol from artifacts or other confounding factors, is introduced in this paper. This new method is based on the fusion between correlation analysis and the discrete wavelet transform, to identify similarities in the time course of the BOLD signal in a group of volunteers. We illustrate the usefulness of this approach by analyzing fMRI data from normal subjects presented with standardized human face pictures expressing different degrees of sadness. The results show that the proposed wavelet correlation analysis has greater statistical power than conventional GLM or time domain intersubject correlation analysis. (C) 2010 Elsevier B.V. All rights reserved.

Gamma ventral capsulotomy for treatment of resistant obsessive-compulsive disorder: A structural MRI pilot prospective study

Objective: The purpose of this study was to investigate regional structural abnormalities in the brains of five patients with refractory obsessive-compulsive disorder (OCD) submitted to gamma ventral capsulotomy. Methods: We acquired morphometric magnetic resonance imaging (MRI) data before and after 1 year of radiosurgery using a 1.5-T MRI scanner. Images were spatially normalized and segmented using optimized voxel-based morphometry (VBM) methods. Voxelwise statistical comparisons between pre- and post-surgery MRI scans were performed using a general linear model. Findings in regions predicted a priori to show volumetric changes (orbitofrontal cortex, anterior cingulate gyrus, basal ganglia and thalamus) were reported as significant if surpassing a statistical threshold of p<0.001 (uncorrected for multiple comparisons). Results: We detected a significant regional postoperative increase in gray matter volume in the right inferior frontal gyri (Brodmann area 47, BA47) when comparing all patients pre and postoperatively. Conclusions: Our results support the current theory of frontal-striatal-thalamic-cortical (FSTC) circuitry involvement in OCD pathogenesis. Gamma ventral capsulotomy is associated with neurobiological changes in the inferior orbitofrontal cortex in refractory OCD patients. (C) 2008 Elsevier Ireland Ltd. All rights reserved.

Regional gray matter abnormalities in panic disorder: A voxel-based morphometry study

Although abnonnalities in brain structures involved in the neurobiology of fear and anxiety have been implicated in the pathophysiology of panic disorder (PD), relatively few studies have made use of voxel-based morphometry (VBM) magnetic resonance imaging (MRI) to determine structural brain abnormalities in PD. We have assessed gray matter volume in 19 PD patients and 20 healthy volunteers using VBM. Images were acquired using a 1.5 T MRI scanner, and were spatially normalized and segmented using optimized VBM. Statistical comparisons were performed using the general linear model. A relative increase in gay matter volume was found in the left insula of PD patients compared with controls. Additional structures showing differential increases were the left superior temporal gyrus, the midbrain, and the pons. A relative gray matter deficit was found in the right anterior cingulate cortex. The insula and anterior cingulate abnormalities may be relevant to the pathophysiology of PD, since these structures participate in the evaluation process that ascribes negative emotional meaning to potentially distressing cognitive and interoceptive sensory information. The abnormal brain stem structures may be involved in the generation of panic attacks. (C) 2007 Elsevier Ireland Ltd. All rights reserved.

Forefoot Deformity in Diabetic Neuropathic Individuals and its Role in Pressure Distribution and Gait

Background. Foot deformities have been related to diabetic neuropathy progression but their influence on plantar distribution during dynamic tasks is not completely understood. The purpose of the present study was to investigate the influence of metatarsal head prominence and claw toes on regional plantar pressures during gait in patients with diabetic neuropathy Methods Seventy-one adults participated in this study categorized into three groups: a control group (CG, n = 32), patients with diabetic neuropathy without any foot deformities (DG, n = 20), and patients with diabetic neuropathy with metatarsal head prominence and/or claw toes (DMHG, n = 19). Plantar pressure variables (contact area, peak pressure, and maximum mean pressure) were evaluated during gait on rearfoot, midfoot, and forefoot using capacitive insoles (Pedar-X System, Novel Inc., Munich, Germany). A general linear model was applied to repeatedly measure and analyze variance relationships between groups and areas. Results. DMHG. presented larger contact areas at the forefoot and midfoot along with higher peak pressure at the rearfoot compared to the other two groups The DG showed higher mean pressure at the midfoot compared to:the other two groups. Conclusion. The coexistence of diabetic neuropathy and metatarsal head prominence in addition to claw toes, resulted in overloading the rearfoot and enhancing the contact area of forefoot and midfoot while walking. This plantar pressure distribution is a result of a different coordination pattern adopted in order to reduce plantar loads at the anterior parts of the foot that were structurally altered. Patients with diabetic neuropathy without any forefoot deformities presented a different plantar pressure distribution than patients with deformities suggesting that both neuropathy and structural foot alterations can influence foot rollover mechanisms.

Role of ankle mobility in foot rollover during gait in individuals with diabetic neuropathy

Background: The purpose of this study was to investigate the ankle range of motion during neuropathic gait and its influence on plantar pressure distribution in two phases during stance: at heel-strike and at push-off. Methods: Thirty-one adults participated in this study (control group, n = 16; diabetic neuropathic group, n = 15). Dynamic ankle range of motion (electrogoniometer) and plantar pressures (PEDAR-X system) were acquired synchronously during walking. Plantar pressures were evaluated at rearfoot. midfoot and forefoot during the two phases of stance. General linear model repeated measures analysis of variance was applied to investigate relationships between groups, areas and stance phases. Findings: Diabetic neuropathy patients walked using a smaller ankle range of motion in stance phase and smaller ankle flexion at heel-strike (P = 0.0005). Peak pressure and pressure-time integral values were higher in the diabetic group in the midfoot at push-off phase when compared to heel-strike phase. On the other hand, the control group showed similar values of peak pressure in midfoot during both stance phases. Interpretation: The ankle mobility reduction observed could be associated to altered plantar pressure distribution observed in neuropathic subjects. Results demonstrated that midfoot and forefoot play a different role in subjects with neuropathy by receiving higher loads at push-off phase that are probably due to smaller ankle flexion at stance phase. This may explain the higher loads in anterior areas of the foot observed in diabetic neuropathy subjects and confirm an inadequate foot rollover associated to the smaller ankle range of motion at the heel-strike phase. (C) 2009 Elsevier Ltd. All rights reserved.

Detection of Auditory Cortex Activity by fMRI Using a Dependent Component Analysis

Functional MRI (fMRI) data often have low signal-to-noise-ratio (SNR) and are contaminated by strong interference from other physiological sources. A promising tool for extracting signals, even under low SNR conditions, is blind source separation (BSS), or independent component analysis (ICA). BSS is based on the assumption that the detected signals are a mixture of a number of independent source signals that are linearly combined via an unknown mixing matrix. BSS seeks to determine the mixing matrix to recover the source signals based on principles of statistical independence. In most cases, extraction of all sources is unnecessary; instead, a priori information can be applied to extract only the signal of interest. Herein we propose an algorithm based on a variation of ICA, called Dependent Component Analysis (DCA), where the signal of interest is extracted using a time delay obtained from an autocorrelation analysis. We applied such method to inspect functional Magnetic Resonance Imaging (fMRI) data, aiming to find the hemodynamic response that follows neuronal activation from an auditory stimulation, in human subjects. The method localized a significant signal modulation in cortical regions corresponding to the primary auditory cortex. The results obtained by DCA were also compared to those of the General Linear Model (GLM), which is the most widely used method to analyze fMRI datasets.

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.)

Coactivation of the shoulder and arm muscles during closed kinetic chain exercises on an unstable surface

Introduction: The purpose of this study was to compare the electromyography index of muscle coactivation of the following muscle pairs: posterior deltoid and pectoralis major (PD/PM); triceps brachii and biceps brachii (TB/BB); and serratus anterior and upper trapezius (SA/UT) during three different closed kinetic chain exercises (wall-press, bench-press and push-up) on an unstable surface at the maximal load. Methods: A total of 20 healthy sedentary men participated in the study. Integral linear values were obtained from three sustained contractions of six seconds each for the three proposed exercises. Mean coactivation index values were compared using the mixed-effects linear model, with a five percent significance level. Results: Electromyography indexes of muscle coactivation showed significant differences for the PD/PM and TB/BB muscle pairs. No differences were found between exercises for the SA/UT muscle pair. Conclusion: Our results seem to differ from those of previous studies, which reported that the similarity in exercises performed is responsible for the comparable muscle activation levels.

Non-linear boundary element formulation with tangent operator to analyse crack propagation in quasi-brittle materials

This work deals with analysis of cracked structures using BEM. Two formulations to analyse the crack growth process in quasi-brittle materials are discussed. They are based on the dual formulation of BEM where two different integral equations are employed along the opposite sides of the crack surface. The first presented formulation uses the concept of constant operator, in which the corrections of the nonlinear process are made only by applying appropriate tractions along the crack surfaces. The second presented BEM formulation to analyse crack growth problems is an implicit technique based on the use of a consistent tangent operator. This formulation is accurate, stable and always requires much less iterations to reach the equilibrium within a given load increment in comparison with the classical approach. Comparison examples of classical problem of crack growth are shown to illustrate the performance of the two formulations. (C) 2009 Elsevier Ltd. All rights reserved.