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.

Oxygen uptake kinetics and middle distance swimming performance.

OBJECTIVE: The aim of this study was to determine whether V˙O(2) kinetics and specifically, the time constant of transitions from rest to heavy (τ(p)H) and severe (τ(p)S) exercise intensities, are related to middle distance swimming performance. DESIGN: Fourteen highly trained male swimmers (mean ± SD: 20.5 ± 3.0 yr; 75.4 ± 12.4 kg; 1.80 ± 0.07 m) performed an discontinuous incremental test, as well as square wave transitions for heavy and severe swimming intensities, to determine V˙O(2) kinetics parameters using two exponential functions. METHODS: All the tests involved front-crawl swimming with breath-by-breath analysis using the Aquatrainer swimming snorkel. Endurance performance was recorded as the time taken to complete a 400 m freestyle swim within an official competition (T400), one month from the date of the other tests. RESULTS: T400 (Mean ± SD) (251.4 ± 12.4 s) was significantly correlated with τ(p)H (15.8 ± 4.8s; r=0.62; p=0.02) and τ(p)S (15.8 ± 4.7s; r=0.61; p=0.02). The best single predictor of 400 m freestyle time, out of the variables that were assessed, was the velocity at V˙O(2max)vV˙O(2max), which accounted for 80% of the variation in performance between swimmers. However, τ(p)H and V˙O(2max) were also found to influence the prediction of T400 when they were included in a regression model that involved respiratory parameters only. CONCLUSIONS: Faster kinetics during the primary phase of the V˙O(2) response is associated with better performance during middle-distance swimming. However, vV˙O(2max) appears to be a better predictor of T400.

The value of serum procalcitonin level for differentiation of infectious from noninfectious causes of fever after orthopaedic surgery.

BACKGROUND: Early diagnosis of postoperative orthopaedic infections is important in order to rapidly initiate adequate antimicrobial therapy. There are currently no reliable diagnostic markers to differentiate infectious from noninfectious causes of postoperative fever. We investigated the value of the serum procalcitonin level in febrile patients after orthopaedic surgery. METHODS: We prospectively evaluated 103 consecutive patients with new onset of fever within ten days after orthopaedic surgery. Fever episodes were classified by two independent investigators who were blinded to procalcitonin results as infectious or noninfectious origin. White blood-cell count, C-reactive protein level, and procalcitonin level were assessed on days 0, 1, and 3 of the postoperative fever. RESULTS: Infection was diagnosed in forty-five (44%) of 103 patients and involved the respiratory tract (eighteen patients), urinary tract (eighteen), joints (four), surgical site (two), bloodstream (two), and soft tissues (one). Unlike C-reactive protein levels and white blood-cell counts, procalcitonin values were significantly higher in patients with infection compared with patients without infection on the day of fever onset (p = 0.04), day 1 (p = 0.07), and day 3 (p = 0.003). Receiver-operating characteristics demonstrated that procalcitonin had the highest diagnostic accuracy, with a value of 0.62, 0.62, and 0.71 on days 0, 1, and 3, respectively. In a multivariate logistic regression analysis, procalcitonin was a significant predictor for postoperative infection on days 0, 1, and 3 of fever with an odds ratio of 2.3 (95% confidence interval, 1.1 to 4.4), 2.3 (95% confidence interval, 1.1 to 5.2), and 3.3 (95% confidence interval, 1.2 to 9.0), respectively. CONCLUSIONS: Serum procalcitonin is a helpful diagnostic marker supporting clinical and microbiological findings for more reliable differentiation of infectious from noninfectious causes of fever after orthopaedic surgery.

Adiponectin and C-reactive protein are positively associated with microalbuminuria in an adult Caucasian population

