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.

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.

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

rrreg fits a linear probability model for randomized response data

Counterbalancing effects of competition for resources and facilitation against grazing in alpine snowbed communities

Alpine snowbeds are habitats where the major limiting factors for plant growth are herbivory and a small time window for growth due to late snowmelt. Despite these limitations, snowbed vegetation usually forms a dense carpet of palatable plants due to favourable abiotic conditions for plant growth within the short growing season. These environmental characteristics make snowbeds particularly interesting to study the interplay of facilitation and competition. We hypothesised an interplay between resource competition and facilitation against herbivory. Further, we investigated whether these predicted neighbour effects were species-specific and/or dependent on ontogeny, and whether the balance of positive and negative plant–plant interactions shifted along a snowmelt gradient. We determined the neighbour effects by means of neighbour removal experiments along the snowmelt gradient, and linear mixed model analyses. The results showed that the effects of neighbour removal were weak but generally consistent among species and snowmelt dates, and depended on whether biomass production or survival was considered. Higher total biomass and increased fruiting in removal plots indicated that plants competed for nutrients, water, and light, thereby supporting the hypothesis of prevailing competition for resources in snowbeds. However, the presence of neighbours reduced herbivory and thereby also facilitated survival. For plant growth the facilitative effects against herbivores in snowbeds counterbalanced competition for resources, leading to a weak negative net effect. Overall the neighbour effects were not species-specific and did not change with snowmelt date. Our finding of counterbalancing effects of competition and facilitation within a plant community is of special theoretical value for species distribution models and can explain the success of models that give primary importance to abiotic factors and tend to overlook interrelations between biotic and abiotic effects on plants.

Coping personality type and environmental enrichment affect aggression at weaning in pigs

This study investigated the effects of different environmental treatments and personality types on aggression at mixing of newly weaned domestic piglets. From birth to weaning, 16 litters were housed with their dams in either barren (B) or larger, substrate-enriched (E) environments. At 15 days old, piglets were classified as 'high' (HR) or low resistant' (LR) in a manual restraint test (backtest), which is thought to identify proactive (HR) and reactive (LR) stress coping strategies that may reflect different personality types. At 30 days old, 128 piglets were weaned, relocated and mixed into 32 pens comprising two HR and two LR unfamiliar pigs, balanced for sex and weaning weight. Eight B and eight E groups changed environmental condition whereas the others remained in the same type of environment. Number and duration of fights. fight outcomes and unilateral fighting were scored for 5 h post-mixing and skin lesions were counted before and 5 h, 1 day and 2 days after mixing. On the day following weaning, fighting and also exploratory and oral manipulative behaviours were measured for 6 h. Generalized Linear Mixed Model analyses suggested interactions between pre-weaning environment, post-weaning environment and personality type. Overall, pre-weaning E pigs had longer fights at weaning and mixing (P=0.01) and fought for longer on the next day (P=0.02) than pre-weaning B pigs, and inflicted more skin lesions (P=0.02). Post-weaning enrichment did not affect fighting at mixing but reduced the time spent fighting the next day (P=0.03). Personality had subtle and environment-dependent effects on fighting, and influenced the "structure" rather than the amount of aggressive behaviour. HR pigs, for instance, bullied (i.e. chased surrendering pigs) more often (P=0.009) and their fighting behaviour was less affected by their relative body weight than that of LR pigs. Post-weaning E pigs showed relatively higher levels of exploratory behaviour (P=0.02) and less oral manipulative behaviour (P=0.04) than post-weaning B pigs. In particular, switching from a good quality environment (E) to a worse quality one (B) at weaning decreased exploratory behaviour on the next day, especially for LR pigs, who also tended to fight with and orally manipulate their pen mates more in that condition, and seemed to be more affected by a deterioration of the environment. Overall, pre-weaning enrichment increased aggression after weaning whereas post-weaning enrichment reduced it, and personality type related to some aspects of fighting behaviour. (C) 2011 Elsevier B.V. All rights reserved.

Clinical and molecular characterization of the potential CF disease modifier syntaxin 1A

