RStan code