https://legendhjc.blogspot.com/2023/02/how-to-add-legend-to-plot-in-r.html