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.

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.

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.

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.

Clinical trial design in non-invasive brain stimulation psychiatric research

Major depressive disorder (MDD) trials - investigating either non-pharmacological or pharmacological interventions - have shown mixed results. Many reasons explain this heterogeneity, but one that stands out is the trial design due to specific challenges in the field. We aimed therefore to review the methodology of non-invasive brain stimulation (NIBS) trials and provide a framework to improve clinical trial design. We performed a systematic review for randomized, controlled MDD trials whose intervention was transcranial magnetic stimulation (rTMS) or transcranial direct current stimulation (tDCS) in MEDLINE and other databases from April 2002 to April 2008. We created an unstructured checklist based on CONSORT guidelines to extract items such as power analysis, sham method, blinding assessment, allocation concealment, operational criteria used for MDD, definition of refractory depression and primary study hypotheses. Thirty-one studies were included. We found that the main methodological issues can be divided in to three groups: (1) issues related to phase II/small trials, (2) issues related to MDD trials and, (3) specific issues of NIBS studies. Taken together, they can threaten study validity and lead to inconclusive results. Feasible solutions include: estimating the sample size a priori; measuring the degree of refractoriness of the subjects; specifying the primary hypothesis and statistical tests; controlling predictor variables through stratification randomization methods or using strict eligibility criteria; adjusting the study design to the target population; using adaptive designs and exploring NIBS efficacy employing biological markers. In conclusion, our study summarizes the main methodological issues of NIBS trials and proposes a number of alternatives to manage them. Copyright (C) 2011 John Wiley & Sons, Ltd.

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.

Evidence-based guidelines for interpreting change scores for the European Organisation for the Research and Treatment of Cancer Quality of Life Questionnaire Core 30

Aim: To use published literature and experts' opinion to investigate the clinical meaning and magnitude of changes in the Quality of Life (QOL) of groups of patients measured with the European Organisation for the Research and Treatment of Cancer Quality of Life Questionnaire Core 30 (EORTC QLQ-C30). Methods: An innovative method combining systematic review of published studies, expert opinions and meta-analysis was used to estimate large, medium, and small mean changes over time for QLQ-C30 scores. Results: Nine hundred and eleven papers were identified, leading to 118 relevant papers. One thousand two hundred and thirty two mean changes in QOL over time were combined in the meta-analysis, with timescales ranging from four days to five years. Guidelines were produced for trivial, small, and medium size classes, for each subscale and for improving and declining scores separately. Estimates for improvements were smaller than respective estimates for declines. Conclusions: These guidelines can be used to aid sample size calculations and interpretation of mean changes over time from groups of patients. Observed mean changes in the QLQ-C30 scores are generally small in most clinical situations, possibly due to response shift. Careful consideration is needed when planning studies where QOL changes over time are of primary interest; the timing of follow up, sample attrition, direction of QOL changes, and subscales of primary interest are key considerations. (C) 2012 Elsevier Ltd. All rights reserved.

Does ragging play a role in medical student depression - Cause or effect?

Background: Medical students experience a lot of stress what may contribute to symptoms of depression. In this study we set out to look at the environmental factors which may be contributing in one medical school in Brazil. Methods: We assessed depressive symptoms using Beck's Depression Inventory in 465 and 267 medical students in 2001 and 2006 respectively. We explored possible social and environmental causes using qualitative data. Results: Nearly 15% scored above the cut off for depression in both the samples. Males in the pre-clinical stage in 2006 showed an increase in depressive symptoms than males in the same cycle in 2001 (aOR = 7.36 [95% CI = 0.85-63.5] p = 0.07). Qualitative data confirmed that factors such as ragging and low social involvement were correlated with depressive symptoms in pre-clinical stage males. Limitations: The sample size was small both for quantitative and qualitative aspects of the study. Conclusions: It appears that ragging plays an important role in the genesis of depressive symptoms in medical students. (C) 2012 Elsevier B.V. All rights reserved.

Social phobia in obsessive-compulsive disorder: Prevalence and correlates

Background: Social Phobia (SP) is an anxiety disorder that frequently co-occurs with obsessive-compulsive disorder (OCD); however, studies that evaluate clinical factors associated with this specific comorbidity are rare. The aim was to estimate the prevalence of SP in a large multicenter sample of OCD patients and compare the characteristics of individuals with and without SP. Method: A cross-sectional study with 1001 patients of the Brazilian Research Consortium on Obsessive-Compulsive Spectrum Disorders using several assessment instruments, including the Dimensional Yale-Brown Obsessive-Compulsive Scale and the Structured Clinical Interview for DSM-IV Axis I Disorders. Univariate analyses were followed by logistic regression. Results: Lifetime prevalence of SP was 34.6% (N=346). The following variables remained associated with SP comorbidity after logistic regression: male sex, lower socioeconomic status, body dysmorphic disorder, specific phobia, dysthymia, generalized anxiety disorder, agoraphobia, Tourette syndrome and binge eating disorder. Limitations: The cross-sectional design does not permit the inference of causal relationships; some retrospective information may have been subject to recall bias; all patients were being treated in tertiary services, therefore generalization of the results to other samples of OCD sufferers should be cautious. Despite the large sample size, some hypotheses may not have been confirmed due to the small number of cases with these characteristics (type 2 error). Conclusion: SP is frequent among OCD patients and co-occurs with other disorders that have common phenomenological features. These findings have important implications for clinical practice, indicating the need for broader treatment approaches for individuals with this profile. (C) 2012 Elsevier B.V. All rights reserved.

