The Endothelial Sample Size Analysis in Corneal Specular Microscopy Clinical Examinations

Purpose: To evaluate endothelial cell sample size and statistical error in corneal specular microscopy (CSM) examinations. Methods: One hundred twenty examinations were conducted with 4 types of corneal specular microscopes: 30 with each BioOptics, CSO, Konan, and Topcon corneal specular microscopes. All endothelial image data were analyzed by respective instrument software and also by the Cells Analyzer software with a method developed in our lab(US Patent). A reliability degree (RD) of 95% and a relative error (RE) of 0.05 were used as cut-off values to analyze images of the counted endothelial cells called samples. The sample size mean was the number of cells evaluated on the images obtained with each device. Only examinations with RE<0.05 were considered statistically correct and suitable for comparisons with future examinations. The Cells Analyzer software was used to calculate the RE and customized sample size for all examinations. Results: Bio-Optics: sample size, 97 +/- 22 cells; RE, 6.52 +/- 0.86; only 10% of the examinations had sufficient endothelial cell quantity (RE<0.05); customized sample size, 162 +/- 34 cells. CSO: sample size, 110 +/- 20 cells; RE, 5.98 +/- 0.98; only 16.6% of the examinations had sufficient endothelial cell quantity (RE<0.05); customized sample size, 157 +/- 45 cells. Konan: sample size, 80 +/- 27 cells; RE, 10.6 +/- 3.67; none of the examinations had sufficient endothelial cell quantity (RE>0.05); customized sample size, 336 +/- 131 cells. Topcon: sample size, 87 +/- 17 cells; RE, 10.1 +/- 2.52; none of the examinations had sufficient endothelial cell quantity (RE>0.05); customized sample size, 382 +/- 159 cells. Conclusions: A very high number of CSM examinations had sample errors based on Cells Analyzer software. The endothelial sample size (examinations) needs to include more cells to be reliable and reproducible. The Cells Analyzer tutorial routine will be useful for CSM examination reliability and reproducibility.

Parameter recovery for a skew-normal IRT model under a Bayesian approach: hierarchical framework, prior and kernel sensitivity and sample size

Item response theory (IRT) comprises a set of statistical models which are useful in many fields, especially when there is an interest in studying latent variables (or latent traits). Usually such latent traits are assumed to be random variables and a convenient distribution is assigned to them. A very common choice for such a distribution has been the standard normal. Recently, Azevedo et al. [Bayesian inference for a skew-normal IRT model under the centred parameterization, Comput. Stat. Data Anal. 55 (2011), pp. 353-365] proposed a skew-normal distribution under the centred parameterization (SNCP) as had been studied in [R. B. Arellano-Valle and A. Azzalini, The centred parametrization for the multivariate skew-normal distribution, J. Multivariate Anal. 99(7) (2008), pp. 1362-1382], to model the latent trait distribution. This approach allows one to represent any asymmetric behaviour concerning the latent trait distribution. Also, they developed a Metropolis-Hastings within the Gibbs sampling (MHWGS) algorithm based on the density of the SNCP. They showed that the algorithm recovers all parameters properly. Their results indicated that, in the presence of asymmetry, the proposed model and the estimation algorithm perform better than the usual model and estimation methods. Our main goal in this paper is to propose another type of MHWGS algorithm based on a stochastic representation (hierarchical structure) of the SNCP studied in [N. Henze, A probabilistic representation of the skew-normal distribution, Scand. J. Statist. 13 (1986), pp. 271-275]. Our algorithm has only one Metropolis-Hastings step, in opposition to the algorithm developed by Azevedo et al., which has two such steps. This not only makes the implementation easier but also reduces the number of proposal densities to be used, which can be a problem in the implementation of MHWGS algorithms, as can be seen in [R.J. Patz and B.W. Junker, A straightforward approach to Markov Chain Monte Carlo methods for item response models, J. Educ. Behav. Stat. 24(2) (1999), pp. 146-178; R. J. Patz and B. W. Junker, The applications and extensions of MCMC in IRT: Multiple item types, missing data, and rated responses, J. Educ. Behav. Stat. 24(4) (1999), pp. 342-366; A. Gelman, G.O. Roberts, and W.R. Gilks, Efficient Metropolis jumping rules, Bayesian Stat. 5 (1996), pp. 599-607]. Moreover, we consider a modified beta prior (which generalizes the one considered in [3]) and a Jeffreys prior for the asymmetry parameter. Furthermore, we study the sensitivity of such priors as well as the use of different kernel densities for this parameter. Finally, we assess the impact of the number of examinees, number of items and the asymmetry level on the parameter recovery. Results of the simulation study indicated that our approach performed equally as well as that in [3], in terms of parameter recovery, mainly using the Jeffreys prior. Also, they indicated that the asymmetry level has the highest impact on parameter recovery, even though it is relatively small. A real data analysis is considered jointly with the development of model fitting assessment tools. The results are compared with the ones obtained by Azevedo et al. The results indicate that using the hierarchical approach allows us to implement MCMC algorithms more easily, it facilitates diagnosis of the convergence and also it can be very useful to fit more complex skew IRT models.

