Inference and local influence assessment in skew-normal null intercept measurement error model

In this article, we discuss inferential aspects of the measurement error regression models with null intercepts when the unknown quantity x (latent variable) follows a skew normal distribution. We examine first the maximum-likelihood approach to estimation via the EM algorithm by exploring statistical properties of the model considered. Then, the marginal likelihood, the score function and the observed information matrix of the observed quantities are presented allowing direct inference implementation. In order to discuss some diagnostics techniques in this type of models, we derive the appropriate matrices to assessing the local influence on the parameter estimates under different perturbation schemes. The results and methods developed in this paper are illustrated considering part of a real data set used by Hadgu and Koch [1999, Application of generalized estimating equations to a dental randomized clinical trial. Journal of Biopharmaceutical Statistics, 9, 161-178].

Epsilon Birnbaum-Saunders distribution family: properties and inference

In this paper we introduce a new extension for the Birnbaum-Saunder distribution based on the family of the epsilon-skew-symmetric distributions studied in Arellano-Valle et al. (J Stat Plan Inference 128(2):427-443, 2005). The extension allows generating Birnbaun-Saunders type distributions able to deal with extreme or outlying observations (Dupuis and Mills, IEEE Trans Reliab 47:88-95, 1998). Basic properties such as moments and Fisher information matrix are also studied. Results of a real data application are reported illustrating good fitting properties of the proposed model.

Inflated beta distributions

This paper considers the issue of modeling fractional data observed on [0,1), (0,1] or [0,1]. Mixed continuous-discrete distributions are proposed. The beta distribution is used to describe the continuous component of the model since its density can have quite different shapes depending on the values of the two parameters that index the distribution. Properties of the proposed distributions are examined. Also, estimation based on maximum likelihood and conditional moments is discussed. Finally, practical applications that employ real data are presented.

Properties and Inference on the Skew-Curved-Symmetric Family of Distributions

In this article, we study some results related to a specific class of distributions, called skew-curved-symmetric family of distributions that depends on a parameter controlling the skewness and kurtosis at the same time. Special elements of this family which are studied include symmetric and well-known asymmetric distributions. General results are given for the score function and the observed information matrix. It is shown that the observed information matrix is always singular for some special cases. We illustrate the flexibility of this class of distributions with an application to a real dataset on characteristics of Australian athletes.

Matrix Representation of the Stationary Measure for the Multispecies TASEP

In this work we construct the stationary measure of the N species totally asymmetric simple exclusion process in a matrix product formulation. We make the connection between the matrix product formulation and the queueing theory picture of Ferrari and Martin. In particular, in the standard representation, the matrices act on the space of queue lengths. For N > 2 the matrices in fact become tensor products of elements of quadratic algebras. This enables us to give a purely algebraic proof of the stationary measure which we present for N=3.

An R implementation for generalized Birnbaum-Saunders distributions

The Birnbaum-Saunders (BS) model is a positively skewed statistical distribution that has received great attention in recent decades. A generalized version of this model was derived based on symmetrical distributions in the real line named the generalized BS (GBS) distribution. The R package named gbs was developed to analyze data from GBS models. This package contains probabilistic and reliability indicators and random number generators from GBS distributions. Parameter estimates for censored and uncensored data can also be obtained by means of likelihood methods from the gbs package. Goodness-of-fit and diagnostic methods were also implemented in this package in order to check the suitability of the GBS models. in this article, the capabilities and features of the gbs package are illustrated by using simulated and real data sets. Shape and reliability analyses for GBS models are presented. A simulation study for evaluating the quality and sensitivity of the estimation method developed in the package is provided and discussed. (C) 2008 Elsevier B.V. All rights reserved.

Testing hypotheses in the Birnbaum-Saunders distribution under type-II censored samples

The two-parameter Birnbaum-Saunders distribution has been used successfully to model fatigue failure times. Although censoring is typical in reliability and survival studies, little work has been published on the analysis of censored data for this distribution. In this paper, we address the issue of performing testing inference on the two parameters of the Birnbaum-Saunders distribution under type-II right censored samples. The likelihood ratio statistic and a recently proposed statistic, the gradient statistic, provide a convenient framework for statistical inference in such a case, since they do not require to obtain, estimate or invert an information matrix, which is an advantage in problems involving censored data. An extensive Monte Carlo simulation study is carried out in order to investigate and compare the finite sample performance of the likelihood ratio and the gradient tests. Our numerical results show evidence that the gradient test should be preferred. Further, we also consider the generalized Birnbaum-Saunders distribution under type-II right censored samples and present some Monte Carlo simulations for testing the parameters in this class of models using the likelihood ratio and gradient tests. Three empirical applications are presented. (C) 2011 Elsevier B.V. All rights reserved.

