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.

Delivery of health care at the end of life in cancer patients of four swiss cantons a retrospective database study (SAKK 89/09)

Background The use of cancer related therapy in cancer patients at the end-of-life has increased over time in many countries. Given a lack of published Swiss data, the objective of this study was to describe delivery of health care during the last month before death of cancer patients. Methods Claims data were used to assess health care utilization of cancer patients (identified by cancer registry data of four participating cantons), deceased between 2006-2008. Primary endpoints were hospitalization rate and delivery of cancer related therapies during the last 30 days before death. Multivariate logistic regression assessed the explanatory value of patient and geographic characteristics. Results 3809 identified cancer patients were included. Hospitalization rate (mean 68.5%, 95%CI 67.0-69.9) and percentage of patients receiving anti-cancer drug therapies (ACDT, mean 14.5%, 95%CI 13.4-15.6) and radiotherapy (mean 7.7%, 95%CI 6.7-8.4) decreased with age. Canton of residence and insurance type status most significantly influenced the odds for hospitalization or receiving ACDT. Conclusions The intensity of cancer specific care showed substantial variation by age, cancer type, place of residence and insurance type status. This may be partially driven by cultural differences within Switzerland and the cantonal organization of the Swiss health care system.

Determinants of Acidobacteria activity inferred from the relative abundances of 16S rRNA transcripts in German grassland and forest soils

16S rRNA genes and transcripts of Acidobacteria were investigated in 57 grassland and forest soils of three different geographic regions. Acidobacteria contributed 9-31% of bacterial 16S rRNA genes whereas the relative abundances of the respective transcripts were 4-16%. The specific cellular 16S rRNA content (determined as molar ratio of rRNA:rRNA genes) ranged between 3 and 80, indicating a low in situ growth rate. Correlations with flagellate numbers, vascular plant diversity and soil respiration suggest that biotic interactions are important determinants of Acidobacteria 16S rRNA transcript abundances in soils. While the phylogenetic composition of Acidobacteria differed significantly between grassland and forest soils, high throughput denaturing gradient gel electrophoresis and terminal restriction fragment length polymorphism fingerprinting detected 16S rRNA transcripts of most phylotypes in situ. Partial least squares regression suggested that chemical soil conditions such as pH, total nitrogen, C:N ratio, ammonia concentrations and total phosphorus affect the composition of this active fraction of Acidobacteria. Transcript abundance for individual Acidobacteria phylotypes was found to correlate with particular physicochemical (pH, temperature, nitrogen or phosphorus) and, most notably, biological parameters (respiration rates, abundances of ciliates or amoebae, vascular plant diversity), providing culture-independent evidence for a distinct niche specialization of different Acidobacteria even from the same subdivision.

Regression of the inflammatory microenvironment of the peritoneal cavity in women with endometriosis by GnRHa treatment.

OBJECTIVE To investigate the effect of gonadotropin-releasing hormone analogues (GnRHa) on the peritoneal fluid microenvironment in women with endometriosis. STUDY DESIGN Peritoneal fluid was collected from 85 women with severe endometriosis (rAFS stage III and IV) during laparoscopic surgery during the proliferative phase. Prior to surgery clinical data were collected. The concentrations of specific markers for endometriosis in the peritoneal fluid were determined using an ELISA and a comparison between peritoneal fluid markers in women using GnRHa and no hormonal treatment was performed using a non-parametric Mann-Whitney U test. RESULTS The study included peritoneal fluid from 39 patients who had been administered GnRHa (Zoladex(®)) in the three months prior to surgery and 46 from women with no hormonal treatment in this period. Concentrations of IL-8, PAPP-A, glycodelin-A and midkine were significantly reduced in the GnRHa treatment group compared to women receiving no hormonal treatment. RANTES, MCP-1, ENA-78, TNF-α, OPG, IP-10 and defensin showed no significant change between the two groups. CONCLUSIONS GnRHa mediate a significant regression in the inflammatory nature of the peritoneal microenvironment in women with endometriosis.

Independent measurement of biogenic silica in sediments by FTIR spectroscopy and PLS regression

