- data <- read.csv("PresidentHeight.csv", header = TRUE)
- data <- na.omit(data)
- HtDiff <- data[, 3] - data[, 5]
- par(mar = c(5, 5, 5, 7))
- plot(data[, 3], data[, 5], xlab = "Winner's height", las = 1, ylab = "Loser's height", xlim = c(160, 200), ylim = c(160, 200), pch = 18, col = "blue", yaxs = "i", xaxs = "i", font.lab = 2)
- rect(par("usr")[1], par("usr")[3], par("usr")[2], par("usr")[4], col = "gray")
- points(data[, 3], data[, 5], pch = 18, col = "blue")
- grid(8, 8, lty = 1, col = "black")
- abline(a = 0, b = 1, col = "red")
- abline(a = 10, b = 1, col = "red")
- abline(a = 20, b = 1, col = "red")
- abline(a = -10, b = 1, col = "red")
- abline(a = -20, b = 1, col = "red")
- mtext("Winner \n>20 cm \nshorter", side = 3, cex = 0.75, line = 0, at = c(170))
- mtext("Winner \n10-20 cm \nshorter", side = 3, cex = 0.75, line = 0, at = c(185))
- mtext("Winner \n0-10 cm \nshorter", side = 3, cex = 0.75, line = 0, at = c(195))
- mtext("Winner \n0-10 cm \ntaller", side = 4, cex = 0.75, line = 1, las = 1, at = c(195))
- mtext("Winner \n10-20 cm \ntaller", side = 4, cex = 0.75, line = 1, las = 1, at = c(185))
- mtext("Winner \n>20 cm \ntaller", side = 4, cex = 0.75, line = 1, las = 1, at = c(170))
- mtext("Same height", side = 4, cex = 0.75, line = 0, las = 1, at = c(202), col = "red")
复制代码
|