The uncertainty is an important attribute about data that can arise from different sources including randomness and fuzziness, therefore in uncertain environments, especially, in modeling, planning, decision-making, and control under uncertainty, most data available contain some degree of fuzziness, randomness, or both, and at the same time, some of this data may be anomalous (outliers). In this regard, the new fuzzy regression approaches by creating a functional relationship between response and explanatory variables can provide efficient tools to explanation, prediction and possibly control of randomness, fuzziness, and outliers in the data obtained from uncertain environments. In the present study, we propose a new two-stage fuzzy linear regression model based on a new interval type-2 (IT2) fuzzy least absolute deviation (FLAD) method so that regression coefficients and dependent variables are trapezoidal IT2 fuzzy numbers and independent variables are crisp. In the first stage, to estimate the IT2 fuzzy regression coefficients and provide an initial model (by original dataset), we introduce two new distance measures for comparison of IT2 fuzzy numbers and propose a novel framework for solving fuzzy mathematical programming problems. In the second stage, we introduce a new procedure to determine the mild and extreme fuzzy outlier cutoffs and apply them to remove the outliers, and then provide the final model based on a clean dataset. Furthermore, to evaluate the performance of the proposed methodology, we introduce and employ suitable goodness of fit indices. Finally, to illustrate the theoretical results of the proposed method and explain how it can be used to derive the regression model with IT2 trapezoidal fuzzy data, as well as compare the performance of the proposed model with some well-known models using training data designed by Tanaka et al. [55], we provide two numerical examples.