Missing data mechanisms and their implications on the analysis of categorical data

We review some issues related to the implications of different missing data mechanisms on statistical inference for contingency tables and consider simulation studies to compare the results obtained under such models to those where the units with missing data are disregarded. We confirm that although, in general, analyses under the correct missing at random and missing completely at random models are more efficient even for small sample sizes, there are exceptions where they may not improve the results obtained by ignoring the partially classified data. We show that under the missing not at random (MNAR) model, estimates on the boundary of the parameter space as well as lack of identifiability of the parameters of saturated models may be associated with undesirable asymptotic properties of maximum likelihood estimators and likelihood ratio tests; even in standard cases the bias of the estimators may be low only for very large samples. We also show that the probability of a boundary solution obtained under the correct MNAR model may be large even for large samples and that, consequently, we may not always conclude that a MNAR model is misspecified because the estimate is on the boundary of the parameter space.

Dynamic Changes in the Mental Rotation Network Revealed by Pattern Recognition Analysis of fMRI Data

We investigated the temporal dynamics and changes in connectivity in the mental rotation network through the application of spatio-temporal support vector machines (SVMs). The spatio-temporal SVM [Mourao-Miranda, J., Friston, K. J., et al. (2007). Dynamic discrimination analysis: A spatial-temporal SVM. Neuroimage, 36, 88-99] is a pattern recognition approach that is suitable for investigating dynamic changes in the brain network during a complex mental task. It does not require a model describing each component of the task and the precise shape of the BOLD impulse response. By defining a time window including a cognitive event, one can use spatio-temporal fMRI observations from two cognitive states to train the SVM. During the training, the SVM finds the discriminating pattern between the two states and produces a discriminating weight vector encompassing both voxels and time (i.e., spatio-temporal maps). We showed that by applying spatio-temporal SVM to an event-related mental rotation experiment, it is possible to discriminate between different degrees of angular disparity (0 degrees vs. 20 degrees, 0 degrees vs. 60 degrees, and 0 degrees vs. 100 degrees), and the discrimination accuracy is correlated with the difference in angular disparity between the conditions. For the comparison with highest accuracy (08 vs. 1008), we evaluated how the most discriminating areas (visual regions, parietal regions, supplementary, and premotor areas) change their behavior over time. The frontal premotor regions became highly discriminating earlier than the superior parietal cortex. There seems to be a parcellation of the parietal regions with an earlier discrimination of the inferior parietal lobe in the mental rotation in relation to the superior parietal. The SVM also identified a network of regions that had a decrease in BOLD responses during the 100 degrees condition in relation to the 0 degrees condition (posterior cingulate, frontal, and superior temporal gyrus). This network was also highly discriminating between the two conditions. In addition, we investigated changes in functional connectivity between the most discriminating areas identified by the spatio-temporal SVM. We observed an increase in functional connectivity between almost all areas activated during the 100 degrees condition (bilateral inferior and superior parietal lobe, bilateral premotor area, and SMA) but not between the areas that showed a decrease in BOLD response during the 100 degrees condition.

Visual analysis of image collections

Multidimensional Visualization techniques are invaluable tools for analysis of structured and unstructured data with variable dimensionality. This paper introduces PEx-Image-Projection Explorer for Images-a tool aimed at supporting analysis of image collections. The tool supports a methodology that employs interactive visualizations to aid user-driven feature detection and classification tasks, thus offering improved analysis and exploration capabilities. The visual mappings employ similarity-based multidimensional projections and point placement to layout the data on a plane for visual exploration. In addition to its application to image databases, we also illustrate how the proposed approach can be successfully employed in simultaneous analysis of different data types, such as text and images, offering a common visual representation for data expressed in different modalities.

Spatial analysis of urban violence based on emergency room data

