Generalized linear and generalized additive models in studies of species distributions: setting the scene

An important statistical development of the last 30 years has been the advance in regression analysis provided by generalized linear models (GLMs) and generalized additive models (GAMs). Here we introduce a series of papers prepared within the framework of an international workshop entitled: Advances in GLMs/GAMs modeling: from species distribution to environmental management, held in Riederalp, Switzerland, 6-11 August 2001.We first discuss some general uses of statistical models in ecology, as well as provide a short review of several key examples of the use of GLMs and GAMs in ecological modeling efforts. We next present an overview of GLMs and GAMs, and discuss some of their related statistics used for predictor selection, model diagnostics, and evaluation. Included is a discussion of several new approaches applicable to GLMs and GAMs, such as ridge regression, an alternative to stepwise selection of predictors, and methods for the identification of interactions by a combined use of regression trees and several other approaches. We close with an overview of the papers and how we feel they advance our understanding of their application to ecological modeling.

Fate of linear and supercoiled multinucleosomic templates during transcription.

Electron microscopy was used to monitor the fate of reconstituted nucleosome cores during in vitro transcription of long linear and supercoiled multinucleosomic templates by the prokaryotic T7 RNA polymerase and the eukaryotic RNA polymerase II. Transcription by T7 RNA polymerase disrupted the nucleosomal configuration in the transcribed region, while nucleosomes were preserved upstream of the transcription initiation site and in front of the polymerase. Nucleosome disruption was independent of the topology of the template, linear or supercoiled, and of the presence or absence of nucleosome positioning sequences in the transcribed region. In contrast, the nucleosomal configuration was preserved during transcription from the vitellogenin B1 promoter with RNA polymerase II in a rat liver total nuclear extract. However, the persistence of nucleosomes on the template was not RNA polymerase II-specific, but was dependent on another activity present in the nuclear extract. This was demonstrated by addition of the extract to the T7 RNA polymerase transcription reaction, which resulted in retention of the nucleosomal configuration. This nuclear activity, also found in HeLa cell nuclei, is heat sensitive and could not be substituted by nucleoplasmin, chromatin assembly factor (CAF-I) or a combination thereof. Altogether, these results identify a novel nuclear activity, called herein transcription-dependent chromatin stabilizing activity I or TCSA-I, which may be involved in a nucleosome transfer mechanism during transcription.

Concomitant chemoradiotherapy versus acceleration of radiotherapy with or without concomitant chemotherapy in locally advanced head and neck carcinoma (GORTEC 99-02): an open-label phase 3 randomised trial.

