Regression and Other Stories: Restaurant

Demonstration of using Stan for optimization. See Appendix B in Regression and Other Stories.

---

Load packages

library("rprojroot")
root<-has_file(".ROS-Examples-root")$make_fix_file()
library("rstan")
rstan_options(auto_write = TRUE)

Show Stan code

writeLines(readLines(root("Restaurant","restaurant.stan")))

parameters {
  real<lower=0,upper=100> x;
}
model {
  target += (5000/x^2)*(x-11);
}

Compile Stan code

resto <- stan_model(root("Restaurant","restaurant.stan"))

Optimize

fit <- optimizing(resto)
print(fit)

$par
       x 
21.99954 

$value
[1] 113.6364

$return_code
[1] 0

$theta_tilde
            x
[1,] 21.99954