Plotting regression coefficients and other estimates in Stata

Graphical presentation of regression results has become increasingly popular in the scientific literature, as graphs are much easier to read than tables in many cases. In Stata such plots can be produced by the -marginsplot- command. However, while -marginsplot- is very versatile and flexible, it has two major limitations: it can only process results left behind by -margins- and it can only handle one set of results at the time. In this article I introduce a new command called -coefplot- that overcomes these limitations. It plots results from any estimation command and combines results from several models into a single graph. The default behavior of -coefplot- is to plot markers for coefficients and horizontal spikes for confidence intervals. However, -coefplot- can also produce various other types of graphs. The capabilities of -coefplot- are illustrated in this article using a series of examples.

COEFPLOT: Stata module to plot regression coefficients and other results

coefplot plots results from estimation commands or Stata matrices. Results from multiple models or matrices can be combined in a single graph. The default behavior of coefplot is to draw markers for coefficients and horizontal spikes for confidence intervals. However, coefplot can also produce various other types of graphs.

Plotting regression coefficients and other estimates

Graphical display of regression results has become increasingly popular in presentations and in scientific literature because graphs are often much easier to read than tables. Such plots can be produced in Stata by the marginsplot command (see [R] marginsplot). However, while marginsplot is versatile and flexible, it has two major limitations: it can only process results left behind by margins (see [R] margins), and it can handle only one set of results at a time. In this article, I introduce a new command called coefplot that overcomes these limitations. It plots results from any estimation command and combines results from several models into one graph. The default behavior of coefplot is to plot markers for coefficients and horizontal spikes for confidence intervals. However, coefplot can also produce other types of graphs. I illustrate the capabilities of coefplot by using a series of examples.

Ebola virus disease outbreak in Nigeria: Transmission dynamics and rapid control.

International air travel has already spread Ebola virus disease (EVD) to major cities as part of the unprecedented epidemic that started in Guinea in December 2013. An infected airline passenger arrived in Nigeria on July 20, 2014 and caused an outbreak in Lagos and then Port Harcourt. After a total of 20 reported cases, including 8 deaths, Nigeria was declared EVD free on October 20, 2014. We quantified the impact of early control measures in preventing further spread of EVD in Nigeria and calculated the risk that a single undetected case will cause a new outbreak. We fitted an EVD transmission model to data from the outbreak in Nigeria and estimated the reproduction number of the index case at 9.0 (95% confidence interval [CI]: 5.2-15.6). We also found that the net reproduction number fell below unity 15 days (95% CI: 11-21 days) after the arrival of the index case. Hence, our study illustrates the time window for successful containment of EVD outbreaks caused by infected air travelers.

Combining spatial-temporal and phylogenetic analysis approaches for improved understanding on global H5N1 transmission

Background Since late 2003, the highly pathogenic influenza A H5N1 had initiated several outbreak waves that swept across the Eurasia and Africa continents. Getting prepared for reassortment or mutation of H5N1 viruses has become a global priority. Although the spreading mechanism of H5N1 has been studied from different perspectives, its main transmission agents and spread route problems remain unsolved. Methodology/Principal Findings Based on a compilation of the time and location of global H5N1 outbreaks from November 2003 to December 2006, we report an interdisciplinary effort that combines the geospatial informatics approach with a bioinformatics approach to form an improved understanding on the transmission mechanisms of H5N1 virus. Through a spherical coordinate based analysis, which is not conventionally done in geographical analyses, we reveal obvious spatial and temporal clusters of global H5N1 cases on different scales, which we consider to be associated with two different transmission modes of H5N1 viruses. Then through an interdisciplinary study of both geographic and phylogenetic analysis, we obtain a H5N1 spreading route map. Our results provide insight on competing hypotheses as to which avian hosts are responsible for the spread of H5N1. Conclusions/Significance We found that although South China and Southeast Asia may be the virus pool of avian flu, East Siberia may be the source of the H5N1 epidemic. The concentration of migratory birds from different places increases the possibility of gene mutation. Special attention should be paid to East Siberia, Middle Siberia and South China for improved surveillance of H5N1 viruses and monitoring of migratory birds.

