The raw data with outliers were transformed to normal distribution,which was necessary for the studied variables. All the results were stored in Microsoft Excel® spreadsheets. The descriptive parameters were calculated using Microsoft Excel® and SPSS® for Windows (version 15.0). Local Indicators of Spatial Association values were measured in the GeoDa software (version 0.95i, Spatial Analysis Laboratory, 2007). All maps were produced using GIS software ArcMap