OBJECTIVE: To estimate the spatial intensity of urban violence events using wavelet-based methods and emergency room data. METHODS: Information on victims attended at the emergency room of a public hospital in the city of São Paulo, Southeastern Brazil, from January 1, 2002 to January 11, 2003 were obtained from hospital records. The spatial distribution of 3,540 events was recorded and a uniform random procedure was used to allocate records with incomplete addresses. Point processes and wavelet analysis technique were used to estimate the spatial intensity, defined as the expected number of events by unit area. RESULTS: Of all georeferenced points, 59% were accidents and 40% were assaults. There is a non-homogeneous spatial distribution of the events with high concentration in two districts and three large avenues in the southern area of the city of São Paulo. CONCLUSIONS: Hospital records combined with methodological tools to estimate intensity of events are useful to study urban violence. The wavelet analysis is useful in the computation of the expected number of events and their respective confidence bands for any sub-region and, consequently, in the specification of risk estimates that could be used in decision-making processes for public policies.

Short report: historical analysis of a near disaster: Anopheles gambiae in Brazil

Attributed to human-mediated dispersal, a species of the Anopheles gambiae complex invaded northeastern Brazil in 1930. This event is considered unique among the intercontinental introductions of disease vectors and the most serious one: "Few threats to the future health of the Americas have equalled that inherent in the invasion of Brazil, in 1930, by Anopheles gambiae." Because it was only in the 1960s that An. gambiae was recognized as a species complex now including seven species, the precise species identity of the Brazilian invader remains a mystery. Here we used historical DNA analysis of museum specimens, collected at the time of invasion from Brazil, and aimed at the identification of the Brazilian invader. Our results identify the arid-adapted Anopheles arabiensis as being the actual invading species. Establishing the identity of the species, in addition to being of intrinsic historical interest, can inform future threats of this sort especially in a changing environment. Furthermore, these results highlight the potential danger of human-mediated range expansions of insect disease vectors and the importance of museum collections in retrieving historical information

Qualitative Analysis of the Interdisciplinary Interaction between Data Analysis Specialists and Novice Clinical Researchers

Background: The inherent complexity of statistical methods and clinical phenomena compel researchers with diverse domains of expertise to work in interdisciplinary teams, where none of them have a complete knowledge in their counterpart's field. As a result, knowledge exchange may often be characterized by miscommunication leading to misinterpretation, ultimately resulting in errors in research and even clinical practice. Though communication has a central role in interdisciplinary collaboration and since miscommunication can have a negative impact on research processes, to the best of our knowledge, no study has yet explored how data analysis specialists and clinical researchers communicate over time. Methods/Principal Findings: We conducted qualitative analysis of encounters between clinical researchers and data analysis specialists (epidemiologist, clinical epidemiologist, and data mining specialist). These encounters were recorded and systematically analyzed using a grounded theory methodology for extraction of emerging themes, followed by data triangulation and analysis of negative cases for validation. A policy analysis was then performed using a system dynamics methodology looking for potential interventions to improve this process. Four major emerging themes were found. Definitions using lay language were frequently employed as a way to bridge the language gap between the specialties. Thought experiments presented a series of ""what if'' situations that helped clarify how the method or information from the other field would behave, if exposed to alternative situations, ultimately aiding in explaining their main objective. Metaphors and analogies were used to translate concepts across fields, from the unfamiliar to the familiar. Prolepsis was used to anticipate study outcomes, thus helping specialists understand the current context based on an understanding of their final goal. Conclusion/Significance: The communication between clinical researchers and data analysis specialists presents multiple challenges that can lead to errors.

Dispersive analysis of KL mu 3 and KLe3 scalar and vector form factors using KTeV data

Using the published KTeV samples of K(L) -> pi(+/-)e(-/+)nu and K(L) -> pi(+/-)mu(-/+)nu decays, we perform a reanalysis of the scalar and vector form factors based on the dispersive parametrization. We obtain phase-space integrals I(K)(e) = 0.15446 +/- 0.00025 and I(K)(mu) = 0.10219 +/- 0.00025. For the scalar form factor parametrization, the only free parameter is the normalized form factor value at the Callan-Treiman point (C); our best-fit results in InC = 0.1915 +/- 0.0122. We also study the sensitivity of C to different parametrizations of the vector form factor. The results for the phase-space integrals and C are then used to make tests of the standard model. Finally, we compare our results with lattice QCD calculations of F(K)/F(pi) and f(+)(0).

Artificial neural networks analysis of Sao Paulo subway tunnel settlement data

