Stochastic search for semiparametric linear regression models

This paper introduces and analyzes a stochastic search method for parameter estimation in linear regression models in the spirit of Beran and Millar [Ann. Statist. 15(3) (1987) 1131–1154]. The idea is to generate a random finite subset of a parameter space which will automatically contain points which are very close to an unknown true parameter. The motivation for this procedure comes from recent work of Dümbgen et al. [Ann. Statist. 39(2) (2011) 702–730] on regression models with log-concave error distributions.

The Blinder-Oaxaca decomposition for linear regression models

The counterfactual decomposition technique popularized by Blinder (1973, Journal of Human Resources, 436–455) and Oaxaca (1973, International Economic Review, 693–709) is widely used to study mean outcome differences between groups. For example, the technique is often used to analyze wage gaps by sex or race. This article summarizes the technique and addresses several complications, such as the identification of effects of categorical predictors in the detailed decomposition or the estimation of standard errors. A new command called oaxaca is introduced, and examples illustrating its usage are given.

On Graphically Checking Goodness-of-fit of Binary Logistic Regression Models

OBJECTIVES: This paper is concerned with checking goodness-of-fit of binary logistic regression models. For the practitioners of data analysis, the broad classes of procedures for checking goodness-of-fit available in the literature are described. The challenges of model checking in the context of binary logistic regression are reviewed. As a viable solution, a simple graphical procedure for checking goodness-of-fit is proposed. METHODS: The graphical procedure proposed relies on pieces of information available from any logistic analysis; the focus is on combining and presenting these in an informative way. RESULTS: The information gained using this approach is presented with three examples. In the discussion, the proposed method is put into context and compared with other graphical procedures for checking goodness-of-fit of binary logistic models available in the literature. CONCLUSION: A simple graphical method can significantly improve the understanding of any logistic regression analysis and help to prevent faulty conclusions.

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.

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.

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.

A longitudinal three-center study of dental arch relationship in patients with bilateral cleft lip and palate

OBJECTIVE: To compare and evaluate longitudinally the dental arch relationships from 4.5 to 13.5 years of age with the Bauru-BCLP Yardstick in a large sample of patients with bilateral cleft lip and palate (BCLP). DESIGN: Retrospective longitudinal intercenter outcome study. PATIENTS: Dental casts of 204 consecutive patients with complete BCLP were evaluated at 6, 9, and 12 years of age. All models were identified only by random identification numbers. SETTING: Three cleft palate centers with different treatment protocols. MAIN OUTCOME MEASURES: Dental arch relationships were categorized with the Bauru-BCLP yardstick. Increments for each interval (from 6 to 9 years, 6 to 12 years, and 9 to 12 years) were analyzed by logistic and linear regression models. RESULTS: There were no significant differences in outcome measures between the centers at age 12 or at age 9. At age 6, center B showed significantly better results (p=.027), but this difference diminished as the yardstick score for this group increased over time (linear regression analysis), the difference with the reference category (center C, boys) for the intervals 6 to 12 and 9 to 12 years being 10.4% (p=.041) and 12.9% (p=.009), respectively. CONCLUSIONS: Despite different treatment protocols, dental arch relationships in the three centers were comparable in final scores at age 9 and 12 years. Delaying hard palate closure and employing infant orthopedics did not appear to be advantageous in the long run. Premaxillary osteotomy employed in center B appeared to be associated with less favorable development of the dental arch relationship between 9 and 12 years.

Sonographic follow-up of patients with carpal tunnel syndrome undergoing surgical or nonsurgical treatment: prospective cohort study