A Single Statistic for Monitoring the Covariance Matrix of Bivariate Processes

In this article, we present a new control chart for monitoring the covariance matrix in a bivariate process. In this method, n observations of the two variables were considered as if they came from a single variable (as a sample of 2n observations), and a sample variance was calculated. This statistic was used to build a new control chart specifically as a VMIX chart. The performance of the new control chart was compared with its main competitors: the generalized sampled variance chart, the likelihood ratio test, Nagao's test, probability integral transformation (v(t)), and the recently proposed VMAX chart. Among these statistics, only the VMAX chart was competitive with the VMIX chart. For shifts in both variances, the VMIX chart outperformed VMAX; however, VMAX showed better performance for large shifts (higher than 10%) in one variance.

Variable patterns of Y chromosome homology in Akodontini rodents (Sigmodontinae): a phylogenetic signal revealed by chromosome painting

The Akodontini is the second most speciose tribe of sigmodontine rodents, one of the most diverse groups of neotropical mammals. Molecular phylogenetic analyses are discordant regarding the interrelationships of genera, with low support for some clades. However, two clades are concordant, one (clade A) with Akodon sensu strictu (excluding Akodon serrensis), "Akodon" serrensis, Bibimys, Deltamys, Juscelinomys, Necromys, Oxymycterus, Podoxymys, Thalpomys and Thaptomys, and another (clade B) with Blarinomys, Brucepattersonius, Kunsia, Lenoxus and Scapteromys. Here, we present chromosome painting using Akodon paranaensis (APA) Y paint, after suppression of simple repetitive sequences, on ten Akodontini genera. Partial Y chromosome homology, in addition to the homology already reported on the Akodon genus, was detected on the Y chromosomes of "A." serrensis, Thaptomys, Deltamys, Necromys and Thalpomys and on Y and X chromosomes in Oxymycterus. In Blarinomys, Brucepattersonius, Scapteromys and Kunsia, no APA Y signal was observed using different hybridization conditions; APA X paint gave positive signals only on the X chromosome in all genera. The Y chromosome homology was variable in size and positioning among the species studied as follow: (1) whole acrocentric Y chromosome in Akodon and "A." serrensis, (2) Yp and pericentromeric region in submetacentric Y of Necromys and Thaptomys, (3) pericentromeric region in acrocentric Y of Deltamys, (4) distal Yq in the acrocentric Y chromosome of Thalpomys and (5) proximal Yq in the acrocentric Y and Xp in the basal clade A genus Oxymycterus. The results suggest that the homology involves pairing (pseudoautosomal) and additional regions that have undergone rearrangement during divergence. The widespread Y homology represents a phylogenetic signal in Akodontini that provides additional evidence supporting the monophyly of clade A. The findings also raise questions about the evolution of the pseudoautosomal region observed in Oxymycterus. The Y chromosomes of these closely related species seem to have undergone dynamic rearrangements, including restructuring and reduction of homologous segments. Furthermore, the changes observed may indicate progressive attrition of the Y chromosome in more distantly related species.

Integrating on-line process control and imperfect corrective maintenance: An economical design

Existing studies of on-line process control are concerned with economic aspects, and the parameters of the processes are optimized with respect to the average cost per item produced. However, an equally important dimension is the adoption of an efficient maintenance policy. In most cases, only the frequency of the corrective adjustment is evaluated because it is assumed that the equipment becomes "as good as new" after corrective maintenance. For this condition to be met, a sophisticated and detailed corrective adjustment system needs to be employed. The aim of this paper is to propose an integrated economic model incorporating the following two dimensions: on-line process control and a corrective maintenance program. Both performances are objects of an average cost per item minimization. Adjustments are based on the location of the measurement of a quality characteristic of interest in a three decision zone. Numerical examples are illustrated in the proposal. (c) 2012 Elsevier B.V. All rights reserved.

