## Maps of some fish species
# Divide the plot window into 4 frames, 2 per row
par(mfrow = c(2,2))
# Plot four species
plot(spa,
asp = 1,
cex.axis = 0.8,
col = "brown",
cex = spe$Satr,
main = "Brown trout",
xlab = "x coordinate (km)",
ylab = "y coordinate (km)"
)
lines(spa, col = "light blue")
plot(spa,
asp = 1,
cex.axis = 0.8,
col = "brown",
cex = spe$Thth,
main = "Grayling",
xlab = "x coordinate (km)",
ylab = "y coordinate (km)"
)
lines(spa, col = "light blue")
plot(spa,
asp = 1,
cex.axis = 0.8,
col = "brown",
cex = spe$Baba,
main = "Barbel",
xlab = "x coordinate (km)",
ylab = "y coordinate (km)"
)
lines(spa, col = "light blue")
plot(spa,
asp = 1,
cex.axis = 0.8,
col = "brown",
cex = spe$Abbr,
main = "Common bream",
xlab = "x coordinate (km)",
ylab = "y coordinate (km)"
)
lines(spa, col = "light blue")
2.2 Data Exploration
17
# Divide the plot window into 4 frames, 2 per row
par(mfrow = c(2,2))
# Plot four species
plot(spa,
asp = 1,
cex.axis = 0.8,
col = "brown",
cex = spe$Satr,
main = "Brown trout",
xlab = "x coordinate (km)",
ylab = "y coordinate (km)"
)
lines(spa, col = "light blue")
plot(spa,
asp = 1,
cex.axis = 0.8,
col = "brown",
cex = spe$Thth,
main = "Grayling",
xlab = "x coordinate (km)",
ylab = "y coordinate (km)"
)
lines(spa, col = "light blue")
plot(spa,
asp = 1,
cex.axis = 0.8,
col = "brown",
cex = spe$Baba,
main = "Barbel",
xlab = "x coordinate (km)",
ylab = "y coordinate (km)"
)
lines(spa, col = "light blue")
plot(spa,
asp = 1,
cex.axis = 0.8,
col = "brown",
cex = spe$Abbr,
main = "Common bream",
xlab = "x coordinate (km)",
ylab = "y coordinate (km)"
)
lines(spa, col = "light blue")
2.2 Data Exploration
17