Purpose To compare changes in the largest cross-sectional area (CSA) of the median nerve in wrists undergoing surgical decompression with changes in wrists undergoing non-surgical treatment of carpal tunnel syndrome (CTS). Methods This study was a prospective cohort study in 55 consecutive patients with 78 wrists with established CTS, including 60 wrists treated with surgical decompression and 18 wrists with non-surgical treatment. A sonographic examination was scheduled before and 4 months after initiation of treatment. We compared changes in CSA of the median nerve between wrists with surgical treatment and wrists with non-surgical treatment using linear regression models. Results Decreases in CSA of the median nerve were more pronounced in wrists with CTS release than in wrists undergoing nonsurgical treatment (difference in means, 1.0 mm2; 95% confidence interval, 0.3–1.8 mm2). Results were robust to the adjustment for age, gender, and neurological severity at baseline. Among wrists with CTS release, those with postoperative CSA of 10 mm2 or less tended to have better clinical outcomes than those with postoperative CSA of greater than 10 mm2 (p=.055). Postoperative sonographic workup in the 3 patients with unfavorable outcome or recurrence identified likely causes for treatment failure in 2 patients. Conclusions In this observational study, surgical decompression was associated with a greater decrease in median nerve CSA than was nonsurgical treatment. Smaller postoperative CSAs may be associated with better clinical outcomes. Additional randomized trials are necessary to determine the optimal treatment strategy in different subgroups of patients with CTS. Type of study/level of evidence Therapeutic III.

Population-based analysis of 4113 patients with acute cholecystitis: defining the optimal time-point for laparoscopic cholecystectomy

Objective: To compare clinical outcomes after laparoscopic cholecystectomy (LC) for acute cholecystitis performed at various time-points after hospital admission. Background: Symptomatic gallstones represent an important public health problem with LC the treatment of choice. LC is increasingly offered for acute cholecystitis, however, the optimal time-point for LC in this setting remains a matter of debate. Methods: Analysis was based on the prospective database of the Swiss Association of Laparoscopic and Thoracoscopic Surgery and included patients undergoing emergency LC for acute cholecystitis between 1995 and 2006, grouped according to the time-points of LC since hospital admission (admission day (d0), d1, d2, d3, d4/5, d ≥6). Linear and generalized linear regression models assessed the effect of timing of LC on intra- or postoperative complications, conversion and reoperation rates and length of postoperative hospital stay. Results: Of 4113 patients, 52.8% were female, median age was 59.8 years. Delaying LC resulted in significantly higher conversion rates (from 11.9% at d0 to 27.9% at d ≥6 days after admission, P < 0.001), surgical postoperative complications (5.7% to 13%, P < 0.001) and re-operation rates (0.9% to 3%, P = 0.007), with a significantly longer postoperative hospital stay (P < 0.001). Conclusions: Delaying LC for acute cholecystitis has no advantages, resulting in significantly increased conversion/re-operation rate, postoperative complications and longer postoperative hospital stay. This investigation—one of the largest in the literature—provides compelling evidence that acute cholecystitis merits surgery within 48 hours of hospital admission if impact on the patient and health care system is to be minimized.

Predictors and overestimation of recalled mobile phone use among children and adolescents

A growing body of literature addresses possible health effects of mobile phone use in children and adolescents by relying on the study participants' retrospective reconstruction of mobile phone use. In this study, we used data from the international case-control study CEFALO to compare self-reported with objectively operator-recorded mobile phone use. The aim of the study was to assess predictors of level of mobile phone use as well as factors that are associated with overestimating own mobile phone use. For cumulative number and duration of calls as well as for time since first subscription we calculated the ratio of self-reported to operator-recorded mobile phone use. We used multiple linear regression models to assess possible predictors of the average number and duration of calls per day and logistic regression models to assess possible predictors of overestimation. The cumulative number and duration of calls as well as the time since first subscription of mobile phones were overestimated on average by the study participants. Likelihood to overestimate number and duration of calls was not significantly different for controls compared to cases (OR=1.1, 95%-CI: 0.5 to 2.5 and OR=1.9, 95%-CI: 0.85 to 4.3, respectively). However, likelihood to overestimate was associated with other health related factors such as age and sex. As a consequence, such factors act as confounders in studies relying solely on self-reported mobile phone use and have to be considered in the analysis.

Do Central Hypersensitivity and Altered Pain Modulation Predict the Course of Chronic Low Back and Neck Pain?

