This is a plot that I will need for my paper. I decided to publish it so that I don’t lose the code in the future.

par(mfrow=c(1,3), mar=c(4,4,2,1))
curve(-0.0936*x^3+0.08999*x^2+0.0098*x, from = -1, to=1, lwd=2, xaxt='n', xlab="a. Overall", ylab="Logged Time Spent on Cooking, Women", yaxt='n', ylim=c(-0.02,0.2))
axis(3)
mtext('low', side = 1, line=0.3, adj=0, cex = 0.8)
mtext('high', side = 1, line=0.3, adj=1, cex = 0.8)
mtext('Relative Economic Contribution', side = 1, line=1.8, cex = 0.8)
mtext('less', side = 2, line=1, adj=0, cex = 0.8)
mtext('more', side = 2, line=1, adj=1, cex = 0.8)
curve(-0.26014*x^3+0.0458*x^2+0.2569*x+0.09, from = -1, to=1, xaxt='n', xlab="b. Weekend", ylab="Logged Time Spent on Cooking, Women", yaxt='n', ylim=c(-0.02,0.2))
axis(3)
mtext('low', side = 1, line=0.3, adj=0, cex = 0.8)
mtext('high', side = 1, line=0.3, adj=1, cex = 0.8)
mtext('Relative Economic Contribution', side = 1, line=1.8, cex = 0.8)
mtext('less', side = 2, line=1, adj=0, cex = 0.8)
mtext('more', side = 2, line=1, adj=1, cex = 0.8)
curve(-0.03753*x^3+0.10597*x^2-0.0807*x, from = -1, to=1, lty='dashed', xaxt='n', xlab="c. Weekday", ylab="Logged Time Spent on Cooking, Women", yaxt='n', ylim=c(-0.02,0.2))
axis(3)
mtext('low', side = 1, line=0.3, adj=0, cex = 0.8)
mtext('high', side = 1, line=0.3, adj=1, cex = 0.8)
mtext('Relative Economic Contribution', side = 1, line=1.8, cex = 0.8)
mtext('less', side = 2, line=1, adj=0, cex = 0.8)
mtext('more', side = 2, line=1, adj=1, cex = 0.8)