Artificial neural networks have been used to analyze a number of engineering problems, including settlement caused by different tunneling methods in various types of ground mass. This paper focuses on settlement over shotcrete- supported tunnels on Sao Paulo subway line 2 (West Extension) that were excavated in Tertiary sediments using the sequential excavation method. The adjusted network is a good tool for predicting settlement above new tunnels to be excavated in similar conditions. The influence of network training parameters on the quality of results is also discussed. (C) 2007 Elsevier Ltd. All rights reserved.

Short report: Historical analysis of a near disaster: Anopheles gambiae in Brazil

Attributed to human-mediated dispersal, a species of the Anopheles gambiae complex invaded northeastern Brazil in 1930. This event is considered unique among the intercontinental introductions of disease vectors and the most serious one: ""Few threats to the future health of the Americas have equalled that inherent in the invasion of Brazil, in 1930, by Anopheles gambiae."" Because it was only in the 1960s that An. gambiae was recognized as a species complex now including seven species, the precise species identity of the Brazilian invader remains a mystery. Here we used historical DNA analysis of museum specimens, collected at the time of invasion from Brazil, and aimed at the identification of the Brazilian invader. Our results identify the arid-adapted Anopheles arabiensis as being the actual invading species. Establishing the identity of the species, in addition to being of intrinsic historical interest, can inform future threats of this sort especially in a changing environment. Furthermore, these results highlight the potential danger of human-mediated range expansions of insect disease vectors and the importance of museum collections in retrieving historical information.

A Comparative Cytogenetic Analysis of 2 Bothriuridae Species and Overview of the Chromosome Data of Scorpiones

The order Scorpiones is one of the most cytogenetically interesting groups within Arachnida by virtue of the combination of chromosome singularities found in the 59 species analyzed so far. In this work, mitotic and meiotic chromosomes of 2 species of the family Bothriuridae were detailed. This family occupies a basal position within the superfamily Scorpionoidea. Furthermore, review of the cytogenetic data of all previously studied scorpions is presented. Light microscopy chromosome analysis showed that Bothriurus araguayae and Bothriurus rochensis possess low diploid numbers compared with those of species belonging to closely related families. Gonadal cells examined under light and in transmission electron microscopy revealed, for the first time, that the Bothriuridae species possess typical monocentric chromosomes, and male meiosis presented chromosomes with synaptic and achiasmatic behavior. Moreover, in the sample of B. araguayae studied, heterozygous translocations were verified. The use of techniques to highlight specific chromosomal regions also revealed additional differences between the 2 Bothriurus species. The results herein recorded and the overview elaborated using the available cytogenetic information of Scorpiones elucidated current understanding regarding the processes of chromosome evolution that have occurred in Bothriuridae and in Scorpiones as a whole.

Parametric analysis of discrepant sets of data

In this paper a new parametric method to deal with discrepant experimental results is developed. The method is based on the fit of a probability density function to the data. This paper also compares the characteristics of different methods used to deduce recommended values and uncertainties from a discrepant set of experimental data. The methods are applied to the (137)Cs and (90)Sr published half-lives and special emphasis is given to the deduced confidence intervals. The obtained results are analyzed considering two fundamental properties expected from an experimental result: the probability content of confidence intervals and the statistical consistency between different recommended values. The recommended values and uncertainties for the (137)Cs and (90)Sr half-lives are 10,984 (24) days and 10,523 (70) days, respectively. (C) 2009 Elsevier B.V. All rights reserved.

Multiple Pathways Analysis of Brain Functional Networks from EEG Signals: An Application to Real Data

In the present study, we propose a theoretical graph procedure to investigate multiple pathways in brain functional networks. By taking into account all the possible paths consisting of h links between the nodes pairs of the network, we measured the global network redundancy R (h) as the number of parallel paths and the global network permeability P (h) as the probability to get connected. We used this procedure to investigate the structural and dynamical changes in the cortical networks estimated from a dataset of high-resolution EEG signals in a group of spinal cord injured (SCI) patients during the attempt of foot movement. In the light of a statistical contrast with a healthy population, the permeability index P (h) of the SCI networks increased significantly (P < 0.01) in the Theta frequency band (3-6 Hz) for distances h ranging from 2 to 4. On the contrary, no significant differences were found between the two populations for the redundancy index R (h) . The most significant changes in the brain functional network of SCI patients occurred mainly in the lower spectral contents. These changes were related to an improved propagation of communication between the closest cortical areas rather than to a different level of redundancy. This evidence strengthens the hypothesis of the need for a higher functional interaction among the closest ROIs as a mechanism to compensate the lack of feedback from the peripheral nerves to the sensomotor areas.