OBJECTIVES:: Widespread central hypersensitivity and altered conditioned pain modulation (CPM) have been documented in chronic pain conditions. Information on their prognostic values is limited. This study tested the hypothesis that widespread central hypersensitivity (WCH) and altered CPM, assessed during the chronic phase of low back and neck pain, predict poor outcome. METHODS:: A total of 169 consecutive patients with chronic low back or neck pain, referred to the pain clinic during 1 year, were analyzed. Pressure pain tolerance threshold at the second toe and tolerance time during cold pressor test at the hand assessed WCH. CPM was measured by the change in pressure pain tolerance threshold (test stimulus) after cold pressor test (conditioning stimulus). A structured telephone interview was performed 12 to 15 months after testing to record outcome parameters. Linear regression models were used, with average and maximum pain intensity of the last 24 hours at follow-up as endpoints. Multivariable analyses included sex, age, catastrophizing scale, Beck Depression Inventory, pain duration, intake of opioids, and type of pain syndrome. RESULTS:: Statistically significant reductions from baseline to follow-up were observed in pain intensity (P<0.001). No evidence for an association between the measures of WCH or CPM and intensity of chronic pain at follow-up was found. DISCUSSION:: A major predictive value of the measures that we used is unlikely. Future studies adopting other assessment modalities and possibly standardized treatments are needed to further elucidate the prognostic value of WCH and altered CPM in chronic pain.

Physical activity intensity and surrogate markers for cardiovascular health in adolescents

We examined the impact of physical activity (PA) on surrogate markers of cardiovascular health in adolescents. 52 healthy students (28 females, mean age 14.5 ± 0.7 years) were investigated. Microvascular endothelial function was assessed by peripheral arterial tonometry to determine reactive hyperemic index (RHI). Vagal activity was measured using 24 h analysis of heart rate variability [root mean square of successive normal-to-normal intervals (rMSSD)]. Exercise testing was performed to determine peak oxygen uptake ([Formula: see text]) and maximum power output. PA was assessed by accelerometry. Linear regression models were performed and adjusted for age, sex, skinfolds, and pubertal status. The cohort was dichotomized into two equally sized activity groups (low vs. high) based on the daily time spent in moderate-to-vigorous PA (MVPA, 3,000-5,200 counts(.)min(-1), model 1) and vigorous PA (VPA, >5,200 counts(.)min(-1), model 2). MVPA was an independent predictor for rMSSD (β = 0.448, P = 0.010), and VPA was associated with maximum power output (β = 0.248, P = 0.016). In model 1, the high MVPA group exhibited a higher vagal tone (rMSSD 49.2 ± 13.6 vs. 38.1 ± 11.7 ms, P = 0.006) and a lower systolic blood pressure (107.3 ± 9.9 vs. 112.9 ± 8.1 mmHg, P = 0.046). In model 2, the high VPA group had higher maximum power output values (3.9 ± 0.5 vs. 3.4 ± 0.5 W kg(-1), P = 0.012). In both models, no significant differences were observed for RHI and [Formula: see text]. In conclusion, in healthy adolescents, PA was associated with beneficial intensity-dependent effects on vagal tone, systolic blood pressure, and exercise capacity, but not on microvascular endothelial function.

Stable virulence levels in the HIV epidemic of Switzerland over two decades

OBJECTIVE: To determine whether the virulence of HIV-1 has been changing since its introduction into Switzerland. DESIGN: A prospective cohort study of HIV-1 infected individuals with well-characterized pre-therapy disease history. METHODS: To minimize the effect of recently imported viruses and ethnicity-associated host factors, the analysis was restricted to the white, north-west-European majority population of the cohort. Virulence was characterized by the decline slope of the CD4 cell count (n = 817 patients), the decline slope of the CD4:CD8 ratio (n = 815 patients) and the viral setpoint (n = 549 patients) in untreated patients with sufficient data points. Linear regression models were used to detect correlations between the date of diagnosis (ranging between 1984 and 2003) and the virulence markers, controlling for gender, exposure category, age and CD4 cell count at entry. RESULTS: We found no correlation between any of the virulence markers and the date of diagnosis. Inspection of short-term trends confirmed that virulence has fluctuated around a stable level over time. CONCLUSIONS: The lack of long-term time trends in the virulence markers indicates that HIV-1 is not evolving towards increasing or decreasing virulence at a perceptible rate. Both highly virulent and attenuated strains have apparently been unable to spread at the population level. This result suggests that either the evolution of virulence may be slow or inhibited due to evolutionary constraints, or HIV-1 may have already evolved to optimal virulence in the human host.