ANALYSIS OF THE BALANCE CONTROL IN SURFERS DURING THE ERECT POSTURE

Background: Surfing is a sport that has become considerably popular, which increased interest in research about the aspects that can influence on the performance of these athletes, such as injuries, aerobic fitness and reaction time. Due to the ever-changing environment and high instability required for surfing, the surfers must develop some neuromuscular skills (agility, balance, muscle strength and flexibility) to acquire better performance in this modality. Nevertheless, there are still few scientific studies concerned about the investigation of these motor skills in surfing. Objective: The aim of this study was to evaluate the balance control in surfers compared to practitioners of other physical activities. Methods: Participants remained on a force platform while performing tasks involving visual deprivation (eyes open or closed) and somatosensory disturbance (steady surface or use of foam), with covariation of experimental conditions. The following variables were analyzed: speed and root mean square (RMS) displacement of the center of pressure in the anteroposterior (AP) and mediolateral (ML) directions. Results: The results showed no difference between groups during the experimental conditions, that is to say, both surfers and the control group varied over the conditions of eyes closed and on foam. Conclusion: Although surfing requires the surfer to have great balance control, the results did not reveal a relationship between this sport and better performance in balance control. However, we must consider the small sample size and the fact that this sport requires dynamic balance, while the study evaluated static balance.

Brazilian multicenter study on prevalence of preterm birth and associated factors

Abstract Background The occurrence of preterm birth remains a complex public health condition. It is considered the main cause of neonatal morbidity and mortality, resulting in a high likelihood of sequelae in surviving children. With variable incidence in several countries, it has grown markedly in the last decades. In Brazil, however, there are still difficulties to estimate its real occurrence. Therefore, it is essential to establish the prevalence and causes of this condition in order to propose prevention actions. This study intend to collect information from hospitals nationwide on the prevalence of preterm births, their associated socioeconomic and environmental factors, diagnostic and treatment methods resulting from causes such as spontaneous preterm labor, prelabor rupture of membranes, and therapeutic preterm birth, as well as neonatal results. Methods/Design This proposal is a multicenter cross-sectional study plus a nested case-control study, to be implemented in 27 reference obstetric centers in several regions of Brazil (North: 1; Northeast: 10; Central-west: 1; Southeast: 13; South: 2). For the cross sectional component, the participating centers should perform, during a period of six months, a prospective surveillance of all patients hospitalized to give birth, in order to identify preterm birth cases and their main causes. In the first three months of the study, an analysis of the factors associated with preterm birth will also be carried out, comparing women who have preterm birth with those who deliver at term. For the prevalence study, 37,000 births will be evaluated (at term and preterm), corresponding to approximately half the deliveries of all participating centers in 12 months. For the case-control study component, the estimated sample size is 1,055 women in each group (cases and controls). The total number of preterm births estimated to be followed in both components of the study is around 3,600. Data will be collected through a questionnaire all patients will answer after delivery. The data will then be encoded in an electronic form and sent online by internet to a central database. The data analysis will be carried out by subgroups according to gestational age at preterm birth, its probable causes, therapeutic management, and neonatal outcomes. Then, the respective rates, ratios and relative risks will be estimated for the possible predictors. Discussion These findings will provide information on preterm births in Brazil and their main social and biological risk factors, supporting health policies and the implementation of clinical trials on preterm birth prevention and treatment strategies, a condition with many physical and emotional consequences to children and their families.

A search for anisotropy in the arrival directions of ultra high energy cosmic rays recorded at the Pierre Auger Observatory

Observations of cosmic rays arrival directions made with the Pierre Auger Observatory have previously provided evidence of anisotropy at the 99% CL using the correlation of ultra high energy cosmic rays (UHECRs) with objects drawn from the Veron-Cetty Veron catalog. In this paper we report on the use of three catalog independent methods to search for anisotropy. The 2pt-L, 2pt+ and 3pt methods, each giving a different measure of self-clustering in arrival directions, were tested on mock cosmic ray data sets to study the impacts of sample size and magnetic smearing on their results, accounting for both angular and energy resolutions. If the sources of UHECRs follow the same large scale structure as ordinary galaxies in the local Universe and if UHECRs are deflected no more than a few degrees, a study of mock maps suggests that these three method can efficiently respond to the resulting anisotropy with a P-value = 1.0% or smaller with data sets as few as 100 events. using data taken from January 1, 2004 to July 31, 2010 we examined the 20, 30, ... , 110 highest energy events with a corresponding minimum energy threshold of about 49.3 EeV. The minimum P-values found were 13.5% using the 2pt-L method, 1.0% using the 2pt+ method and 1.1% using the 3pt method for the highest 100 energy events. In view of the multiple (correlated) scans performed on the data set, these catalog-independent methods do not yield strong evidence of anisotropy in the highest energy cosmic rays.