Objective: Microalbuminuria (MAU) is a marker of early kidney injury and cardiovascular risk. We assessed the association of MAU with plasma adiponectin, leptin, and hsCRP as inflammatory marker, accounting for hypertension, diabetes and obesity. Design and Methods: Population based, cross-sectional study in Caucasian subjects aged 35 to 75 years in Lausanne, Switzerland. MAU, measured by quantitative immunonephelometry on spot morning urine, was used either as a continuous (MAU) or dichotomized variable (MA defined as MAU > 2.5 and >3.5 mg/mmol creatinine in men and women, respectively). Results: The 2955 women (age 53.3_10.7, mean_SD years) had mean body mass index (BMI) 24.9_4.5 kg/m. The 2479 men (age 53.1_10.8 years) hadmean BMI 27.0_3.9 kg/m2.Median hsCRP was 1.3 and 1.3 mg/L, median adiponectin 6.2 and 10.6mg/mL in men and women, respectively. MA prevalence was 4.9% in women and 9.8% in men. In multivariate regression analysis adjusting for potential confounders (age, sex, hypertension, diabetes, eGFR, BMI, percent fat mass, insulin and smoking), logtransformed MAU was positively associated with hsCRP (P<0.001) and adiponectin (P¼0.002), but not with leptin. The association of adiponectin with MAU was stronger in subjects with low hsCRP, and vice versa (P interaction<0.001). Conclusion: Adiponectin and hsCRP are significant positive determinants of MAU, independently of diabetes, hypertension and fat mass. A negative interaction between hsCRP and adiponectin was found for their effect on MAU. Whether hyperadiponectinemia represents an adequate protective response to vascular stress or has negative causal impact on the development of MAU should be assessed in further studies.

Condom use with steady partners among heterosexual people living with HIV in Europe: testing the information-motivation-behavioral skills model.

Guided by a modified information-motivation-behavioral skills model, this study identified predictors of condom use among heterosexual people living with HIV with their steady partners. Consecutive patients at 14 European HIV outpatient clinics received an anonymous, standardized, self-administered questionnaire between March and December 2007. Data were analyzed using descriptive statistics and two-step backward elimination regression analyses stratified by gender. The survey included 651 participants (n = 364, 56% women; n = 287, 44%). Mean age was 39 years for women and 43 years for men. Most had acquired HIV sexually and more than half were in a serodiscordant relationship. Sixty-three percent (n = 229) of women and 59% of men (n = 169) reported at least one sexual encounter with a steady partner 6 months prior to the survey. Fifty-one percent (n = 116) of women and 59% of men (n = 99) used condoms consistently with that partner. In both genders, condom use was positively associated with subjective norm conducive to condom use, and self-efficacy to use condoms. Having a partner whose HIV status was positive or unknown reduced condom use. In men, higher education and knowledge about condom use additionally increased condom use, while the use of erectile-enhancing medication decreased it. For women, HIV disclosure to partners additionally reduced the likelihood of condom use. Positive attitudes to condom use and subjective norm increased self-efficacy in both genders, however, a number of gender-related differences appeared to influence self-efficacy. Service providers should pay attention to the identified predictors of condom use and adopt comprehensive and gender-related approaches for preventive interventions with people living with HIV.

Body surface area and body weight predict total liver volume in Western adults.

Computed tomography (CT) is used increasingly to measure liver volume in patients undergoing evaluation for transplantation or resection. This study is designed to determine a formula predicting total liver volume (TLV) based on body surface area (BSA) or body weight in Western adults. TLV was measured in 292 patients from four Western centers. Liver volumes were calculated from helical computed tomographic scans obtained for conditions unrelated to the hepatobiliary system. BSA was calculated based on height and weight. Each center used a different established method of three-dimensional volume reconstruction. Using regression analysis, measurements were compared, and formulas correlating BSA or body weight to TLV were established. A linear regression formula to estimate TLV based on BSA was obtained: TLV = -794.41 + 1,267.28 x BSA (square meters; r(2) = 0.46; P &lt;.0001). A formula based on patient weight also was derived: TLV = 191.80 + 18.51 x weight (kilograms; r(2) = 0.49; P &lt;.0001). The newly derived TLV formula based on BSA was compared with previously reported formulas. The application of a formula obtained from healthy Japanese individuals underestimated TLV. Two formulas derived from autopsy data for Western populations were similar to the newly derived BSA formula, with a slight overestimation of TLV. In conclusion, hepatic three-dimensional volume reconstruction based on helical CT predicts TLV based on BSA or body weight. The new formulas derived from this correlation should contribute to the estimation of TLV before liver transplantation or major hepatic resection.

