PSPLINE: Stata module providing a penalized spline scatterplot smoother based on linear mixed model technology

Pspline uses xtmixed to fit a penalized spline regression and plots the smoothed function. Additional covariates can be specified to adjust the smooth and plot partial residuals.

Two-dimensional linear-combination model fitting of magnetic resonance spectra to define the macromolecule baseline using FiTAID, a Fitting Tool for Arrays of Interrelated Datasets

To propose the determination of the macromolecular baseline (MMBL) in clinical 1H MR spectra based on T(1) and T(2) differentiation using 2D fitting in FiTAID, a general Fitting Tool for Arrays of Interrelated Datasets.

A prediction model for assessing residential radon concentration in Switzerland

Indoor radon is regularly measured in Switzerland. However, a nationwide model to predict residential radon levels has not been developed. The aim of this study was to develop a prediction model to assess indoor radon concentrations in Switzerland. The model was based on 44,631 measurements from the nationwide Swiss radon database collected between 1994 and 2004. Of these, 80% randomly selected measurements were used for model development and the remaining 20% for an independent model validation. A multivariable log-linear regression model was fitted and relevant predictors selected according to evidence from the literature, the adjusted R², the Akaike's information criterion (AIC), and the Bayesian information criterion (BIC). The prediction model was evaluated by calculating Spearman rank correlation between measured and predicted values. Additionally, the predicted values were categorised into three categories (50th, 50th-90th and 90th percentile) and compared with measured categories using a weighted Kappa statistic. The most relevant predictors for indoor radon levels were tectonic units and year of construction of the building, followed by soil texture, degree of urbanisation, floor of the building where the measurement was taken and housing type (P-values <0.001 for all). Mean predicted radon values (geometric mean) were 66 Bq/m³ (interquartile range 40-111 Bq/m³) in the lowest exposure category, 126 Bq/m³ (69-215 Bq/m³) in the medium category, and 219 Bq/m³ (108-427 Bq/m³) in the highest category. Spearman correlation between predictions and measurements was 0.45 (95%-CI: 0.44; 0.46) for the development dataset and 0.44 (95%-CI: 0.42; 0.46) for the validation dataset. Kappa coefficients were 0.31 for the development and 0.30 for the validation dataset, respectively. The model explained 20% overall variability (adjusted R²). In conclusion, this residential radon prediction model, based on a large number of measurements, was demonstrated to be robust through validation with an independent dataset. The model is appropriate for predicting radon level exposure of the Swiss population in epidemiological research. Nevertheless, some exposure misclassification and regression to the mean is unavoidable and should be taken into account in future applications of the model.

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.

Quantitative 1H-magnetic resonance spectroscopy of human brain: Influence of composition and parameterization of the basis set in linear combination model-fitting

Localized short-echo-time (1)H-MR spectra of human brain contain contributions of many low-molecular-weight metabolites and baseline contributions of macromolecules. Two approaches to model such spectra are compared and the data acquisition sequence, optimized for reproducibility, is presented. Modeling relies on prior knowledge constraints and linear combination of metabolite spectra. Investigated was what can be gained by basis parameterization, i.e., description of basis spectra as sums of parametric lineshapes. Effects of basis composition and addition of experimentally measured macromolecular baselines were investigated also. Both fitting methods yielded quantitatively similar values, model deviations, error estimates, and reproducibility in the evaluation of 64 spectra of human gray and white matter from 40 subjects. Major advantages of parameterized basis functions are the possibilities to evaluate fitting parameters separately, to treat subgroup spectra as independent moieties, and to incorporate deviations from straightforward metabolite models. It was found that most of the 22 basis metabolites used may provide meaningful data when comparing patient cohorts. In individual spectra, sums of closely related metabolites are often more meaningful. Inclusion of a macromolecular basis component leads to relatively small, but significantly different tissue content for most metabolites. It provides a means to quantitate baseline contributions that may contain crucial clinical information.

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.

RRREG: Stata module to estimate linear probability model for randomized response data

rrreg fits a linear probability model for randomized response data

Interrelation of peri-operative morbidity and ASA class assignment in patients undergoing gynaecological surgery

