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.

Improved testing inference in mixed linear models

Mixed linear models are commonly used in repeated measures studies. They account for the dependence amongst observations obtained from the same experimental unit. Often, the number of observations is small, and it is thus important to use inference strategies that incorporate small sample corrections. In this paper, we develop modified versions of the likelihood ratio test for fixed effects inference in mixed linear models. In particular, we derive a Bartlett correction to such a test, and also to a test obtained from a modified profile likelihood function. Our results generalize those in [Zucker, D.M., Lieberman, O., Manor, O., 2000. Improved small sample inference in the mixed linear model: Bartlett correction and adjusted likelihood. Journal of the Royal Statistical Society B, 62,827-838] by allowing the parameter of interest to be vector-valued. Additionally, our Bartlett corrections allow for random effects nonlinear covariance matrix structure. We report simulation results which show that the proposed tests display superior finite sample behavior relative to the standard likelihood ratio test. An application is also presented and discussed. (C) 2008 Elsevier B.V. All rights reserved.

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.

Efetividade de programa de suplementação alimentar no ganho ponderal de crianças

OBJETIVO: Avaliar a efetividade de programa governamental de suplementação alimentar no ganho ponderal de crianças. MÉTODOS: Estudo de coorte com dados secundários de 25.433 crianças de baixa renda com idade entre seis e 24 meses que ingressaram em programa de distribuição de leite fortificado Projeto Vivaleite, realizado no Estado de São Paulo de 2003 a 2008. O ganho ponderal foi medido por meio dos valores de escores z de peso para idade, calculados pelo padrão da Organização Mundial da Saúde (2007), obtidos, na rotina do programa, ao ingressar e a cada quatro meses durante a permanência. As crianças foram divididas em três grupos de escore z ao entrar: sem comprometimento de peso (z > -1); risco de baixo peso (-2 < z < -1) e baixo peso (z < -2). Utilizou-se regressão linear multinível (modelo misto), permitindo a comparação, em cada idade, das médias ajustadas do escore z dos ingressantes e participantes há pelo menos quatro meses, ajustadas para correlação entre medidas repetidas. RESULTADOS: Verificou-se efeito positivo do programa no ganho de peso das crianças, variando em função do estado nutricional ao ingressar; para as que entraram sem comprometimento de peso, o ganho médio ajustado foi 0,183 escore z;entre as que entraram com risco de baixo peso, foi 0,566; e entre as ingressantes com baixo peso, foi 1,005 escore z. CONCLUSÕES: O programa é efetivo para o ganho ponderal de crianças menores de dois anos, com efeito mais pronunciado entre as crianças que entram no programa em condições menos favoráveis de peso.

Matter-wave two-dimensional solitons in crossed linear and nonlinear optical lattices

The existence of multidimensional matter-wave solitons in a crossed optical lattice (OL) with a linear optical lattice (LOL) in the x direction and a nonlinear optical lattice (NOL) in the y direction, where the NOL can be generated by a periodic spatial modulation of the scattering length using an optically induced Feshbach resonance is demonstrated. In particular, we show that such crossed LOLs and NOLs allow for stabilizing two-dimensional solitons against decay or collapse for both attractive and repulsive interactions. The solutions for the soliton stability are investigated analytically, by using a multi-Gaussian variational approach, with the Vakhitov-Kolokolov necessary criterion for stability; and numerically, by using the relaxation method and direct numerical time integrations of the Gross-Pitaevskii equation. Very good agreement of the results corresponding to both treatments is observed.

A Fault-Tolerant Manipulator Robot Based on H(2), H(infinity), and Mixed H(2)/H(infinity) Markovian Controls

This paper develops a Markovian jump model to describe the fault occurrence in a manipulator robot of three joints. This model includes the changes of operation points and the probability that a fault occurs in an actuator. After a fault, the robot works as a manipulator with free joints. Based on the developed model, a comparative study among three Markovian controllers, H(2), H(infinity), and mixed H(2)/H(infinity) is presented, applied in an actual manipulator robot subject to one and two consecutive faults.

A solid-like FEM for geometrically non-linear 3D frames

This study presents a solid-like finite element formulation to solve geometric non-linear three-dimensional inhomogeneous frames. To achieve the desired representation, unconstrained vectors are used instead of the classic rigid director triad; as a consequence, the resulting formulation does not use finite rotation schemes. High order curved elements with any cross section are developed using a full three-dimensional constitutive elastic relation. Warping and variable thickness strain modes are introduced to avoid locking. The warping mode is solved numerically in FEM pre-processing computational code, which is coupled to the main program. The extra calculations are relatively small when the number of finite elements. with the same cross section, increases. The warping mode is based on a 2D free torsion (Saint-Venant) problem that considers inhomogeneous material. A scheme that automatically generates shape functions and its derivatives allow the use of any degree of approximation for the developed frame element. General examples are solved to check the objectivity, path independence, locking free behavior, generality and accuracy of the proposed formulation. (C) 2009 Elsevier B.V. All rights reserved.

