Robust linear mixed models with skew-normal independent distributions from a Bayesian perspective

Linear mixed models were developed to handle clustered data and have been a topic of increasing interest in statistics for the past 50 years. Generally. the normality (or symmetry) of the random effects is a common assumption in linear mixed models but it may, sometimes, be unrealistic, obscuring important features of among-subjects variation. In this article, we utilize skew-normal/independent distributions as a tool for robust modeling of linear mixed models under a Bayesian paradigm. The skew-normal/independent distributions is an attractive class of asymmetric heavy-tailed distributions that includes the skew-normal distribution, skew-t, skew-slash and the skew-contaminated normal distributions as special cases, providing an appealing robust alternative to the routine use of symmetric distributions in this type of models. The methods developed are illustrated using a real data set from Framingham cholesterol study. (C) 2009 Elsevier B.V. All rights reserved.

Leverage analysis for linear mixed models

We consider a generalized leverage matrix useful for the identification of influential units and observations in linear mixed models and show how a decomposition of this matrix may be employed to identify high leverage points for both the marginal fitted values and the random effect component of the conditional fitted values. We illustrate the different uses of the two components of the decomposition with a simulated example as well as with a real data set.

Likelihood ratio tests for variance components in linear mixed models

Although the asymptotic distributions of the likelihood ratio for testing hypotheses of null variance components in linear mixed models derived by Stram and Lee [1994. Variance components testing in longitudinal mixed effects model. Biometrics 50, 1171-1177] are valid, their proof is based on the work of Self and Liang [1987. Asymptotic properties of maximum likelihood estimators and likelihood tests under nonstandard conditions. J. Amer. Statist. Assoc. 82, 605-610] which requires identically distributed random variables, an assumption not always valid in longitudinal data problems. We use the less restrictive results of Vu and Zhou [1997. Generalization of likelihood ratio tests under nonstandard conditions. Ann. Statist. 25, 897-916] to prove that the proposed mixture of chi-squared distributions is the actual asymptotic distribution of such likelihood ratios used as test statistics for null variance components in models with one or two random effects. We also consider a limited simulation study to evaluate the appropriateness of the asymptotic distribution of such likelihood ratios in moderately sized samples. (C) 2008 Elsevier B.V. All rights reserved.

Local Influence Analysis for Skew-Normal Linear Mixed Models