The impact of the CACNA1C risk allele on limbic structures and facial emotions recognition in bipolar disorder subjects and healthy controls

Introduction: Impairments in facial emotion recognition (PER) have been reported in bipolar disorder (BD) during all mood states. FER has been the focus of functional magnetic resonance imaging studies evaluating differential activation of limbic regions. Recently, the alpha 1-C subunit of the L-type voltage-gated calcium channel (CACNA1C) gene has been described as a risk gene for BD and its Met allele found to increase CACNA1C mRNA expression. In healthy controls, the CACNA1C risk (Met) allele has been reported to increase limbic system activation during emotional stimuli and also to impact on cognitive function. The aim of this study was to investigate the impact of CACNA1C genotype on FER scores and limbic system morphology in subjects with BD and healthy controls. Material and methods: Thirty-nine euthymic BD I subjects and 40 healthy controls were submitted to a PER recognition test battery and genotyped for CACNA1C. Subjects were also examined with a 3D 3-Tesla structural imaging protocol. Results: The CACNA1C risk allele for BD was associated to FER impairment in BD, while in controls nothing was observed. The CACNA1C genotype did not impact on amygdala or hippocampus volume neither in BD nor controls. Limitations: Sample size. Conclusion: The present findings suggest that a polymorphism in calcium channels interferes FER phenotype exclusively in BD and doesn't interfere on limbic structures morphology. (C) 2012 Elsevier B.V. All rights reserved.

The effect of noise on the voice of preschool institution educators

OBJECTIVE: To analyze the association between noise levels present in preschool institutions and vocal disorders among educators. METHODS: Cross-sectional study conducted in 2009 with 28 teachers from three preschool institutions located in the city of Sao Paulo (Southeastern Brazil). Sound pressure levels were measured according to Brazilian Technical Standards Association, with the use of a sound level meter. The averages were classified according to the levels of comfort, discomfort, and auditory damage proposed by the Pan American Health Organization. The educators underwent voice evaluation: self-assessment with visual analogue scale, auditory perceptual evaluation using the GRBAS scale, and acoustic analysis utilizing the Praat program. To analyze the association between noise and voice evaluation, descriptive statistics and the chi-square test were employed, with significance of 10% due to sample size. RESULTS: The teachers' age ranged between 21 and 56 years. The noise average was 72.7 dB, considered as damage 2. The professionals' vocal self-assessment ranked an average of 5.1 on the scale, being considered as moderate alteration. In the auditory-perceptual assessment, 74% presented vocal alteration, especially hoarseness; of these, 52% were considered mild alterations. In the acoustic assessment the majority presented fundamental frequency below the expected level. Averages for jitter, shimmer and harmonic-noise ratio showed alterations. An association between the presence of noise between the harmonics and vocal disorders was observed. CONCLUSIONS: There is an association between presence of noise between the harmonics and vocal alteration, with high noise levels. Although most teachers presented mild voice alteration, the self-evaluation showed moderate alteration, probably due to the difficulty in projection.

ABSENCE OF STABILIZATION POINT OF THE SPECIES ACCUMULATION CURVE IN TROPICAL FORESTS

The definition of the sample size is a major problem in studies of phytosociology. The species accumulation curve is used to define the sampling sufficiency, but this method presents some limitations such as the absence of a stabilization point that can be objectively determined and the arbitrariness of the order of sampling units in the curve. A solution to this problem is the use of randomization procedures, e. g. permutation, for obtaining a mean species accumulation curve and empiric confidence intervals. However, the randomization process emphasizes the asymptotical character of the curve. Moreover, the inexistence of an inflection point in the curve makes it impossible to define objectively the point of optimum sample size.

On Wald Residuals in Generalized Linear Models

A rigorous asymptotic theory for Wald residuals in generalized linear models is not yet available. The authors provide matrix formulae of order O(n(-1)), where n is the sample size, for the first two moments of these residuals. The formulae can be applied to many regression models widely used in practice. The authors suggest adjusted Wald residuals to these models with approximately zero mean and unit variance. The expressions were used to analyze a real dataset. Some simulation results indicate that the adjusted Wald residuals are better approximated by the standard normal distribution than the Wald residuals.