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.

Oscillometric measurement of ankle-brachial index in patients with suspected peripheral disease: comparison with Doppler method

QUESTION UNDER STUDY: Purpose was to validate accuracy and reliability of automated oscillometric ankle-brachial (ABI) measurement prospectively against the current gold standard of Doppler-assisted ABI determination. METHODS: Oscillometric ABI was measured in 50 consecutive patients with peripheral arterial disease (n = 100 limbs, mean age 65 +/- 6 years, 31 men, 19 diabetics) after both high and low ABI had been determined conventionally by Doppler under standardised conditions. Correlation was assessed by linear regression and Pearson product moment correlation. Degree of inter-modality agreement was quantified by use of Bland and Altman method. RESULTS: Oscillometry was performed significantly faster than Doppler-assisted ABI (3.9 +/- 1.3 vs 11.4 +/- 3.8 minutes, P <0.001). Mean readings were 0.62 +/- 0.25, 0.70 +/- 0.22 and 0.63 +/- 0.39 for low, high and oscillometric ABI, respectively. Correlation between oscillometry and Doppler ABI was good overall (r = 0.76 for both low and high ABI) and excellent in oligo-symptomatic, non-diabetic patients (r = 0.81; 0.07 +/- 0.23); it was, however, limited in diabetic patients and in patients with critical limb ischaemia. In general, oscillometric ABI readings were slightly higher (+0.06), but linear regression analysis showed that correlation was sustained over the whole range of measurements. CONCLUSIONS: Results of automated oscillometric ABI determination correlated well with Doppler-assisted measurements and could be obtained in shorter time. Agreement was particularly high in oligo-symptomatic non-diabetic patients.

Tidal volume single breath washout of two tracer gases--a practical and promising lung function test

Background Small airway disease frequently occurs in chronic lung diseases and may cause ventilation inhomogeneity (VI), which can be assessed by washout tests of inert tracer gas. Using two tracer gases with unequal molar mass (MM) and diffusivity increases specificity for VI in different lung zones. Currently washout tests are underutilised due to the time and effort required for measurements. The aim of this study was to develop and validate a simple technique for a new tidal single breath washout test (SBW) of sulfur hexafluoride (SF6) and helium (He) using an ultrasonic flowmeter (USFM). Methods The tracer gas mixture contained 5% SF6 and 26.3% He, had similar total MM as air, and was applied for a single tidal breath in 13 healthy adults. The USFM measured MM, which was then plotted against expired volume. USFM and mass spectrometer signals were compared in six subjects performing three SBW. Repeatability and reproducibility of SBW, i.e., area under the MM curve (AUC), were determined in seven subjects performing three SBW 24 hours apart. Results USFM reliably measured MM during all SBW tests (n = 60). MM from USFM reflected SF6 and He washout patterns measured by mass spectrometer. USFM signals were highly associated with mass spectrometer signals, e.g., for MM, linear regression r-squared was 0.98. Intra-subject coefficient of variation of AUC was 6.8%, and coefficient of repeatability was 11.8%. Conclusion The USFM accurately measured relative changes in SF6 and He washout. SBW tests were repeatable and reproducible in healthy adults. We have developed a fast, reliable, and straightforward USFM based SBW method, which provides valid information on SF6 and He washout patterns during tidal breathing.

Estimation of pelvic tilt on anteroposterior X-rays--a comparison of six parameters

OBJECTIVE: To compare six different parameters described in literature for estimation of pelvic tilt on an anteroposterior pelvic radiograph and to create a simple nomogram for tilt correction of prosthetic cup version in total hip arthroplasty. DESIGN: Simultaneous anteroposterior and lateral pelvic radiographs are taken routinely in our institution and were analyzed prospectively. The different parameters (including three distances and three ratios) were measured and compared to the actual pelvic tilt on the lateral radiograph using simple linear regression analysis. PATIENTS: One hundred and four consecutive patients (41 men, 63 women with a mean age of 31.7 years, SD 9.2 years, range 15.7-59.1 years) were studied. RESULTS: The strongest correlation between pelvic tilt and one of the six parameters for both men and women was the distance between the upper border of the symphysis and the sacrococcygeal joint. The correlation coefficient was 0.68 for men (P<0.001) and 0.61 for women (P<0.001). Based on this linear correlation, a nomogram was created that enables fast, tilt-corrected cup version measurements in clinical routine use. CONCLUSION: This simple method for correcting variations in pelvic tilt on plain radiographs can potentially improve the radiologist's ability to diagnose and interpret malformations of the acetabulum (particularly acetabular retroversion and excessive acetabular overcoverage) and post-operative orientation of the prosthetic acetabulum.