In this article, we consider local influence analysis for the skew-normal linear mixed model (SN-LMM). As the observed data log-likelihood associated with the SN-LMM is intractable, Cook`s well-known approach cannot be applied to obtain measures of local influence. Instead, we develop local influence measures following the approach of Zhu and Lee (2001). This approach is based on the use of an EM-type algorithm and is measurement invariant under reparametrizations. Four specific perturbation schemes are discussed. Results obtained for a simulated data set and a real data set are reported, illustrating the usefulness of the proposed methodology.

Optimal synthesis of Fenton reactor networks for phenol degradation

There is an increasing need to treat effluents contaminated with phenol with advanced oxidation processes (AOPs) to minimize their impact on the environment as well as on bacteriological populations of other wastewater treatment systems. One of the most promising AOPs is the Fenton process that relies on the Fenton reaction. Nevertheless, there are no systematic studies on Fenton reactor networks. The objective of this paper is to develop a strategy for the optimal synthesis of Fenton reactor networks. The strategy is based on a superstructure optimization approach that is represented as a mixed integer non-linear programming (MINLP) model. Network superstructures with multiple Fenton reactors are optimized with the objective of minimizing the sum of capital, operation and depreciation costs of the effluent treatment system. The optimal solutions obtained provide the reactor volumes and network configuration, as well as the quantities of the reactants used in the Fenton process. Examples based on a case study show that multi-reactor networks yield decrease of up to 45% in overall costs for the treatment plant. (C) 2010 The Institution of Chemical Engineers. Published by Elsevier B.V. All rights reserved.

Model for facilities or vendors location in a global scale considering several echelons in the Chain

The facilities location problem for companies with global operations is very complex and not well explored in the literature. This work proposes a MILP model that solves the problem through minimization of the total logistic cost. Main contributions of the model are the pioneer carrying cost calculation, the treatment given to the take-or-pay costs and to the international tax benefits such as drawback and added value taxes in Brazil. The model was successfully applied to a real case of a chemical industry with industrial plants and sales all over the world. The model application recommended a totally new sourcing model for the company.

Job rotation in assembly lines employing disabled workers

In this paper we consider the programming of job rotation in the assembly line worker assignment and balancing problem. The motivation for this study comes from the designing of assembly lines in sheltered work centers for the disabled, where workers have different task execution times. In this context, the well-known training aspects associated with job rotation are particularly desired. We propose a metric along with a mixed integer linear model and a heuristic decomposition method to solve this new job rotation problem. Computational results show the efficacy of the proposed heuristics. (C) 2009 Elsevier B.V. All rights reserved.

Orthogonal packing of identical rectangles within isotropic convex regions

A mixed integer continuous nonlinear model and a solution method for the problem of orthogonally packing identical rectangles within an arbitrary convex region are introduced in the present work. The convex region is assumed to be made of an isotropic material in such a way that arbitrary rotations of the items, preserving the orthogonality constraint, are allowed. The solution method is based on a combination of branch and bound and active-set strategies for bound-constrained minimization of smooth functions. Numerical results show the reliability of the presented approach. (C) 2010 Elsevier Ltd. All rights reserved.

Multicriteria analysis in decision making under information uncertainty

This paper presents results of research related to multicriteria decision making under information uncertainty. The Bell-man-Zadeh approach to decision making in a fuzzy environment is utilized for analyzing multicriteria optimization models (< X, M > models) under deterministic information. Its application conforms to the principle of guaranteed result and provides constructive lines in obtaining harmonious solutions on the basis of analyzing associated maxmin problems. This circumstance permits one to generalize the classic approach to considering the uncertainty of quantitative information (based on constructing and analyzing payoff matrices reflecting effects which can be obtained for different combinations of solution alternatives and the so-called states of nature) in monocriteria decision making to multicriteria problems. Considering that the uncertainty of information can produce considerable decision uncertainty regions, the resolving capacity of this generalization does not always permit one to obtain unique solutions. Taking this into account, a proposed general scheme of multicriteria decision making under information uncertainty also includes the construction and analysis of the so-called < X, R > models (which contain fuzzy preference relations as criteria of optimality) as a means for the subsequent contraction of the decision uncertainty regions. The paper results are of a universal character and are illustrated by a simple example. (c) 2007 Elsevier Inc. All rights reserved.

Optimized harvest scheduling in eucalyptus plantations under operational, spatial and climatic constraints

The representation of sustainability concerns in industrial forests management plans, in relation to environmental, social and economic aspects, involve a great amount of details when analyzing and understanding the interaction among these aspects to reduce possible future impacts. At the tactical and operational planning levels, methods based on generic assumptions usually provide non-realistic solutions, impairing the decision making process. This study is aimed at improving current operational harvesting planning techniques, through the development of a mixed integer goal programming model. This allows the evaluation of different scenarios, subject to environmental and supply constraints, increase of operational capacity, and the spatial consequences of dispatching harvest crews to certain distances over the evaluation period. As a result, a set of performance indicators was selected to evaluate all optimal solutions provided to different possible scenarios and combinations of these scenarios, and to compare these outcomes with the real results observed by the mill in the study case area. Results showed that it is possible to elaborate a linear programming model that adequately represents harvesting limitations, production aspects and environmental and supply constraints. The comparison involving the evaluated scenarios and the real observed results showed the advantage of using more holistic approaches and that it is possible to improve the quality of the planning recommendations using linear programming techniques.

A Lagrangian relaxation approach to a coupled lot-sizing and cutting stock problem

Industrial production processes involving both lot-sizing and cutting stock problems are common in many industrial settings. However, they are usually treated in a separate way, which could lead to costly production plans. In this paper, a coupled mathematical model is formulated and a heuristic method based on Lagrangian relaxation is proposed. Computational results prove its effectiveness. (C) 2009 Elsevier B.V. All rights reserved.

Sampling, WLS, and Mixed Models

Mixed models may be defined with or without reference to sampling, and can be used to predict realized random effects, as when estimating the latent values of study subjects measured with response error. When the model is specified without reference to sampling, a simple mixed model includes two random variables, one stemming from an exchangeable distribution of latent values of study subjects and the other, from the study subjects` response error distributions. Positive probabilities are assigned to both potentially realizable responses and artificial responses that are not potentially realizable, resulting in artificial latent values. In contrast, finite population mixed models represent the two-stage process of sampling subjects and measuring their responses, where positive probabilities are only assigned to potentially realizable responses. A comparison of the estimators over the same potentially realizable responses indicates that the optimal linear mixed model estimator (the usual best linear unbiased predictor, BLUP) is often (but not always) more accurate than the comparable finite population mixed model estimator (the FPMM BLUP). We examine a simple example and provide the basis for a broader discussion of the role of conditioning, sampling, and model assumptions in developing inference.

Changes in soil organic carbon storage under different agricultural management systems in the Southwest Amazon Region of Brazil

