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.

Internal Residual Stress Measurements in a Bioactive Glass-Ceramic Using Vickers Indentation

The residual stress distribution that arises in the glass matrix during cooling of a partially crystallized 17.2Na(2)O-32.1CaO-48.1SiO(2)-2.5P(2)O(5) (mol%) bioactive glass-ceramic was measured using the Vickers indentation method proposed by Zeng and Rowcliffe (ZR). The magnitude of the determined residual stress at the crystal/glass boundary was 1/4-1/3 of the values measured using X-ray diffraction (within the crystals) and calculated using Selsing`s model. A correction for the crack geometry factor, assuming a semi-elliptical shape, is proposed and then good agreement between experimental and theoretical values is found. Thus, if the actual crack geometry is taken into account, the indentation technique of ZR can be successfully used. In addition, a numerical model for the calculation of residual stresses that takes into account the hemispherical shape of the crystalline precipitates at a free surface was developed. The result is that near the sample surface, the radial component of the residual stress is increased by 70% in comparison with the residual stress calculated by Selsing`s model.

Micro Far-Infrared Reflectivity of CaNb(2)O(6) Single Crystal Fibers Grown by the Laser-Heated Pedestal Growth Technique

CaNb(2)O(6) single crystal fibers were grown by the laser-heated pedestal growth technique, directly from the starting reagents. Optically transparent fibers were obtained in the form of rods with elliptical cross-section, free from cracks, impurities, and secondary phases, with an average diameter of 0.4 mm and about 20 mm of length. The fibers grew within the orthorhombic Pbcn columbite structure, with the growth axis nearly parallel to the crystallographic a-direction. The parameters b and c were parallel to the shorter and larger ellipsis axes. A special setup using a microscope was developed to obtain the far-infrared reflectivity spectra of these micrometer-sized fibers, allowing the identification and assignment of 34 of the 38 polar phonons foreseen for the material. From these phonons, the intrinsic dielectric constant ( = 18.2) and quality factor ( of 185 THz) could be estimated, showing the potential of the material for applications in microwave circuitry. These results, along with previous polarized Raman data (Cryst. Growth Des. 2010, 10, 1569), allow us to present a comprehensive set of optical phonon modes and to discuss the potential use of designed CaNb(2)O(6) microcrystals in compact optical devices.

Two- and three-dimensional shape fabric analysis by the intercept method in grey levels

The count intercept is a robust method for the numerical analysis of fabrics Launeau and Robin (1996). It counts the number of intersections between a set of parallel scan lines and a mineral phase, which must be identified on a digital image. However, the method is only sensitive to boundaries and therefore supposes the user has some knowledge about their significance. The aim of this paper is to show that a proper grey level detection of boundaries along scan lines is sufficient to calculate the two-dimensional anisotropy of grain or crystal distributions without any particular image processing. Populations of grains and crystals usually display elliptical anisotropies in rocks. When confirmed by the intercept analysis, a combination of a minimum of 3 mean length intercept roses, taken on 3 more or less perpendicular sections, allows the calculation of 3-dimensional ellipsoids and the determination of their standard deviation with direction and intensity in 3 dimensions as well. The feasibility of this quick method is attested by numerous examples on theoretical objects deformed by active and passive deformation, on BSE images of synthetic magma flow, on drawing or direct analysis of thin section pictures of sandstones and on digital images of granites directly taken and measured in the field. (C) 2010 Elsevier B.V. All rights reserved.

Random number generators for the generalized Birnbaum-Saunders distribution

The generalized Birnbaum-Saunders distribution pertains to a class of lifetime models including both lighter and heavier tailed distributions. This model adapts well to lifetime data, even when outliers exist, and has other good theoretical properties and application perspectives. However, statistical inference tools may not exist in closed form for this model. Hence, simulation and numerical studies are needed, which require a random number generator. Three different ways to generate observations from this model are considered here. These generators are compared by utilizing a goodness-of-fit procedure as well as their effectiveness in predicting the true parameter values by using Monte Carlo simulations. This goodness-of-fit procedure may also be used as an estimation method. The quality of this estimation method is studied here. Finally, through a real data set, the generalized and classical Birnbaum-Saunders models are compared by using this estimation method.

