ROBREG: Stata module providing robust regression estimators

robreg provides a number of robust estimators for linear regression models. Among them are the high breakdown-point and high efficiency MM-estimator, the Huber and bisquare M-estimator, and the S-estimator, each supporting classic or robust standard errors. Furthermore, basic versions of the LMS/LQS (least median of squares) and LTS (least trimmed squares) estimators are provided. Note that the moremata package, also available from SSC, is required.

Prognostic factors and outcomes in primary urethral cancer: results from the international collaboration on primary urethral carcinoma.

PURPOSE To evaluate risk factors for survival in a large international cohort of patients with primary urethral cancer (PUC). METHODS A series of 154 patients (109 men, 45 women) were diagnosed with PUC in ten referral centers between 1993 and 2012. Kaplan-Meier analysis with log-rank test was used to investigate various potential prognostic factors for recurrence-free (RFS) and overall survival (OS). Multivariate models were constructed to evaluate independent risk factors for recurrence and death. RESULTS Median age at definitive treatment was 66 years (IQR 58-76). Histology was urothelial carcinoma in 72 (47 %), squamous cell carcinoma in 46 (30 %), adenocarcinoma in 17 (11 %), and mixed and other histology in 11 (7 %) and nine (6 %), respectively. A high degree of concordance between clinical and pathologic nodal staging (cN+/cN0 vs. pN+/pN0; p < 0.001) was noted. For clinical nodal staging, the corresponding sensitivity, specificity, and overall accuracy for predicting pathologic nodal stage were 92.8, 92.3, and 92.4 %, respectively. In multivariable Cox-regression analysis for patients staged cM0 at initial diagnosis, RFS was significantly associated with clinical nodal stage (p < 0.001), tumor location (p < 0.001), and age (p = 0.001), whereas clinical nodal stage was the only independent predictor for OS (p = 0.026). CONCLUSIONS These data suggest that clinical nodal stage is a critical parameter for outcomes in PUC.

Increasingly Rank: The Use and Misuse of Rankings in Economic Development

Debates over the merits of competing schemes for ranking metropolitan areas as hightech centers shed little light on the important policy questions that should be the core of economic development policy. There are no strong theoretical reasons for preferring one ranking system to others. Rankings often conflate different industries and ignore history, obscuring the varied and often idiosyncratic processes that drive growth in different regions. Although an occupational perspective is a useful one for examining economic activity, it is a supplement to, not a replacement for, a careful understanding of metropolitan industrial specialization. Practitioners should not put too much weight on any ranking system but instead should work to develop detailed knowledge of their region’s special economic niche and to develop relationships and strategies that build on established strengths.

Modeling absolute differences in life expectancy with a censored skew-normal regression approach.

Parameter estimates from commonly used multivariable parametric survival regression models do not directly quantify differences in years of life expectancy. Gaussian linear regression models give results in terms of absolute mean differences, but are not appropriate in modeling life expectancy, because in many situations time to death has a negative skewed distribution. A regression approach using a skew-normal distribution would be an alternative to parametric survival models in the modeling of life expectancy, because parameter estimates can be interpreted in terms of survival time differences while allowing for skewness of the distribution. In this paper we show how to use the skew-normal regression so that censored and left-truncated observations are accounted for. With this we model differences in life expectancy using data from the Swiss National Cohort Study and from official life expectancy estimates and compare the results with those derived from commonly used survival regression models. We conclude that a censored skew-normal survival regression approach for left-truncated observations can be used to model differences in life expectancy across covariates of interest.

Consistency of Concave Regression with an Application to Current-Status Data

We consider the problem of nonparametric estimation of a concave regression function F. We show that the supremum distance between the least square s estimatorand F on a compact interval is typically of order(log(n)/n)2/5. This entails rates of convergence for the estimator’s derivative. Moreover, we discuss the impact of additional constraints on F such as monotonicity and pointwise bounds. Then we apply these results to the analysis of current status data, where the distribution function of the event times is assumed to be concave.

Interlaboratory variability of MIB1 staining in well-differentiated pancreatic neuroendocrine tumors

