The linear time series model refers to the class of models for which fixed correlation parameters can fully explain the dependency between two random variables, but many real-life circumstances, such as monthly unemployment results, supplies and demands, interest rate, exchange rate, share prices, rainfall, etc., violate the assumption of linearity. For fitting and forecasting of nonlinear time series data, the self-exiting threshold autoregressive (SETAR) model was suggested. Using R to generate random nonlinear autoregressive data, a Monte Carlo simulation was performed, the SETAR model was fitted to the simulated data and Lafia rainfall data, Nasarawa State, Nigeria to determine the best regime orders and/or scheme number to make future forecast. Using Mean Square Error (MSE) and Akaike Information Criteria (AIC), the relative performance of models was examined. At a specific autoregressive order, regime order, sample size and step ahead, the model with minimum criteria was considered as the best. The results show that the best autoregressive and regime orders to be chosen are 3rd and 2nd [SETAR (3, 2)] respectively for fitting and forecasting nonlinear autoregressive time series data with small and moderate sample sizes. As the sample size increases, the output of the four models increases. Finally, it is shown that when sample size and number of steps forward are increased, the efficiency and forecasting capacity of the four models improves.