I'd like to make a density plot of wages by marrital status (married/not married). This is what I have, but it's only providing a density plot for one variable.
nls80 = read.csv('https://raw.githubusercontent.com/bandcar/Examples/main/nls80(1).csv')
library(ggplot2)
ggplot(nls80, aes(x = wage, colour = married)) +
geom_density()
Suggestions on a fix?