Examining the Agreement of Structural Loss Measured With Heidelberg Retinal Tomograhy and Functional Loss With Standard Automated Perimetry (SAP; Octopus), Pulsar Perimetry (PP), and Moorfields Motion Displacement Test (MDT) Perimetry

Purpose: To examine the relationship of functional measurements with structural measures. Methods: 146 eyes of 83 test subjects underwent Heidelberg Retinal Tomography (HRTIII) (disc area<2.43, mphsd<40), and perimetry testing with Octopus (SAP; Dynamic), Pulsar (PP; TOP) and Moorfields MDT (ESTA). Glaucoma was defined as progressive structural or functional loss (20 eyes). Perimetry test points were grouped into 6 sectors based on the estimated optic nerve head angle into which the associated nerve fiber bundle enters (Garway-Heath map). Perimetry summary measures (PSM) (MD SAP/ MD PP/ PTD MDT) were calculated from the average total deviation of each measured threshold from the normal for each sector. We calculated the 95% significance level of the sectorial PSM from the respective normative data. We calculated the percentage agreement with group1 (G1), healthy on HRT and within normal perimetric limits, and group 2 (G2), abnormal on HRT and outside normal perimetric limits. We also examined the relationship of PSM and rim area (RA) in those sectors classified as abnormal by MRA (Moorfields Regression Analysis) of HRT. Results: The mean age was 65 (range= [37, 89]). The global sensitivity versus specificity of each instrument in detecting glaucomatous eyes was: MDT 80% vs. 88%, SAP 80% vs. 80%, PP 70% vs. 89% and HRT 80% vs. 79%. Highest percentage agreement of HRT (respectively G1, G2, sector) with PSM were MDT (89%, 57%, nasal superior), SAP (83%, 74%, temporal superior), PP (74%, 63%, nasal superior). Globally percentage agreement (respectively G1, G2) was MDT (92%, 28%), SAP (87%, 40%) and PP (77%, 49%). Linear regression showed there was no significant trend globally associating RA and PSM. However, sectorally the supero-nasal sector had a statistically significant (p<0.001) trend with each instrument, the associated r2 coefficients are (MDT 0.38 SAP 0.56 and PP 0.39). Conclusions: There were no significant differences in global sensitivity or specificity between instruments. Structure-function relationships varied significantly between instruments and were consistently strongest supero-nasally. Further studies are required to investigate these relationships in detail.

Effect of age on intraoperative cerebrovascular autoregulation and near-infrared spectroscopy-derived cerebral oxygenation.

Background. Age is an important risk factor for perioperative cerebral complications such as stroke, postoperative cognitive dysfunction, and delirium. We explored the hypothesis that intraoperative cerebrovascular autoregulation is less efficient and brain tissue oxygenation lower in elderly patients, thus, increasing the vulnerability of elderly brains to systemic insults such as hypotension.Methods. We monitored intraoperative cerebral perfusion in 50 patients aged 18-40 and 77 patients >65 yr at two Swiss university hospitals. Mean arterial pressure (MAP) was measured continuously using a plethysmographic method. An index of cerebrovascular autoregulation (Mx) was calculated based on changes in transcranial Doppler flow velocity due to changes in MAP. Cerebral oxygenation was assessed by the tissue oxygenation index (TOI) using near-infrared spectroscopy. End-tidal CO(2), O(2), and sevoflurane concentrations and peripheral oxygen saturation were recorded continuously. Standardized anaesthesia was administered in all patients (thiopental, sevoflurane, fentanyl, atracurium).Results. Autoregulation was less efficient in patients aged >65 yr [by 0.10 (SE 0.04; P=0.020)] in a multivariable linear regression analysis. This difference was not attributable to differences in MAP, end-tidal CO2, or higher doses of sevoflurane. TOI was not significantly associated with age, sevoflurane dose, or Mx but increased with increasing flow velocity [by 0.09 (SE 0.04; P=0.028)] and increasing MAP [by 0.11 (SE 0.05; P=0.043)].Conclusions. Our results do not support the hypothesis that older patients' brains are more vulnerable to systemic insults. The difference of autoregulation between the two groups was small and most likely clinically insignificant.

