curve(-0.2764464*x^3+0.0267897*x^2-0.0627775*x, -1, 1, xaxt='n', xlab="Cubic", ylab="Logged Time (min.)", yaxt='n', ylim=c(-0.3,0.4))
curve(0.0203028*x^3+0.1884012*x^2+0.1527854*x, -1, 1, add=T, xaxt='n', xlab="", ylab="", yaxt='n', ylim=c(-0.3,0.4))
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)
mtext('Cooking', side = 3, line=2.9, cex = 1.2)
mtext('American women', side = 3, line=2, cex = 0.8)
x1 <- c(seq(-1, 0, 0.01), seq(0, -1, -0.01))
y1 <- c(0.0203028*x1[1:101]^3+0.1884012*x1[1:101]^2+0.1527854*x1[1:101], -0.2764464*x1[102:202]^3+0.0267897*x1[102:202]^2-0.0627775*x1[102:202])
polygon(x1,y1, col=gray(0.9), border = NA)
x2 <- c(seq(0, 1, 0.01), seq(1, 0, -0.01))
y2 <- c(0.0203028*x2[1:101]^3+0.1884012*x2[1:101]^2+0.1527854*x2[1:101], -0.2764464*x2[102:202]^3+0.0267897*x2[102:202]^2-0.0627775*x2[102:202])
polygon(x2,y2, col=gray(0.9), border = NA)
curve(-0.1280718*x^3+0.1075954*x^2+0.045004*x, -1, 1, add=T, lwd=2, xaxt='n', xlab="", ylab="", yaxt='n', ylim=c(-0.3,0.4))