A multivariate ultrastructural errors-in-variables model with equation error

This paper deals with asymptotic results on a multivariate ultrastructural errors-in-variables regression model with equation errors Sufficient conditions for attaining consistent estimators for model parameters are presented Asymptotic distributions for the line regression estimators are derived Applications to the elliptical class of distributions with two error assumptions are presented The model generalizes previous results aimed at univariate scenarios (C) 2010 Elsevier Inc All rights reserved

Semiparametric Bayesian measurement error modeling

This work presents a Bayesian semiparametric approach for dealing with regression models where the covariate is measured with error. Given that (1) the error normality assumption is very restrictive, and (2) assuming a specific elliptical distribution for errors (Student-t for example), may be somewhat presumptuous; there is need for more flexible methods, in terms of assuming only symmetry of errors (admitting unknown kurtosis). In this sense, the main advantage of this extended Bayesian approach is the possibility of considering generalizations of the elliptical family of models by using Dirichlet process priors in dependent and independent situations. Conditional posterior distributions are implemented, allowing the use of Markov Chain Monte Carlo (MCMC), to generate the posterior distributions. An interesting result shown is that the Dirichlet process prior is not updated in the case of the dependent elliptical model. Furthermore, an analysis of a real data set is reported to illustrate the usefulness of our approach, in dealing with outliers. Finally, semiparametric proposed models and parametric normal model are compared, graphically with the posterior distribution density of the coefficients. (C) 2009 Elsevier Inc. All rights reserved.

Multivariate measurement error models based on scale mixtures of the skew-normal distribution

Scale mixtures of the skew-normal (SMSN) distribution is a class of asymmetric thick-tailed distributions that includes the skew-normal (SN) distribution as a special case. The main advantage of these classes of distributions is that they are easy to simulate and have a nice hierarchical representation facilitating easy implementation of the expectation-maximization algorithm for the maximum-likelihood estimation. In this paper, we assume an SMSN distribution for the unobserved value of the covariates and a symmetric scale mixtures of the normal distribution for the error term of the model. This provides a robust alternative to parameter estimation in multivariate measurement error models. Specific distributions examined include univariate and multivariate versions of the SN, skew-t, skew-slash and skew-contaminated normal distributions. The results and methods are applied to a real data set.

An extension of the generalized Birnbaum-Saunders distribution

In this paper we present an extension of the generalized Birnbaum-Saunders distribution family introduced in [Diaz-Garcia, J.A., Leiva-Sanchez, V., 2005. A new family of life distributions based on the contoured elliptically distributions. Journal of Statistical Planning and Inference 128 (2), 445-457] with a view to make it even more flexible in terms of its kurtosis coefficient. Properties involving moments and asymmetry and kurtosis indexes are studied for some special members of this family such as the slash Birnbaum-Saunders and slash-t Birnbaum-Saunders. Simulation studies for some particular cases and a real data analysis are also reported, illustrating the usefulness of the extension considered. (C) 2008 Elsevier B.V. All rights reserved.

Efficient solutions to robust, semi-implicit discretizations of the immersed boundary method