Interaction between effects of genes coding for dopamine and glutamate transmission on striatal and parahippocampal function

The genes for the dopamine transporter (DAT) and the D-Amino acid oxidase activator (DAOA or G72) have been independently implicated in the risk for schizophrenia and in bipolar disorder and/or their related intermediate phenotypes. DAT and G72 respectively modulate central dopamine and glutamate transmission, the two systems most robustly implicated in these disorders. Contemporary studies have demonstrated that elevated dopamine function is associated with glutamatergic dysfunction in psychotic disorders. Using functional magnetic resonance imaging we examined whether there was an interaction between the effects of genes that influence dopamine and glutamate transmission (DAT and G72) on regional brain activation during verbal fluency, which is known to be abnormal in psychosis, in 80 healthy volunteers. Significant interactions between the effects of G72 and DAT polymorphisms on activation were evident in the striatum, parahippocampal gyrus, and supramarginal/angular gyri bilaterally, the right insula, in the right pre-/postcentral and the left posterior cingulate/retrosplenial gyri (P < 0.05, FDR-corrected across the whole brain). This provides evidence that interactions between the dopamine and the glutamate system, thought to be altered in psychosis, have an impact in executive processing which can be modulated by common genetic variation.

Impact of viral load and the duration of primary infection on HIV transmission

Abstract Objectives: HIV 'treatment as prevention' (TasP) describes early treatment of HIV-infected patients intended to reduce viral load and transmission. Crucial assumptions for estimating TasP's effectiveness are the underlying estimates of transmission risk. We aimed to determine transmission risk during primary infection, and of the relation of HIV transmission risk to viral load. Design: A systematic review and meta-analysis. Methods: We searched PubMed and Embase databases for studies that established a relationship between viral load and transmission risk, or primary infection and transmission risk, in serodiscordant couples. We analysed assumptions about the relationship between viral load and transmission risk, and between duration of primary infection and transmission risk. Results: We found 36 eligible articles, based on six different study populations. Studies consistently found that larger viral loads lead to higher HIV transmission rates, but assumptions about the shape of this increase varied from exponential increase to saturation. The assumed duration of primary infection ranged from 1.5 to 12 months; for each additional month, the log10 transmission rate ratio between primary and asymptomatic infection decreased by 0.40. Conclusion: Assumptions and estimates of the relationship between viral load and transmission risk, and the relationship between primary infection and transmission risk, vary substantially and predictions of TasP's effectiveness should take this uncertainty into account.

The Transmission of Manuscripts and the Genetics of Texts. A new critical edition of Wolfram von Eschenbach’s Parzival

The discussion on the New Philology triggered by French and North American scholars in the last decade of the 20th century emphasized the material character of textual transmission inside and outside the written evidences of medieval manuscripts by downgrading the active role of the historical author. However, the reception of the ideas propagated by the New Philology adherents was rather divided. Some researchers questioned its innovative status (K. Stackmann: “Neue Philologie?”), others saw a new era of the “powers of philology” evoked (H.-U. Gumbrecht). Besides the debates on the New Philology another concept of textual materiality strengthened in the last decade, maintaining that textual alterations somewhat relate to biogenetic mutations. In a matter of fact, phenomena such as genetic and textual variation, gene recombination and ‘contamination’ (the mixing of different exemplars in one manuscript text) share common features. The paper discusses to what extent the biogenetic concepts can be used for evaluating manifestations of textual production (as the approach of ‘critique génétique’ does) and of textual transmission (as the phylogenetic analysis of manuscript variation does). In this context yet the genealogical concept of stemmatology – the treelike representation of textual development abhorred by the New Philology adepts – might prove to be useful for describing the history of texts. The textual material to be analyzed will be drawn from the Parzival Project, which is currently preparing a new electronic edition of Wolfram von Eschenbach’s Parzival novel written shortly after 1200 and transmitted in numerous manuscripts up to the age of printing (www.parzival.unibe.ch). Researches of the project have actually resulted in suggesting that the advanced knowledge of the manuscript transmission yields a more precise idea on the author’s own writing process.