Neuroendocrine tumors (NET) are routinely graded and staged to judge prognosis. Proliferation index using MIB1 staining has been introduced to assess grading. There are vivid discussions on cutoff definitions, automated counting, and interobserver variability. However, no data exist regarding interlaboratory reproducibility for low proliferation indices which are of importance to discriminate between G1 and G2 NET. We performed MIB1 staining in three different university hospital-based pathology laboratories on a tissue micro array (TMA) of a well-characterized patient cohort, containing pancreatic NET of 61 patients. To calculate the proliferation index, number of positive tumor nuclei was divided by the total number of tumor nuclei. Labeling index was compared to mitotic counts in whole tissue sections and to clinical outcome. Linear regression analysis, intraclass comparison, and log-rank analysis were performed. Intraclass correlation showed moderate-to-fair agreement. Especially low proliferating tumors were affected by interlaboratory differences. Log-rank analysis was performed for each lab and resulted in three different cutoffs (5.0, 3.0, and 0.5 %). Every calculated cutoff stratified the patient cohort to a significant extent for the underlying stain (p < 0.001, <0.001, and <0.001) but showed no or lesser significance when applied to the other stains. Significant and relevant interlab differences for MIB1 exist. Since the MIB1 proliferation index influences grading, local cutoffs or external standardization should urgently be introduced to achieve reliability and reproducibility.

Detecting outlying studies in meta-regression models using a forward search algorithm.

When considering data from many trials, it is likely that some of them present a markedly different intervention effect or exert an undue influence on the summary results. We develop a forward search algorithm for identifying outlying and influential studies in meta-analysis models. The forward search algorithm starts by fitting the hypothesized model to a small subset of likely outlier-free studies and proceeds by adding studies into the set one-by-one that are determined to be closest to the fitted model of the existing set. As each study is added to the set, plots of estimated parameters and measures of fit are monitored to identify outliers by sharp changes in the forward plots. We apply the proposed outlier detection method to two real data sets; a meta-analysis of 26 studies that examines the effect of writing-to-learn interventions on academic achievement adjusting for three possible effect modifiers, and a meta-analysis of 70 studies that compares a fluoride toothpaste treatment to placebo for preventing dental caries in children. A simple simulated example is used to illustrate the steps of the proposed methodology, and a small-scale simulation study is conducted to evaluate the performance of the proposed method. Copyright © 2016 John Wiley & Sons, Ltd.

Regression of coronary atherosclerosis: Current evidence and future perspectives

Coronary atherosclerosis has been considered a chronic disease characterized by ongoing progression in response to systemic risk factors and local pro-atherogenic stimuli. As our understanding of the pathobiological mechanisms implicated in atherogenesis and plaque progression is evolving, effective treatment strategies have been developed that led to substantial reduction of the clinical manifestations and acute complications of coronary atherosclerotic disease. More recently, intracoronary imaging modalities have enabled detailed in vivo quantification and characterization of coronary atherosclerotic plaque, serial evaluation of atherosclerotic changes over time, and assessment of vascular responses to effective anti-atherosclerotic medications. The use of intracoronary imaging modalities has demonstrated that intensive lipid lowering can halt plaque progression and may even result in regression of coronary atheroma when the highest doses of the most potent statins are used. While current evidence indicates the feasibility of atheroma regression and of reversal of presumed high-risk plaque characteristics in response to intensive anti-atherosclerotic therapies, these changes of plaque size and composition are modest and their clinical implications remain largely elusive. Growing interest has focused on achieving more pronounced regression of coronary plaque using novel anti-atherosclerotic medications, and more importantly on elucidating ways toward clinical translation of favorable changes of plaque anatomy into more favorable clinical outcomes for our patients.

Eosinophilic oesophagitis: relationship of quality of life with clinical, endoscopic and histological activity.

BACKGROUND Knowledge about determinants of quality of life (QoL) in eosinophilic oesophagitis (EoO) patients helps to identify patients at risk of experiencing poor QoL and to tailor therapeutic interventions accordingly. AIM To evaluate the impact of symptom severity, endoscopic and histological activity on EoE-specific QoL in adult EoE patients. METHODS Ninety-eight adult EoE patients were prospectively included (64% male, median age 39 years). Patients completed two validated instruments to assess EoE-specific QoL (EoO-QoL-A) and symptom severity (adult EoE activity index patient-reported outcome) and then underwent esophagogastroduodenoscopy with biopsy sampling. Physicians reported standardised information on EoE-associated endoscopic and histological alterations. The Spearman's rank correlation coefficient was calculated to determine the relationship between QoL and symptom severity. Linear regression and analysis of variance was used to quantify the extent to which variations in severity of EoE symptoms, endoscopic and histological findings explain variations in QoL. RESULTS Quality of life strongly correlated with symptom severity (r = 0.610, P < 0.001). While the variation in severity of symptoms, endoscopic and histological findings alone explained 38%, 35% and 22% of the variability in EoE-related QoL, respectively, these together explained 60% of variation. Symptom severity explained 18-35% of the variation in each of the five QoL subscale scores. CONCLUSIONS Eosinophilic oesophagitis symptom severity and biological disease activity determine QoL in adult patients with eosinophilic oesophagitis. Therefore, reduction in both eosinophilic oesophagitis symptoms as well as biological disease activity is essential for improvement of QoL in adult patients. Clinicaltrials.gov number, NCT00939263.

