Modeling and forecasting exchange rate volatility in Bangladesh using GARCH models: a comparison based on normal and Student’s t-error distribution

Modeling exchange rate volatility has remained crucially important because of its diverse implications. This study aimed to address the issue of error distribution assumption in modeling and forecasting exchange rate volatility between the Bangladeshi taka (BDT) and the US dollar ($). Using daily exchange rates for 7 years (January 1, 2008, to April 30, 2015), this study attempted to model dynamics following generalized autoregressive conditional heteroscedastic (GARCH), asymmetric power ARCH (APARCH), exponential generalized autoregressive conditional heteroscedstic (EGARCH), threshold generalized autoregressive conditional heteroscedstic (TGARCH), and integrated generalized autoregressive conditional heteroscedstic (IGARCH) processes under both normal and Student’s t-distribution assumptions for errors. It was found that, in contrast with the normal distribution, the application of Student’s t-distribution for errors helped the models satisfy the diagnostic tests and show improved forecasting accuracy. With such error distribution for out-of-sample volatility forecasting, AR(2)–GARCH(1, 1) is considered the best.

can help to design policies to mitigate the adverse effects of exchange rate volatility on important economic indicators.
Accordingly, beginning with Engle's (1982) autoregressive conditional heteroscedasticity (ARCH) model, several models have already been developed to model volatility. Different models have aimed to capture different features of volatility. For example, while some models are only used to model "volatility clustering," 1 others are used to capture the "leverage effect." 2 However, it is widely recognized among researchers that because of increased kurtosis, along with an increase in data frequency, the rate of return in financial or macroeconomic variables as exchange rates might have a levy distribution, or a fat tail (Mandelbrot, 1963). Thus, volatility models should not be considered under the general normality assumption for errors, because if the errors are not thus, then volatility forecasting based on such models will be misleading. The current study makes an effort to appropriately model and forecast the volatility of the exchange rate return of the taka against the US dollar, taking into account the issue of error distribution assumption.
The rest of this paper is organized as follows. Literature review presents a review of previous literature that has attempted to model and forecast exchange rate volatility using various ARCH and generalized autoregressive conditional heteroscedasticity (GARCH) models. Methods describes the data and the theoretical methodology. Estimation results and findings analyzes and compares the results of different ARCH and GARCH models. Forecasting exchange rate return volatility is discussed in Volatility Forecasting, and Conclusion concludes the paper.