QUASISTATIONARY DISTRIBUTIONS AND FLEMING-VIOT PROCESSES IN FINITE SPACES

Consider a continuous-time Markov process with transition rates matrix Q in the state space Lambda boolean OR {0}. In In the associated Fleming-Viot process N particles evolve independently in A with transition rates matrix Q until one of them attempts to jump to state 0. At this moment the particle jumps to one of the positions of the other particles, chosen uniformly at random. When Lambda is finite, we show that the empirical distribution of the particles at a fixed time converges as N -> infinity to the distribution of a single particle at the same time conditioned on not touching {0}. Furthermore, the empirical profile of the unique invariant measure for the Fleming-Viot process with N particles converges as N -> infinity to the unique quasistationary distribution of the one-particle motion. A key element of the approach is to show that the two-particle correlations are of order 1/N.

A fast and robust statistical test based on likelihood ratio with Bartlett correction to identify Granger causality between gene sets

We propose a likelihood ratio test ( LRT) with Bartlett correction in order to identify Granger causality between sets of time series gene expression data. The performance of the proposed test is compared to a previously published bootstrapbased approach. LRT is shown to be significantly faster and statistically powerful even within non- Normal distributions. An R package named gGranger containing an implementation for both Granger causality identification tests is also provided.

Evaluation of in vitro human gingival fibroblast seeding on acellular dermal matrix

The acellular dermal matrix (ADM) was introduced in periodontology as a substitute for the autogenous grafts, which became restricted because of the limited source of donor's tissue. The aim of this study was to investigate, in vitro, the distribution, proliferation and viability of human gingival fibroblasts seeded onto ADM. ADM was seeded with human gingival fibroblasts for up to 21 days. The following parameters were evaluated: cell distribution, proliferation and viability. Results revealed that, at day 7, fibroblasts were adherent and spread on ADM surface, and were unevenly distributed, forming a discontinuous single cell layer; at day 14, a confluent fibroblastic monolayer lining ADM surface was noticed. At day 21, the cell monolayer exhibited a reduction in cell density. At 7 days, about to 90% of adherent cells on ADM surface were cycling while at 14 and 21 days this proportion was significantly reduced. A high proportion of viable cell was detected on AMD surface both on 14 and 21 days. The results suggest that fibroblast seeding onto ADM for 14 days can allow good conditions for cell adhesion and spreading on the matrix; however, migration inside the matrix was limited.

Production and body composition traits of broilers in relation to breast weight evaluated by path analysis

Breast weight has great economic importance in poultry industry, and may be associated with other variables. This work aimed to estimate phenotypic correlations between performance (live body weight at 7 and 28 days, and at slaughter, and depth of the breast muscle measured by ultrasonography), carcass (eviscerated body weight and leg weight) and body composition (heart, liver and abdominal fat weight) traits in a broiler line, and quantify the direct and indirect influence of these traits on breast weight. Path analysis was used by expanding the matrix of partial correlation in coefficients which give the direct influence of one trait on another, regardless the effect of the other traits. The simultaneous maintenance of live body weight at slaughter and eviscerated body weight in the matrix of correlations might be harmful for statistical analysis involving systems of normal equations, like path analysis, due to the observed multicollinearity. The live body weight at slaughter and the depth of the breast muscle as measured by ultrasonography directly affected breast weight and were identified as the most responsible factors for the magnitude of the correlation coefficients obtained between the studied traits and breast weight. Individual pre-selection for these traits could favor an increased breast weight in the future reproducer candidates of this line if the broilers' environmental conditions and housing are maintained, since the live body weight at slaughter and the depth of breast muscle measured by ultrasonography were directly related to breast weight.

Metronidazole release using natural rubber latex as matrix