BACKGROUND: Concomitant chemoradiotherapy and accelerated radiotherapy independently improve outcomes for patients with locally advanced head and neck squamous-cell carcinoma (HNSCC). We aimed to assess the efficacy and safety of a combination of these approaches. METHODS: In our open-label phase 3 randomised trial, we enrolled patients with locally advanced, stage III and IV (non-metastatic) HNSCC and an Eastern Cooperative Oncology Group performance status of 0-2. We randomly allocated patients centrally with a computer program (with centre, T stage, N stage, and localisation as minimisation factors) in a 1:1:1 ratio to receive conventional chemoradiotherapy (70 Gy in 7 weeks plus three cycles of 4 days' concomitant carboplatin-fluorouracil), accelerated radiotherapy-chemotherapy (70 Gy in 6 weeks plus two cycles of 5 days' concomitant carboplatin-fluorouracil), or very accelerated radiotherapy alone (64·8 Gy [1·8 Gy twice daily] in 3·5 weeks). The primary endpoint, progression-free survival (PFS), was assessed in all enrolled patients. This trial is completed. The trial is registered with ClinicalTrials.gov, number NCT00828386. FINDINGS: Between Feb 29, 2000, and May 9, 2007, we randomly allocated 279 patients to receive conventional chemoradiotherapy, 280 to accelerated radiotherapy-chemotherapy, and 281 to very accelerated radiotherapy. Median follow-up was 5·2 years (IQR 4·9-6·2); rates of chemotherapy and radiotherapy compliance were good in all groups. Accelerated radiotherapy-chemotherapy offered no PFS benefit compared with conventional chemoradiotherapy (HR 1·02, 95% CI 0·84-1·23; p=0·88) or very accelerated radiotherapy (0·83, 0·69-1·01; p=0·060); conventional chemoradiotherapy improved PFS compared with very accelerated radiotherapy (0·82, 0·67-0·99; p=0·041). 3-year PFS was 37·6% (95% CI 32·1-43·4) after conventional chemoradiotherapy, 34·1% (28·7-39·8) after accelerated radiotherapy-chemotherapy, and 32·2% (27·0-37·9) after very accelerated radiotherapy. More patients in the very accelerated radiotherapy group had RTOG grade 3-4 acute mucosal toxicity (226 [84%] of 268 patients) compared with accelerated radiotherapy-chemotherapy (205 [76%] of 271 patients) or conventional chemoradiotherapy (180 [69%] of 262; p=0·0001). 158 (60%) of 265 patients in the conventional chemoradiotherapy group, 176 (64%) of 276 patients in the accelerated radiotherapy-chemotherapy group, and 190 (70%) of 272 patients in the very accelerated radiotherapy group were intubated with feeding tubes during treatment (p=0·045). INTERPRETATION: Chemotherapy has a substantial treatment effect given concomitantly with radiotherapy and acceleration of radiotherapy cannot compensate for the absence of chemotherapy. We noted the most favourable outcomes for conventional chemoradiotherapy, suggesting that acceleration of radiotherapy is probably not beneficial in concomitant chemoradiotherapy schedules. FUNDING: French Ministry of Health.

Neonatal steroids induce a down-regulation of tenascin-C and elastin and cause a deceleration of the first phase and an acceleration of the second phase of lung alveolarization.

Pre- and postnatal corticosteroids are often used in perinatal medicine to improve pulmonary function in preterm infants. To mimic this clinical situation, newborn rats were treated systemically with dexamethasone (Dex), 0.1-0.01 mg/kg/day on days P1-P4. We hypothesized that postnatal Dex may have an impact on alveolarization by interfering with extracellular matrix proteins and cellular differentiation. Morphological alterations were observed on 3D images obtained by high-resolution synchrotron radiation X-ray tomographic microscopy. Alveolarization was quantified stereologically by estimating the formation of new septa between days P4 and P60. The parenchymal expression of tenascin-C (TNC), smooth muscle actin (SMA), and elastin was measured by immunofluorescence and gene expression for TNC by qRT-PCR. After Dex treatment, the first phase of alveolarization was significantly delayed between days P6 and P10, whereas the second phase was accelerated. Elastin and SMA expressions were delayed by Dex treatment, whereas TNC expression was delayed and prolonged. A short course of neonatal steroids impairs the first phase of alveolarization, most likely by altering the TNC and elastin expression. Due to an overshooting catch-up during the second phase of alveolarization, the differences disappear when the animals reach adulthood.

Drug-induced linear IgA bullous dermatosis probably induced by furosemide.

Linear IgA bullous dermatosis (LABD) is an autoimmune disease, characterized by linear deposition of IgA along the basement membrane zone. Drug-induced LABD is rare but increasing in frequency. A new case of drug-induced LABD associated with the administration of furosemide is described.

Resting-State Functional Connectivity Emerges from Structurally and Dynamically Shaped Slow Linear Fluctuations.

