summary(wffc)
with(wffc, table(water, session))
# Obtain some simple plots
waihou = subset(wffc, water == "Waihou")
waimak = subset(wffc, water == "Waimakariri")
whang = subset(wffc, water == "Whanganui")
otam = subset(wffc, water == "Otamangakau")
roto = subset(wffc, water == "Rotoaira")
minlength = min(wffc[,"length"])
maxlength = max(wffc[,"length"])
nwater = c("Waihou"=nrow(waihou), "Waimakariri"=nrow(waimak),
"Whanganui"=nrow(whang), "Otamangakau"=nrow(otam),
"Rotoaira"=nrow(roto))
par(mfrow=c(2,3), las=1)
# Overall distribution of length
with(wffc, boxplot(length/10 ~ water, ylim=c(minlength, maxlength)/10,
border="blue", main="Length (cm)", cex.axis=0.5))
# Overall distribution of LOG length
with(wffc, boxplot(length/10 ~ water, ylim=c(minlength, maxlength)/10,
border="blue", log="y", cex.axis=0.5,
main="Length (cm) on a log scale"))
# Overall distribution of number of captures
pie(nwater, border="blue", main="Proportion of captures",
labels=names(nwater), density=10, col=1:length(nwater),
angle=85+30* 1:length(nwater))
# Overall distribution of number of captures
with(wffc, barplot(nwater, main="Number of captures", cex.names=0.5,
col="lightblue"))
# Overall distribution of proportion of number of captures
with(wffc, barplot(nwater / sum(nwater), cex.names=0.5, col="lightblue",
main="Proportion of captures"))
# An interesting lake
with(roto, hist(length/10, xlab="Fish length (cm)", col="lightblue",
breaks=seq(18, 70, by=3), prob=TRUE, ylim=c(0, 0.08),
border="blue", ylab="", main="Lake Rotoaira", lwd=2))
Run the code above in your browser using DataLab