Natural Rubber Latex (NRL) can be used successfully in controlled release drug delivery due to their excellent matrix forming properties. Recently, NRL has shown to stimulate angiogenesis, cellular adhesion and the formation of extracellular matrix, promoting the replacement and regeneration of tissue. A dermatological delivery system comprising a topically acceptable, inert support impregnated with a metronidazole (MET) solution was developed. MET 2-(2- methyl- 5-nitro- 1H- imidazol- 1-yl) ethanol, has been widely used for the treatment of protozoa and anaerobic bacterial infections. MET is a nitroimidazole anti-infective medication used mainly in the treatment of infections caused by susceptible organisms, particularly anaerobic bacteria and protozoa. In a previous study, we have tested NRL as an occlusive membrane for GBR with promising results. One possible way to decrease the inflammatory process, it was incorporated the MET in NRL. MET was incorporated into the NRL, by mixing it in solution for in vitro protein delivery experiments. The solutions of latex and MET were polymerized at different temperatures, from -100 to 40 °C, in order to control the membrane morphology. SEM microscopy analysis showed that the number, size and distribution of pores in NRL membranes varied depending on polymerization temperature, as well as its overall morphology. Results demonstrated that the best drug-delivery system was the membrane polymerized at -100 °C, which does release 77,1% of its MET content for up 310 hours.

Neural networks and statistical analysis for classification of corneal videokeratography maps based on Zernike coefficients: a quantitative comparison

PURPOSE: The main goal of this study was to develop and compare two different techniques for classification of specific types of corneal shapes when Zernike coefficients are used as inputs. A feed-forward artificial Neural Network (NN) and discriminant analysis (DA) techniques were used. METHODS: The inputs both for the NN and DA were the first 15 standard Zernike coefficients for 80 previously classified corneal elevation data files from an Eyesys System 2000 Videokeratograph (VK), installed at the Departamento de Oftalmologia of the Escola Paulista de Medicina, São Paulo. The NN had 5 output neurons which were associated with 5 typical corneal shapes: keratoconus, with-the-rule astigmatism, against-the-rule astigmatism, "regular" or "normal" shape and post-PRK. RESULTS: The NN and DA responses were statistically analyzed in terms of precision ([true positive+true negative]/total number of cases). Mean overall results for all cases for the NN and DA techniques were, respectively, 94% and 84.8%. CONCLUSION: Although we used a relatively small database, results obtained in the present study indicate that Zernike polynomials as descriptors of corneal shape may be a reliable parameter as input data for diagnostic automation of VK maps, using either NN or DA.

Dose de exposição radiométrica de granitos do estado de Rondônia, Brasil

Este trabalho avaliou a concentração dos radioelementos K, eU e eTh em amostras de granitos do Estado de Rondônia, Brasil. A análise estatística dos dados obtidos indicou que eles seguem distribuições lognormais. Os valores modais encontrados correspondem a cerca de 11% para K, 29 ppm para eU e 85 ppm para eTh. Correlações diretas significativas foram determinadas entre as concentrações dos três radioelementos, isto é, r = 0,71 (entre K e eU), r = 0,72 (entre K e eTh)e r = 0,72 (entre eU e eTh), sugerindo que são congruentes os processos de seu acúmulo nos minerais das rochas analisadas. Os dados de concentração permitiram estimar a taxa de dose absorvida de radiação no ar acima de 1 m do nível do terreno, a qual também segue uma distribuição lognormal, com valor modal de 2,7 mSv/ano, que é ligeiramente superior à média global de 2,4 mSv/ano. Os resultados obtidos também permitiram avaliar, do ponto de vista radiométrico, se os granitos analisados são adequados para emprego como revestimento em construção civil.

Geographical and seasonal distributions of the seedeaters Sporophila bouvreuil and Sporophila pileata (aves: Emberizidae)

Many species in the genus Sporophila are migratory. Migration patterns, while poorly studied, may be influenced by seed production which can be very seasonal in some regions. The distribution of S. bouvreuil extends from the Amazon and Suriname south through a large part of the open regions of Brazil. Sporophila pileata, on the other hand, is found in southeastern and southern Brazil as well as Argentina and Paraguay. Both of these species migrate, but their movement patterns are poorly known. To better understand the geographical and the seasonal distributions of S. bouvreuil and S. pileata, we grouped the records into two categories: the breeding season (September to March) and the putative migration season (April to August). We found two areas of sympatry between S. bouvreuil and S. pileata in the Brazilian states of Minas Gerais and São Paulo. For S. bouvreuil we suggest that populations that breed in the Amazon migrate to the Cerrado or Caatinga, where they will encounter resident populations of the same species. These resident populations may take part in short distance migrations. Sporophila pileata, on the other hand, occur in the Cerrado and open areas within the Atlantic Forest and it is not yet possible to determine migratory tendencies or destinations in the non-breeding season.