Fit the data

fit = stan_glm(obesity ~ strengthTraining, data = exercise, refresh= 0)
print(fit, detail=FALSE)

##                  Median MAD_SD
## (Intercept)      61.8    3.8  
## strengthTraining -1.1    0.1  
## 
## Auxiliary parameter(s):
##       Median MAD_SD
## sigma 2.5    0.2

Note: refresh=0 suppresses the usual stan chain progress reports.

Chapter 8 will discuss fitting procedure, but given the fit, what can we learn from this?