Validation of Logistic Regression Models for Landslide Susceptibility Maps

A wide range of numerical models and tools have been developed over the last decades to support the decision making process in environmental applications, ranging from physical models to a variety of statistically-based methods. In this study, a landslide susceptibility map of a part of Three Gorges Reservoir region of China was produced, employing binary logistic regression analyses. The available information includes the digital elevation model of the region, geological map and different GIS layers including land cover data obtained from satellite imagery. The landslides were observed and documented during the field studies. The validation analysis is exploited to investigate the quality of mapping.

Imputation in data fusion of heterogeneous data sets a model-based numerical experiment

Given the very large amount of data obtained everyday through population surveys, much of the new research again could use this information instead of collecting new samples. Unfortunately, relevant data are often disseminated into different files obtained through different sampling designs. Data fusion is a set of methods used to combine information from different sources into a single dataset. In this article, we are interested in a specific problem: the fusion of two data files, one of which being quite small. We propose a model-based procedure combining a logistic regression with an Expectation-Maximization algorithm. Results show that despite the lack of data, this procedure can perform better than standard matching procedures.

Parenteral is more efficient than mucosal immunization to induce regression of human papillomavirus-associated genital tumors.

Cervical cancer is a public health concern as it represents the second cause of cancer death in women worldwide. High-risk human papillomaviruses (HPV) are the etiologic agents, and HPV E6 and/or E7 oncogene-specific therapeutic vaccines are under development to treat HPV-related lesions in women. Whether the use of mucosal routes of immunization may be preferable for inducing cell-mediated immune responses able to eradicate genital tumors is still debated because of the uniqueness of the female genital mucosa (GM) and the limited experimentation. Here, we compared the protective activity resulting from immunization of mice via intranasal (i.n.), intravaginal (IVAG) or subcutaneous (s.c.) routes with an adjuvanted HPV type 16 E7 polypeptide vaccine. Our data show that s.c. and i.n. immunizations elicited similar frequencies and avidity of TetE71CD81 and E7-specific Interferon-gamma-secreting cells in the GM, whereas slightly lower immune responses were induced by IVAG immunization. In a novel orthotopic murine model, both s.c. and i.n. immunizations allowed for complete long-term protection against genital E7-expressing tumor challenge. However, only s.c. immunization induced complete regression of already established genital tumors. This suggests that the higher E7-specific systemic response observed after s.c. immunization may contribute to the regression of growing genital tumors, whereas local immune responses may be sufficient to impede genital challenges. Thus, our data show that for an efficiently adjuvanted protein-based vaccine, parenteral vaccination route is superior to mucosal vaccination route for inducing regression of established genital tumors in a murine model of HPV-associated genital cancer.

Intravaginal and subcutaneous immunization induced vaccine specific CD8 T cells and tumor regression in the bladder.

PURPOSE: Vaccines targeting tumor associated antigens are in development for bladder cancer. Most of these cancers are nonmuscle invasive at diagnosis and confined in the mucosa and submucosa. However, to our knowledge how vaccination may induce the regression of tumors at such mucosal sites has not been examined previously. We compared different immunization routes for the ability to induce vaccine specific antitumor CD8 T cells in the bladder and bladder tumor regression in mice. MATERIALS AND METHODS: In the absence of a murine bladder tumor model expressing a tumor antigen relevant for human use we established an orthotopic model expressing the HPV-16 tumor antigen E7 as a model. We used an adjuvant E7 polypeptide to induce CD8 T cell mediated tumor regression. RESULTS: Subcutaneous and intravaginal but not intranasal vaccination induced a high number of TetE7(+)CD8(+) T cells in the bladder as well as bladder tumor regression. The entry of vaccine specific T cells in the bladder was not the only key since persistent regression of established bladder tumors by intravaginal or subcutaneous immunization was associated with tumor infiltration of total CD4 and CD8 T cells. This resulted in an increase in TetE7(+)CD8(+) T cells and a decrease in T regulatory cells, leading to an increased number of effector interferon-γ secreting vaccine specific CD8 T cells in the regressing bladder tumor. CONCLUSIONS: These data show that immunization routes should be tailored to each mucosal tumor site. Subcutaneous or intravaginal vaccination may be of additional value to treat patients with bladder cancer.

Evaluating thermal treeline indicators based on air and soil temperature using an air-to-soil temperature transfer model