OBJECTIVE: The aim of this study was to estimate intra- and post-operative risk using the American Society of Anaesthesiologists (ASA) classification which is an important predictor of an intervention and of the entire operating programme. STUDY DESIGN: In this retrospective study, 4435 consecutive patients undergoing elective and emergency surgery at the Gynaecological Clinic of the University Hospital of Zurich were included. The ASA classification for pre-operative risk assessment was determined by an anaesthesiologist after a thorough physical examination. We observed several pre-, intra- and post-operative parameters, such as age, body-mass-index, duration of anaesthesia, duration of surgery, blood loss, duration of post-operative stay, complicated post-operative course, morbidity and mortality. The investigation of different risk factors was achieved by a multiple linear regression model for log-transformed duration of hospitalisation. RESULTS: Age and obesity were responsible for a higher ASA classification. ASA grade correlates with the duration of anaesthesia and the duration of the surgery itself. There was a significant difference in blood loss between ASA grades I (113+/-195 ml) and III (222+/-470 ml) and between classes II (176+/-432 ml) and III. The duration of post-operative hospitalisation could also be correlated with ASA class. ASA class I=1.7+/-3.0 days, ASA class II=3.6+/-4.3 days, ASA class III=6.8+/-8.2 days, and ASA class IV=6.2+/-3.9 days. The mean post-operative in-hospital stay was 2.5+/-4.0 days without complications, and 8.7+/-6.7 days with post-operative complications. Multiple linear regression model showed that not only the ASA classification contained an important information for the duration of hospitalisation. Parameters such as age, class of diagnosis, post-operative complications, etc. also have an influence on the duration of hospitalisation. CONCLUSION: This study shows that the ASA classification can be used as a good and early available predictor for the planning of an intervention in gynaecological surgery. The ASA classification helps the surgeon to assess the peri-operative risk profile of which important information can be derived for the planning of the operation programme.

The influence of a continuous rate infusion of dexmedetomidine on the nociceptive withdrawal reflex and temporal summation during isoflurane anaesthesia in dogs

