Assessing influence in survival data with a cure fraction and covariates

Diagnostic methods have been an important tool in regression analysis to detect anomalies, such as departures from error assumptions and the presence of outliers and influential observations with the fitted models. Assuming censored data, we considered a classical analysis and Bayesian analysis assuming no informative priors for the parameters of the model with a cure fraction. A Bayesian approach was considered by using Markov Chain Monte Carlo Methods with Metropolis-Hasting algorithms steps to obtain the posterior summaries of interest. Some influence methods, such as the local influence, total local influence of an individual, local influence on predictions and generalized leverage were derived, analyzed and discussed in survival data with a cure fraction and covariates. The relevance of the approach was illustrated with a real data set, where it is shown that, by removing the most influential observations, the decision about which model best fits the data is changed.

Data collapse, scaling functions, and analytical solutions of generalized growth models

We consider a nontrivial one-species population dynamics model with finite and infinite carrying capacities. Time-dependent intrinsic and extrinsic growth rates are considered in these models. Through the model per capita growth rate we obtain a heuristic general procedure to generate scaling functions to collapse data into a simple linear behavior even if an extrinsic growth rate is included. With this data collapse, all the models studied become independent from the parameters and initial condition. Analytical solutions are found when time-dependent coefficients are considered. These solutions allow us to perceive nontrivial transitions between species extinction and survival and to calculate the transition's critical exponents. Considering an extrinsic growth rate as a cancer treatment, we show that the relevant quantity depends not only on the intensity of the treatment, but also on when the cancerous cell growth is maximum.

Indications of Conical Emission of Charged Hadrons at the BNL Relativistic Heavy Ion Collider

Three-particle azimuthal correlation measurements with a high transverse momentum trigger particle are reported for pp, d + Au, and Au + Au collisions at root(S)NN = 200 GeV by the STAR experiment. Dijet structures are observed in pp, d + Au and peripheral Au + Au collisions. An additional structure is observed in central Au + Au data, signaling conical emission of correlated charged hadrons. The conical emission angle is found to be theta = 1.37 +/- 0.02(stat)(-0.07)(+0.06)(syst), independent of p perpendicular to.

Urban air pollution and chronic obstructive pulmonary disease-related emergency department visits

Background: Patients with chronic obstructive pulmonary disease (COPD) can have recurrent disease exacerbations triggered by several factors, including air pollution. Visits to the emergency respiratory department can be a direct result of short-term exposure to air pollution. The aim of this study was to investigate the relationship between the daily number of COPD emergency department visits and the daily environmental air concentrations of PM(10), SO(2), NO(2), CO and O(3) in the City of Sao Paulo, Brazil. Methods: The sample data were collected between 2001 and 2003 and are categorised by gender and age. Generalised linear Poisson regression models were adopted to control for both short-and long-term seasonal changes as well as for temperature and relative humidity. The non-linear dependencies were controlled using a natural cubic spline function. Third-degree polynomial distributed lag models were adopted to estimate both lag structures and the cumulative effects of air pollutants. Results: PM(10) and SO(2) readings showed both acute and lagged effects on COPD emergency department visits. Interquartile range increases in their concentration (28.3 mg/m(3) and 7.8 mg/m(3), respectively) were associated with a cumulative 6-day increase of 19% and 16% in COPD admissions, respectively. An effect on women was observed at lag 0, and among the elderly the lag period was noted to be longer. Increases in CO concentration showed impacts in the female and elderly groups. NO(2) and O(3) presented mild effects on the elderly and in women, respectively. Conclusion: These results indicate that air pollution affects health in a gender-and age-specific manner and should be considered a relevant risk factor that exacerbates COPD in urban environments.

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.

Work disability remains a major problem in rheumatoid arthritis in the 2000s: data from 32 countries in the QUEST-RA Study