Literature review
Modeling exchange rate volatility has remained crucially important because of its diverse implications. Bala and Asemota (2013) examined exchange rate volatility using GARCH models. They used monthly exchange rate return series for the naira (Nigerian currency) against the US dollar ($), British pound, and euro. To compare the estimates, various GARCH models were estimated with and without volatility breaks. It was revealed that most of the models rejected the existence of a leverage effect, except for those with volatility breaks. Since it was observed that results improved when the volatility models considered breaks, incorporating significant events in the GARCH models was suggested. Clement and Samuel (2011) also aimed to model Nigerian exchange rate volatility. They used the monthly exchange rate of the naira against the US dollar and British pound for the period from 2007 to 2010. They found that the exchange rate return series was nonstationary and that the series residuals were asymmetric. Since return volatility was found to be persistent, the study recommended further investigation of the impact of government policies on foreign exchange rates. Rofael and Hosni (2015) aimed to forecast and estimate exchange rate volatility in Egypt using ARCH and state space (SS) models. Using daily exchange rate data covering about 10 years, they found volatility clustering in Egypt's exchange rate returns, as well as a risk of mismatch between exchange rates and the stock market. Similar results were obtained by Choo, Loo, and Ahmad (2002), who used GARCH model variants to capture the exchange rate volatility dynamics of the Malaysian ringgit (RM) against the pound sterling. They used daily data for the period from 1990 to 1997 and concluded that the volatility of the RM-sterling exchange rate was persistent. For within-sample modeling, they found the GARCH models to be the best, while for forecasting, almost all the GARCHin -mean models outperformed ordinary GARCH models. Dhamija and Bhalla (2010) argued that conditionally heteroscedastic models can be used to model exchange rate volatility. They found that integrated generalized autoregressive conditional heteroscedsticity (IGARCH) and threshold generalized autoregressive conditional heteroscedsticity (TGARCH) models performed better than others when forecasting the volatility of five daily currencies: the British pound, German mark, Japanese yen, Indian rupee and Euro. Later, Ramasamy and Munisamy (2012) concluded that GARCH models were efficient for predicting exchange rate volatility. They examined the daily exchange rates of four currencies-the Australian dollar, Singapore dollar, Thai bhat, and Philippine peso-using GARCH, Glostern -Jagannathan -Runkle GARCH (GJR-GARCH), and exponential generalized autoregressive conditional heteroscedsticity (EGARCH) models. They argued that the improvements made by leveraging in EGARCH and GJR-GARCH models did not improve forecasting accuracy.
Meanwhile, Brooks and Burke (1998) used modified information criteria to select models from the GARCH family. Using weekly exchange rate returns for the Canadian dollar, German mark, and Japanese yen against the US dollar for the period from March 1973 to September 1989, they compared the performance of different out-ofsample models. They found that the out-of-sample forecasting accuracy of the models compared favorably on mean absolute errors but less favorably on mean squared errors with those generated by commonly used GARCH(1, 1) models. However, Hansen and Lunde (2005) did an out-of-sample comparison of 330 different GARCH family models using two types of assets: daily exchange rates (Deutsche Mark (DM)-US dollar) and IBM stock prices. They found that different models were the best for forecasting the volatility of the two types of assets. In terms of forecasting accuracy, GARCH(1, 1) outperformed the other models. For IBM stock prices, however, the models with leverage effects performed better than GARCH(1, 1). Herwartz and Reimers (2002) analyzed daily exchange rate changes between the Deutsche mark (DM) and the US dollar and the DM and the Japanese yen (JPY) for the period from 1975 to 1998. They used a GARCH(1, 1) model with leptokurtic innovations to capture volatility clustering, and they found that the identified points of structural change were subject to changes in monetary policies in the US and Japan. Similarly, Çağlayan, Ün, and Dayıoğlu (2013) modeled the exchange rate volatility of MIST (Mexico, Indonesia, South Korea, and Turkey) countries against the US dollar using asymmetric GARCH models. They used monthly exchange rate data for the period from 1993 to 2013 to investigate leverage effects and fat-tailed features. They identified the existence of asymmetrical and leveraging effects in the exchange rates of MIST countries against the US dollar. Vee, Gonpot, and Sookia (2011) also examined the forecasting accuracy of GARCH(1, 1) using Student's t and generalized error distribution (GED). Using daily data for exchange rates between the US dollar and the Mauritian rupee, they compared the mean absolute error (MAE) and root mean squared error (RMSE) of the models based on forecasting estimates. They found that GARCH(1, 1) with GED had better forecasting accuracy compared to that using Student's t-distribution.
In contrast, Tse (1998) examined the conditional heteroscedasticity of yen-US dollar exchange rates using daily observations for the period from 1978 to 1994. Extending APARCH models to a process that was fractionally integrated; they found that, unlike the stock market, the appreciation and depreciation shock of the yen against the US dollar affected future volatility in a similar fashion. They argued that there is no substantial difference between fractionally integrated models and stable models. More recently, Pelinescu (2014) analyzed exchange rates between the Romanian leu and the euro considering the influence of other macroeconomic variables. Using daily observations for the period from 2000 to 2013, the study found that the exchange rates consisted of ARCH processes, and exchange returns were correlated with volatility.
To the best of our knowledge, relatively few studies have aimed to model BDT-US dollar exchange rates. In one recent example, Alam and Rahman (2012) modeled the aforementioned exchange rates using daily data for the period from July 2006 to April 2012. That study has limitations since excess skewness and kurtosis issues were detected but not addressed. Further, the lag specification of the mean equation was not properly addressed. Such serious flaws can result in misleading conclusions. Therefore, finding a proper model with an appropriate assumption for the distribution of errors would be a valuable contribution to the field of exchange rate modeling in Bangladesh regarding the BDT against the US dollar.

Data and variable construction
The data used for this study covered January 1, 2008 to April 30, 2015. These were daily observations of the nominal exchange rate of the Bangladeshi taka (BDT) against the US dollar supplied by Bangladesh Bank, the country's central bank. Since the nominal exchange rate series is usually nonstationary, and hence is not appropriate for analysis, we converted the series into the rate of return on exchange rate by following logarithmic transformation. In particular, we used the following formula to calculate the rate of return on exchange rate: Here, r t stands for exchange rate return at period t; f x t and f x t-1 stand for the nominal exchange rate of the BDT-US dollar at period t and (t-1). The statistical software EViews 9 (Econometric Views Version 9, IHS Markit Ltd.) was used to perform the quantitative exercise.

