“…As mentioned earlier, in the GLM, the error variable follows a distribution of the exponential family, which includes the Normal, Poisson, Binomial, and Gamma distributions. The linear coefficients are estimated using the maximum likelihood estimation (MLE) method if the residuals are non-Normal or the ordinary least squares (OLS), if Normal(Yuan and Yang 2005, Yan and Su 2009, Mahmoud 2019). However, if there is a large number of dummy variables; as a result…”