Introduction: Work disability is a major consequence of rheumatoid arthritis (RA), associated not only with traditional disease activity variables, but also more significantly with demographic, functional, occupational, and societal variables. Recent reports suggest that the use of biologic agents offers potential for reduced work disability rates, but the conclusions are based on surrogate disease activity measures derived from studies primarily from Western countries. Methods: The Quantitative Standard Monitoring of Patients with RA (QUEST-RA) multinational database of 8,039 patients in 86 sites in 32 countries, 16 with high gross domestic product (GDP) (>24K US dollars (USD) per capita) and 16 low-GDP countries (<11K USD), was analyzed for work and disability status at onset and over the course of RA and clinical status of patients who continued working or had stopped working in high-GDP versus low-GDP countries according to all RA Core Data Set measures. Associations of work disability status with RA Core Data Set variables and indices were analyzed using descriptive statistics and regression analyses. Results: At the time of first symptoms, 86% of men (range 57%-100% among countries) and 64% (19%-87%) of women <65 years were working. More than one third (37%) of these patients reported subsequent work disability because of RA. Among 1,756 patients whose symptoms had begun during the 2000s, the probabilities of continuing to work were 80% (95% confidence interval (CI) 78%-82%) at 2 years and 68% (95% CI 65%-71%) at 5 years, with similar patterns in high-GDP and low-GDP countries. Patients who continued working versus stopped working had significantly better clinical status for all clinical status measures and patient self-report scores, with similar patterns in high-GDP and low-GDP countries. However, patients who had stopped working in high-GDP countries had better clinical status than patients who continued working in low-GDP countries. The most significant identifier of work disability in all subgroups was Health Assessment Questionnaire (HAQ) functional disability score. Conclusions: Work disability rates remain high among people with RA during this millennium. In low-GDP countries, people remain working with high levels of disability and disease activity. Cultural and economic differences between societies affect work disability as an outcome measure for RA.

Integration and visualization of systems biology data in context of the genome

Background: High-density tiling arrays and new sequencing technologies are generating rapidly increasing volumes of transcriptome and protein-DNA interaction data. Visualization and exploration of this data is critical to understanding the regulatory logic encoded in the genome by which the cell dynamically affects its physiology and interacts with its environment. Results: The Gaggle Genome Browser is a cross-platform desktop program for interactively visualizing high-throughput data in the context of the genome. Important features include dynamic panning and zooming, keyword search and open interoperability through the Gaggle framework. Users may bookmark locations on the genome with descriptive annotations and share these bookmarks with other users. The program handles large sets of user-generated data using an in-process database and leverages the facilities of SQL and the R environment for importing and manipulating data. A key aspect of the Gaggle Genome Browser is interoperability. By connecting to the Gaggle framework, the genome browser joins a suite of interconnected bioinformatics tools for analysis and visualization with connectivity to major public repositories of sequences, interactions and pathways. To this flexible environment for exploring and combining data, the Gaggle Genome Browser adds the ability to visualize diverse types of data in relation to its coordinates on the genome. Conclusions: Genomic coordinates function as a common key by which disparate biological data types can be related to one another. In the Gaggle Genome Browser, heterogeneous data are joined by their location on the genome to create information-rich visualizations yielding insight into genome organization, transcription and its regulation and, ultimately, a better understanding of the mechanisms that enable the cell to dynamically respond to its environment.

New Epidemiological Data on Brazilian Spotted Fever in an Endemic Area of the State of Sao Paulo, Brazil

The present work evaluated rickettsial infection in dogs and their ticks in an area endemic for Brazilian spotted fever (BSF) in the metropolitan area of Sao Paulo, Brazil, where the tick Amblyomma aureolatum was presumed to be the vector of the disease. Ticks were collected on dogs from 185 houses, encompassing single infestations by Rhipicephalus sanguineus, Amblyomma aureolatum, Amblyomma longirostre, or Amblyomma sp. in dogs from 60 (32.4%), 77 (41.6%), 2 (1.1%), and 25 (13.5%) houses, respectively; 19 (10.3%) houses had dogs with mixed infestations by R. sanguineus and A. aureolatum; 1 (0.5%) house had dogs with infestations by A. aureolatum and A. longirostre; and 1 (0.5%) house had dogs with infestations by R. sanguineus and Amblyomma sp. Overall, A. aureolatum was present in dogs from 97 (52.4%) houses, and R. sanguineus in dogs from 80 (43.2%) houses. A total of 287 ticks (130 A. aureolatum and 157 R. sanguineus) infesting dogs from 98 houses were selected for testing by polymerase chain reaction (PCR) targeting rickettsial genes. Overall, 3.1% of the A. aureolatum ticks were infected by Rickettsia bellii, and 1.3% of the R. sanguineus were infected by Ricketttsii rickettsii. For serology, we selected 23 dogs living in and in the vicinity of the house where the R. rickettsii-infected ticks were collected. The indirect fluorescent antibody (IFA) test detected antibodies reactive with R. rickettsii in sera from 16 (69.6%) dogs, with titers ranging from 256 to 32,768. It is established that Amblyomma aureolatum is a vector of R. rickettsii in the Sao Paulo metropolitan area, but our results highlight for the first time in Brazil, a possible role of R. sanguineus in the epidemiology of R. rickettsii, corroborating previous findings in Mexico and the United States, where R. sanguineus has been implicated in the transmission of R. rickettsii to humans.