Specification of different models
It is well established that when modeling volatility using GARCH family models, the appropriate specification of the mean equation is vitally important. The misspecification of that equation could fail to address the autocorrelation problem that could arise in the volatility model. Therefore, we estimated all of the GARCH family models using three different specifications for the mean equations. The first is with a constant only, the second is with a constant and one autoregressive variable (i.e., AR(1)), and the third is constant and has two autoregressive variables (i.e., AR(2)). They can be expressed as follows: To specify the variance equation to model volatility presence in differenced logarithmic exchange rates (i.e., exchange rate returns), we used five different GARCH family models, each of which has a different purpose. In particular, we modeled the variance for the above three mean equations using GARCH, APARCH, EGARCH, TGARCH, and IGARCH models. The sensitivity of the estimation results of the models was checked by changing the distribution assumptions. More specifically, the sensitivity and appropriateness of the estimation results were observed by changing the distribution assumption from normal to Student's t-distribution. The reason for this is that in the literature on financial asset returns, it is established that the return variable is more likely to follow a "levy distribution" with "fat tails," and kurtosis is likely to increase with data frequency (Alexander, 1961;Andersen & Bollerslev, 1998). Engle (1982) pioneered modeling volatility using conditionally heteroscedastic regression with the autoregressive conditional heteroscedasticity (ARCH) model. One major problem with such modeling is that the required lag length usually remains large, which means it needs to estimate a large number of parameters to predict volatility. In contrast, the generalized autoregressive conditional heteroscedasticity (GARCH) models proposed by Bollerslev (1986) allow conditional variance to depend upon its own lag. This typically reduces the number of required ARCH lags when predicting volatility.
Nevertheless, we assumed the following general form for the variance equation: Depending on the specification of h t in Eq. (4), we can have several possible models within the GARCH family. In particular, we estimated the following specification of GARCH family models for estimating volatility in logarithmic exchange rates: Here, the GARCH(1, 1) model consists of one ARCH term denoted as ε 2 t−1 and one GARCH term denoted as h t−1 . For the variance to remain well behaved, some restrictions needed to be imposed: η > 0, α ≥ 0, and β ≥ 0. The sum of the ARCH coefficient and GARCH coefficient governs the persistency of volatility shocks. Their sum should be less than the unit (α + β < 1) to ensure that series ε t is stationary and the variance is positive.
To allow for possible nonlinearities in the parameters of the variance equation, Ding, Granger and Engel (1993) developed the asymmetric power autoregressive conditional heteroscedasticity (APARCH) model. In the above APARCH model, δ denotes the power parameter that requires the condition δ > 0, and γ is the parameter capturing asymmetry or leverage effect, which requires the condition |γ| ≤ 1.
The standard GARCH(1, 1) model suggests that the shock in ε t−1 has the same effect regardless of whether ε t−1 > 0 or ε t−1 < 0. Nevertheless, a typical feature of financial data is that negative shocks generate more volatility compared to positive shocks. Aiming to incorporate this asymmetrical effect in financial data, Nelson (1991) proposed the exponential generalized autoregressive conditional heteroscedasticity (EGARCH) model. In the EGARCH specification, γ is the asymmetry parameter measuring leverage effect, α is the size parameter measuring the magnitude of shocks, and persistency is captured through β. An important feature of the EGARCH specification is that conditional variance is an exponential function, thus there is no need for nonnegetivity restrictions, as in earlier GARCH specifications.
Introduced by Zakoian (1994) and Glosten, Jagannathan, and Runkle (1993), threshold generalized autoregressive conditional heteroscedasticity (TGARCH) is another model developed to analyze leverage affects. In the above TGARCH(1, 1) model, ε t−1 > 0 (good news) and ε t−1 < 0 (bad news) produce a differential effect on volatility. Here, good news has an impact of α while bad news has an impact of (α + γ). Thus, when γ > 0, the increase in volatility due to bad news is greater than that for good news, and we can conclude that there is a leverage effect. Similar to standard GARCH models, we also need nonnegative restrictions here for α, γ, and β. A nonstationary GARCH model (i.e., GARCH model with unit roots) can be regarded as an integrated generalized autoregressive conditional heteroscedasticity (IGARCH) model. This model was originally developed by Engle and Bollerslev (1986). In practice, when the parameters of the GARCH model are restricted to a sum equal to one and ignore the constant term, a standard GARCH model is transformed into an IGARCH model. In the above IGARCH model, the additional constraints are {α + (1 − α)} = 1 and 0 < α < 1. Table 1 shows the estimation results of the different conditional mean models for the differenced logarithmic exchange rate using Ordinary Least Squares (OLS) method. OLS was applied to estimate regression since the series was observed to be mean reverting ( Table 9 in Appendix). In particular, three different conditional mean equations were estimated: one with the constant only, one with the AR(1) term, and one with AR(2) terms. The models were not augmented with further AR terms since they were not significant ( Fig. 1 in Appendix). As the results show, in Model 1, the constant term is significant, and the F-statistic, testing the null of the absence of the ARCH effect in the differenced logarithmic exchange rate, is statistically significant at the 1% level. In Model 2, the conditional mean model was augmented with the AR(1) term, which was found to be statistically significant along with the test for heteroscedastic ARCH effects. The conditional mean model was augmented further with AR(2) terms in Model 3. The estimation results reveal that both AR(1) and AR(2) terms are significant along with the F-statistic for the heteroscedastic ARCH effect. Thus, it is established that the differenced logarithmic exchange rate of the taka against the US dollar is conditionally heteroscedastic, and this remains crucially important for modeling volatility clustering. The existence of volatility clustering in the taka-US dollar exchange rate is also evident in the residual plot in Fig. 2 (Appendix).

