if (FALSE) fishlength <- seq(0.0, 0.72, by = 0.001)
plot(fishlength, wffc.P2star(fishlength), type = "l", col = "blue",
las = 1, lty = "dashed", lwd = 2, las = 1, cex.main = 0.8,
xlab = "Fish length (m)", ylab = "Competition points",
main = "Current (red) & proposed (blue&green) WFFC point system")
lines(fishlength, wffc.P1star(fishlength), col = "red", lwd = 2)
lines(fishlength, wffc.P3star(fishlength), col = "darkgreen",
lwd = 2, lty = "dashed")
lines(fishlength, wffc.P4star(fishlength), col = "orange",
lwd = 2, lty = "dashed")
abline(v = (1:4) * 0.18, lty = "dotted")
abline(h = (1:13) * wffc.P1star(0.18), lty = "dotted")
# Successive slopes:
(wffc.P1star((2:8)*0.18) - wffc.P1star((1:7)*0.18)) / (0.18 * 2000)
(wffc.P3star((2:8)*0.18) - wffc.P3star((1:7)*0.18)) / (0.18 * 2000)
(wffc.P4star((2:8)*0.18) - wffc.P4star((1:7)*0.18)) / (0.18 * 2000)
# More successive slopes:
MM2 <- 0.18 / 2
ind1 <- 2:12
(wffc.P4star((ind1)*MM2) - wffc.P4star((ind1-1)*MM2)) / (MM2 * 2000)
# About 50 percent more points:
wffc.P2 (2 * 0.18) / wffc.P1 (2 * 0.18)
wffc.P2star(2 * 0.18) / wffc.P1star(2 * 0.18)
Run the code above in your browser using DataLab