The BATS model is related to the TBATS model. It can be fitted using the bats() command in the forecast package for R.