Conservation of threatened habitat types under future climate change – Lessons from plant-distribution models and current extinction trends in southern Germany

A higher risk of future range losses as a result of climate change is expected to be one of the main drivers of extinction trends in vascular plants occurring in habitat types of high conservation value. Nevertheless, the impact of the climate changes of the last 60 years on the current distribution and extinction patterns of plants is still largely unclear. We applied species distribution models to study the impact of environmental variables (climate, soil conditions, land cover, topography), on the current distribution of 18 vascular plant species characteristic of three threatened habitat types in southern Germany: (i) xero-thermophilous vegetation, (ii) mesophilous mountain grasslands (mountain hay meadows and matgrass communities), and (iii) wetland habitats (bogs, fens, and wet meadows). Climate and soil variables were the most important variables affecting plant distributions at a spatial level of 10 × 10 km. Extinction trends in our study area revealed that plant species which occur in wetland habitats faced higher extinction risks than those in xero-thermophilous vegetation, with the risk for species in mesophilous mountain grasslands being intermediary. For three plant species characteristic either of mesophilous mountain grasslands or wetland habitats we showed exemplarily that extinctions from 1950 to the present day have occurred at the edge of the species’ current climatic niche, indicating that climate change has likely been the main driver of extinction. This is largely consistent with current extinction trends reported in other studies. Our study indicates that the analysis of past extinctions is an appropriate means to assess the impact of climate change on species and that vulnerability to climate change is both species- and habitat-specific.

SPECTRAL DOMAIN-OPTICAL COHERENCE TOMOGRAPHY IMAGE CONTRAST AND BACKGROUND COLOR SETTINGS INFLUENCE IDENTIFICATION OF RETINAL STRUCTURES.

PURPOSE To evaluate image contrast and color setting on assessment of retinal structures and morphology in spectral-domain optical coherence tomography. METHODS Two hundred and forty-eight Spectralis spectral-domain optical coherence tomography B-scans of 62 patients were analyzed by 4 readers. B-scans were extracted in 4 settings: W + N = white background with black image at normal contrast 9; W + H = white background with black image at maximum contrast 16; B + N = black background with white image at normal contrast 12; B + H = black background with white image at maximum contrast 16. Readers analyzed the images to identify morphologic features. Interreader correlation was calculated. Differences between Fleiss-kappa correlation coefficients were examined using bootstrap method. Any setting with significantly higher correlation coefficient was deemed superior for evaluating specific features. RESULTS Correlation coefficients differed among settings. No single setting was superior for all respective spectral-domain optical coherence tomography parameters (P = 0.3773). Some variables showed no differences among settings. Hard exudates and subretinal fluid were best seen with B + H (κ = 0.46, P = 0.0237 and κ = 0.78, P = 0.002). Microaneurysms were best seen with W + N (κ = 0.56, P = 0.025). Vitreomacular interface, enhanced transmission signal, and epiretinal membrane were best identified using all color/contrast settings together (κ = 0.44, P = 0.042, κ = 0.57, P = 0.01, and κ = 0.62, P ≤ 0.0001). CONCLUSION Contrast and background affect the evaluation of retinal structures on spectral-domain optical coherence tomography images. No single setting was superior for all features, though certain changes were best seen with specific settings.

Sexual Activity Without Condoms and Risk of HIV Transmission in Serodifferent Couples When the HIV-Positive Partner Is Using Suppressive Antiretroviral Therapy