DESCRIPTION OF ADULTS AND NYMPH, AND REDESCRIPTION OF THE LARVA, OF ORNITHODOROS MARINKELLEI (ACARI: ARGASIDAE), WITH DATA ON ITS PHYLOGENETIC POSITION

The argasid tick Ornithodoros marinkellei Kohls, Clifford, and Jones, 1969 was described 4 decades ago based on larval specimens collected from bats (Pteronotus spp.) in Colombia and Panama. Thereafter, larval O. marinkellei parasitizing bats were reported from Venezuela, Guyana, and Brazil. Herein, we describe the adults and nymph, and redescribe the larva of O. marinkellei based on specimens recently collected in the western Brazilian Amazon region. In contrast to all other known adult argasids, the idiosoma of both males and females of O. marinkellei is covered with sclerotized plaques. The idiosoma of the nymph of O. marinkellei is entirely micromamillated, and differs from the adults by the absence of plaques. The larva of O. marinkellei is morphologically similar to the larvae of the 2 other species belonging to the subgenus Subparmatus, i.e., Ornithodoros viguerasi Cooley and Kohls, 1941 and Ornithodoros mormoops Kohls, Clifford, and Jones, 1969. Because of the long and narrow dorsal plate, the larva of O. marinkellei is readily distinguished from O. viguerasi and O. mormoops. Comparison of our larvae from Brazil with O. marinkellei paratype specimens from Colombia confirmed their taxonomic identification. However, a few morphological differences, particularly in the size of the gnathosoma, were observed. Further studies are necessary to clarify whether O. marinkellei is a complex of different species, or a single species represented by morphologically polymorphic, and geographically distinct populations. Partial mitochondrial 16S rDNA gene sequences were generated for O. marinkellei specimens from Brazil, and compared with available homologous sequences in GenBank. Phylogenetic analyses revealed O. marinkellei to be distinct from the remaining argasid species available in GenBank, including other bat-associated tick species that are found in sympatry with O. marinkellei in the Neotropical region.

SPITZER OBSERVATIONS OF A1763. II. CONSTRAINING THE NATURE OF ACTIVITY IN THE CLUSTER-FEEDING FILAMENT WITH VLA AND XMM-NEWTON DATA

The A1763 superstructure at z = 0.23 contains the first galaxy filament to be directly detected using mid-infrared observations. Our previous work has shown that the frequency of starbursting galaxies, as characterized by 24 mu m emission is much higher within the filament than at either the center of the rich galaxy cluster, or the field surrounding the system. New Very Large Array and XMM-Newton data are presented here. We use the radio and X-ray data to examine the fraction and location of active galaxies, both active galactic nuclei (AGNs) and starbursts (SBs). The radio far-infrared correlation, X-ray point source location, IRAC colors, and quasar positions are all used to gain an understanding of the presence of dominant AGNs. We find very few MIPS-selected galaxies that are clearly dominated by AGN activity. Most radio-selected members within the filament are SBs. Within the supercluster, three of eight spectroscopic members detected both in the radio and in the mid-infrared are radio-bright AGNs. They are found at or near the core of A1763. The five SBs are located further along the filament. We calculate the physical properties of the known wide angle tail (WAT) source which is the brightest cluster galaxy of A1763. A second double lobe source is found along the filament well outside of the virial radius of either cluster. The velocity offset of the WAT from the X-ray centroid and the bend of the WAT in the intracluster medium are both consistent with ram pressure stripping, indicative of streaming motions along the direction of the filament. We consider this as further evidence of the cluster-feeding nature of the galaxy filament.

Spiral-like structure at the centre of nearby clusters of galaxies

Context. X-ray data analysis have found that fairly complex structures at cluster centres are more common than expected. Many of these structures have similar morphologies, which exhibit spiral-like substructure. Aims. It is not yet well known how these structures are formed or maintained. Understanding the origin of these spiral-like features at the centre of some clusters is the major motivation behind this work. Methods. We analyse deep Chandra observations of 15 nearby galaxy clusters ( 0.01 < z < 0.06), and use X-ray temperature and substructure maps to detect small features at the cores of the clusters. Results. We detect spiral-like features at the centre of 7 clusters: A85, A426, A496, Hydra A cluster, Centaurus, Ophiuchus, and A4059. These patterns are similar to those found in numerical hydrodynamic simulations of cluster mergers with non-zero impact parameter. In some clusters of our sample, a strong radio source also occupies the inner region of the cluster, which indicates a possible connection between the two. Our investigation implies that these spiral-like structures may be caused by off-axis minor mergers. Since these features occur in regions of high density, they may confine radio emission from the central galaxy producing, in some cases, unusual radio morphology.