Cystic fibrosis (CF) is caused by mutations in the CF transmembrane conductance regulator gene (CFTR). Disease severity in CF varies greatly, and sibling studies strongly indicate that genes other than CFTR modify disease outcome. Syntaxin 1A (STX1A) has been reported as a negative regulator of CFTR and other ion channels. We hypothesized that STX1A variants act as a CF modifier by influencing the remaining function of mutated CFTR. We identified STX1A variants by genomic resequencing patients from the Bernese CF Patient Data Registry and applied linear mixed model analysis to establish genotype-phenotype correlations, revealing STX1A rs4363087 (c.467-38A>G) to significantly influence lung function. The same STX1A risk allele was recognized in the European CF Twin and Sibling Study (P=0.0027), demonstrating that the genotype-phenotype association of STX1A to CF disease severity is robust enough to allow replication in two independent CF populations. rs4363087 is in linkage disequilibrium to the exonic variant rs2228607 (c.204C>T). Considering that neither rs4363087 nor rs2228607 changes the amino-acid sequence of STX1A, we investigated their effects on mRNA level. We show that rs2228607 reinforces aberrant splicing of STX1A mRNA, leading to nonsense-mediated mRNA decay. In conclusion, we demonstrate the clinical relevance of STX1A variants in CF, and evidence the functional relevance of STX1A variant rs2228607 at molecular level. Our findings show that genes interacting with CFTR can modify CF disease progression.European Journal of Human Genetics advance online publication, 10 April 2013; doi:10.1038/ejhg.2013.57.

Long-term gas exchange characteristics as marker of deterioration in patients with cystic fibrosis

Background and Aim In patients with cystic fibrosis (CF) the architecture of the developing lungs and the ventilation of lung units are progressively affected, influencing intrapulmonary gas mixing and gas exchange. We examined the long-term course of blood gas measurements in relation to characteristics of lung function and the influence of different CFTR genotype upon this process. Methods Serial annual measurements of PaO2 and PaCO2 assessed in relation to lung function, providing functional residual capacity (FRCpleth), lung clearance index (LCI), trapped gas (VTG), airway resistance (sReff), and forced expiratory indices (FEV1, FEF50), were collected in 178 children (88 males; 90 females) with CF, over an age range of 5 to 18 years. Linear mixed model analysis and binary logistic regression analysis were used to define predominant lung function parameters influencing oxygenation and carbon dioxide elimination. Results PaO2 decreased linearly from age 5 to 18 years, and was mainly associated with FRCpleth, (p < 0.0001), FEV1 (p < 0.001), FEF50 (p < 0.002), and LCI (p < 0.002), indicating that oxygenation was associated with the degree of pulmonary hyperinflation, ventilation inhomogeneities and impeded airway function. PaCO2 showed a transitory phase of low PaCO2 values, mainly during the age range of 5 to 12 years. Both PaO2 and PaCO2 presented with different progression slopes within specific CFTR genotypes. Conclusion In the long-term evaluation of gas exchange characteristics, an association with different lung function patterns was found and was closely related to specific genotypes. Early examination of blood gases may reveal hypocarbia, presumably reflecting compensatory mechanisms to improve oxygenation.

Only Neochordoplasty Achieves Physiological Trans-valvular Pressure Gradients After Repair Of Acute Posterior Leaflet Prolapse In Porcine Mitral Valves

Objective: Minimizing resection and preserving leaflet tissue has been previously shown to be beneficial for mitral valve function and leaflet kinematics after repair of acute posterior leaflet prolapse in porcine valves. We examined the effects of different additional methods of mitral valve repair (neochordoplasty, ring annuloplasty, edge-to-edge repair and triangular resection) on hemodynamics at different heart rates in an experimental model. Methods: Severe acute P2 prolapse was created in eight porcine mitral valves by resecting the posterior marginal chordae. Valve hemodynamics was quantified under pulsatile conditions in an in vitro heart simulator before and after surgical manipulation. Mitral regurgitation was corrected using four different methods of repair on the same valve: neochordoplasty with expanded polytetrafluoroethylene sutures alone and together with ring annuloplasty, edge-to-edge repair and triangular resection, both with non-restrictive annuloplasty. Residual mitral valve leak, trans-valvular pressure gradients, flow and cardiac output were measured at 60 and 80 beats/min. A validated statistical linear mixed model was used to analyze the effect of treatment. The p values were calculated using a two-sided Wald test. Results: Only neochordoplasty with expanded polytetrafluoroethylene sutures but without ring annuloplasty achieved similar hemodynamics compared to those of the native mitral valve (p range 0.071-0.901). Trans-valvular diastolic pressure gradients were within a physiologic range but significantly higher than those of the native valve following neochordoplasty with ring annuloplasty (p=0.000), triangular resection (p=0.000) and edge-to-edge repair (p=0.000). Neochordoplasty alone was significantly better in terms of hemodynamic than neochordoplasty with a ring annuloplasty (p=0.000). These values were stable regardless of heart rate or ring size. Conclusions: Neochordoplasty without ring annuloplasty is the only repair technique able to achieve almost native physiological hemodynamics after correction of leaflet prolapse in a porcine experimental model of acute chordal rupture.