Influence of support material on the immobilization of biomass for the degradation of linear alkylbenzene sulfonate in anaerobic reactors

Two horizontal-flow anaerobic immobilized biomass reactors (HAIB) were used to study the degradation of the LAS surfactant: one filled with charcoal (HAIB1) and the other with a mixed bed of expanded clay and polyurethane foam (HAIB2). The reactors were fed with synthetic substrate supplemented with 14 mg l(-1) of LAS, kept at 30 +/- 2 degrees C and operated with a hydraulic retention time (HRT) of 12 h. The surfactant was quantified by HPLC. Spatial variation analyses were done to quantify organic matter and LAS consumption along the reactor length. The presence of the surfactant in the load did not affect the removal of organic matter (COD), which was close to 90% in both reactors for an influent COD of 550 ring l(-1). The results of a mass balance indicated that 28% of all LAS added to HAIB1 was removed by degradation. HAIB2 presented 27% degradation. Molecular biology techniques revealed microorgan isms belonging the uncultured Holophaga sp., uncultured delta Proteobacterium, uncultured Verrucomicrobium sp., Bacteroides sp. and uncultured gamma Proteobacterium sp. The reactor with biomass immobilized on charcoal presented lower adsorption and a higher kinetic degradation coefficient. So, it was the most suitable support for LAS anaerobic treatment. (c) 2008 Elsevier Ltd. All rights reserved.

Verification of a mixed high-order accurate DNS code for laminar turbulent transition by the method of manufactured solutions

This paper presents results on a verification test of a Direct Numerical Simulation code of mixed high-order of accuracy using the method of manufactured solutions (MMS). This test is based on the formulation of an analytical solution for the Navier-Stokes equations modified by the addition of a source term. The present numerical code was aimed at simulating the temporal evolution of instability waves in a plane Poiseuille flow. The governing equations were solved in a vorticity-velocity formulation for a two-dimensional incompressible flow. The code employed two different numerical schemes. One used mixed high-order compact and non-compact finite-differences from fourth-order to sixth-order of accuracy. The other scheme used spectral methods instead of finite-difference methods for the streamwise direction, which was periodic. In the present test, particular attention was paid to the boundary conditions of the physical problem of interest. Indeed, the verification procedure using MMS can be more demanding than the often used comparison with Linear Stability Theory. That is particularly because in the latter test no attention is paid to the nonlinear terms. For the present verification test, it was possible to manufacture an analytical solution that reproduced some aspects of an instability wave in a nonlinear stage. Although the results of the verification by MMS for this mixed-order numerical scheme had to be interpreted with care, the test was very useful as it gave confidence that the code was free of programming errors. Copyright (C) 2009 John Wiley & Sons, Ltd.

A hybrid-mixed finite element formulation for the geometrically exact analysis of three-dimensional framed structures

This paper addresses the development of a hybrid-mixed finite element formulation for the quasi-static geometrically exact analysis of three-dimensional framed structures with linear elastic behavior. The formulation is based on a modified principle of stationary total complementary energy, involving, as independent variables, the generalized vectors of stress-resultants and displacements and, in addition, a set of Lagrange multipliers defined on the element boundaries. The finite element discretization scheme adopted within the framework of the proposed formulation leads to numerical solutions that strongly satisfy the equilibrium differential equations in the elements, as well as the equilibrium boundary conditions. This formulation consists, therefore, in a true equilibrium formulation for large displacements and rotations in space. Furthermore, this formulation is objective, as it ensures invariance of the strain measures under superposed rigid body rotations, and is not affected by the so-called shear-locking phenomenon. Also, the proposed formulation produces numerical solutions which are independent of the path of deformation. To validate and assess the accuracy of the proposed formulation, some benchmark problems are analyzed and their solutions compared with those obtained using the standard two-node displacement/ rotation-based formulation.

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.

A linear optimization approach to the combined production planning model

Two fundamental processes usually arise in the production planning of many industries. The first one consists of deciding how many final products of each type have to be produced in each period of a planning horizon, the well-known lot sizing problem. The other process consists of cutting raw materials in stock in order to produce smaller parts used in the assembly of final products, the well-studied cutting stock problem. In this paper the decision variables of these two problems are dependent of each other in order to obtain a global optimum solution. Setups that are typically present in lot sizing problems are relaxed together with integer frequencies of cutting patterns in the cutting problem. Therefore, a large scale linear optimizations problem arises, which is exactly solved by a column generated technique. It is worth noting that this new combined problem still takes the trade-off between storage costs (for final products and the parts) and trim losses (in the cutting process). We present some sets of computational tests, analyzed over three different scenarios. These results show that, by combining the problems and using an exact method, it is possible to obtain significant gains when compared to the usual industrial practice, which solve them in sequence. (C) 2010 The Franklin Institute. Published by Elsevier Ltd. 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.