data(CreelMN)
str(CreelMN)
head(CreelMN)
levels(CreelMN$species)
## ONLY RUN IN INTERACTIVE MODE
if (FALSE) {
require(FSA)
wae <- filterD(CreelMN,species=="WAE")
( waetbl <- table(wae$harvest) )
( waeptbl <- prop.table(waetbl)*100 )
( waercum <- rcumsum(waeptbl) )
op <- par(mfrow=c(1,2),mar=c(3.5,3.5,0.5,0.5),mgp=c(2,0.5,0))
barplot(waeptbl,xlab="Number of Walleye Harvested",ylab="Percentage of Individuals")
barplot(waercum,xlab="Minimum Number of Walleye Harvested",ylab="Percentage of Individuals")
par(op)
( ttlwae <- sum(wae$harvest) )
wae$svd5 <- wae$harvest-5 # Used to determine the number of fish
xtabs(~svd5,data=wae) # saved if limit was reduced to 5.
( svd5 <- sum(wae$svd5[wae$svd5>0]) ) # Only sum positive numbers (i.e., saved)
round(svd5/ttlwae*100,1) # Show as a percentage
wae$svd1 <- wae$harvest-1 # same but if reduced to 1 bag.
xtabs(~svd1,data=wae)
( svd1 <- sum(wae$svd1[wae$svd1>0]) )
round(svd1/ttlwae*100,1)
}
Run the code above in your browser using DataLab