pgr <- gray.colors(n = 25, start = 0.95, end = 0.0)
with (deepCmin,
filled.contour(x = multser, y = kseries, z = outcost,
ylab = "k (/day)", xlab = "multiplication factor (-)",
main = "Model cost landscape", col = pgr, nlevels = 25,
plot.axes = {
axis(1); axis(2);
points(optpar20$poppar[,2], optpar20$poppar[,1], pch = "o", cex = .5);
points(optpar25$poppar[,2], optpar25$poppar[,1], pch = "+", cex = 1);
points(optpar$par[2], optpar$par[1], pch = 16, cex = 2)
}
)
)
Run the code above in your browser using DataLab