The Cerrado and Amazon regions of Brazil are probably the largest agricultural frontier of the world, and Could be a sink or source for C depending on the net effect of land use change and subsequent management on soil organic C pools. We evaluated the effects of agricultural management systems on soil organic C (SOC) stocks in the Brazilian states of Rondonia and Mato Grosso, and derived regional specific factors for soil C stock change associated with different management systems. We used 50 observations (data points) in this study, including 42 dealing with annual cropping practices and 8 dealing with perennial cropping, and analyzed the data in linear mixed-effect models. No tillage (NT) systems in Cerrado areas increased SOC Storage by 1.08 +/- 0.06 relative to SOC stocks under native conditions, while SOC storage increased by a modest factor of 1.01 +/- 0.17 in Cerradao and Amazon Forest conditions. Full tillage (FT) had negative effect on SOC storage relative to NT, decreasing SOC stocks by a factor of 0.94 +/- 0.04. but did not significantly reduce SOC stocks relative to native levels when adopted in the Cerrado region. Perennial cropping had a minimal impact on SOC stocks, estimated at a factor Value of 0.98 +/- 0.14, suggesting these systems maintain about 98% of the SOC stock found under native vegetation. The results Suggest that NT adoption may be increasing SOC with land use change from native vegetation to cropland management in the Cerrado region of Brazil. (C) 2009 Elsevier B.V. All rights reserved.

The association between air pollution and blood pressure in traffic controllers in Santo Andre, Sao Paulo, Brazil

Background: Urban air pollutants are associated with cardiovascular events. Traffic controllers are at high risk for pollution exposure during outdoor work shifts. Objective: The purpose of this study was to evaluate the relationship between air pollution and systemic blood pressure in traffic controllers during their work shifts. Methods: This cross-sectional study enrolled 19 male traffic controllers from Santo Andre city (Sao Paulo, Brazil) who were 30-60 years old and exposed to ambient air during outdoor work shifts. Systolic and diastolic blood pressure readings were measured every 15 min by an Ambulatory Arterial Blood Pressure Monitoring device. Hourly measurements (lags of 0-5 h) and the moving averages (2-5 h) of particulate matter (PM(10)), ozone (O(3)) ambient concentrations and the acquired daily minimum temperature and humidity means from the Sao Paulo State Environmental Agency were correlated with both systolic and diastolic blood pressures. Statistical methods included descriptive analysis and linear mixed effect models adjusted for temperature, humidity, work periods and time of day. Results: Interquartile increases of PM(10) (33 mu g/m(3)) and O(3) (49 mu g/m(3)) levels were associated with increases in all arterial pressure parameters, ranging from 1.06 to 2.53 mmHg. PM(10) concentration was associated with early effects (lag 0), mainly on systolic blood pressure. However, O(3) was weakly associated most consistently with diastolic blood pressure and with late cumulative effects. Conclusions: Santo Andre traffic controllers presented higher blood pressure readings while working their outdoor shifts during periods of exposure to ambient pollutant fluctuations. However, PM(10) and O(3) induced cardiovascular effects demonstrated different time courses and end-point behaviors and probably acted through different mechanisms. (C) 2011 Elsevier Inc. All rights reserved.

ELECTROMYOGRAPHIC AMPLITUDE RATIO OF SERRATUS ANTERIOR AND UPPER TRAPEZIUS MUSCLES DURING MODIFIED PUSH-UPS AND BENCH PRESS EXERCISES

Imbalance and weakness of the serratus anterior and upper trapezius force couple have been described in patients with shoulder dysfunction. There is interest in identifying exercises that selectively activate these muscles and including it in rehabilitation protocols. This study aims to verify the UT/SA electromyographic (EMG) amplitude ratio, performed in different upper limb exercises and on two bases of support. Twelve healthy men were tested (average age = 22.8 +/- 3.1 years), and surface EMG was recorded from the upper trapezius and serratus anterior using single differential surface electrodes. Volunteers performed isometric contractions over a stable base of support and on a Swiss ball during the wall push-up (WP), bench press (BP), and push-up (PU) exercises. All SEMG data are reported as a percentage of root mean square or integral of linear envelope from the maximal value obtained in one of three maximal voluntary contractions for each muscle studied. A linear mixed-effect model was performed to compare UT/SA ratio values. The WP, BP, and PU exercises showed UT/SA ratio mean +/- SD values of 0.69 +/- 0.72, 0.14 +/- 0.12, and 0.39 +/- 0.37 for stable surfaces, respectively, whereas for unstable surfaces, the values were 0.73 +/- 0.67, 0.43 +/- 0.39, and 0.32 +/- 0.30. The results demonstrate that UT/SA ratio was influenced by the exercises and by the upper limb base of support. The practical application is to show that BP on a stable surface is the exercise preferred over WP and PU on either surfaces for serratus anterior muscle training in patients with imbalance between the UT/SA force couple or serratus anterior weakness.