Pulse Rates

fit_obj <- stan_glm(data = x, 
                    formula = pulse ~ gender + age, 
                    family = gaussian, 
                    refresh = 0,
                    seed = 9)

fit_obj

stan_glm
 family:       gaussian [identity]
 formula:      pulse ~ gender + age
 observations: 8563
 predictors:   3
------
            Median MAD_SD
(Intercept) 81.4    0.3  
genderMale  -4.0    0.3  
age         -0.1    0.0  

Auxiliary parameter(s):
      Median MAD_SD
sigma 11.7    0.1  

------
* For help interpreting the printed output see ?print.stanreg
* For info on the priors used see ?prior_summary.stanreg