Systematic analysis of factors associated with progression and regression of ulcerative colitis in 918 patients.

BACKGROUND Studies that systematically assess change in ulcerative colitis (UC) extent over time in adult patients are scarce. AIM To assess changes in disease extent over time and to evaluate clinical parameters associated with this change. METHODS Data from the Swiss IBD cohort study were analysed. We used logistic regression modelling to identify factors associated with a change in disease extent. RESULTS A total of 918 UC patients (45.3% females) were included. At diagnosis, UC patients presented with the following disease extent: proctitis [199 patients (21.7%)], left-sided colitis [338 patients (36.8%)] and extensive colitis/pancolitis [381 (41.5%)]. During a median disease duration of 9 [4-16] years, progression and regression was documented in 145 patients (15.8%) and 149 patients (16.2%) respectively. In addition, 624 patients (68.0%) had a stable disease extent. The following factors were identified to be associated with disease progression: treatment with systemic glucocorticoids [odds ratio (OR) 1.704, P = 0.025] and calcineurin inhibitors (OR: 2.716, P = 0.005). No specific factors were found to be associated with disease regression. CONCLUSIONS Over a median disease duration of 9 [4-16] years, about two-thirds of UC patients maintained the initial disease extent; the remaining one-third had experienced either progression or regression of the disease extent.

The Risk of Seizure After Surgery for Unruptured Intracranial Aneurysms: A Prospective Cohort Study.

BACKGROUND We aimed to identify a group of patients with a low risk of seizure after surgery for unruptured intracranial aneurysms (UIA). OBJECTIVE To determine the risk of seizure after discharge from surgery for UIA. METHODS A consecutive prospectively collected cohort database was interrogated for all surgical UIA cases. There were 726 cases of UIA (excluding cases proximal to the superior cerebellar artery on the vertebrobasilar system) identified and analyzed. Cox proportional hazards regression models and Kaplan-Meier life table analyses were generated assessing risk factors. RESULTS Preoperative seizure history and complication of aneurysm repair were the only risk factors found to be significant. The risk of first seizure after discharge from hospital following surgery for patients with neither preoperative seizure, treated middle cerebral artery aneurysm, nor postoperative complications (leading to a modified Rankin Scale score >1) was <0.1% and 1.1% at 12 months and 7 years, respectively. The risk for those with preoperative seizures was 17.3% and 66% at 12 months and 7 years, respectively. The risk for seizures with either complications (leading to a modified Rankin Scale score >1) from surgery or treated middle cerebral artery aneurysm was 1.4% and 6.8% at 12 months and 7 years, respectively. These differences in the 3 Kaplan-Meier curves were significant (log-rank P < .001). CONCLUSION The risk of seizures after discharge from hospital following surgery for UIA is very low when there is no preexisting history of seizures. If this result can be supported by other series, guidelines that restrict returning to driving because of the risk of postoperative seizures should be reconsidered. ABBREVIATIONS MCA, middle cerebral arterymRS, modified Rankin ScaleUIA, unruptured intracranial aneurysms.

Older patients with low Charlson score and high-risk prostate cancer benefit from radical prostatectomy

INTRODUCTION The aim of the study was to identify the appropriate level of Charlson comorbidity index (CCI) in older patients (>70 years) with high-risk prostate cancer (PCa) to achieve survival benefit following radical prostatectomy (RP). METHODS We retrospectively analyzed 1008 older patients (>70 years) who underwent RP with pelvic lymph node dissection for high-risk prostate cancer (preoperative prostate-specific antigen >20 ng/mL or clinical stage ≥T2c or Gleason ≥8) from 14 tertiary institutions between 1988 and 2014. The study population was further grouped into CCI < 2 and ≥2 for analysis. Survival rate for each group was estimated with Kaplan-Meier method and competitive risk Fine-Gray regression to estimate the best explanatory multivariable model. Area under the curve (AUC) and Akaike information criterion were used to identify ideal 'Cut off' for CCI. RESULTS The clinical and cancer characteristics were similar between the two groups. Comparison of the survival analysis using the Kaplan-Meier curve between two groups for non-cancer death and survival estimations for 5 and 10 years shows significant worst outcomes for patients with CCI ≥ 2. In multivariate model to decide the appropriate CCI cut-off point, we found CCI 2 has better AUC and p value in log rank test. CONCLUSION Older patients with fewer comorbidities harboring high-risk PCa appears to benefit from RP. Sicker patients are more likely to die due to non-prostate cancer-related causes and are less likely to benefit from RP.

Environmental correlates of species rank − abundance distributions in global drylands

