The regression was calculated in the open-source statistics software R (R Core Team 2018) using the glm.bn function from the MASS package. Results were expressed as rate ratios. Unpaired t tests of log-transformed data using the Holm-Šidák method for multiplicity adjustment were used for ...