Role of Chlamydia trachomatis in miscarriage.

To determine the role of Chlamydia trachomatis in miscarriage, we prospectively collected serum, cervicovaginal swab specimens, and placental samples from 386 women with and without miscarriage. Prevalence of immunoglobulin G against C. trachomatis was higher in the miscarriage group than in the control group (15.2% vs. 7.3%; p = 0.018). Association between C. trachomatis-positive serologic results and miscarriage remained significant after adjustment for age, origin, education, and number of sex partners (odds ratio 2.3, 95% confidence interval 1.1-4.9). C. trachomatis DNA was more frequently amplified from products of conception or placenta from women who had a miscarriage (4%) than from controls (0.7%; p = 0.026). Immunohistochemical analysis confirmed C. trachomatis in placenta from 5 of 7 patients with positive PCR results, whereas results of immunohistochemical analysis were negative in placenta samples from all 8 negative controls tested. Associations between miscarriage and serologic/molecular evidence of C. trachomatis infection support its role in miscarriage.

Comparison of neoadjuvant cisplatin-based chemotherapy versus radiochemotherapy followed by resection for stage III (N2) NSCLC.

OBJECTIVE: Comparison of prospectively treated patients with neoadjuvant cisplatin-based chemotherapy vs radiochemotherapy followed by resection for mediastinoscopically proven stage III N2 non-small cell lung cancer with respect to postoperative morbidity, pathological nodal downstaging, overall and disease-free survival, and site of recurrence. METHODS: Eighty-two patients were enrolled between January 1994 to June 2003, 36 had cisplatin and doxetacel-based chemotherapy (group I) and 46 cisplatin-based radiochemotherapy up to 44 Gy (group II), either as sequential (25 patients) or concomitant (21 patients) treatment. All patients had evaluation of absence of distant metastases by bone scintigraphy, thoracoabdominal CT scan or PET scan, and brain MRI, and all underwent pre-induction mediastinoscopy, resection and mediastinal lymph node dissection by the same surgeon. RESULTS: Group I and II comprised T1/2 tumors in 47 and 28%, T3 tumors in 45 and 41%, and T4 tumors in 8 and 31% of the patients, respectively (P=0.03). There was a similar distribution of the extent of resection (lobectomy, sleeve lobectomy, left and right pneumonectomy) in both groups (P=0.9). Group I and II revealed a postoperative 90-d mortality of 3 and 4% (P=0.6), a R0-resection rate of 92 and 94% (P=0.9), and a pathological mediastinal downstaging in 61 and 78% of the patients (P<0.01), respectively. 5y-overall survival and disease-free survival of all patients were 40 and 36%, respectively, without significant difference between T1-3 and T4 tumors. There was no significant difference in overall survival rate in either induction regimens, however, radiochemotherapy was associated with a longer disease-free survival than chemotherapy (P=0.04). There was no significant difference between concurrent vs sequential radiochemotherapy with respect to postoperative morbidity, resectability, pathological nodal downstaging, survival and disease-free survival. CONCLUSIONS: Neoadjuvant cisplatin-based radiochemotherapy was associated with a similar postoperative mortality, an increased pathological nodal downstaging and a better disease-free survival as compared to cisplatin doxetacel-based chemotherapy in patients with stage III (N2) NSCLC although a higher number of T4 tumors were admitted to radiochemotherapy.

Prediction of myocardial infarction size using the SYNTAX score in patients treated with primary percutaneous coronary intervention for acute ST- segment elevation myocardial infarction