Estimation Results and Discussion
Since the ARCH effect was detected in the model, GARCH estimation was performed with different specifications. In particular, the GARCH(1, 1) and APARCH(1, 1) models were initially estimated with normal error distribution assumption. While the GARCH model was applied for the purpose of capturing variance dynamics, the APARCH model was applied to test for the presence of asymmetric volatility effects in the variance dynamics of exchange rate returns. Table 2 shows the estimation results. Here, it is evident that the autoregressive coefficients of the lagged dependent variables for mean equations ρ 1 and ρ 2 are statistically significant. In all of the specifications, the coefficients of the GARCH components (i.e., α and β) are positive and statistically significant at the 1% level. However, as their sum exceeds 1, the residuals of the regressions would be nonstationary, and hence variance would be indefinite. Since the residual of the GARCH model should be white noise, a diagnostic test was performed in the form of the Ljung-Box Q-test under the null hypothesis (H 0 : No Serial Correlation in the Error Term). We calculated Q-statistics for the standardized residuals (Q1) and for their squared values (Q2) (especially testing for the fourth and eighth lag, inspired by Tse (1998)). All of the Q1-statistics were significant at the 5% level. Thus, there was not enough evidence to reject the null hypothesis of no serial correlation in the error term. Also, the F-statistic confirmed that the models still had ARCH effects.
From the APARCH specification in Table 2, it is evident that the autoregressive coefficients of the lagged dependent variables for mean equations ρ 1 and ρ 2 are Standard Errors are in Parenthesis. *** indicates significant at 10% level, ** significant at 5% level and * indicates that at 1% level statistically significant. The coefficients α and β were again found to be positive and statistically significant, having a sum exceeding 1 for all of the specifications. The coefficient δ is positive and significant, and not equal to 2, establishing that it is not a standard GARCH model (Ding, Granger, & Engle 1993). The significance and sign of the coefficient γ determine the leverage effect. The positive significant value of γ indicates the existence of a leverage effect where negative past values of ε t increase volatility more than positive past values of the same magnitude. Meanwhile, the negative significant value of γ indicates the existence of a leverage effect in the opposite direction: positive past values of ε t increase volatility more than negative past values of the same magnitude. As the table shows, the leverage coefficient is negative and significant. Thus, an asymmetric volatility effect exists for exchange rate returns. More specifically, the appreciation and depreciation of the taka against the US dollar does not necessarily cause symmetric variation in the exchange rate return. However, the test results for detecting serial correlation establish that all of the specifications have a serial correlation problem, and the F-statistic shows that an ARCH effect still exists. Figure 3 (Appendix) shows that none of the above models have normally distributed errors on the basis of the Jarque-Bera test statistic. Since it was established that the distribution of error terms in different models was not normal, we checked the estimation results using the assumption of Student's t-distribution for the error terms. To justify this, we also checked the skewness and kurtosis of the exchange rate return series. Table 10 (Appendix) shows that the exchange rate return is highly skewed. The distribution is also leptokurtic, implying that its central peak is higher and sharper with longer and fatter tails. This finding is in line with Mandelbrot (1963), who found fat-tailed and excess kurtosis for the rate of return, which has been further clarified for different financial assets in other studies. Thus, following Bollerslev (1987); Vee, Gonpot, and Sookia (2011); and Çağlayan, Ün, and Dayıoğlu (2013), we examined the estimation results using Student's t-distribution for the error terms. Table 3 shows the results. The autoregressive coefficients in both GARCH and APARCH specifications were found significant, as before. Here, the ARCH parameter denoted by β measures the reaction of conditional volatility to market shocks (i.e., it measures the variance response in exchange rate returns against appreciation or depreciation). On the other hand, the GARCH parameter denoted by β measures the persistence of conditional volatility, regardless of shocks to the market. Both coefficients were found to be positive and significant, implying that variance in exchange rate returns Robust Standard Errors are in Parenthesis. *** indicates significant at 10% level, ** indicates significant at 5% level and *indicates that at 1% level for the taka against the US dollar is responsive and persistent against shocks defined as appreciation or depreciation. However, in all of the models, since their sum exceeds a value of 1, the variance is arguably infinite and nonstationary. The existence of an asymmetric volatility effect is again established in the APARCH model. Nonetheless, the diagnostic indicators reveal that there is no ARCH effect in either the GARCH and APARCH specifications; also, all models successfully passed the noautocorrelation test in residuals. Therefore, it is more appropriate to perform further analysis using Student's t-distribution. It is worth noting that in all of the types of models, regardless of different distribution assumptions for the residuals, the coefficients for the ARCH and GARCH components sum to more than 1. Also, the constant terms in the variance equations of the GARCH and APARCH models using Student's t-distribution were not significantly different from 0. Thus, as the necessary nonnegativity restrictions have not been satisfied, the variance is not well behaved.
Thus, to incorporate the asymmetric volatility effect and eliminate the problem of nonnegativity restrictions, we estimated an EGARCH model. Table 4 shows the estimation results of the EGARCH model with different specifications for the mean equations, using both normal and Student's t-distribution for the residual. The autoregressive coefficients were significant when we estimated EGARCH using normally distributed errors. Here, in the variance equation, γ represents what is Robust Standard Errors are in Parenthesis. *** indicates significant at 10% level, ** indicates significant at 5% level and * indicates that at 1% level popularly known as the "asymmetry parameter," and α represents the "size parameter." The former measures the asymmetric effect on volatility while the latter measures the effect of the magnitude of shocks about their mean. It can be observed that though the size parameter is significant, the asymmetric parameter is insignificant, establishing the possible absence of an asymmetric effect on volatility.
A look at the diagnostic indicators reveals that the model still contains an ARCH effect and has a serial correlation problem. Furthermore, when the model is reestimated using Student's t-distribution for residuals, both the asymmetry and size parameters become statistically insignificant. Therefore, according to the EGARCH specification, the taka-US dollar exchange rate return exhibits symmetric volatility; appreciation and depreciation could possibly have a similar effect on future volatility. This is in line with Diebold and Nerlove (1989);Bollerslev, Chou, and Kroner (1992); and Tse (1998) since the diagnostic indicators showed no ARCH effect, and the model overcame the autocorrelation problem based on the Ljung-Box Q-test using squared residuals. However, the autocorrelation problem remained when regular residuals were used for the test.
Since EGARCH disregarded the possible existence of an asymmetric volatility effect, we tried another parameterization-namely, TGARCH. Table 5 shows the estimation results. All of the nonnegativity restrictions required for model validity were satisfied. Also, the parameter λ, which captures the asymmetric response of volatility to shocks, Robust Standard Errors are in Parenthesis. *** indicates significant at 10% level, ** indicates significant at 5% level and * indicates that at 1% level was found to be negative and consistently significant for all of the models, indicating the possible existence of an asymmetric volatility effect. However, the findings should be regarded with caution since this estimation is valid only under normal error distribution. 3 Furthermore, this model also shows an autocorrelation problem, and there is still an ARCH effect. Finally, since it has been observed that the sum of the persistence parameters exceeds the unitary value in earlier GARCH estimations, it can be deduced that the variance might not be well behaved in such models. Therefore, it would be interesting to model volatility clustering with such models while imposing restrictions on the persistence parameters. One popular restriction on the persistence parameters of GARCH models is referred to as "persistence parameters sum up to unit." The estimation of GARCH models with this restriction leads to IGARCH specifications. Table 6 shows IGARCH estimation with normally distributed errors as well as errors following Student's t-distribution.
It can be observed from the estimation results of the IGARCH models that the restriction that was applied is valid. Moreover, the IGARCH specifications successfully overcame all of the diagnostic tests when Student's t-distribution, rather than normal distribution, was used as the error distribution. There was no ARCH effect and no autocorrelation detected in the regular and squared residuals.
Finally, the evidence indicates that all of the models (GARCH, APARCH, EGARCH, and IGARCH) satisfy the required diagnostic standard under Student's t-distribution as the assumption for residuals. Further, the log-likelihood for all of the models was improved when the aforementioned distribution was used.

