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.

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.

Clustering and diffusion in a symplectic map lattice with non-local coupling

We study a symplectic chain with a non-local form of coupling by means of a standard map lattice where the interaction strength decreases with the lattice distance as a power-law, in Such a way that one can pass continuously from a local (nearest-neighbor) to a global (mean-field) type of coupling. We investigate the formation of map clusters, or spatially coherent structures generated by the system dynamics. Such clusters are found to be related to stickiness of chaotic phase-space trajectories near periodic island remnants, and also to the behavior of the diffusion coefficient. An approximate two-dimensional map is derived to explain some of the features of this connection. (C) 2008 Elsevier Ltd. All rights reserved.

ASSOCIATED A-TYPE SUBALKALINE AND HIGH-K CALC-ALKALINE GRANITES IN THE ITU GRANITE PROVINCE, SOUTHEASTERN BRAZIL: PETROLOGICAL AND TECTONIC SIGNIFICANCE

The 590-580 Ma Itu Granite Province (IGP) is a roughly linear belt of post-orogenic granite plutons similar to 60 km wide extending for some 350 km along the southern edge of the Apia-Guaxupe Terrane in southeastern Brazil. Typical components are subalkaline A-type granites (some with rapakivi texture) that crystallized at varied, but mostly strongly oxidizing conditions, and contrast with a coeval association of also oxidized high-K calc-alkaline granites in terms of major (e. g., lower Ca/Fe) and trace elements (higher Nb, Y, Zr). Mantle-derived magmas (such as those forming the LILE-rich Piracaia Monzodiorite, with epsilon(Nd(t)) = -7 to -10, (87)Sr/(86)Sr((t)) = 0.7045-0.7055) are inferred to derive from enriched subcontinental lithosphere modified during previous subduction, and may have played a role in the generation of the A-type granites, adding melts or fluids or both to the lower crust from which the latter were generated. The IGP is interpreted as a reflection of crust uplift and increased heat flux during ascent of hot, less dense asthenosphere after continental collision, probably reflecting breakoff of an oceanic slab coeval to the right-lateral accretion of a terrane related to the Mantiqueira Orogenic System.

Groundwater favorability map using GIS multicriteria data analysis on crystalline terrain, Sao Paulo State, Brazil

This paper presents the groundwater favorability mapping on a fractured terrain in the eastern portion of Sao Paulo State, Brazil. Remote sensing, airborne geophysical data, photogeologic interpretation, geologic and geomorphologic maps and geographic information system (GIS) techniques have been used. The results of cross-tabulation between these maps and well yield data allowed groundwater prospective parameters in a fractured-bedrock aquifer. These prospective parameters are the base for the favorability analysis whose principle is based on the knowledge-driven method. The mutticriteria analysis (weighted linear combination) was carried out to give a groundwater favorabitity map, because the prospective parameters have different weights of importance and different classes of each parameter. The groundwater favorability map was tested by cross-tabulation with new well yield data and spring occurrence. The wells with the highest values of productivity, as well as all the springs occurrence are situated in the excellent and good favorabitity mapped areas. It shows good coherence between the prospective parameters and the well yield and the importance of GIS techniques for definition of target areas for detail study and wells location. (c) 2008 Elsevier B.V. All rights reserved.

Improved score tests in symmetric linear regression models

