Plot.plot({
marks: [
Plot.linearRegressionY(data2, {x: "culmen_length_mm", y: "culmen_depth_mm", stroke: "species"}),
Plot.dot(data2, {x: "culmen_length_mm", y: "culmen_depth_mm", stroke: "species"}),
Plot.density(data2, {x: "culmen_length_mm", y: "culmen_depth_mm", stroke: "species", opacity: 0.2})
],
color: {legend: true}
})