Volatility forecasting
In-sample estimation accuracy To check whether the accuracy of volatility forecasting among the different models varied with distribution assumptions, we compared the log-likelihood, Schwarz Bayesian information criterion (SBC), and Akaike information criterion (AIC) for all of the models, estimating for whole-sample observations under normal and Student's t-error distribution. Table 7 shows the results. It is clear that the performance and goodness of fit of each model improved when Student's t-distribution was used for the residuals. This is because likelihood increased while both SBC and AIC decreased compared to when normal distribution was used for the residuals. Considering the Student's t-distribution for the residuals, a comparison of indicators reveals that among all of the models used for in-sample estimation, AR(2)-APARCH(1, 1) is the best since it has maximum likelihood and minimal SBC and AIC.

Out-of-sample forecasting accuracy
To check the forecasting accuracy of the models, we created a pseudo sample using the period from June 1, 2008 to December 31, 2012. All of the models were estimated for the pseudo sample period and then employed to estimate the variance in exchange rate returns for the period from January 1, 2013 to April 30, 2015. The forecasting performance of the models was compared on the basis of four different indicators under normal distribution and Student's t-distribution: root mean square error (RMSE), mean absolute error (MAE), mean absolute percent error (MAPE), and Theil inequality (TI). Table 8 shows the comparative forecasting accuracy of the different models under normal and Student's t-distribution for the residuals. For the AR(2)-GARCH(1, 1) model, RMSE and TI decreased while MAE remained constant and MAPE increased when changing the distribution assumption. Thus, this model's forecasting accuracy was improved under Student's t-distribution. Likewise, the forecasting accuracy of AR(2)-IGARCH(1, 1) was improved with the above distribution (all of the four indicators decreased when Student's t-distribution was used). Meanwhile, AR(2)-APARCH(1, 1), which had the best in-sample estimation accuracy under Student's t-distribution, showed lower accuracy when such distribution was used for out-of-sample forecasting. Similarly, the forecasting accuracy of AR(2)-EGARCH(1, 1) also became poor when Student's t-distribution was used;

Conclusion
Bangladesh, which the World Bank recently promoted from a low-income country to a lower-middle-income one, is moderately dependent on trade, but remittance holds high importance. In 2014, trade accounted for about 45% of Bangladesh's GDP (about 55% for other low-middle-income countries) while remittance accounted for 8.7% of the GDP (only 1.5% for other low-middle-income countries). Moreover, most exports from Bangladesh are destined for the US. Therefore, exchange rate volatility will have a significant effect on trade and remittance, and consequently on the whole economy. For this reason, it is extremely important to properly model and forecast exchange rate volatility. This study aimed to model the volatility of the taka-US dollar exchange rate return using daily observations over a span of 7 years. The leptokurtic fat-tailed nature of the exchange rate return series usually establishes a rationale for using skewed distribution-such as Student's t-rather than normal distribution to estimate volatility models. The main focus was on whether the same nature exists in the taka-US dollar exchange rate return and whether the results were improved with Student's t-distribution. In particular, this study tried to model the volatility dynamics of the taka-US dollar exchange rate return using GARCH, APARCH, EGARCH, TGARCH, and IGARCH models. The findings from the models were compared under the regular normal distribution assumption for the residuals against Student's t-distribution. All of the models successfully passed the diagnostic tests when Student's t-distribution was used for the residuals. Further, in-sample estimation accuracy was observed to be improved when such a distribution was used. For modeling in-sample volatility dynamics, AR(2)-IGARCH(1, 1) was found to be the most accurate. We also tried to find the appropriate model for volatility forecasting, and we found that the forecasting accuracy of AR(2)-GARCH(1, 1) and AR(2)-IGARCH(1, 1) was improved with Student's t-distribution against normal distribution. However, in terms of out-of-sample forecasting accuracy, AR(2)-GARCH(1, 1) is considered the best compared to AR(2)-IGARCH(1, 1) since MAE, MAPE, and TI were all observed to be lower for the former model. That said, the issue of structural change presents a potential limitation of this exercise. During the data span, however, there was no regime shift in exchange rates, which minimizes its likelihood.   I  II  III  IV  I  II  III  IV  I  II   2013 2014 2015 RF ± 2 S.E.