Theoretical models predict lognormal species abundance distributions (SADs) in stable and productive environments, with log-series SADs in less stable, dispersal driven communities. We studied patterns of relative species abundances of perennial vascular plants in global dryland communities to: (i) assess the influence of climatic and soil characteristics on the observed SADs, (ii) infer how environmental variability influences relative abundances, and (iii) evaluate how colonisation dynamics and environmental filters shape abundance distributions. We fitted lognormal and log-series SADs to 91 sites containing at least 15 species of perennial vascular plants. The dependence of species relative abundances on soil and climate variables was assessed using general linear models. Irrespective of habitat type and latitude, the majority of the SADs (70.3%) were best described by a lognormal distribution. Lognormal SADs were associated with low annual precipitation, higher aridity, high soil carbon content, and higher variability of climate variables and soil nitrate. Our results do not corroborate models predicting the prevalence of log-series SADs in dryland communities. As lognormal SADs were particularly associated with sites with drier conditions and a higher environmental variability, we reject models linking lognormality to environmental stability and high productivity conditions. Instead our results point to the prevalence of lognormal SADs in heterogeneous environments, allowing for more evenly distributed plant communities, or in stressful ecosystems, which are generally shaped by strong habitat filters and limited colonisation. This suggests that drylands may be resilient to environmental changes because the many species with intermediate relative abundances could take over ecosystem functioning if the environment becomes suboptimal for dominant species.

ALTERATIONS IN HYPOTHALAMIC CONTENT OF LUTEINIZING HORMONE RELEASING HORMONE AND PITUITARY RESPONSIVENESS ASSOCIATED WITH TESTICULAR REGRESSION INDUCED BY PINEAL RELATED TREATMENTS IN THE GOLDEN HAMSTER

The adult male golden hamster, when exposed to blinding (BL), short photoperiod (SP), or daily melatonin injections (MEL) demonstrates dramatic reproductive collapse. This collapse can be blocked by removal of the pineal gland prior to treatment. Reproductive collapse is characterized by a dramatic decrease in both testicular weight and serum gonadotropin titers. The present study was designed to examine the interactions of the hypothalamus and pituitary gland during testicular regression, and to specifically compare and contrast changes caused by the three commonly employed methods of inducing testicular regression (BL,SP,MEL). Hypothalamic LHRH content was altered by all three treatments. There was an initial increase in content of LHRH that occurred concomitantly with the decreased serum gonadotropin titers, followed by a precipitous decline in LHRH content which reflected the rapid increases in both serum LH and FSH which occur during spontaneous testicular recrudescence. In vitro pituitary responsiveness was altered by all three treatments: there was a decline in basal and maximally stimulatable release of both LH and FSH which paralleled the fall of serum gonadotropins. During recrudescence both basal and maximal release dramatically increased in a manner comparable to serum hormone levels. While all three treatments were equally effective in their ability to induce changes at all levels of the endocrine system, there were important temporal differences in the effects of the various treatments. Melatonin injections induced the most rapid changes in endocrine parameters, followed by exposure to short photoperiod. Blinding required the most time to induce the same changes. This study has demonstrated that pineal-mediated testicular regression is a process which involves dynamic changes in multiply-dependent endocrine relationships, and proper evaluation of these changes must be performed with specific temporal events in mind. ^

Physiogenomic Analysis of Weight Loss Induced by Dietary Carbohydrate Restriction

ABSTRACT : BACKGROUND : Diets that restrict carbohydrate (CHO) have proven to be a successful dietary treatment of obesity for many people, but the degree of weight loss varies across individuals. The extent to which genetic factors associate with the magnitude of weight loss induced by CHO restriction is unknown. We examined associations among polymorphisms in candidate genes and weight loss in order to understand the physiological factors influencing body weight responses to CHO restriction. METHODS : We screened for genetic associations with weight loss in 86 healthy adults who were instructed to restrict CHO to a level that induced a small level of ketosis (CHO ~10% of total energy). A total of 27 single nucleotide polymorphisms (SNPs) were selected from 15 candidate genes involved in fat digestion/metabolism, intracellular glucose metabolism, lipoprotein remodeling, and appetite regulation. Multiple linear regression was used to rank the SNPs according to probability of association, and the most significant associations were analyzed in greater detail. RESULTS : Mean weight loss was 6.4 kg. SNPs in the gastric lipase (LIPF), hepatic glycogen synthase (GYS2), cholesteryl ester transfer protein (CETP) and galanin (GAL) genes were significantly associated with weight loss. CONCLUSION : A strong association between weight loss induced by dietary CHO restriction and variability in genes regulating fat digestion, hepatic glucose metabolism, intravascular lipoprotein remodeling, and appetite were detected. These discoveries could provide clues to important physiologic adaptations underlying the body mass response to CHO restriction.