Objective To examine the influence of a low dose dexmedetomidine infusion on the nociceptive withdrawal reflex and temporal summation in dogs during isoflurane anaesthesia. Study design Prospective experimental blinded cross-over study. Animals Eight healthy mixed breed dogs, body weight Mean +/- SD 26.5 +/- 8.4 kg and age 25 +/- 16 months. Methods Anaesthesia was induced with propofol and maintained with isoflurane (Fe'ISO 1.3%) delivered in oxygen and air. After stabilization, baseline recordings (time 0) were obtained, then a dexmedetomidine bolus (1 mug kg(-1) IV) followed by a continuous rate infusion (1 mug kg(-1) hour(-1) ) or saline placebo were administered. At times 10, 30 and 60 minutes after the initial bolus, electrical stimulations of increasing intensity were applied over the lateral plantar digital nerve, and administered both as single and as repeated stimuli. The resulting reflex responses were recorded using electromyography. Data were analysed using a multivariable linear regression model and a Kruskal Wallis test for single stimulation data, and repeated measures anova and paired t-test for repeated stimulation data. Results The AUC for the stimulus-response curves after single stimulation were similar for both treatments at time 0. At times 10, 30 and 60 the AUCs for the stimulus-response curves were significantly lower with dexmedetomidine treatment than with placebo. Temporal summation was evident in both treatments at times 0, 10, 30 and 60 starting from a stimulation intensity of 10 mA. The magnitude of temporal summation was smaller in dexmedetomidine than in placebo treated dogs at time 10, 30 and 60, but not at time 0. Conclusions During isoflurane anaesthesia, low dose dexmedetomidine suppresses the nociceptive reflex responses after single and repeated electrical stimulation. Clinical relevance This experimental study confirms previous reports on its peri-operative efficacy under clinical conditions, and further indicates that dexmedetomidine might reduce the risk of post-operative chronic pain development.

Dental erosive wear assessment among adolescents and adults utilizing the basic erosive wear examination (BEWE) scoring system.

OBJECTIVES To investigate erosive tooth wear and related variables among adolescents and adults in Israel, utilizing the new basic erosive wear examination (BEWE) scoring system, in an attempt to contribute to the ongoing review, evaluation, and further development of an international standardized index. MATERIAL AND METHODS A cross-sectional, descriptive, and analytic survey was conducted among 500 subjects of five age groups. Dental erosion was measured according to the new BEWE scoring system. Independent variables included gender, age, origin, education, employment status, and diet. A backward stepwise linear regression model was applied to identify significantly associated variables. RESULTS Fifty percent of the survey subjects demonstrated erosive tooth wear; among them, 10 % had distinct erosion of over 50 % of the dental surface. Total BEWE score differences by age groups were statistically significant; as the age increased, the mean total BEWE scores increased (p < 0.001). The association between acidic foods and erosion was evident among the younger population (p = 0.038). In a multiple regression model, age (p < 0.001) and diet (p = 0.044) achieved statistical significance as variables associated with dental erosive wear. CONCLUSIONS Our study is one of the first to use the BEWE scoring system in an epidemiological survey among adolescents and adults. It was found that the BEWE index is straightforward, easy to conduct, and comfortably accepted by the examinees. CLINICAL RELEVANCE The present findings, together with further international research, should contribute toward continued evaluation of the BEWE system as an international standard and thereby, toward more optimal understanding, evidence-based treatment, and prevention of dental erosive wear.

New anthropometry-based age- and sex-specific reference values for urinary 24-hour creatinine excretion based on the adult Swiss population.

BACKGROUND Urinary creatinine excretion is used as a marker of completeness of timed urine collections, which are a keystone of several metabolic evaluations in clinical investigations and epidemiological surveys. The current reference values for 24-hour urinary creatinine excretion rely on observations performed in the 1960s and 1970s in relatively small and mostly selected groups, and may thus poorly fit to the present-day general European population. The aim of this study was to establish and validate anthropometry-based age- and sex-specific reference values of the 24-hour urinary creatinine excretion on adult populations with preserved renal function. METHODS We used data from two independent Swiss cross-sectional population-based studies with standardised 24-hour urinary collection and measured anthropometric variables. Only data from adults of European descent, with estimated glomerular filtration rate (eGFR) ≥60 ml/min/1.73 m(2) and reported completeness of the urinary collection were retained. A linear regression model was developed to predict centiles of the 24-hour urinary creatinine excretion in 1,137 participants from the Swiss Survey on Salt and validated in 994 participants from the Swiss Kidney Project on Genes in Hypertension. RESULTS The mean urinary creatinine excretion was 193 ± 41 μmol/kg/24 hours in men and 151 ± 38 μmol/kg/24 hours in women in the Swiss Survey on Salt. The values were inversely correlated with age and body mass index (BMI). Based on current reference values (177 to 221 μmol/kg/24 hours in men and 133 to 177 μmol/kg/24 hours in women), 56% of the urinary collections in the whole population and 67% in people >60 years old would have been considered as inaccurate. A linear regression model with sex, BMI and age as predictor variables was found to provide the best prediction of the observed values and showed a good fit when applied to the validation population. CONCLUSIONS We propose a validated prediction equation for 24-hour urinary creatinine excretion in the general European population, based on readily available variables such as age, sex and BMI, and a few derived normograms to ease its clinical application. This should help healthcare providers to interpret the completeness of a 24-hour urine collection in daily clinical practice and in epidemiological population studies.

Retrievability of implant-supported crowns when using three different cements: a controlled clinical trial

PURPOSE The purpose of this study was to analyze the removal of implant-supported crowns retained by three different cements using an air-accelerated crown remover and to evaluate the patients' response to the procedure. MATERIALS AND METHODS This controlled clinical trial was conducted with 21 patients (10 women, 11 men; mean age: 51 ± 10.2 years) who had received a total of 74 implants (all placed in the posterior zone of the mandible). Four months after implant surgery, the crowns were cemented on standard titanium abutments of different heights. Three different cements (two temporary: Harvard TEMP and Improv; and one definitive: Durelon) were used and randomly assigned to the patients. Eight months later, one blinded investigator removed all crowns. The number of activations of the instrument (CORONAflex, KaVo) required for crown removal was recorded. The patients completed a questionnaire retrospectively to determine the impact of the procedure and to gauge their subjective perception. A linear regression model and descriptive statistics were used for data analysis. RESULTS All crowns could be retrieved without any technical complications or damage. Both abutment height (P = .019) and cement type (P = .004) had a significant effect on the number of activations, but the type of cement was more important. An increased total number of activations had no or only a weak correlation to the patients' perception of concussion, noise, pain, and unwillingness to use the device. CONCLUSIONS Cemented implant crowns can be removed, and the application of an air-accelerated device is a practicable method. A type of cement with appropriate retention force has to be selected. The impact on the patients' subjective perception should be taken into account.