Objectives The relevance of the SYNTAX score for the particular case of patients with acute ST- segment elevation myocardial infarction (STEMI) undergoing primary percutaneous coronary intervention (PPCI)  has previously only been studied in the setting of post hoc analysis of large prospective randomized clinical trials. A "real-life" population approach has never been explored before. The aim of this study was to evaluate the impact of the SYNTAX score for the prediction of the myocardial infarction size, estimated by the creatin-kinase (CK) peak value, using the SYNTAX score in patients treated with primary coronary intervention for acute ST-segment elevation myocardial infarction. Methods The primary endpoint of the study was myocardial infarction size as measured by the CK peak value. The SYNTAX score was calculated retrospectively in 253 consecutive patients with acute ST-segment elevation myocardial infarction (STEMI) undergoing primary percutaneous coronary intervention (PPCI) in a large tertiary referral center in Switzerland, between January 2009 and June 2010. Linear regression analysis was performed to compare myocardial infarction size with the SYNTAX score. This same endpoint was then stratified according to SYNTAX score tertiles: low <22 (n=178), intermediate [22-32] (n=60), and high >=33 (n=15). Results There were no significant differences in terms of clinical characteristics between the three groups. When stratified according to the SYNTAX score tertiles, average CK peak values of 1985 (low<22), 3336 (intermediate [22-32]) and 3684 (high>=33) were obtained with a p-value <0.0001. Bartlett's test for equal variances between the three groups was 9.999 (p-value <0.0067). A moderate Pearson product-moment correlation coefficient (r=0.4074) with a high statistical significance level (p-value <0.0001) was found. The coefficient of determination (R^2=0.1660) showed that approximately 17% of the variation of CK peak value (myocardial infarction size) could be explained by the SYNTAX score, i.e. by the coronary disease complexity. Conclusion In an all-comers population, the SYNTAX score is an additional tool in predicting myocardial infarction size in patients treated with primary percutaneous coronary intervention (PPCI). The stratification of patients in different risk groups according to SYNTAX enables to identify a high-risk population that may warrant particular patient care.

Quetiapine affects neuropeptide Y and corticotropin-releasing hormone in cerebrospinal fluid from schizophrenia patients: relationship to depression and anxiety symptoms and to treatment response.

Cumulative evidence indicates that neuropeptides play a role in the pathophysiology of schizophrenia. Early data showed increased neuropeptide Y (NPY) in cerebrospinal fluid (CSF) from schizophrenia patients and data from rodents show that antipsychotic drugs modulate NPY levels in and release from selected rat brain regions. In view of these findings we investigated whether the atypical antipsychotic quetiapine, originally used as an antipsychotic but subsequently shown to be efficient also in major depressive disorder and in both poles of bipolar disorder, would affect NPY-like immunoreactivity (-LI), and corticotropin-releasing hormone (CRH)-LI levels in CSF of schizophrenia patients. NPY-LI and CRH-LI in CSF were determined in 22 patients with schizophrenia. Lumbar puncture was performed at baseline and again after 4 wk of quetiapine treatment (600 mg/d). Patients were assessed with the Positive and Negative Syndrome Scale (PANSS) at baseline and at weekly intervals. Quetiapine treatment was associated with a significant increase in NPY-LI (p<0.001) and decrease in CRH-LI (p<0.01). Stepwise multiple regression analysis revealed that ΔNPY-LI and ΔCRH-LI levels predicted 63% (p<0.001) of the variability of the ΔPANSS total score, ΔNPY-LI 42% (p<0.05) of the ΔPANSS anxiety items (G2) and ΔCRH-LI 40% (p=0.05) of the ΔPANSS depression items (G6). These results suggest that while quetiapine's effects on monoamines are probably related to its antipsychotic properties, the modulation of NPY and CRH accounts for its antidepressant and anxiolytic effects and can be markers of response.

Off-resonance angiography: a new method to depict vessels--phantom and rabbit studies.

