Transforming dataframe to plot interaction effects in R

Question

I have a dataframe that looks like this:

ppl <- data.frame(
  ID = sample(1:50),
  X = rnorm(50, 0.5),
  Y1 = rnorm(50),
  Y2 = rnorm(50, 1)
)