In this article, we calculate response surface models for a large range of quantiles of the Kapetanios, Shin, and Snell (2003, Journal of Econometrics 112: 359–379) and Kapetanios and Shin (2008, Economics Letters 100: 377–380) tests for the null hypothesis of a unit root against the alternative—that the series of interest follows a globally stationary exponential smooth transition autoregressive process. The response surface models allow estimation of finite-sample critical values and approximate p-values for different combinations of the number of ob- servations, T, and the lag order in the test regression, p. The latter can be either specified by the user or optimally selected using a data-dependent procedure. We present the new commands kssur and ksur and illustrate their use with an empirical example.