In recent research, both soil (root-zone) and air temperature have been used as predictors for the treeline position worldwide. In this study, we intended to (a) test the proposed temperature limitation at the treeline, and (b) investigate effects of season length for both heat sum and mean temperature variables in the Swiss Alps. As soil temperature data are available for a limited number of sites only, we developed an air-to-soil transfer model (ASTRAMO). The air-to-soil transfer model predicts daily mean root-zone temperatures (10cm below the surface) at the treeline exclusively from daily mean air temperatures. The model using calibrated air and root-zone temperature measurements at nine treeline sites in the Swiss Alps incorporates time lags to account for the damping effect between air and soil temperatures as well as the temporal autocorrelations typical for such chronological data sets. Based on the measured and modeled root-zone temperatures we analyzed. the suitability of the thermal treeline indicators seasonal mean and degree-days to describe the Alpine treeline position. The root-zone indicators were then compared to the respective indicators based on measured air temperatures, with all indicators calculated for two different indicator period lengths. For both temperature types (root-zone and air) and both indicator periods, seasonal mean temperature was the indicator with the lowest variation across all treeline sites. The resulting indicator values were 7.0 degrees C +/- 0.4 SD (short indicator period), respectively 7.1 degrees C +/- 0.5 SD (long indicator period) for root-zone temperature, and 8.0 degrees C +/- 0.6 SD (short indicator period), respectively 8.8 degrees C +/- 0.8 SD (long indicator period) for air temperature. Generally, a higher variation was found for all air based treeline indicators when compared to the root-zone temperature indicators. Despite this, we showed that treeline indicators calculated from both air and root-zone temperatures can be used to describe the Alpine treeline position.

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.

Radioimmunotherapy of human colon carcinoma by 131I-labelled monoclonal anti-CEA antibodies in a nude mouse model.