Anemia is associated with abdominal aortic aneurysm (AAA) size and decreased long-term survival after endovascular AAA repair

OBJECTIVE: Anemia is a common comorbid condition in various inflammatory states and an established predictor of mortality in patients with chronic heart failure, ischemic heart disease, and end-stage renal disease. The present study of patients with abdominal aortic aneurysm (AAA) undergoing endovascular repair (EVAR) assessed the relationships between baseline hemoglobin concentration and AAA size, as well as anemia and long-term survival. METHODS: Between March 1994 and November 2006, 711 patients (65 women, mean age 75.8 +/- 7.8 years) underwent elective EVAR. Anemia was defined as a hemoglobin level <13 g/dL in men and <12 g/dL in women. Post-EVAR mean follow-up was 48.3 +/- 32.0 months. Association of hemoglobin level with AAA size was assessed with multiple linear regression. Mortality was determined with use of the internet-based Social Security Death Index and the electronic hospital record. Kaplan-Meier survival curves of anemic and nonanemic patient groups were compared by the log-rank method. Multivariable logistic regression models were used to determine the influence of anemia on vital status after EVAR. RESULTS: A total of 218/711 (30.7%) of AAA patients undergoing EVAR had anemia at baseline. After adjustment for various risk factors, hemoglobin level was inversely related to maximum AAA diameter (beta: - .144, 95%-CI: -1.482 - .322, P = .002). Post-EVAR survival was 65.5% at 5 years and 44.4% at 10 years. In long-term follow-up, survival was significantly lower in patients with anemia as compared to patients without anemia (P < .0001 by log-rank). Baseline hemoglobin levels were independently related to long-term mortality in multivariable Cox regression analysis adjusted for various risk factors (adjusted HR: 0.866, 95% CI: .783 to .958, P = .005). Within this model, statin use (adjusted HR: .517, 95% CI: .308 to .868, P = .013) was independently related to long-term survival, whereas baseline AAA diameter (adjusted HR: 1.022, 95% CI: 1.009 to 1.036, P = .001) was an independently associated with increased mortality. CONCLUSIONS: Baseline hemoglobin concentration is independently associated with AAA size and reduced long-term survival following EVAR. Thus, the presence or absence of anemia offers a potential refinement of existing risk stratification instruments.

Bloodstream infections in pediatric ECLS: usefulness of daily blood culture monitoring and predictive value of biological markers. The British Columbia experience

INTRODUCTION: The incidence of bloodstream infection (BSI) in extracorporeal life support (ECLS) is reported between 0.9 and 19.5%. In January 2006, the Extracorporeal Life Support Organization (ELSO) reported an overall incidence of 8.78% distributed as follows: respiratory: 6.5% (neonatal), 20.8% (pediatric); cardiac: 8.2% (neonatal) and 12.6% (pediatric). METHOD: At BC Children's Hospital (BCCH) daily surveillance blood cultures (BC) are performed and antibiotic prophylaxis is not routinely recommended. Positive BC (BC+) were reviewed, including resistance profiles, collection time of BC+, time to positivity and mortality. White blood cell count, absolute neutrophile count, immature/total ratio, platelet count, fibrinogen and lactate were analyzed 48, 24 and 0 h prior to BSI. A univariate linear regression analysis was performed. RESULTS: From 1999 to 2005, 89 patients underwent ECLS. After exclusion, 84 patients were reviewed. The attack rate was 22.6% (19 BSI) and 13.1% after exclusion of coagulase-negative staphylococci (n = 8). BSI patients were significantly longer on ECLS (157 h) compared to the no-BSI group (127 h, 95% CI: 106-148). Six BSI patients died on ECLS (35%; 4 congenital diaphragmatic hernias, 1 hypoplastic left heart syndrome and 1 after a tetralogy repair). BCCH survival on ECLS was 71 and 58% at discharge, which is comparable to previous reports. No patient died primarily because of BSI. No BSI predictor was identified, although lactate may show a decreasing trend before BSI (P = 0.102). CONCLUSION: Compared with ELSO, the studied BSI incidence was higher with a comparable mortality. We speculate that our BSI rate is explained by underreporting of "contaminants" in the literature, the use of broad-spectrum antibiotic prophylaxis and a higher yield with daily monitoring BC. We support daily surveillance blood cultures as an alternative to antibiotic prophylaxis in the management of patients on ECLS.

Sampling-Design Effects on Properties of Species-Area Relationships - A Case Study from Estonian Dry Grassland Communities