Testing allele homogeneity: the problem of nested hypotheses

Background: The evaluation of associations between genotypes and diseases in a case-control framework plays an important role in genetic epidemiology. This paper focuses on the evaluation of the homogeneity of both genotypic and allelic frequencies. The traditional test that is used to check allelic homogeneity is known to be valid only under Hardy-Weinberg equilibrium, a property that may not hold in practice. Results: We first describe the flaws of the traditional (chi-squared) tests for both allelic and genotypic homogeneity. Besides the known problem of the allelic procedure, we show that whenever these tests are used, an incoherence may arise: sometimes the genotypic homogeneity hypothesis is not rejected, but the allelic hypothesis is. As we argue, this is logically impossible. Some methods that were recently proposed implicitly rely on the idea that this does not happen. In an attempt to correct this incoherence, we describe an alternative frequentist approach that is appropriate even when Hardy-Weinberg equilibrium does not hold. It is then shown that the problem remains and is intrinsic of frequentist procedures. Finally, we introduce the Full Bayesian Significance Test to test both hypotheses and prove that the incoherence cannot happen with these new tests. To illustrate this, all five tests are applied to real and simulated datasets. Using the celebrated power analysis, we show that the Bayesian method is comparable to the frequentist one and has the advantage of being coherent. Conclusions: Contrary to more traditional approaches, the Full Bayesian Significance Test for association studies provides a simple, coherent and powerful tool for detecting associations.

CONTEXT TREE SELECTION AND LINGUISTIC RHYTHM RETRIEVAL FROM WRITTEN TEXTS

The starting point of this article is the question "How to retrieve fingerprints of rhythm in written texts?" We address this problem in the case of Brazilian and European Portuguese. These two dialects of Modern Portuguese share the same lexicon and most of the sentences they produce are superficially identical. Yet they are conjectured, on linguistic grounds, to implement different rhythms. We show that this linguistic question can be formulated as a problem of model selection in the class of variable length Markov chains. To carry on this approach, we compare texts from European and Brazilian Portuguese. These texts are previously encoded according to some basic rhythmic features of the sentences which can be automatically retrieved. This is an entirely new approach from the linguistic point of view. Our statistical contribution is the introduction of the smallest maximizer criterion which is a constant free procedure for model selection. As a by-product, this provides a solution for the problem of optimal choice of the penalty constant when using the BIC to select a variable length Markov chain. Besides proving the consistency of the smallest maximizer criterion when the sample size diverges, we also make a simulation study comparing our approach with both the standard BIC selection and the Peres-Shields order estimation. Applied to the linguistic sample constituted for our case study, the smallest maximizer criterion assigns different context-tree models to the two dialects of Portuguese. The features of the selected models are compatible with current conjectures discussed in the linguistic literature.

Decisions on Temporal Sampling Protocol Influence the Detection of Ecological Patterns

Within-site variability in species detectability is a problem common to many biodiversity assessments and can strongly bias the results. Such variability can be caused by many factors, including simple counting inaccuracies, which can be solved by increasing sample size, or by temporal changes in species behavior, meaning that the way the temporal sampling protocol is designed is also very important. Here we use the example of mist-netted tropical birds to determine how design decisions in the temporal sampling protocol can alter the data collected and how these changes might affect the detection of ecological patterns, such as the species-area relationship (SAR). Using data from almost 3400 birds captured from 21,000 net-hours at 31 sites in the Brazilian Atlantic Forest, we found that the magnitude of ecological trends remained fairly stable, but the probability of detecting statistically significant ecological patterns varied depending on sampling effort, time of day and season in which sampling was conducted. For example, more species were detected in the wet season, but the SAR was strongest in the dry season. We found that the temporal distribution of sampling effort was more important than its total amount, discovering that similar ecological results could have been obtained with one-third of the total effort, as long as each site had been equally sampled over 2 yr. We discuss that projects with the same sampling effort and spatial design, but with different temporal sampling protocol are likely to report different ecological patterns, which may ultimately lead to inappropriate conservation strategies.