A nonrgodic ground-motion model explicitly takes systematic local source, path and site effects on the predicted ground-motion into account. With an increasing number of ground-motion records, it s possible to estimate these effects. Landwehr et al. (2016) proposed a varying coefficient model as a tool to estimate nonergodic ground-motion models, based on Gaussian processes. Gaussian processes are computationally expen- sive, so for large data sets, some approximations have to be made. Here, we compare different Bayesian implementations of varying coefficient models, using the probabilistic programming language Stan (Carpenter et al., 2017), and the integrated nested Laplace approximation (INLA) (Rue et al., 2009). The models are used to fit nonergodic models on the California subset of the NGA-West2 data set (Ancheta et al., 2014). We find that both implementations lead to very similar results, both in the estimated parameters and the predicted ground motions.