Effect of an intervertebral disk spacer on stiffness after monocortical screw/polymethylmethacrylate fixation in simulated and cadaveric canine cervical vertebral columns.

OBJECTIVE To determine the biomechanical effect of an intervertebral spacer on construct stiffness in a PVC model and cadaveric canine cervical vertebral columns stabilized with monocortical screws/polymethylmethacrylate (PMMA). STUDY DESIGN Biomechanical study. SAMPLE POPULATION PVC pipe; cadaveric canine vertebral columns. METHODS PVC model-PVC pipe was used to create a gap model mimicking vertebral endplate orientation and disk space width of large-breed canine cervical vertebrae; 6 models had a 4-mm gap with no spacer (PVC group 1); 6 had a PVC pipe ring spacer filling the gap (PCV group 2). Animals-large breed cadaveric canine cervical vertebral columns (C2-C7) from skeletally mature dogs without (cadaveric group 1, n = 6, historical data) and with an intervertebral disk spacer (cadaveric group 2, n = 6) were used. All PVC models and cadaver specimens were instrumented with monocortical titanium screws/PMMA. Stiffness of the 2 PVC groups was compared in extension, flexion, and lateral bending using non-destructive 4-point bend testing. Stiffness testing in all 3 directions was performed of the unaltered C4-C5 vertebral motion unit in cadaveric spines and repeated after placement of an intervertebral cortical allograft ring and instrumentation. Data were compared using a linear mixed model approach that also incorporated data from previously tested spines with the same screw/PMMA construct but without disk spacer (cadaveric group 1). RESULTS Addition of a spacer increased construct stiffness in both the PVC model (P < .001) and cadaveric vertebral columns (P < .001) compared to fixation without a spacer. CONCLUSIONS Addition of an intervertebral spacer significantly increased construct stiffness of monocortical screw/PMMA fixation.

A continuous-time MILP model for short-term scheduling of make-and-pack production processes

In process industries, make-and-pack production is used to produce food and beverages, chemicals, and metal products, among others. This type of production process allows the fabrication of a wide range of products in relatively small amounts using the same equipment. In this article, we consider a real-world production process (cf. Honkomp et al. 2000. The curse of reality – why process scheduling optimization problems are diffcult in practice. Computers & Chemical Engineering, 24, 323–328.) comprising sequence-dependent changeover times, multipurpose storage units with limited capacities, quarantine times, batch splitting, partial equipment connectivity, and transfer times. The planning problem consists of computing a production schedule such that a given demand of packed products is fulfilled, all technological constraints are satisfied, and the production makespan is minimised. None of the models in the literature covers all of the technological constraints that occur in such make-and-pack production processes. To close this gap, we develop an efficient mixed-integer linear programming model that is based on a continuous time domain and general-precedence variables. We propose novel types of symmetry-breaking constraints and a preprocessing procedure to improve the model performance. In an experimental analysis, we show that small- and moderate-sized instances can be solved to optimality within short CPU times.

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.

Modeling of red blood cell life-spans in hematologically normal populations

Despite the impact of red blood cell (RBC) Life-spans in some disease areas such as diabetes or anemia of chronic kidney disease, there is no consensus on how to quantitatively best describe the process. Several models have been proposed to explain the elimination process of RBCs: random destruction process, homogeneous life-span model, or a series of 4-transit compartment model. The aim of this work was to explore the different models that have been proposed in literature, and modifications to those. The impact of choosing the right model on future outcomes prediction--in the above mentioned areas--was also investigated. Both data from indirect (clinical data) and direct life-span measurement (biotin-labeled data) methods were analyzed using non-linear mixed effects models. Analysis showed that: (1) predictions from non-steady state data will depend on the RBC model chosen; (2) the transit compartment model, which considers variation in life-span in the RBC population, better describes RBC survival data than the random destruction or homogenous life-span models; and (3) the additional incorporation of random destruction patterns, although improving the description of the RBC survival data, does not appear to provide a marked improvement when describing clinical data.