The immersed boundary method is a versatile tool for the investigation of flow-structure interaction. In a large number of applications, the immersed boundaries or structures are very stiff and strong tangential forces on these interfaces induce a well-known, severe time-step restriction for explicit discretizations. This excessive stability constraint can be removed with fully implicit or suitable semi-implicit schemes but at a seemingly prohibitive computational cost. While economical alternatives have been proposed recently for some special cases, there is a practical need for a computationally efficient approach that can be applied more broadly. In this context, we revisit a robust semi-implicit discretization introduced by Peskin in the late 1970s which has received renewed attention recently. This discretization, in which the spreading and interpolation operators are lagged. leads to a linear system of equations for the inter-face configuration at the future time, when the interfacial force is linear. However, this linear system is large and dense and thus it is challenging to streamline its solution. Moreover, while the same linear system or one of similar structure could potentially be used in Newton-type iterations, nonlinear and highly stiff immersed structures pose additional challenges to iterative methods. In this work, we address these problems and propose cost-effective computational strategies for solving Peskin`s lagged-operators type of discretization. We do this by first constructing a sufficiently accurate approximation to the system`s matrix and we obtain a rigorous estimate for this approximation. This matrix is expeditiously computed by using a combination of pre-calculated values and interpolation. The availability of a matrix allows for more efficient matrix-vector products and facilitates the design of effective iterative schemes. We propose efficient iterative approaches to deal with both linear and nonlinear interfacial forces and simple or complex immersed structures with tethered or untethered points. One of these iterative approaches employs a splitting in which we first solve a linear problem for the interfacial force and then we use a nonlinear iteration to find the interface configuration corresponding to this force. We demonstrate that the proposed approach is several orders of magnitude more efficient than the standard explicit method. In addition to considering the standard elliptical drop test case, we show both the robustness and efficacy of the proposed methodology with a 2D model of a heart valve. (C) 2009 Elsevier Inc. All rights reserved.

A robust Bayesian approach to null intercept measurement error model with application to dental data

Measurement error models often arise in epidemiological and clinical research. Usually, in this set up it is assumed that the latent variable has a normal distribution. However, the normality assumption may not be always correct. Skew-normal/independent distribution is a class of asymmetric thick-tailed distributions which includes the Skew-normal distribution as a special case. In this paper, we explore the use of skew-normal/independent distribution as a robust alternative to null intercept measurement error model under a Bayesian paradigm. We assume that the random errors and the unobserved value of the covariate (latent variable) follows jointly a skew-normal/independent distribution, providing an appealing robust alternative to the routine use of symmetric normal distribution in this type of model. Specific distributions examined include univariate and multivariate versions of the skew-normal distribution, the skew-t distributions, the skew-slash distributions and the skew contaminated normal distributions. The methods developed is illustrated using a real data set from a dental clinical trial. (C) 2008 Elsevier B.V. All rights reserved.

Könsroller i Film : En kvalitativ analys av hur unga tjejer skildras i tre moderna svenska ungdomsfilmer

Min utgångspunkt i denna uppsats var könsroller i film. Utifrån detta skapade jag mitt syfte som var att undersöka hur valda kvinnliga karaktärer i huvudrollen i tre svenska ungdomsfilmer skildrades, samt analysera huruvida skildringarna av karaktärerna bevarar eller bryter mot etablerade könsroller utifrån aktuell forskning. ¨Mina frågeställningar var:•Vilka utmärkande könsrollsmönster går att utröna i de valda karaktärernas agerande/handlande gentemot andra?•Vilken reaktion bemöts de valda karaktärerna med, utifrån detta agerade/handlande?•Vilka budskap landar filmen i, när det gäller huvudrollsinnehavarnas inordning eller revoltering utifrån de etablerade könsroller?De tre karaktärerna som jag valde att analysera var Sofie i Hipp Hipp Hora, Emma i Fjorton Suger och Ida i Sandor /slash/ Ida. Jag skapade ett analysverktyg bestående utav några nyckelkategorier kring könsroller. Utifrån dessa kategorier kunde jag analysera hur karaktärernas agerande gentemot andra bemöttes, och om karaktärerna bröt eller bevarade etablerade könsrollerna utifrån aktuell forskning. Mitt resultat visar karaktärerna vågar bryta emot de etablerade könsrollerna. Filmerna förmedlar även en bild utav unga smala vackra blonda tjejer som är tuffa och vågar gå sina egna vägar trots att det inte alltid är den lätta eller den rätta vägen.

