The stochastic response surface method (SRSM) and the response surface method (RSM) are often used for structural reliability analysis, especially for reliability problems with implicit performance functions. This paper aims to compare these two methods in terms of fitting the performance function, accuracy and efficiency in estimating probability of failure as well as statistical moments of system output response. The computational procedures of two response surface methods are briefly introduced first. Then their capabilities are demonstrated and compared in detail through two examples. The results indicate that the probability of failure mainly reflects the accuracy of the response surface function (RSF) fitting the performance function in the vicinity of the design point, while the statistical moments of system output response reflect the accuracy of the RSF fitting the performance function in the entire space. In addition, the performance function can be well fitted by the SRSM with an optimal order polynomial chaos expansion both in the entire physical and in the independent standard normal spaces. However, it can be only well fitted by the RSM in the vicinity of the design point. For reliability problems involving random variables with approximate normal distributions, such as normal, lognormal, and Gumbel Max distributions, both the probability of failure and statistical moments of system output response can be accurately estimated by the SRSM, whereas the RSM can only produce the probability of failure with a reasonable accuracy.