Despite widespread use of species-area relationships (SARs), dispute remains over the most representative SAR model. Using data of small-scale SARs of Estonian dry grassland communities, we address three questions: (1) Which model describes these SARs best when known artifacts are excluded? (2) How do deviating sampling procedures (marginal instead of central position of the smaller plots in relation to the largest plot; single values instead of average values; randomly located subplots instead of nested subplots) influence the properties of the SARs? (3) Are those effects likely to bias the selection of the best model? Our general dataset consisted of 16 series of nested-plots (1 cm(2)-100 m(2), any-part system), each of which comprised five series of subplots located in the four corners and the centre of the 100-m(2) plot. Data for the three pairs of compared sampling designs were generated from this dataset by subsampling. Five function types (power, quadratic power, logarithmic, Michaelis-Menten, Lomolino) were fitted with non-linear regression. In some of the communities, we found extremely high species densities (including bryophytes and lichens), namely up to eight species in 1 cm(2) and up to 140 species in 100 m(2), which appear to be the highest documented values on these scales. For SARs constructed from nested-plot average-value data, the regular power function generally was the best model, closely followed by the quadratic power function, while the logarithmic and Michaelis-Menten functions performed poorly throughout. However, the relative fit of the latter two models increased significantly relative to the respective best model when the single-value or random-sampling method was applied, however, the power function normally remained far superior. These results confirm the hypothesis that both single-value and random-sampling approaches cause artifacts by increasing stochasticity in the data, which can lead to the selection of inappropriate models.

Development and validation of energy expenditure prediction models based on GT3X accelerometer data in 5- to 9-year-old children

Background: Accelerometry has been established as an objective method that can be used to assess physical activity behavior in large groups. The purpose of the current study was to provide a validated equation to translate accelerometer counts of the triaxial GT3X into energy expenditure in young children. Methods: Thirty-two children aged 5–9 years performed locomotor and play activities that are typical for their age group. Children wore a GT3X accelerometer and their energy expenditure was measured with indirect calorimetry. Twenty-one children were randomly selected to serve as development group. A cubic 2-regression model involving separate equations for locomotor and play activities was developed on the basis of model fit. It was then validated using data of the remaining children and compared with a linear 2-regression model and a linear 1-regression model. Results: All 3 regression models produced strong correlations between predicted and measured MET values. Agreement was acceptable for the cubic model and good for both linear regression approaches. Conclusions: The current linear 1-regression model provides valid estimates of energy expenditure for ActiGraph GT3X data for 5- to 9-year-old children and shows equal or better predictive validity than a cubic or a linear 2-regression model.

Predictors of depressed mood 12 months after injury. Contribution of self-efficacy and social support

Abstract Purpose: There is evidence that depressed mood and perception of pain are related in patients with chronic illness. However, how individual resources such as self-efficacy and social support play a role in this association remains unclear. The aim of this study was to investigate the influence of both variables as either moderator or mediator. Method: In a longitudinal study, 274 injured workers (M = 43.24 years) were investigated. Data were collected on sociodemographics, depressed mood, pain, social support, and self-efficacy at three months post-injury, and depressed mood one year post-injury. Results: Hierarchical multiple linear regression analyses revealed that pain (β = 0.14; p < 0.01) and social support (β = -0.18; p < 0.001) were significant predictors of depressed mood. Self-efficacy moderated the relationship of pain (β = -0.12; p < 0.05) and depressed mood after one year. Lower self-efficacy in combination with pain had a stronger impact than higher self-efficacy and pain on depressed mood. Social support did not moderate the association. Conclusions: Self-efficacy for managing pain is important in the development of depressed mood. According to the results of this study, we suggest that the detection of low social support and low self-efficacy might be important in long-term rehabilitation process. Implications for Rehabilitation Risk for depressed mood one year after an accident is high: One in five workers report depressed mood. Protective factors for depressed mood in injured workers needs to be considered in the rehabilitation. Focusing on resources like social support and self-efficacy could be protective against depressed mood. The early detection of low social support and low self-efficacy might be important in long-term rehabilitation processes

Living Renal Allograft Transplantation: Diffusion-weighted MR Imaging in Longitudinal Follow-up of the Donated and the Remaining Kidney.