Importance A key factor in assessing the effectiveness and cost-effectiveness of antiretroviral therapy (ART) as a prevention strategy is the absolute risk of HIV transmission through condomless sex with suppressed HIV-1 RNA viral load for both anal and vaginal sex. Objective To evaluate the rate of within-couple HIV transmission (heterosexual and men who have sex with men [MSM]) during periods of sex without condoms and when the HIV-positive partner had HIV-1 RNA load less than 200 copies/mL. Design, Setting, and Participants The prospective, observational PARTNER (Partners of People on ART-A New Evaluation of the Risks) study was conducted at 75 clinical sites in 14 European countries and enrolled 1166 HIV serodifferent couples (HIV-positive partner taking suppressive ART) who reported condomless sex (September 2010 to May 2014). Eligibility criteria for inclusion of couple-years of follow-up were condomless sex and HIV-1 RNA load less than 200 copies/mL. Anonymized phylogenetic analysis compared couples' HIV-1 polymerase and envelope sequences if an HIV-negative partner became infected to determine phylogenetically linked transmissions. Exposures Condomless sexual activity with an HIV-positive partner taking virally suppressive ART. Main Outcomes and Measures Risk of within-couple HIV transmission to the HIV-negative partner. Results Among 1166 enrolled couples, 888 (mean age, 42 years [IQR, 35-48]; 548 heterosexual [61.7%] and 340 MSM [38.3%]) provided 1238 eligible couple-years of follow-up (median follow-up, 1.3 years [IQR, 0.8-2.0]). At baseline, couples reported condomless sex for a median of 2 years (IQR, 0.5-6.3). Condomless sex with other partners was reported by 108 HIV-negative MSM (33%) and 21 heterosexuals (4%). During follow-up, couples reported condomless sex a median of 37 times per year (IQR, 15-71), with MSM couples reporting approximately 22 000 condomless sex acts and heterosexuals approximately 36 000. Although 11 HIV-negative partners became HIV-positive (10 MSM; 1 heterosexual; 8 reported condomless sex with other partners), no phylogenetically linked transmissions occurred over eligible couple-years of follow-up, giving a rate of within-couple HIV transmission of zero, with an upper 95% confidence limit of 0.30/100 couple-years of follow-up. The upper 95% confidence limit for condomless anal sex was 0.71 per 100 couple-years of follow-up. Conclusions and Relevance Among serodifferent heterosexual and MSM couples in which the HIV-positive partner was using suppressive ART and who reported condomless sex, during median follow-up of 1.3 years per couple, there were no documented cases of within-couple HIV transmission (upper 95% confidence limit, 0.30/100 couple-years of follow-up). Additional longer-term follow-up is necessary to provide more precise estimates of risk.

Audiological results with Baha in conductive and mixed hearing loss

The level of improvement in the audiological results of Baha(®) users mainly depends on the patient's preoperative hearing thresholds and the type of Baha sound processor used. This investigation shows correlations between the preoperative hearing threshold and postoperative aided thresholds and audiological results in speech understanding in quiet of 84 Baha users with unilateral conductive hearing loss, bilateral conductive hearing loss and bilateral mixed hearing loss. Secondly, speech understanding in noise of 26 Baha users with different Baha sound processors (Compact, Divino, and BP100) is investigated. Linear regression between aided sound field thresholds and bone conduction (BC) thresholds of the better ear shows highest correlation coefficients and the steepest slope. Differences between better BC thresholds and aided sound field thresholds are smallest for mid-frequencies (1 and 2 kHz) and become larger at 0.5 and 4 kHz. For Baha users, the gain in speech recognition in quiet can be expected to lie in the order of magnitude of the gain in their hearing threshold. Compared to its predecessor sound processors Baha(®) Compact and Baha(®) Divino, Baha(®) BP100 improves speech understanding in noise significantly by +0.9 to +4.6 dB signal-to-noise ratio, depending on the setting and the use of directional microphone. For Baha users with unilateral and bilateral conductive hearing loss and bilateral mixed hearing loss, audiological results in aided sound field thresholds can be estimated with the better BC hearing threshold. The benefit in speech understanding in quiet can be expected to be similar to the gain in their sound field hearing threshold. The most recent technology of Baha sound processor improves speech understanding in noise by an order of magnitude that is well perceived by users and which can be very useful in everyday life.