The main purpose of this study is to apply three bivariate statistical models, namely weight of evidence (WoE), evidence belief function (EBF) and index of entropy (IoE), and their ensembles with logistic regression (LR) for landslide susceptibility mapping in Muchuan County, China. First, a landslide inventory map contained 279 landslides was obtained through the field investigation and interpretation of aerial photographs. Next, the landslides were randomly divided into two parts for training and validation with the ratio of 70/30. In addition, according to the regional geological environment characteristics, twelve landslide conditioning factors were selected, including altitude, plan curvature, profile curvature, slope angle, distance to roads, distance to rivers, topographic wetness index (TWI), normalized different vegetation index (NDVI), land use, soil, and lithology. Subsequently, the landslide susceptibility mapping was carried out by the above models. Eventually, the accuracy of this research was validated by the area under the receiver operating characteristic (ROC) curve and the results indicated that the landslide susceptibility map produced by EBF-LR model has the highest accuracy (0.826), followed by IoE-LR model (0.825), WoE-LR model (0.792), EBF model (0.791), IoE model (0.778), and WoE model (0.753). The results of this study can provide references of landslide prevention and land use planning for local government.Symmetry 2019, 11, 762 2 of 24years, machine learning method has been gradually applied in landslide susceptibility mapping researches, such as artificial neural network (ANN) [17][18][19], support vector machine (SVM) [20][21][22], logistic model tree (LMT) [23,24], rotation forest (RF) [25,26], classification and regression tree (CART) [27,28], adaptive neuro-fuzzy inference systems (ANFIS) [29,30], and genetic algorithm (GA) [31,32]. Furthermore, statistical approach is another widely-used model which can also be divided into two types: bivariate and multivariate. In statistical approaches, the weight of each class of every factor was calculated by overlying the landslide inventory map and landslide conditioning factors map [33,34]. The frequently-used statistical approaches are frequency ratio (FR) [35][36][37], logistic regression (LR), evidential belief function (EBF) [38,39], weights of evidence (WoE) [40,41], certainty factor (CF) [42,43], and information value (IV) [44,45].In this paper, three bivariate statistical methods, namely weight of evidence (WoE), evidence belief function (EBF), index of entropy (IoE), coupled with logistic regression (LR) were introduced to carry out the landslide susceptibility mapping in Muchuan County, which lies in the southeast of Sichuan Province, China. All the approaches were first applied in Muchuan County and the main purpose of this article is to obtain the landslide susceptibility maps through the above models. The main difference between this study and literatures referenced above is to achieve accurate landslide susceptibility maps for...