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.

Knowledge and valuation of Andean agroforestry species: the role of sex, age, and migration among members of a rural community in Bolivia

Background Agroforestry is a sustainable land use method with a long tradition in the Bolivian Andes. A better understanding of people’s knowledge and valuation of woody species can help to adjust actor-oriented agroforestry systems. In this case study, carried out in a peasant community of the Bolivian Andes, we aimed at calculating the cultural importance of selected agroforestry species, and at analysing the intracultural variation in the cultural importance and knowledge of plants according to peasants’ sex, age, and migration. Methods Data collection was based on semi-structured interviews and freelisting exercises. Two ethnobotanical indices (Composite Salience, Cultural Importance) were used for calculating the cultural importance of plants. Intracultural variation in the cultural importance and knowledge of plants was detected by using linear and generalised linear (mixed) models. Results and discussion The culturally most important woody species were mainly trees and exotic species (e.g. Schinus molle, Prosopis laevigata, Eucalyptus globulus). We found that knowledge and valuation of plants increased with age but that they were lower for migrants; sex, by contrast, played a minor role. The age effects possibly result from decreasing ecological apparency of valuable native species, and their substitution by exotic marketable trees, loss of traditional plant uses or the use of other materials (e.g. plastic) instead of wood. Decreasing dedication to traditional farming may have led to successive abandonment of traditional tool uses, and the overall transformation of woody plant use is possibly related to diminishing medicinal knowledge. Conclusions Age and migration affect how people value woody species and what they know about their uses. For this reason, we recommend paying particular attention to the potential of native species, which could open promising perspectives especially for the young migrating peasant generation and draw their interest in agroforestry. These native species should be ecologically sound and selected on their potential to provide subsistence and promising commercial uses. In addition to offering socio-economic and environmental services, agroforestry initiatives using native trees and shrubs can play a crucial role in recovering elements of the lost ancient landscape that still forms part of local people’s collective identity.

Factors that trigger emergency physicians to contact a poison centre: findings from a Swiss study.

OBJECTIVE Poison centres offer rapid and comprehensive support for emergency physicians managing poisoned patients. This study investigates institutional, case-specific and poisoning-specific factors which influence the decision of emergency physicians to contact a poison centre. METHODS Retrospective, consecutive review of all poisoning-related admissions to the emergency departments (EDs) of a primary care hospital and a university hospital-based tertiary referral centre during 2007. Corresponding poison centre consultations were extracted from the poison centre database. Data were matched and analysed by logistic regression and generalised linear mixed models. RESULTS 545 poisonings were treated in the participating EDs (350 (64.2%) in the tertiary care centre, 195 (35.8%) in the primary care hospital). The poison centre was consulted in 62 (11.4%) cases (38 (61.3%) by the tertiary care centre and 24 (38.7%) by the primary care hospital). Factors significantly associated with poison centre consultation included gender (female vs male) (OR 2.99; 95% CI 1.69 to 5.29; p<0.001), number of ingested substances (>1 vs 1) (OR 2.84; 95% CI 1.65 to 4.9; p<0.001) and situation (accidental vs intentional) (OR 2.76; 95% CI 1.05 to 7.25; p=0.039). In contrast, age, medical history and hospital size did not influence poison centre consultation. Poison centre consultation was significantly higher during the week, and significantly less during night shifts. The poison centre was consulted significantly more when patients were admitted to intensive care units (OR 5.81; 95% CI 3.25 to 10.37; p<0.001). Asymptomatic and severe versus mild cases were associated with more frequent consultation (OR 4.48; 95% CI 1.78 to 11.26; p=0.001 and OR 2.76; 95% CI 1.42 to 5.38; p=0.003). CONCLUSIONS We found low rates of poison centre consultation by emergency physicians. It appears that intensive care unit admission and other factors reflecting either complexity or uncertainty of the clinical situation are the strongest predictors for poison centre consultation. Hospital size did not influence referral behaviour.

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.