For dose - response analysis, we used the method proposed by Greenland and Longnecker [13] to estimate study - specific slopes from
the correlated natural logarithm of the RR across categories of coffee consumption, assigning to each class the dose corresponding to the midpoint of upper and lower boundaries.