REVISITING THE FOSSIL GROUP CANDIDATES UGC 842 AND NGC 6034

We present a new insight on NGC 6034 and UGC 842, two groups of galaxies previously reported in the literature as being fossil groups. The study is based on optical photometry and spectroscopy obtained with the CTIO Blanco telescope and Sloan Digital Sky Survey archival data. We find that NGC 6034 is embedded in a large structure, dominated by three rich clusters and other small groups. Its first and next four ranked galaxies have magnitude differences in the r band and projected distances which violate the optical criteria to classify it as a fossil group. We confirm that the UGC 842 group is a fossil group, but with about half the velocity dispersion that is reported in previous works. The velocity distribution of its galaxies reveals the existence of two structures in its line of sight, one with sigma(nu) similar to 223 km s(-1) and another with sigma(nu) similar to 235 km s(-1), with a difference in velocity of similar to 820 km s(-1). The main structure is dominated by passive galaxies, while these represent similar to 60% of the second structure. The X-ray temperature for the intragroup medium of a group with such a velocity dispersion is expected to be kT similar to 0.5-1 keV, against the observed value of kT similar to 1.9 keV reported in the literature. This result makes UGC 842 a special case among fossil groups because (1) it represents more likely the interaction between two small groups, which warms the intragroup medium and/or (2) it could constitute evidence that member galaxies lost energy in the process of spiraling toward the group center, and decreased the velocity dispersion of the system. As far as we know, UGC 842 is the first low-mass fossil group studied in detail.

The B0.5IVe CoRoT target HD 49330 I. Photometric analysis from CoRoT data

Context. Be stars undergo outbursts producing a circumstellar disk from the ejected material. The beating of non-radial pulsations has been put forward as a possible mechanism of ejection. Aims. We analyze the pulsational behavior of the early B0.5IVe star HD 49330 observed during the first CoRoT long run towards the Galactical anticenter (LRA1). This Be star is located close to the lower edge of the beta Cephei instability strip in the HR diagram and showed a 0.03 mag outburst during the CoRoT observations. It is thus an ideal case for testing the aforementioned hypothesis. Methods. We analyze the CoRoT light curve of HD 49330 using Fourier methods and non-linear least square fitting. Results. In this star, we find pulsation modes typical of beta Cep stars (p modes) and SPB stars (g modes) with amplitude variations along the run directly correlated with the outburst. These results provide new clues about the origin of the Be phenomenon as well as strong constraints on the seismic modelling of Be stars.

The BRAzilian Seismographic Integrated Systems (BRASIS): infrastructure and data management

In geophysics and seismology, raw data need to be processed to generate useful information that can be turned into knowledge by researchers. The number of sensors that are acquiring raw data is increasing rapidly. Without good data management systems, more time can be spent in querying and preparing datasets for analyses than in acquiring raw data. Also, a lot of good quality data acquired at great effort can be lost forever if they are not correctly stored. Local and international cooperation will probably be reduced, and a lot of data will never become scientific knowledge. For this reason, the Seismological Laboratory of the Institute of Astronomy, Geophysics and Atmospheric Sciences at the University of Sao Paulo (IAG-USP) has concentrated fully on its data management system. This report describes the efforts of the IAG-USP to set up a seismology data management system to facilitate local and international cooperation.

Kinematic constraints to the transition redshift from supernovae type Ia union data

The kinematic approach to cosmological tests provides direct evidence to the present accelerating stage of the Universe that does not depend on the validity of general relativity, as well as on the matter-energy content of the Universe. In this context, we consider here a linear two-parameter expansion for the decelerating parameter, q(z)=q(0)+q(1)z, where q(0) and q(1) are arbitrary constants to be constrained by the union supernovae data. By assuming a flat Universe we find that the best fit to the pair of free parameters is (q(0),q(1))=(-0.73,1.5) whereas the transition redshift is z(t)=0.49(-0.07)(+0.14)(1 sigma) +0.54-0.12(2 sigma). This kinematic result is in agreement with some independent analyses and more easily accommodates many dynamical flat models (like Lambda CDM).