A mixture of 3 MAbs directed against 3 different CEA epitopes was radiolabelled with 131I and used for the treatment of a human colon carcinoma transplanted s.c. into nude mice. Intact MAbs and F(ab')2 fragments were mixed because it had been shown by autoradiography that these 2 antibody forms can penetrate into different areas of the tumor nodule. Ten days after transplantation of colon tumor T380 a single dose of 600 microCi of 131I MAbs was injected i.v. The tumor grafts were well established (as evidenced by exponential growth in untreated mice) and their size continued to increase up to 6 days after radiolabelled antibody injection. Tumor shrinking was then observed lasting for 4-12 weeks. In a control group injected with 600 microCi of 131I coupled to irrelevant monoclonal IgG, tumor growth was delayed, but no regression was observed. Tumors of mice injected with the corresponding amount of unlabelled antibodies grew like those of untreated mice. Based on measurements of the effective whole-body half-life of injected 131I, the mean radiation dose received by the animals was calculated to be 382 rads for the antibody group and 478 rads for the normal IgG controls. The genetically immunodeficient animals exhibited no increase in mortality, and only limited bone-marrow toxicity was observed. Direct measurement of radioactivity in mice dissected 1, 3 and 7 days after 131I-MAb injection showed that 25, 7.2 and 2.2% of injected dose were recovered per gram of tumor, the mean radiation dose delivered to the tumor being thus more than 5,000 rads. These experiments show that therapeutic doses of radioactivity can be selectively directed to human colon carcinoma by i.v. injection of 131I-labelled anti-CEA MAbs.

In vivo targeting of an anti-tumor antibody coupled to antigenic MHC class I complexes induces specific growth inhibition and regression of established syngeneic tumor grafts.

The concept of antibody-mediated targeting of antigenic MHC/peptide complexes on tumor cells in order to sensitize them to T-lymphocyte cytotoxicity represents an attractive new immunotherapy strategy. In vitro experiments have shown that an antibody chemically conjugated or fused to monomeric MHC/peptide can be oligomerized on the surface of tumor cells, rendering them susceptible to efficient lysis by MHC-peptide restricted specific T-cell clones. However, this strategy has not yet been tested entirely in vivo in immunocompetent animals. To this aim, we took advantage of OT-1 mice which have a transgenic T-cell receptor specific for the ovalbumin (ova) immunodominant peptide (257-264) expressed in the context of the MHC class I H-2K(b). We prepared and characterized conjugates between the Fab' fragment from a high-affinity monoclonal antibody to carcinoembryonic antigen (CEA) and the H-2K(b) /ova peptide complex. First, we showed in OT-1 mice that the grafting and growth of a syngeneic colon carcinoma line transfected with CEA could be specifically inhibited by systemic injections of the conjugate. Next, using CEA transgenic C57BL/6 mice adoptively transferred with OT-1 spleen cells and immunized with ovalbumin, we demonstrated that systemic injections of the anti-CEA-H-2K(b) /ova conjugate could induce specific growth inhibition and regression of well-established, palpable subcutaneous grafts from the syngeneic CEA-transfected colon carcinoma line. These results, obtained in a well-characterized syngeneic carcinoma model, demonstrate that the antibody-MHC/peptide strategy can function in vivo. Further preclinical experimental studies, using an anti-viral T-cell response, will be performed before this new form of immunotherapy can be considered for clinical use.

Improving generalized regression analysis for the spatial prediction of forest communities

Aim This study used data from temperate forest communities to assess: (1) five different stepwise selection methods with generalized additive models, (2) the effect of weighting absences to ensure a prevalence of 0.5, (3) the effect of limiting absences beyond the environmental envelope defined by presences, (4) four different methods for incorporating spatial autocorrelation, and (5) the effect of integrating an interaction factor defined by a regression tree on the residuals of an initial environmental model. Location State of Vaud, western Switzerland. Methods Generalized additive models (GAMs) were fitted using the grasp package (generalized regression analysis and spatial predictions, http://www.cscf.ch/grasp). Results Model selection based on cross-validation appeared to be the best compromise between model stability and performance (parsimony) among the five methods tested. Weighting absences returned models that perform better than models fitted with the original sample prevalence. This appeared to be mainly due to the impact of very low prevalence values on evaluation statistics. Removing zeroes beyond the range of presences on main environmental gradients changed the set of selected predictors, and potentially their response curve shape. Moreover, removing zeroes slightly improved model performance and stability when compared with the baseline model on the same data set. Incorporating a spatial trend predictor improved model performance and stability significantly. Even better models were obtained when including local spatial autocorrelation. A novel approach to include interactions proved to be an efficient way to account for interactions between all predictors at once. Main conclusions Models and spatial predictions of 18 forest communities were significantly improved by using either: (1) cross-validation as a model selection method, (2) weighted absences, (3) limited absences, (4) predictors accounting for spatial autocorrelation, or (5) a factor variable accounting for interactions between all predictors. The final choice of model strategy should depend on the nature of the available data and the specific study aims. Statistical evaluation is useful in searching for the best modelling practice. However, one should not neglect to consider the shapes and interpretability of response curves, as well as the resulting spatial predictions in the final assessment.

Geomodelling of a fluvial system with semi-supervised support vector regression

Fluvial deposits are a challenge for modelling flow in sub-surface reservoirs. Connectivity and continuity of permeable bodies have a major impact on fluid flow in porous media. Contemporary object-based and multipoint statistics methods face a problem of robust representation of connected structures. An alternative approach to model petrophysical properties is based on machine learning algorithm ? Support Vector Regression (SVR). Semi-supervised SVR is able to establish spatial connectivity taking into account the prior knowledge on natural similarities. SVR as a learning algorithm is robust to noise and captures dependencies from all available data. Semi-supervised SVR applied to a synthetic fluvial reservoir demonstrated robust results, which are well matched to the flow performance

Regression analysis of spatial data.

Many of the most interesting questions ecologists ask lead to analyses of spatial data. Yet, perhaps confused by the large number of statistical models and fitting methods available, many ecologists seem to believe this is best left to specialists. Here, we describe the issues that need consideration when analysing spatial data and illustrate these using simulation studies. Our comparative analysis involves using methods including generalized least squares, spatial filters, wavelet revised models, conditional autoregressive models and generalized additive mixed models to estimate regression coefficients from synthetic but realistic data sets, including some which violate standard regression assumptions. We assess the performance of each method using two measures and using statistical error rates for model selection. Methods that performed well included generalized least squares family of models and a Bayesian implementation of the conditional auto-regressive model. Ordinary least squares also performed adequately in the absence of model selection, but had poorly controlled Type I error rates and so did not show the improvements in performance under model selection when using the above methods. Removing large-scale spatial trends in the response led to poor performance. These are empirical results; hence extrapolation of these findings to other situations should be performed cautiously. Nevertheless, our simulation-based approach provides much stronger evidence for comparative analysis than assessments based on single or small numbers of data sets, and should be considered a necessary foundation for statements of this type in future.

A novel mucosal orthotopic murine model of human papillomavirus-associated genital cancers.

Cervical cancer results from infection with high-risk type human papillomaviruses (HPV). Therapeutic vaccines aiming at controlling existing genital HPV infections and associated lesions are usually tested in mice with HPV-expressing tumor cells subcutaneously implanted into their flank. However, effective vaccine-induced regression of these ectopic tumors strongly contrasts with the poor clinical results of these vaccines produced in patients with HPV-associated genital neoplasia. To assess HPV therapeutic vaccines in a more relevant setting, we have, here, established an orthotopic mouse model where tumors in the genital mucosa (GM) develop after an intravaginal instillation of HPV16 E6/E7-expressing tumor cells transduced with a luciferase-encoding lentiviral vector for in vivo imaging of tumor growth. Tumor take was 80-90% after nonoxynol-9 induced damage of the epithelium. Tumors remained localized in the genital tract, and histological analysis showed that most tumors grew within the squamous epithelium of the vaginal wall. Those tumors induced (i) E7-specific CD8 T cells restricted to the GM and draining lymph nodes, in agreement with their mucosal location and (ii) high Foxp3+ CD4+ infiltrates, similarly to those found in natural non-regressing HPV lesions. This novel genital HPV-tumor model by requiring GM homing of vaccine-induced immune responses able to overcome local immuno-suppression may be more representative of the situation occurring in patients upon therapeutic vaccination.

MGMT methylation analysis of glioblastoma on the Infinium methylation BeadChip identifies two distinct CpG regions associated with gene silencing and outcome, yielding a prediction model for comparisons across datasets, tumor grades, and CIMP-status.

The methylation status of the O(6)-methylguanine-DNA methyltransferase (MGMT) gene is an important predictive biomarker for benefit from alkylating agent therapy in glioblastoma. Recent studies in anaplastic glioma suggest a prognostic value for MGMT methylation. Investigation of pathogenetic and epigenetic features of this intriguingly distinct behavior requires accurate MGMT classification to assess high throughput molecular databases. Promoter methylation-mediated gene silencing is strongly dependent on the location of the methylated CpGs, complicating classification. Using the HumanMethylation450 (HM-450K) BeadChip interrogating 176 CpGs annotated for the MGMT gene, with 14 located in the promoter, two distinct regions in the CpG island of the promoter were identified with high importance for gene silencing and outcome prediction. A logistic regression model (MGMT-STP27) comprising probes cg1243587 and cg12981137 provided good classification properties and prognostic value (kappa = 0.85; log-rank p < 0.001) using a training-set of 63 glioblastomas from homogenously treated patients, for whom MGMT methylation was previously shown to be predictive for outcome based on classification by methylation-specific PCR. MGMT-STP27 was successfully validated in an independent cohort of chemo-radiotherapy-treated glioblastoma patients (n = 50; kappa = 0.88; outcome, log-rank p < 0.001). Lower prevalence of MGMT methylation among CpG island methylator phenotype (CIMP) positive tumors was found in glioblastomas from The Cancer Genome Atlas than in low grade and anaplastic glioma cohorts, while in CIMP-negative gliomas MGMT was classified as methylated in approximately 50 % regardless of tumor grade. The proposed MGMT-STP27 prediction model allows mining of datasets derived on the HM-450K or HM-27K BeadChip to explore effects of distinct epigenetic context of MGMT methylation suspected to modulate treatment resistance in different tumor types.

Robust accelerated failure time regression

Robust estimators for accelerated failure time models with asymmetric (or symmetric) error distribution and censored observations are proposed. It is assumed that the error model belongs to a log-location-scale family of distributions and that the mean response is the parameter of interest. Since scale is a main component of mean, scale is not treated as a nuisance parameter. A three steps procedure is proposed. In the first step, an initial high breakdown point S estimate is computed. In the second step, observations that are unlikely under the estimated model are rejected or down weighted. Finally, a weighted maximum likelihood estimate is computed. To define the estimates, functions of censored residuals are replaced by their estimated conditional expectation given that the response is larger than the observed censored value. The rejection rule in the second step is based on an adaptive cut-off that, asymptotically, does not reject any observation when the data are generat ed according to the model. Therefore, the final estimate attains full efficiency at the model, with respect to the maximum likelihood estimate, while maintaining the breakdown point of the initial estimator. Asymptotic results are provided. The new procedure is evaluated with the help of Monte Carlo simulations. Two examples with real data are discussed.

Model development for capercaillie (Tetrao urogallus) habitat in the Jura Mountains (Western Switzerland)

Capercaillie, Tetrao urogallus, is a threatened species in central Europe, and Swiss populations declined 40 to 50 % between 1970 and 1985. Capercaillie are sensitive to forest structure, and loss of habitat is a major cause of their decline. Knowledge of habitat characteristics is therefore essential for capercaillie conservation. Here, we present models predicting capercaillie probability of occurrence, based on relevant structural habitat variables. Models were built using multiple logistic regression analyses on capercaillie presence/absence data. Vegetation survey was carried out in July 1999 in a 170-km2 forested area (Jura mountains, canton de Vaud, western Switzerland) inhabited by capercaillie and presence/absence of the species was assessed according to dropping presence/absence. The survey was based on 10-m-radius sample plots each in a 1-km2 forest patch (n = 76 with capercaillie droppings, n = 80 without). A first model included seven out of 27 measured habitat variables and a second model only four. The latter model best represents practical needs. It includes three variables which had a negative impact on capercaillie presence: tree and shrub covers and spruce, Picea excelsa, shrub cover, and one which had a positive effect: bilberry, Vaccinium myrtillus, cover, highlighting that capaercaillie selected open forest with high bilberry abundance. The model can be used to map potential capercaillie habitat distribution and to manage the habitat in favour of capercaillie (protection and adapted forestry practices) in the Swiss Jura mountains.