PURPOSE: To evaluate the utility of inversion recovery with on-resonant water suppression (IRON) in combination with injection of the long-circulating monocrystalline iron oxide nanoparticle (MION)-47 for contrast material-enhanced magnetic resonance (MR) angiography. MATERIALS AND METhods: Experiments were approved by the institutional animal care committee. Eleven rabbits were imaged at baseline before injection of a contrast agent and then serially 5-30 minutes, 2 hours, 1 day, and 3 days after a single intravenous bolus injection of 80 micromol of MION-47 per kilogram of body weight (n = 6) or 250 micromol/kg MION-47 (n = 5). Conventional T1-weighted MR angiography and IRON MR angiography were performed on a clinical 3.0-T imager. Signal-to-noise and contrast-to-noise ratios were measured in the aorta of rabbits in vivo. Venous blood was obtained from the rabbits before and after MION-47 injection for use in phantom studies. RESULTS: In vitro blood that contained MION-47 appeared signal attenuated on T1-weighted angiograms, while characteristic signal-enhanced dipolar fields were observed on IRON angiograms. In vivo, the vessel lumen was signal attenuated on T1-weighted MR angiograms after MION-47 injection, while IRON supported high intravascular contrast by simultaneously providing positive signal within the vessels and suppressing background tissue (mean contrast-to-noise ratio, 61.9 +/- 12.4 [standard deviation] after injection vs 1.1 +/- 0.4 at baseline, P < .001). Contrast-to-noise ratio was higher on IRON MR angiograms than on conventional T1-weighted MR angiograms (9.0 +/- 2.5, P < .001 vs IRON MR angiography) and persisted up to 24 hours after MION-47 injection (76.2 +/- 15.9, P < .001 vs baseline). CONCLUSION: IRON MR angiography in conjunction with superparamagnetic nanoparticle administration provides high intravascular contrast over a long time and without the need for image subtraction.

Heritability of left ventricular structure and function in Caucasian families.

AIMS: The aim of this study was to investigate the heritability as well as genetic and environmental correlations of left ventricular (LV) structural and functional traits in complex pedigrees of a Caucasian population. METHODS AND RESULTS: We randomly recruited 459 white European subjects from 52 families (50% women; mean age 45 years). LV structure was measured by M-mode and 2D echocardiography and LV function was measured by conventional Doppler and tissue Doppler imaging (TDI). Other measurements included blood pressure, anthropometric, and biochemical measurements. We estimated the heritability of LV traits while adjusting for covariables, including sex, age, body height and weight, systolic and diastolic blood pressures, and heart rate. With full adjustment, heritability of LV mass was 0.23 (P= 0.025). The TDI-derived mitral annular velocities Ea and Aa showed moderate heritability (h(2)= 0.36 and 0.53, respectively), whereas the mitral inflow A peak had weak heritability (h(2) = 0.25) and the E peak was not heritable (h(2) = 0.11). We partitioned the total phenotypic correlation when it reached significance, into a genetic and an environmental component. The genetic correlations were 0.61 between the E and Ea peaks and 0.90 between the A and Aa peaks. CONCLUSION: Our study demonstrated moderate heritability for LV mass as well as the mitral annular Ea and Aa peaks. We also found significant genetic correlations between the E and Ea peaks and between the A and Aa peaks. Our current findings support the ongoing research to map and detect genetic variants that contribute to the variation in LV mass and other LV structural and functional phenotypes.

Nestmate recognition and the genetic relatedness of nests in the ant Formica pratensis

Genetic relatedness of the mound-building ant Formica pratensis was determined by means of microsatellite DNA polymorphism, and its impact on nestmate recognition was tested in a population in Southern Sweden (Oeland). Recognition between nests was measured by testing aggression levels between single pairs of workers. The genetic distances of nests (Nei's genetic distance) and the spatial distance of nests were correlated and both showed a strong relation to the aggression behavior. Multiple regression analysis revealed a stronger impact of genetic relatedness rather than spatial distances on aggression behavior. Neighbouring nests were more closely related than distant nests, which may reflect budding as a possible spreading mechanism. The genetic distance data showed that nestmate recognition was strongly genetically influenced in F. pratensis.