Brain fluctuations at rest are not random but are structured in spatial patterns of correlated activity across different brain areas. The question of how resting-state functional connectivity (FC) emerges from the brain's anatomical connections has motivated several experimental and computational studies to understand structure-function relationships. However, the mechanistic origin of resting state is obscured by large-scale models' complexity, and a close structure-function relation is still an open problem. Thus, a realistic but simple enough description of relevant brain dynamics is needed. Here, we derived a dynamic mean field model that consistently summarizes the realistic dynamics of a detailed spiking and conductance-based synaptic large-scale network, in which connectivity is constrained by diffusion imaging data from human subjects. The dynamic mean field approximates the ensemble dynamics, whose temporal evolution is dominated by the longest time scale of the system. With this reduction, we demonstrated that FC emerges as structured linear fluctuations around a stable low firing activity state close to destabilization. Moreover, the model can be further and crucially simplified into a set of motion equations for statistical moments, providing a direct analytical link between anatomical structure, neural network dynamics, and FC. Our study suggests that FC arises from noise propagation and dynamical slowing down of fluctuations in an anatomically constrained dynamical system. Altogether, the reduction from spiking models to statistical moments presented here provides a new framework to explicitly understand the building up of FC through neuronal dynamics underpinned by anatomical connections and to drive hypotheses in task-evoked studies and for clinical applications.

Linear relations between writhe and minimal crossing number in Conway families of ideal knots and links

We showed earlier how to predict the writhe of any rational knot or link in its ideal geometric configuration, or equivalently the average of the 3D writhe over statistical ensembles of random configurations of a given knot or link (Cerf and Stasiak 2000 Proc. Natl Acad. Sci. USA 97 3795). There is no general relation between the minimal crossing number of a knot and the writhe of its ideal geometric configuration. However, within individual families of knots linear relations between minimal crossing number and writhe were observed (Katritch et al 1996 Nature 384 142). Here we present a method that allows us to express the writhe as a linear function of the minimal crossing number within Conway families of knots and links in their ideal configuration. The slope of the lines and the shift between any two lines with the same

Linear random knots and their scaling behavior

We present here a nonbiased probabilistic method that allows us to consistently analyze knottedness of linear random walks with up to several hundred noncorrelated steps. The method consists of analyzing the spectrum of knots formed by multiple closures of the same open walk through random points on a sphere enclosing the walk. Knottedness of individual "frozen" configurations of linear chains is therefore defined by a characteristic spectrum of realizable knots. We show that in the great majority of cases this method clearly defines the dominant knot type of a walk, i.e., the strongest component of the spectrum. In such cases, direct end-to-end closure creates a knot that usually coincides with the knot type that dominates the random closure spectrum. Interestingly, in a very small proportion of linear random walks, the knot type is not clearly defined. Such walks can be considered as residing in a border zone of the configuration space of two or more knot types. We also characterize the scaling behavior of linear random knots.

Comparison between a linear ion trap and a triple quadruple MS in the sensitive detection of large peptides at femtomole amounts on column.

In addition to the importance of sample preparation and extract separation, MS detection is a key factor in the sensitive quantification of large undigested peptides. In this article, a linear ion trap MS (LIT-MS) and a triple quadrupole MS (TQ-MS) have been compared in the detection of large peptides at subnanomolar concentrations. Natural brain natriuretic peptide, C-peptide, substance P and D-Junk-inhibitor peptide, a full D-amino acid therapeutic peptide, were chosen. They were detected by ESI and simultaneous MS(1) and MS(2) acquisitions. With direct peptide infusion, MS(2) spectra revealed that fragmentation was peptide dependent, milder on the LIT-MS and required high collision energies on the TQ-MS to obtain high-intensity product ions. Peptide adsorption on surfaces was overcome and peptide dilutions ranging from 0.1 to 25 nM were injected onto an ultra high-pressure LC system with a 1 mm id analytical column and coupled with the MS instruments. No difference was observed between the two instruments when recording in LC-MS(1) acquisitions. However, in LC-MS(2) acquisitions, a better sensitivity in the detection of large peptides was observed with the LIT-MS. Indeed, with the three longer peptides, the typical fragmentation in the TQ-MS resulted in a dramatic loss of sensitivity (> or = 10x).

A general linear relaxometry model of R1 using imaging data.