Morphological analysis of glass, carbon and glass/carbon fiber posts and bonding to self or dual-cured resin luting agents

OBJECTIVE: The aim of this study was to evaluate the morphology of glass (GF), carbon (CF) and glass/carbon (G/CF) fiber posts and their bond strength to self or dual-cured resin luting agents. MATERIAL AND METHODS: Morphological analysis of each post type was conducted under scanning electron microscopy (SEM). Bond strength was evaluated by microtensile test after bisecting the posts and re-bonding the two halves with the luting agents. Data were subjected to two-way ANOVA and Tukey's test (α=0.05). Failure modes were evaluated under optical microscopy and SEM. RESULTS: GF presented wider fibers and higher amount of matrix than CF, and G/CF presented carbon fibers surrounded by glass fibers, and both involved by matrix. For CF and GF, the dual-cured material presented significantly higher (p<0.05) bond strength than the self-cured agent. For the dual agent, CF presented similar bond strength to GF (p>0.05), but higher than that of G/CF (p<0.05). For the self-cured agent, no significant differences (p>0.05) were detected, irrespective of the post type. For GF and G/CF, all failures were considered mixed, while a predominance of adhesive failures was detected for CF. CONCLUSION: The bonding between fiber posts and luting agents was affected by the type of fibers and polymerization mode of the cement. When no surface treatment of the post is performed, the bonding between glass fiber post and dual-cured agent seems to be more reliable.

Analysis of contamination of endodontic absorbent paper points

PURPOSE: The aim of this study was to assess the contamination status of endodontic absorbent paper points from sterilized or not sterilized commercial packs, as well as paper points exposed to the dental office environment. METHODS: Twenty absorbent paper points were evaluated for contamination status packed under different conditions: commercial/sterilized pack, commercial/non-sterilized pack, exposed to the clinical environment, and intentionally contaminated (positive control). Contamination was determined qualitatively and quantitatively by aerobiosis, capnophilic growth, and pour plate. The Petri dishes were analyzed with a colony counter, and the results were expressed as colony-forming units. The data were analyzed by Kruskal-Wallis test (α=0.05). RESULTS: No difference in colony-forming units was found among the groups of endodontic absorbent paper points. All groups were contaminated by fungi and bacteria. CONCLUSION: It can be concluded that the sterilization of absorbent endodontic paper points before clinical use should be recommended regardless of commercial presentation

In vitro analysis of thermocompaction time and gutta-percha type on quality of main canal and lateral canals filling

The aim of this study was to evaluate the quality of filling in main and lateral root canals performed with the McSpadden technique, regarding the time spent on the procedure and the type of gutta-percha employed. Fifty simulated root canals, made with six lateral canals placed two apiece in the cervical, middle and apical thirds of the root, were divided into 5 groups. Group A: McSpadden technique with conventional gutta-percha, performed with sufficient time for canal filling; Group B: McSpadden technique with conventional gutta-percha, performed in twice the mean time used in Group A; Group C: McSpadden technique with TP gutta-percha, performed with sufficient time for canal filling; Group D: McSpadden technique with TP gutta-percha, performed in twice the mean time used in Group C; Group E: lateral condensation technique. Images of the filled root canals were taken using a stereomicroscope and analyzed using the Leica QWIN Pro software for filling material flow, gutta-percha filling extension and sealer flow. Data were analyzed by analysis of variance (ANOVA) and Tukey test (p < 0.05). The best values of penetration in lateral canals in the middle third occurred in the groups where TP gutta-percha was used. However, in the apical third, group B showed the best values. Although a longer time of compactor use allows greater penetration of the filling material into the lateral canals, the presence of voids resulted in bad quality radiographic images, suggesting porosity. The best quality of filling material was observed in Group A (McSpadden technique with conventional Gutta-Percha, performed with sufficient time for root canal filling).