We present an independent calibration model for the determination of biogenic silica (BSi) in sediments, developed from analysis of synthetic sediment mixtures and application of Fourier transform infrared spectroscopy (FTIRS) and partial least squares regression (PLSR) modeling. In contrast to current FTIRS applications for quantifying BSi, this new calibration is independent from conventional wet-chemical techniques and their associated measurement uncertainties. This approach also removes the need for developing internal calibrations between the two methods for individual sediments records. For the independent calibration, we produced six series of different synthetic sediment mixtures using two purified diatom extracts, with one extract mixed with quartz sand, calcite, 60/40 quartz/calcite and two different natural sediments, and a second extract mixed with one of the natural sediments. A total of 306 samples—51 samples per series—yielded BSi contents ranging from 0 to 100 %. The resulting PLSR calibration model between the FTIR spectral information and the defined BSi concentration of the synthetic sediment mixtures exhibits a strong cross-validated correlation ( R2cv = 0.97) and a low root-mean square error of cross-validation (RMSECV = 4.7 %). Application of the independent calibration to natural lacustrine and marine sediments yields robust BSi reconstructions. At present, the synthetic mixtures do not include the variation in organic matter that occurs in natural samples, which may explain the somewhat lower prediction accuracy of the calibration model for organic-rich samples.

Discontinuity in the responses of ecosystem processes and multifunctionality to altered soil community composition

Ecosystem management policies increasingly emphasize provision of multiple, as opposed to single, ecosystem services. Management for such "multifunctionality" has stimulated research into the role that biodiversity plays in providing desired rates of multiple ecosystem processes. Positive effects of biodiversity on indices of multifunctionality are consistently found, primarily because species that are redundant for one ecosystem process under a given set of environmental conditions play a distinct role under different conditions or in the provision of another ecosystem process. Here we show that the positive effects of diversity (specifically community composition) on multifunctionality indices can also arise from a statistical fallacy analogous to Simpson's paradox (where aggregating data obscures causal relationships). We manipulated soil faunal community composition in combination with nitrogen fertilization of model grassland ecosystems and repeatedly measured five ecosystem processes related to plant productivity, carbon storage, and nutrient turnover. We calculated three common multifunctionality indices based on these processes and found that the functional complexity of the soil communities had a consistent positive effect on the indices. However, only two of the five ecosystem processes also responded positively to increasing complexity, whereas the other three responded neutrally or negatively. Furthermore, none of the individual processes responded to both the complexity and the nitrogen manipulations in a manner consistent with the indices. Our data show that multifunctionality indices can obscure relationships that exist between communities and key ecosystem processes, leading us to question their use in advancing theoretical understanding-and in management decisions-about how biodiversity is related to the provision of multiple ecosystem services.

Incidence and Imaging Outcomes of Acute Scaffold Disruption and Late Structural Discontinuity After Implantation of the Absorb Everolimus-Eluting Fully Bioresorbable Vascular Scaffold: Optical Coherence Tomography Assessment in the ABSORB Cohort B Trial (A Clinical Evaluation of the Bioabsorbable Everolimus Eluting Coronary Stent System in the Treatment of Patients With De Novo Native Coronary Artery Lesions).

OBJECTIVES This study sought to describe the frequency and clinical impact of acute scaffold disruption and late strut discontinuity of the second-generation Absorb bioresorbable polymeric vascular scaffolds (Absorb BVS, Abbott Vascular, Santa Clara, California) in the ABSORB (A Clinical Evaluation of the Bioabsorbable Everolimus Eluting Coronary Stent System in the Treatment of Patients With De Novo Native Coronary Artery Lesions) cohort B study by optical coherence tomography (OCT) post-procedure and at 6, 12, 24, and 36 months. BACKGROUND Fully bioresorbable scaffolds are a novel approach to treatment for coronary narrowing that provides transient vessel support with drug delivery capability without the long-term limitations of metallic drug-eluting stents. However, a potential drawback of the bioresorbable scaffold is the potential for disruption of the strut network when overexpanded. Conversely, the structural discontinuity of the polymeric struts at a late stage is a biologically programmed fate of the scaffold during the course of bioresorption. METHODS The ABSORB cohort B trial is a multicenter single-arm trial assessing the safety and performance of the Absorb BVS in the treatment of 101 patients with de novo native coronary artery lesions. The current analysis included 51 patients with 143 OCT pullbacks who underwent OCT at baseline and follow-up. The presence of acute disruption or late discontinuities was diagnosed by the presence on OCT of stacked, overhung struts or isolated intraluminal struts disconnected from the expected circularity of the device. RESULTS Of 51 patients with OCT imaging post-procedure, acute scaffold disruption was observed in 2 patients (3.9%), which could be related to overexpansion of the scaffold at the time of implantation. One patient had a target lesion revascularization that was presumably related to the disruption. Of 49 patients without acute disruption, late discontinuities were observed in 21 patients. There were no major adverse cardiac events associated with this finding except for 1 patient who had a non-ischemia-driven target lesion revascularization. CONCLUSIONS Acute scaffold disruption is a rare iatrogenic phenomenon that has been anecdotally associated with anginal symptoms, whereas late strut discontinuity is observed in approximately 40% of patients and could be viewed as a serendipitous OCT finding of a normal bioresorption process without clinical implications. (ABSORB Clinical Investigation, Cohort B [ABSORB B]; NCT00856856).

