Sorry this is kinda incomplete answer and for python, nor R, but still. Scatter overlay with grouped barplot should be possible, see second example on this showcase page: https://plotly.com/python/graphing-multiple-chart-types/