Purpose To determine whether diffusion-weighted (DW) magnetic resonance (MR) imaging in living renal allograft donation allows monitoring of potential changes in the nontransplanted remaining kidney of the donor because of unilateral nephrectomy and changes in the transplanted kidney before and after transplantation in donor and recipient, respectively, and whether DW MR parameters are correlated in the same kidney before and after transplantation. Materials and Methods The study protocol was approved by the local ethics committee; written informed consent was obtained. Thirteen healthy kidney donors and their corresponding recipients prospectively underwent DW MR imaging (multiple b values) in donors before donation and in donors and recipients at day 8 and months 3 and 12 after donation. Total apparent diffusion coefficient (ADCT) values were determined; contribution of microcirculation was quantified in perfusion fraction (FP). Longitudinal changes of diffusion parameters were compared (repeated-measures one-way analysis of variance with post hoc pairwise comparisons). Correlations were tested (linear regression). Results ADCT values in nontransplanted kidney of donors increased from a preexplantation value of (188 ± 9 [standard deviation]) to (202 ± 11) × 10(-5) mm(2)/sec in medulla and from (199 ± 11) to (210 ± 13) × 10(-5) mm(2)/sec in cortex 1 week after donation (P < .004). Medullary, but not cortical, ADCT values stayed increased up to 1 year. ADCT values in allografts in recipients were stable. Compared with values obtained before transplantation in donors, the corticomedullary difference was reduced in allografts (P < .03). Cortical ADCT values correlated with estimated glomerular filtration rate in recipients (R = 0.56, P < .001) but not donors. Cortical ADCT values in the same kidney before transplantation in donors correlated with those in recipients on day 8 after transplantation (R = 0.77, P = .006). FP did not show significant changes. Conclusion DW MR imaging depicts early adaptations in the remaining nontransplanted kidney of donors after nephrectomy. All diffusion parameters remained constant in allograft recipients after transplantation. This method has potential monitoring utility, although assessment of clinical relevance is needed. © RSNA, 2013 Online supplemental material is available for this article.

Aptitude–Treatment interaction effects in psychooncological interventions

Objective: To examine aptitude–treatment interaction (ATI) effects in cancer patients receiving psychoonco- logical interventions (POIs). Method: N=36 cancer patients were treated with POI. Hierarchical linear regression was used to test two interaction effects between patient baseline characteristics (aptitudes) and process analyses of therapy sessions (treatment) on change in mental health during POI. Results: Patients with high emotional distress did best when their therapy reduced arousal, and patients with lower emotional distress benefited most if therapists emphasized arousal induction. The interaction between the coping style of the patient (internalizing vs. externalizing) and the focus of the treatment (emotion vs. behavior) did not predict POI outcomes. Conclusions: The ATI effect of patient's distress and therapist's arousal induction/reduction may help therapists to make differential treatment decisions in POI. Tailoring treatments to cancer patients based on their personal characteristics may enhance the effectiveness of POI.

Determination of respiratory gas flow by electrical impedance tomography in an animal model of mechanical ventilation

Background A recent method determines regional gas flow of the lung by electrical impedance tomography (EIT). The aim of this study is to show the applicability of this method in a porcine model of mechanical ventilation in healthy and diseased lungs. Our primary hypothesis is that global gas flow measured by EIT can be correlated with spirometry. Our secondary hypothesis is that regional analysis of respiratory gas flow delivers physiologically meaningful results. Methods In two sets of experiments n = 7 healthy pigs and n = 6 pigs before and after induction of lavage lung injury were investigated. EIT of the lung and spirometry were registered synchronously during ongoing mechanical ventilation. In-vivo aeration of the lung was analysed in four regions-of-interest (ROI) by EIT: 1) global, 2) ventral (non-dependent), 3) middle and 4) dorsal (dependent) ROI. Respiratory gas flow was calculated by the first derivative of the regional aeration curve. Four phases of the respiratory cycle were discriminated. They delivered peak and late inspiratory and expiratory gas flow (PIF, LIF, PEF, LEF) characterizing early or late inspiration or expiration. Results Linear regression analysis of EIT and spirometry in healthy pigs revealed a very good correlation measuring peak flow and a good correlation detecting late flow. PIFEIT = 0.702 · PIFspiro + 117.4, r2 = 0.809; PEFEIT = 0.690 · PEFspiro-124.2, r2 = 0.760; LIFEIT = 0.909 · LIFspiro + 27.32, r2 = 0.572 and LEFEIT = 0.858 · LEFspiro-10.94, r2 = 0.647. EIT derived absolute gas flow was generally smaller than data from spirometry. Regional gas flow was distributed heterogeneously during different phases of the respiratory cycle. But, the regional distribution of gas flow stayed stable during different ventilator settings. Moderate lung injury changed the regional pattern of gas flow. Conclusions We conclude that the presented method is able to determine global respiratory gas flow of the lung in different phases of the respiratory cycle. Additionally, it delivers meaningful insight into regional pulmonary characteristics, i.e. the regional ability of the lung to take up and to release air.

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.