PURPOSE: The longitudinal relaxation rate (R1 ) measured in vivo depends on the local microstructural properties of the tissue, such as macromolecular, iron, and water content. Here, we use whole brain multiparametric in vivo data and a general linear relaxometry model to describe the dependence of R1 on these components. We explore a) the validity of having a single fixed set of model coefficients for the whole brain and b) the stability of the model coefficients in a large cohort. METHODS: Maps of magnetization transfer (MT) and effective transverse relaxation rate (R2 *) were used as surrogates for macromolecular and iron content, respectively. Spatial variations in these parameters reflected variations in underlying tissue microstructure. A linear model was applied to the whole brain, including gray/white matter and deep brain structures, to determine the global model coefficients. Synthetic R1 values were then calculated using these coefficients and compared with the measured R1 maps. RESULTS: The model's validity was demonstrated by correspondence between the synthetic and measured R1 values and by high stability of the model coefficients across a large cohort. CONCLUSION: A single set of global coefficients can be used to relate R1 , MT, and R2 * across the whole brain. Our population study demonstrates the robustness and stability of the model. Magn Reson Med, 2014. © 2014 The Authors. Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. Magn Reson Med 73:1309-1314, 2015. © 2014 Wiley Periodicals, Inc.

A non-linear Bayesian approach for upscaling local-scale hydraulic conductivity measurements based on local- and regional-scale geophysical data

Significant progress has been made with regard to the quantitative integration of geophysical and hydrological data at the local scale. However, extending the corresponding approaches to the regional scale represents a major, and as-of-yet largely unresolved, challenge. To address this problem, we have developed an upscaling procedure based on a Bayesian sequential simulation approach. This method is then applied to the stochastic integration of low-resolution, regional-scale electrical resistivity tomography (ERT) data in combination with high-resolution, local-scale downhole measurements of the hydraulic and electrical conductivities. Finally, the overall viability of this upscaling approach is tested and verified by performing and comparing flow and transport simulation through the original and the upscaled hydraulic conductivity fields. Our results indicate that the proposed procedure does indeed allow for obtaining remarkably faithful estimates of the regional-scale hydraulic conductivity structure and correspondingly reliable predictions of the transport characteristics over relatively long distances.

Improvement of walking speed prediction by accelerometry and altimetry, validated by satellite positioning.

Activity monitors based on accelerometry are used to predict the speed and energy cost of walking at 0% slope, but not at other inclinations. Parallel measurements of body accelerations and altitude variation were studied to determine whether walking speed prediction could be improved. Fourteen subjects walked twice along a 1.3 km circuit with substantial slope variations (-17% to +17%). The parameters recorded were body acceleration using a uni-axial accelerometer, altitude variation using differential barometry, and walking speed using satellite positioning (DGPS). Linear regressions were calculated between acceleration and walking speed, and between acceleration/altitude and walking speed. These predictive models, calculated using the data from the first circuit run, were used to predict speed during the second circuit. Finally the predicted velocity was compared with the measured one. The result was that acceleration alone failed to predict speed (mean r = 0.4). Adding altitude variation improved the prediction (mean r = 0.7). With regard to the altitude/acceleration-speed relationship, substantial inter-individual variation was found. It is concluded that accelerometry, combined with altitude measurement, can assess position variations of humans provided inter-individual variation is taken into account. It is also confirmed that DGPS can be used for outdoor walking speed measurements, opening up new perspectives in the field of biomechanics.

Short-term stability of testosterone and epitestosterone conjugates in urine samples: quantification by liquid chromatography-linear ion trap mass spectrometry.

A simple method using liquid chromatography-linear ion trap mass spectrometry for simultaneous determination of testosterone glucuronide (TG), testosterone sulfate (TS), epitestosterone glucuronide (EG) and epitestosterone sulfate (ES) in urine samples was developed. For validation purposes, a urine containing no detectable amount of TG, TS and EG was selected and fortified with steroid conjugate standards. Quantification was performed using deuterated testosterone conjugates to correct for ion suppression/enhancement during ESI. Assay validation was performed in terms of lower limit of detection (1-3ng/mL), recovery (89-101%), intraday precision (2.0-6.8%), interday precision (3.4-9.6%) and accuracy (101-103%). Application of the method to short-term stability testing of urine samples at temperature ranging from 4 to 37 degrees C during a time-storage of a week lead to the conclusion that addition of sodium azide (10mg/mL) is required for preservation of the analytes.