The class of symmetric linear regression models has the normal linear regression model as a special case and includes several models that assume that the errors follow a symmetric distribution with longer-than-normal tails. An important member of this class is the t linear regression model, which is commonly used as an alternative to the usual normal regression model when the data contain extreme or outlying observations. In this article, we develop second-order asymptotic theory for score tests in this class of models. We obtain Bartlett-corrected score statistics for testing hypotheses on the regression and the dispersion parameters. The corrected statistics have chi-squared distributions with errors of order O(n(-3/2)), n being the sample size. The corrections represent an improvement over the corresponding original Rao`s score statistics, which are chi-squared distributed up to errors of order O(n(-1)). Simulation results show that the corrected score tests perform much better than their uncorrected counterparts in samples of small or moderate size.

Three Bartlett-type corrections for score statistics in symmetric nonlinear regression models

We present simple matrix formulae for corrected score statistics in symmetric nonlinear regression models. The corrected score statistics follow more closely a chi (2) distribution than the classical score statistic. Our simulation results indicate that the corrected score tests display smaller size distortions than the original score test. We also compare the sizes and the powers of the corrected score tests with bootstrap-based score tests.

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.

Testing linear cointegration against smooth-transition cointegration

This paper studies a smooth-transition (ST) type cointegration. The proposed ST cointegration allows for regime switching structure in a cointegrated system. It nests the linear cointegration developed by Engle and Granger (1987) and the threshold cointegration studied by Balke and Fomby (1997). We develop F-type tests to examine linear cointegration against ST cointegration in ST-type cointegrating regression models with or without time trends. The null asymptotic distributions of the tests are derived with stationary transition variables in ST cointegrating regression models. And it is shown that our tests have nonstandard limiting distributions expressed in terms of standard Brownian motion when regressors are pure random walks, while have standard asymptotic distributions when regressors contain random walks with nonzero drift. Finite-sample distributions of those tests are studied by Monto Carlo simulations. The small-sample performance of the tests states that our F-type tests have a better power when the system contains ST cointegration than when the system is linearly cointegrated. An empirical example for the purchasing power parity (PPP) data (monthly US dollar, Italy lira and dollar-lira exchange rate from 1973:01 to 1989:10) is illustrated by applying the testing procedures in this paper. It is found that there is no linear cointegration in the system, but there exits the ST-type cointegration in the PPP data.

A Map of Impoundments and Their Primary Purpose on Maine Rivers

http://digitalcommons.colby.edu/atlasofmaine2005/1007/thumbnail.jpg

DETAILED MAP AND RECREATIONAL SUITABILITY ASSESSMENT OF THE COLBY TRAIL SYSTEM

http://digitalcommons.colby.edu/atlasofmaine2005/1026/thumbnail.jpg

Joint dynamic probabilistic constraints with projected linear decision rules

We consider multistage stochastic linear optimization problems combining joint dynamic probabilistic constraints with hard constraints. We develop a method for projecting decision rules onto hard constraints of wait-and-see type. We establish the relation between the original (in nite dimensional) problem and approximating problems working with projections from di erent subclasses of decision policies. Considering the subclass of linear decision rules and a generalized linear model for the underlying stochastic process with noises that are Gaussian or truncated Gaussian, we show that the value and gradient of the objective and constraint functions of the approximating problems can be computed analytically.

Voltammetric determination of L-dopa using an electrode modified with trinuclear ruthenium ammine complex (Ru-red) supported on Y-type zeolite

The L-dopa is the immediate precursor of the neurotransmitter dopamine. Unlike dopamine, L-dopa easily enters the central nervous system and is used in the treatment of Parkinson's disease. A sensitive and selective method is presented for the voltammetric determination of L-dopa in pharmaceutical formulations using a carbon paste electrode modified with trinuclear ruthenium ammine complex [(NH3)(5)Ru-III-O-Ru-IV(NH3)(4)-O-Ru-III(NH3)(5)](6+) (Ru-red) incorporated in NaY zeolite. The parameters which influence on the electrode response (paste composition, potential scan rate, pH and interference) were also investigated. The optimum conditions were found to an electrode composition (m/m) of 25% zeolite containing 6.7% Ru, 50% graphite and 25% mineral oil in acetate buffer at pH 4.8. Voltammetric peak currents showed a linear response for L-dopa concentration in the range between 1.2 x 10(-4) and 1.0 x 10(-2) Mol l(-1) (r = 0.9988) with a detection limit of 8.5 x 10(-5) mol l(-1). The variation coefficient for a 1.0 x 10(-3) mol l(-1) L-dopa (n = 10) was 5.5%. The results obtained for L-dopa in pharmaceutical formulations (tablet) was in agreement with compared official method. In conclusion, this study has illustrated that the proposed electrode modified with Ru-red incorporated zeolite is suitable valuable for selective measurements of L-dopa. (C) 2004 Elsevier B.V. All rights reserved.

Comparison of non-linear growth models to describe the growth curve in West African Dwarf sheep

The objectives of this study were to compare the goodness of fit of four non-linear growth models, i.e. Brody, Gompertz, Logistic and Von Bertalanffy, in West African Dwarf (WAD) sheep. A total of 5274 monthly weight records from birth up to 180 days of age from 889 lambs, collected during 2001 to 2004 in Betecoucou breeding farm in Benin were used. In the preliminary analysis, the General Linear Model Procedure of the Statistical Analysis Systems Institute was applied to the dataset to identify the significant effects of the sex of lamb (male and female), type of birth (single and twin), season of birth (rainy season and dry season), parity of dam (1, 2 and 3) and year of birth (2001, 2002, 2003 and 2004) on the observed birth weight and monthly weight up to 6 months of age. The models parameters (A, B and k), coefficient of determination (112), mean square error (MSE) were calculated using language of technical computing package Matlab(R), 2006. The mean values of A, B and k were substituted into each model to calculate the corresponding Akaike's Information Criterion (AIC). Among the four growth functions, the Brody model has been selected for its accuracy of fit according to the higher R(2), lower MSE and A/C Finally, the parameters A, B and k were adjusted in Matlab(R) 2006 for the sex of lamb, year of birth, season of birth, birth type and the parity of ewe, providing a specific slope of the Brody growth curve. The results of this study suggest that Brody model can be useful for WAD sheep breeding in Betecoucou farm conditions through growth monitoring.

Voltammetric performance and application of a sensor for sodium ions constructed with layered birnessite-type manganese oxide

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)