Assessment of tumor regression of esophageal adenocarcinomas after neoadjuvant chemotherapy: comparison of 2 commonly used scoring approaches.

Histopathologic determination of tumor regression provides important prognostic information for locally advanced gastroesophageal carcinomas after neoadjuvant treatment. Regression grading systems mostly refer to the amount of therapy-induced fibrosis in relation to residual tumor or the estimated percentage of residual tumor in relation to the former tumor site. Although these methods are generally accepted, currently there is no common standard for reporting tumor regression in gastroesophageal cancers. We compared the application of these 2 major principles for assessment of tumor regression: hematoxylin and eosin-stained slides from 89 resection specimens of esophageal adenocarcinomas following neoadjuvant chemotherapy were independently reviewed by 3 pathologists from different institutions. Tumor regression was determined by the 5-tiered Mandard system (fibrosis/tumor relation) and the 4-tiered Becker system (residual tumor in %). Interobserver agreement for the Becker system showed better weighted κ values compared with the Mandard system (0.78 vs. 0.62). Evaluation of the whole embedded tumor site showed improved results (Becker: 0.83; Mandard: 0.73) as compared with only 1 representative slide (Becker: 0.68; Mandard: 0.71). Modification into simplified 3-tiered systems showed comparable interobserver agreement but better prognostic stratification for both systems (log rank Becker: P=0.015; Mandard P=0.03), with independent prognostic impact for overall survival (modified Becker: P=0.011, hazard ratio=3.07; modified Mandard: P=0.023, hazard ratio=2.72). In conclusion, both systems provide substantial to excellent interobserver agreement for estimation of tumor regression after neoadjuvant chemotherapy in esophageal adenocarcinomas. A simple 3-tiered system with the estimation of residual tumor in % (complete regression/1% to 50% residual tumor/>50% residual tumor) maintains the highest reproducibility and prognostic value.

Prognostic value of histopathological regression in 850 neoadjuvantly treated oesophagogastric adenocarcinomas.

BACKGROUND Recently, histopathological tumour regression, prevalence of signet ring cells, and localisation were reported as prognostic factors in neoadjuvantly treated oesophagogastric (junctional and gastric) cancer. This exploratory retrospective study analyses independent prognostic factors within a large patient cohort after preoperative chemotherapy including clinical and histopathological factors. METHODS In all, 850 patients presenting with oesophagogastric cancer staged cT3/4 Nany cM0/x were treated with neoadjuvant chemotherapy followed by resection in two academic centres. Patient data were documented in a prospective database and retrospectively analysed. RESULTS Of all factors prognostic on univariate analysis, only clinical response, complications, ypTNM stage, and R category were independently prognostic (P<0.01) on multivariate analysis. Tumour localisation and signet ring cells were independently prognostic only when investigator-dependent clinical response evaluation was excluded from the multivariate model. Histopathological tumour regression correlates with tumour grading, Laurén classification, clinical response, ypT, ypN, and R categories but was not identified as an independent prognostic factor. Within R0-resected patients only surgical complications and ypTNM stage were independent prognostic factors. CONCLUSIONS Only established prognostic factors like ypTNM stage, R category, and complications were identified as independent prognostic factors in resected patients after neoadjuvant chemotherapy. In contrast, histopathological tumour regression was not found as an independent prognostic marker.

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.