Fanfiction : la anonimicidad, el texto y el armario

This essay looks at slash, a genre within fanfiction, from the perspective of Sedgwick’s theory of the closet, which reflects on the concealing mechanisms associated with homosexuality. While the real author stays in the closet, disguised behind a pseudonym, slash texts present homosexual themes in a very explicit way, often relying on humor or subversive elements. Between these two spheres, the real author and the text, we can find what we call the author’s voice, conscious about the existing homophobic structures, a voice that uses different strategies to shield itself against them. Internet, with the possibility to stay anonymous, serves as a social closet where the masked authors create texts that subvert heteronormativity.

Monitoring Surface Displacement of the Colby Green Retaining Pond Dams

The Colby Green is a campus expansion project which began in October of 2003. The construction would result in three new buildings, additional parking, and an elliptical 75,000-squarefoot green southeast of Mayflower Hill Drive. There were also plans for the construction of three run-off management and sediment ponds below the green, to manage flooding of the green. Three drains in the green transport water to the three retaining ponds which slowly disperse water into the surrounding environment. The ponds were created by constructing earthen dams around the drain outlets. The dams are composed of soil, cobbles, and boulders procured from the surrounding excavation site. Unfortunately, earthen dams are susceptible to many types of erosion which result in their failure. In this case the potential for clay and silt from the underlying Presumpscot Formation to mix with the soil in the earthen dams raised concerns with regards to frost action. In order to monitor the surface displacement of the dams I drove 92 poles into the ground in 8 straight lines across the faces of the dams in the fall of 2005. I returned to the sites during and after the spring thaw of 2006, to check for any signs of movement resulting from frost-heave, surface creep, or any other form of mass wasting. Fortunately, there was no recordable sign of movement in the stakes across any of the retaining ponds. The dams appear to be functioning as designed.

Influência do resíduo de colheita de cana-de-açúcar sem queima sobre a eficiência dos herbicidas imazapic e imazapic + pendimethalin

Os resíduos deixados sobre o solo por ocasião da colheita mecanizada da cana-de-açúcar podem constituir-se em uma barreira física para a ação dos herbicidas no controle de plantas daninhas, quando aplicados em pré-emergência destas plantas sobre a palha da cana. em virtude disso, o presente trabalho teve por objetivos analisar e quantificar a interferência dessa camada de palha sobre o solo na ação dos herbicidas imazapic e imazapic + pendimethalin no controle de plantas daninhas em áreas onde a cana-de-açúcar foi colhida mecanicamente sem a queima da palhada previamente à colheita. Foram realizados dois ensaios simultâneos: um com a retirada da palha dois dias após a aplicação dos herbicidas e o outro com a manutenção desta, ambos conduzidos em casa de vegetação. O imazapic isolado foi aplicado nas dosagens de 0, 122,5 e 147 g i.a.ha-1 e em mistura com pendimethalin na dosagem de 75 + 1500 g i.a.ha-1, com simulação de chuvas nas intensidades de 30, 60 e 90 mm. Após análise dos resultados de biomassa seca, altura e número de folhas das plantas de Sorghum bicolor e Cyperus rotundus, além de nota visual e biomassa seca de Panicum maximum, Brachiaria plantaginea, Digitaria horizontalis, Amaranthus viridis, Ipomoea grandifolia e Brachiaria decumbens, constatou-se eficiência proporcional dos herbicidas à dosagem utilizada, independentemente da presença da palha, à exceção de Ipomoea grandifolia e Brachiaria decumbens, além de haver menor controle nos tratamentos submetidos à chuva de 90 mm. Esses resultados indicam boas perspectivas quanto à aplicação destes herbicidas em áreas de colheita mecanizada de cana-de-açúcar sem queima, para controle de plantas daninhas em condições de pré-emergência.