Dear all,

I would like o to ask about the choice of the family to use to compute RR and ARR,

normally for RR I should choose these families in these lines

```
> str(mydata)
'data.frame': 190 obs. of 2 variables:
$ treat: int 1 1 1 1 1 1 1 1 1 1 ...
$ death: int 1 0 0 0 0 0 0 1 0 0 ...
> model_OR = stan_glm(death~treat, data = mydata,
+ family = binomial(link = "logit"),
+ prior = normal(location =0, scale = 10),
prior_intercept = normal(location =0, scale = 10) ,
QR=FALSE,
+ seed = 12345,
+ refresh = 0)
> OR = exp(coef(model_OR)[2])
> OR
treat
0.5749312
> model_RR = stan_glm(death~treat, data = mydata,
+ family = binomial(link = "log"),
+ prior = normal(location =0, scale = 10),
prior_intercept = normal(location =0, scale = 10) ,
QR=FALSE,
+ seed = 12345,
+ refresh = 0)
> RR = exp(coef(model_RR)[2])
> RR
treat
0.6721806
> model_ARR = stan_glm(death~treat, data = mydata,
+ family = binomial(link = "identity"),
+ prior = normal(location =0, scale = 10),
prior_intercept = normal(location =0, scale = 10) ,
QR=FALSE,
+ seed = 12345,
+ refresh = 0)
Error in stan_glm.fit(x = X, y = Y, weights = weights, offset = offset, :
'link' must be one of logit, probit, cauchit, log, cloglog
ARR = coef(model_ARR)[2]
```

Is it the correct families? if so, what is the correction for the last one?

and thanks in advance.